作者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/m.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