Zhbev Zhbevd Zhbevx Zheev Zheevd Zheevr Zheevx Zhpev Zhpevd Zhpevx

Zhbevx() Sub Zhbevx ( Jobz As  String, Range As  String, Uplo As  String, N As  Long, Kd As  Long, Ab() As  Complex, Q() As  Complex, Vl As  Double, Vu As  Double, Il As  Long, Iu As  Long, AbsTol As  Double, M As  Long, W() As  Double, Z() As  Complex, 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()は三重対角形への変換中に生成される値で上書きされる. [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の固有値・固有ベクトルを求める. ただし, ( 2.20 -0.32-0.81i 0 ) A = ( -0.32+0.81i 2.11 0.37+0.80i ) ( 0 0.37-0.80i 2.93 ) とする. Sub Ex_Zhbevx() Const N = 3, Kd = 1 Dim Ab(Kd, N - 1) As Complex, W(N - 1) As Double, Q(N - 1, N - 1) As Complex 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 Complex, IFail(N - 1) As Long, Info As Long Ab(0, 0) = Cmplx(2.2, 0): Ab(0, 1) = Cmplx(2.11, 0): Ab(0, 2) = Cmplx(2.93, 0) Ab(1, 0) = Cmplx(-0.32, 0.81): Ab(1, 1) = Cmplx(0.37, -0.8) Call Zhbevx("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 Creal(Z(0, 0)), Cimag(Z(0, 0)), Creal(Z(0, 1)), Cimag(Z(0, 1)) Debug.Print Creal(Z(1, 0)), Cimag(Z(1, 0)), Creal(Z(1, 1)), Cimag(Z(1, 1)) Debug.Print Creal(Z(2, 0)), Cimag(Z(2, 0)), Creal(Z(2, 1)), Cimag(Z(2, 1)) Debug.Print Creal(Z(0, 2)), Cimag(Z(0, 2)) Debug.Print Creal(Z(1, 2)), Cimag(Z(1, 2)) Debug.Print Creal(Z(2, 2)), Cimag(Z(2, 2)) Debug.Print "M =", M, "Info =", Info End Sub 実行結果 Eigenvalues = 1.04283948355918 2.52504447979701 3.67211603664382 Eigenvectors = -0.563395195782213 0 0.745522345299636 0 -0.275043343741226 0.696203463844978 -0.102234588437164 0.258781301981572 -0.241207215775279 -0.253094504921304 -0.417819329339414 -0.438410500970378 -0.356065002532476 -0 0.221139195914327 -0.559758589658141 -0.493164614651588 -0.51746898859873 M = 3 Info = 0