[問題] Lapack Block Diagonization

[goldeagleroy (茶葉、咖啡、星巴克)]

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
推 PTT0207445: A是實數對稱矩陣嗎?
推 PTT0207445: 不是所有的矩陣都能拆成塊對角,有特殊原因才可以,你應
→ PTT0207445: 該按照那個特殊原因來寫,把C和D矩陣產生出來再zheev。