Ddisna Dopgtr Dopmtr Dorgtr Dormtr Dpteqr Dsbtrd Dsptrd Dstebz Dstedc Dstein Dstemr Dsteqr Dsterf Dsytrd

Dsterf() Sub Dsterf ( N As  Long, D() As  Double, E() As  Double, Info As  Long  ) 対称三重対角行列の固有値 (QL/QR法) 目的 本ルーチンは対称三重対角行列のすべての固有値を Pal-Walker-Kahan による QL または QR法により求める. 引数 [in] N 行列の行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in,out] D() 配列 D(LD - 1) (LD >= N) [in] 三重対角行列の対角要素. [out] Info = 0 であれば固有値が昇順に入る. [in,out] E() 配列 E(LE - 1) (LE >= N - 1) [in] 三重対角行列の N - 1 個の副対角要素. [out] E() は破壊される. [out] Info = 0: 正常終了. = -1: パラメータ N の誤り. (N < 0) = -2: パラメータ D() の誤り. = -3: パラメータ E() の誤り. = i > 0: 30*N回以内の反復ですべての固有値を見つけることはできなかった. E() の i 個の要素が 0 に収束しなかった. 出典 LAPACK 使用例 対称行列 A の固有値を求める. ただし, ( 2.20 -0.11 -0.32 ) A = ( -0.11 2.93 0.81 ) ( -0.32 0.81 2.37 ) とする. Dsytrdで対称三重対角形に変換したのち, Dsterfを適用する. Sub Ex_Dsytrd_Dsterf() Const N = 3 Dim A(N - 1, N - 1) As Double, Info As Long Dim D(N - 1) As Double, E(N - 2) As Double, Tau(N - 2) As Double A(0, 0) = 2.2 A(1, 0) = -0.11: A(1, 1) = 2.93 A(2, 0) = -0.32: A(2, 1) = 0.81: A(2, 2) = 2.37 Call Dsytrd("L", N, A(), D(), E(), Tau(), Info) If Info <> 0 Then Debug.Print "Error in Dsytrd: Info =", Info Exit Sub End If Call Dsterf(N, D(), E(), Info) If Info <> 0 Then Debug.Print "Error in Dsterf: Info =", Info Exit Sub End If Debug.Print "Eigenvalues =", D(0), D(1), D(2) End Sub 実行結果 Eigenvalues = 1.70705954911046 2.22943643244226 3.56350401844728