Zgsev Zgsgv Znaupd Zneupd

Zneupd() Sub Zneupd ( Rvec As  Long, Howmny As  String, Selct() As  Long, D() As  Complex, Z() As  Complex, Sigma As  Complex, Workev() As  Complex, Bmat As  String, N As  Long, ByVal Which As  String, Nev As  Long, Resid() As  Complex, Ncv As  Long, V() As  Complex, IParam() As  Long, Ipntr() As  Long, Workd() As  Complex, Workl() As  Complex, RWork() As  Double, Info As  Long, Optional Tol As  Double = 0  ) アーノルディ分解の近似固有値・固有ベクトル (複素疎行列) (Arpack) 目的 本ルーチンは固有値問題 A*z = λ*B*z の近似固有値および(必要であれば)対応する近似固有ベクトルを求める. これらの値は Znaupd が返す複線形演算子 OP および行列 B のアーノルディ分解(詳細は Znaupd のマニュアルの「目的」の項を参照)より計算される. 本ルーチン呼び出しの前に Znaupd を呼び出しておかなければならない. 近似固有値は昇順に格納される. 引数 Bmat, N, Which, Nev, Tol, Resid, Ncv, V, Iparam, Ipntr, Workd, Workl および RWork は Zsaupd の最後の呼び出し後, 変更することなく本ルーチンに渡さなければならない. これらのうち V, Ipntr および Workl は本ルーチンによって変更される. 引数 [in] Rvec 固有値問題 A*z = λ*B*z の近似リッツ値に対応する不変部分空間の基底を求めるかどうか指定する. = 0 (False): リッツ値のみ求める. = 1 (True): リッツベクトルまたはシュールベクトルも求める. [in] Howmny 求めるリッツ値に対応する不変部分空間の基底の形式を指定する. = "A": Nev 個のリッツベクトルを求める. = "P": Nev 個のシュールベクトルを求める. = "S": 配列 Select() で指定されたリッツベクトルを求める. [out] Selct() 配列 Selct(LSelct - 1) (LSelct >= Ncv) Howmny = "S" の場合, Selct() は求めるリッツベクトルを指定する. リッツ値 D(j) に対応するリッツベクトルを選択するには, Selct(j) を 1 (True) に設定する必要がある. Howmny = "A" または "P" の場合, Selct() は内部ワークスペースとして使用される. [out] D() 配列 D(LD - 1) (LD >= Nev) A*z = λ*B*z の固有値のリッツ値による近似が入る. [out] Z() 配列 Z(Ldz - 1, LZ - 1) (Ldz >= N, LZ >= Nev) Rvec = 1 (Truw) かつ Howmny = "A" の場合, Z() の列は 固有値問題 A*z = λ*B*z の nconv (Iparam(4)) 個のリッツ値に対応する近似固有ベクトル(リッツベクトル)を表す. Rvec = 0 (False) または Howmny = "P" の場合, Z() は参照されない. [in] Sigma mode = 3 の場合, シフトの実数部を表す. mode = 1 または 2 の場合, 参照されない. [out] Workev() 配列 Workev[LWorkev] (LWorkev >= 2*Ncv) ワークスペース. [in] Bmat Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] N Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] Which Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] Nev Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] Tol Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] Resid() 配列 Resid(LResid - 1) (LResid >= N) Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in] Ncv Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in,out] V() 配列 V(Ldv - 1, LV - 1) (Ldv >= N, LV >= Ncv) [in] Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [out] Rvec = 1 (True) であれば, はじめの nconv 列に求められる不変部分空間を張る近似シュールベクトルが入る. [in] Iparam() 配列 Iparam(LIparam - 1) (LIparam >= 11) Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in,out] Ipntr() 配列 Ipntr(LIpntr - 1) (LIpntr >= 14) [in] Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [out] Ipntr(8) 〜 Ipntr(12) には Workl へのポインターが入る. Ipntr(8): 元の式の Ncv 個のリッツ値へのポインター. Ipntr(9): 使用されない. Ipntr(10): Ncv 個のそれぞれの近似誤差へのポインター. Ipntr(11): H の Ncv ｘ Ncv 上三角シュール行列へのポインター. Ipntr(12): 上ヘッセンベルグ行列 H の固有ベクトルからなる Ncv ｘ Ncv 行列へのポインター. Rvec = 1 (True) のときにのみ参照される. [in] Workd() 配列 Workd(LWorkd - 1) (LWorkd >= 3*N) Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [in,out] Workl() 配列 Workl(LWorkl - 1) (LWorkl >= 3*Ncv^2 + 5*Ncv) [in] Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. Workl(0) 〜 Workl(Ncv^2 + 2*Ncv - 1) には Znaupd で得られた情報が入っている. これらは本ルーチンによって変更されることはない. [out] Workl(Ncv^2 + 2*Ncv) 〜 Workl(3*Ncv^2 + 4*Ncv - 1) には変換していないリッツ値, H の上疑三角行列, および H の関連する不変部分空間の行列表現の実数部および虚数部が入る. [in] RWork() 配列 RWork[LRWork] (LRWork >= Ncv) Znaupd と同じ. Znaupd の最終呼び出し後変更せずに渡すこと. [out] Info リターンコード. = 0: 正常終了. < 0: (-Info)番目の入力パラメータの誤り. = 11: 初期残差ベクトルが 0 である. = 12: Znaupd は十分な精度の固有値を得られなかった. = 13: Zneupd が得た変換後のリッツ値の数が Znaupd が得た数と異なる. = 14: LAPACK ルーチンで固有値の計算に失敗した. = 15: 固有ベクトルの計算に失敗した (LAPACK ルーチン ZTREVC). 出典 ARPACK-NG (https://github.com/opencollab/arpack-ng)