[问题] Lapack Block Diagonization goldeagleroy PTT批踢踢实业坊

[问题] Lapack Block Diagonization

楼主: goldeagleroy (茶叶、咖啡、星巴克) 2017-09-15 04:08:00

Fortran 生手想请问，在LAPACK，有没有办法直接进行Block diagonization?
我有一个分块矩阵:
A B
B* -A*
矩阵是Hermitian，所以我直接用Zheev求解，但是Eigenvalue 却混在一起，
因为LAPACK的输出已经将Eigenvalue由小排到大了。
所以我想将矩阵先对角化:
C 0
0 D
再分别求解
有没有什么subroutine是我用得到的呢?

作者: PTT0207445 (我是乡民) 2017-11-12 06:24:00

A是实数对称矩阵吗?

继续阅读

有人会几A几B fortran95的写法吗？key000079 [问题] Fortran做扑克牌游戏willy7777 [问题] 简单的sum数值积分问题what357159 [问题] Fortran如何抓取指定间隔列的资料fermmy [问题] 字串输出入问题terrylove12 Re: [问题] Fortran 和Tecplot的CFD来绘图MHG Re: [问题] Fortran 和Tecplot的CFD来绘图howermath [问题] Fortran 和Tecplot的CFD来绘图becomesonumb [问题] 如何开一个每列长度不同的矩阵?daumda [问题] intel fortan里的libyooo4017