Dsbev Dsbevd Dsbevx Dspev Dspevd Dspevx Dstev Dstevd Dstevr Dstevx Dsyev Dsyevd Dsyevr Dsyevx

Dsyevd() Sub Dsyevd ( Jobz As  String, Uplo As  String, N As  Long, A() As  Double, W() As  Double, Info As  Long  ) (分割統治法ドライバ) 固有値・固有ベクトル (対称行列) 目的 本ルーチンは実対称行列 A のすべての固有値, および, 必要により固有ベクトルを求める. 固有値のみ求める場合, QL法またはQR法を使用する. 固有ベクトルも求める場合, 分割統治法を使用する. 引数 [in] Jobz = "N": 固有値のみ求める. = "V": 固有値と固有ベクトルを求める. [in] Uplo = "U": A の上三角部分を格納する. = "L": A の下三角部分を格納する. [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in,out] A() 配列 A(LA1 - 1, LA2 - 1) (LA1 >= N, LA2 >= N) [in] 対称行列 A. Uplo = "U" の場合, A() の N x N 上三角部分に行列 A の上三角部分を格納する. Uplo = "L" の場合, A() の N x N 下三角部分に行列 A の下三角部分を格納する. [out] Jobz = "V": Info = 0 の場合, 行列 A の正規直交固有ベクトルを A() に返す.   Jobz = "N": A() の下三角部分(Uplo = "L" の場合)あるいは上三角部分(Uplo = "U" の場合)は対角部分を含め壊される. [out] W() 配列 W(LW - 1) (LW >= N) Info = 0 の場合, 求められた固有値(昇順). [out] Info = 0: 正常終了. = -1: パラメータ Jobz の誤り. (Jobz <> "V" および "N") = -2: パラメータ Uplo の誤り. (Uplo <> "U" および "L") = -3: パラメータ N の誤り. (N < 0) = -4: パラメータ A() の誤り. = -5: パラメータ W() の誤り. = i > 0: Jobz = "N" の場合: 収束しなかった. 中間結果の3重対角形の非対角要素のうち i 個が 0 に収束しなかった.   Jobz = "V" の場合: 行および列 i/(N+1) 〜 mod(i, N+1) の小行列においてアルゴリズムが固有値計算に失敗した. 出典 LAPACK 使用例 対称行列 A の固有値・固有ベクトルを求める. ただし, ( 2.20 -0.11 -0.32 ) A = ( -0.11 2.93 0.81 ) ( -0.32 0.81 2.37 ) とする. Sub Ex_Dsyevd() Const N = 3 Dim A(N - 1, N - 1) As Double, W(N - 1) As Double, Info As Long 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 Dsyevd("V", "L", N, A(), W(), Info) Debug.Print "Eigenvalues =", W(0), W(1), W(2) Debug.Print "Eigenvectors =" Debug.Print A(0, 0), A(0, 1), A(0, 2) Debug.Print A(1, 0), A(1, 1), A(1, 2) Debug.Print A(2, 0), A(2, 1), A(2, 2) Debug.Print "Info =", Info End Sub 実行結果 Eigenvalues = 1.70705954911046 2.22943643244226 3.56350401844728 Eigenvectors = 0.399322077213382 0.894521385341403 0.200931256445799 -0.481026444340548 0.390994588677271 -0.78468897753835 0.780484105216166 -0.216690384632057 -0.586421212707156 Info = 0