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

Dsbevx() Sub Dsbevx ( Jobz As  String, Range As  String, Uplo As  String, N As  Long, Kd As  Long, Ab() As  Double, Q() As  Double, Vl As  Double, Vu As  Double, Il As  Long, Iu As  Long, AbsTol As  Double, M As  Long, W() As  Double, Z() As  Double, IFail() As  Long, Info As  Long  ) (エキスパートドライバ) 固有値・固有ベクトル (対称帯行列) 目的 本ルーチンは実対称帯行列 A の選択された固有値, および, 必要により固有ベクトルを求める. 必要な固有値の範囲あるいは番号の範囲を指定することにより, 求める固有値・固有ベクトルを選択することができる. 引数 [in] Jobz = "N": 固有値のみ求める. = "V": 固有値と固有ベクトルを求める. [in] Range = "A": すべての固有値を求める. = "V": 半開区間(vl, vu]のすべての固有値を求める. = "I": il番目からiu番目までの固有値を求める. [in] Uplo = "U": Aの上三角部分を格納する. = "L": Aの下三角部分を格納する. [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in] Kd 上帯幅(Uplo = "U"の場合)あるいは下帯幅(Uplo = "L"の場合). (Kd >= 0) [in,out] Ab() 配列 Ab(LAb1 - 1, LAb2 - 1) (LAb1 >= Kd + 1, LAb2 >= N) [in] Kd+1×N対称帯行列形式のN×N対称帯行列 A. Uploに従って上または下三角部分を格納する. [out] 配列Ab()は三重対角形への変換中に生成される値で上書きされる.   Uplo = "U": 三重対角行列Tの上副対角要素および対角要素をAb()のKdおよびKd+1番目の行に返す.   Uplo = "L": Tの対角要素および下副対角要素をAb()の最初の2行に返す. [out] Q() 配列 Q(LQ1 - 1, LQ2 - 1) (LQ1 >= N, LQ2 >= N) Jobz = "V": 3重対角行列に変換する際に使われたN×N直交行列. Jobz = "N": 配列q[][]は参照されない. [in] Vl Range = "V": 固有値を求める区間の下端. (Vl < Vu) Range = "A"または"I": 参照されない. [in] Vu Range = "V": 固有値を求める区間の上端. (Vl < Vu) Range = "A"または"I": 参照されない. [in] Il Range = "I": 求める最小固有値の番号. (1 <= Il <= Iu <= N (N > 0 の場合), Il = 1, Iu = 0 (N = 0 の場合)) Range = "A"または"V": 参照されない. [in] Iu Range = "I": 求める最大固有値の番号. (1 <= Il <= Iu <= N (N > 0 の場合), Il = 1, Iu = 0 (N = 0 の場合)) Range = "A"または"V": 参照されない. [in] AbsTol 固有値の絶対誤差限界. 固有値の近似値は区間[a, b]に入っているときに収束したものとみなされる. この区間の幅は AbsTol + eps*max(|a|, |b|) に等しいかこれより小さい. ここでepsはマシンイプシロンである. AbsTol <= 0 の場合, eps*|T|が代わりに使用される. ここで, |T|は行列Aを変換して得られた3重対角行列の1ノルムである. AbsTolを0ではなくアンダフロー限界の2倍(2*Dlamch("S"))に設定したときに固有値が最も正確に求められる. Info > 0 (固有ベクトルのいくつかが収束しなかったことを示す)で戻ったときには, AbsTolを2*Dlamch("S")に設定してみるとよい. "Computing Small Singular Values of Bidiagonal Matrices with Guaranteed High Relative Accuracy," by Demmel and Kahan, LAPACK Working Note #3. を参照のこと. [out] M 求められた固有値の数 (0 <= M <= N) Range = "A"の場合 M = N, Range = "I"の場合 M = Iu - Il + 1 となる. [out] W() 配列 W(LW - 1) (LW >= N) 正常終了時, M個の求められた固有値が先頭から昇順に入る. [out] Z() 配列 Z(LZ1 - 1, LZ2 - 1) (LZ1 >= N, LZ2 >= M) Jobz = "V": Info = 0 の場合, 求められた固有値に対応してZ()の最初のM列に行列Aの正規直交固有ベクトルが入る. W(i)に関連する固有ベクトルがZ()のi列に入る. 固有ベクトルの収束に失敗した場合, Z()のその列には固有ベクトルの最終近似が入り, Ifail()に固有ベクトルの番号が入る. Jobz = "N": Z()は参照されない. 注: 配列Z()は少なくてもmax(1, M)列を割り当てること. Range = "V"の場合, Mの値をあらかじめ知ることはできないが上限値を使用すること. [out] IFail() 配列 IFail(LIFail - 1) (LIFail >= N) Jobz = "V": Info = 0の場合, IFail()の最初のM要素が0に設定される. Info > 0 の場合, 収束しなかった固有ベクトルの番号がIFail()に入る. Jobz = "N": IFail()は参照されない. [out] Info = 0: 正常終了. = -1: パラメータ Jobz の誤り. (Jobz <> "V"および"N") = -2: パラメータ Range の誤り. (Range <> "A", "V"および"I") = -3: パラメータ Uplo の誤り. (Uplo <> "U"および"L") = -4: パラメータ N の誤り. (N < 0) = -5: パラメータ Kd の誤り. (Kd < 0) = -6: パラメータ Ab() の誤り. = -7: パラメータ Q() の誤り. = -9: パラメータ Vu の誤り. (Vu < Vl) = -10: パラメータ Il の誤り. (Il < 1 または Il > N) = -11: パラメータ Iu の誤り. (Iu < min(N, Il) または Iu > N) = -14: パラメータ W() の誤り. = -15: パラメータ Z() の誤り. = -16: パラメータ IFail() の誤り. = i > 0: i個の固有ベクトルが収束しなかった. それらの番号は配列IFail()に入る. 出典 LAPACK 使用例 対称帯行列Aの固有値・固有ベクトルを求める. ただし, ( 0.61 0.79 0 ) A = ( 0.79 2.23 0.25 ) ( 0 0.25 2.87 ) とする. Sub Ex_Dsbevx() Const N = 3, Kd = 1 Dim Ab(Kd, N - 1) As Double, W(N - 1) As Double, Q(N - 1, N - 1) As Double Dim Vl As Double, Vu As Double, Il As Long, Iu As Long, AbsTol As Double Dim M As Long, Z(N - 1, N - 1) As Double, IFail(N - 1) As Long, Info As Long Ab(0, 0) = 0.61: Ab(0, 1) = 2.23: Ab(0, 2) = 2.87 Ab(1, 0) = 0.79: Ab(1, 1) = 0.25 Call Dsbevx("V", "A", "L", N, Kd, Ab(), Q(), Vl, Vu, Il, Iu, AbsTol, M, W(), Z(), IFail(), Info) Debug.Print "Eigenvalues =", W(0), W(1), W(2) Debug.Print "Eigenvectors =" Debug.Print Z(0, 0), Z(0, 1), Z(0, 2) Debug.Print Z(1, 0), Z(1, 1), Z(1, 2) Debug.Print Z(2, 0), Z(2, 1), Z(2, 2) Debug.Print "M =", M, "Info =", Info End Sub 実行結果 Eigenvalues = 0.285074336232239 2.43057331973532 2.99435234403244 Eigenvectors = 0.924204695519984 -0.35289670385953 0.145978070900713 -0.380123828036851 -0.813258637555945 0.440586272822087 3.67635163909111E-02 0.462681645244035 0.885761897474061 M = 3 Info = 0