ddisna dopgtr dopmtr dorgtr dormtr dpteqr dsbtrd dsptrd dstebz dstedc dstein dstemr dsteqr dsterf dsytrd

dpteqr() void dpteqr ( char  compz, int  n, double  d[], double  e[], int  ldz, double  z[], double  work[], int *  info  ) 正定値対称三重対角行列の固有値・固有ベクトル 目的 本ルーチンは正定値対称三重対角行列のすべての固有値, および, 必要により固有ベクトルを求める. 最初に dpttrf を使って行列を分解し, 次に二重対角形の特異値を求めるために dbdsqr を呼び出すことにより, 固有値・固有ベクトルを求める. 本ルーチンは正定値三重対角行列の固有値を高い精度で計算する. すなわち, 固有値の範囲が何桁にもおよぶ場合に小さな固有値とその固有ベクトルが, 例えば通常のQR法よりも, 正確に計算されることを意味する. dsytrd, dsptrd または dsbtrd を使用して三重対角形式に変換した場合には, 元の正定値対称行列または正定値対称帯行列の固有ベクトルを求めることもできる. (ただし, 元の行列の小さな固有値の範囲が何桁にもおよぶ場合, 三重対角行列に変換すると, これらの固有値について高い精度が得られなくなる可能性がある.) 引数 [in] compz = 'N': 固有値のみ求める. = 'V': 元の対称行列の固有値・固有ベクトルを求める. 元の行列から三重対角形への変換を行った直交行列を配列 z[][] に入れておく. = 'I': 三重対角行列の固有値・固有ベクトルを求める. [in] n 行列の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in,out] d[] 配列 d[ld] (ld >= n) [in] 三重対角行列の対角要素. [out] 固有値が降順に入る. [in,out] e[] 配列 e[le] (le >= n - 1) [in] 三重対角行列の n - 1 個の副対角要素. [out] e[] は破壊される. [in] ldz 二次元配列 z[][] の整合寸法. (ldz >= 1 (compz = 'N' の場合), ldz >= max(1, n) (compz = 'V' または 'I' の場合)) [in,out] z[][] 配列 z[lz][ldz] (lz >= n) [in] compz = 'V' の場合, 三重対角形への変換を行った直交行列. [out] compz = 'V' の場合, 元の対称行列の正規直交化した固有ベクトル. compz = 'I' の場合, 三重対角行列の正規直交化した固有ベクトル. info > 0 で終了した場合, 求められた固有値に対応する固有ベクトルのみが z[][] に入る. compz = 'N' の場合, z[][] は参照されない. [out] work[] 配列 work[lwork] (lwork >= max(1, 4*n)) 作業領域. [out] info = 0: 正常終了. = -1: 入力パラメータ compz の誤り. (compz != 'N', 'V' および 'I'') = -2: 入力パラメータ n の誤り. (n < 0) = -5: 入力パラメータ ldz の誤り. (ldz < 1 または (compz = 'V' または 'I' かつ ldz < max(1, n))) = i > 0: i <= n の場合, i 番目の首座小行列式が正定値でなかったため行列のコレスキー分解ができなかった. i > n の場合, SVD アルゴリズムが収束しなかった. 二重対角形の i - n 個の非対角要素が 0 に収束しなかった. 出典 LAPACK