作者goldeagleroy (茶叶、咖啡、星巴克)
看板Fortran
标题[问题] Lapack Block Diagonization
时间Fri Sep 15 04:08:00 2017
Fortran 生手想请问，在LAPACK，有没有办法直接进行Block diagonization?
我有一个分块矩阵:
A B
B* -A*
矩阵是Hermitian，所以我直接用Zheev求解，但是Eigenvalue 却混在一起，
因为LAPACK的输出已经将Eigenvalue由小排到大了。
所以我想将矩阵先对角化:
C 0
0 D
再分别求解
有没有什麽subroutine是我用得到的呢?
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 192.108.69.15
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Fortran/M.1505419682.A.E94.html
1F：推 PTT0207445: A是实数对称矩阵吗? 11/12 06:24
2F：推 PTT0207445: 不是所有的矩阵都能拆成块对角,有特殊原因才可以,你应 11/29 07:06
3F：→ PTT0207445: 该按照那个特殊原因来写,把C和D矩阵产生出来再zheev。 11/29 07:08