[問題] 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