Zhbtrd Zhetrd Zhptrd Zpteqr Zstedc Zstein Zstemr Zsteqr Zungtr Zunmtr Zupgtr Zupmtr

Zstedc() Sub Zstedc ( Compz As  String, N As  Long, D() As  Double, E() As  Double, Z() As  Complex, Info As  Long  ) エルミート行列を変換した対称三重対角行列の固有値・固有ベクトル (分割統治法) 目的 本ルーチンはエルミート行列を変換した対称三重対角行列のすべての固有値, および, 必要により固有ベクトルを分割統治法により求める. Zhetrd, Zhptrd または Zhbtrd を使用して三重対角形に変換した場合は, 元のエルミート行列またはエルミート帯行列の固有ベクトルを求めることもできる. 引数 [in] Compz = "N": 固有値のみ求める. = "I": 三重対角行列の固有値・固有ベクトルを求める. = "V": 元のエルミート行列の固有値・固有ベクトルを求める. 元の行列から三重対角形への変換を行ったユニタリ行列を Z() に入れておかなければならない. [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() は破壊される. [in,out] Z() 配列 Z(LZ1 - 1, LZ2 - 1) (LZ1 >= N, LZ2 >= N) [in] Compz = "V" の場合, 三重対角形への変換を行ったユニタリ行列を Z() に入れる. [out] Info = 0 であれば, Compz = "V" の場合, 元のエルミート行列の正規直交化した固有ベクトルが Z() に入る. Compz = "I" の場合, 対称三重対角行列の正規直交化した固有ベクトルが Z() に入る. Compz = "N" の場合, Z() は参照されない. [out] Info = 0: 正常終了. = -1: パラメータ Compz の誤り. (Compz <> "N", "I" および "V") = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ D() の誤り. = -4: パラメータ E() の誤り. = -5: パラメータ Z() の誤り. = i > 0: 行および列 i/(N + 1) 〜 mod(i, N + 1) にある小行列の計算において固有値を求めることができなかった. 出典 LAPACK