[問題] Lapack Block Diagonization
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://www.ptt.cc/bbs/Fortran/M.1505419682.A.E94.html
推
11/12 06:24, , 1F
11/12 06:24, 1F
推
11/29 07:06, , 2F
11/29 07:06, 2F
→
11/29 07:08, , 3F
11/29 07:08, 3F
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
