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

dsteqr() void dsteqr ( char  compz, int  n, double  d[], double  e[], int  ldz, double  z[], double  work[], int *  info  ) 対称三重対角行列の固有値・固有ベクトル (QL/QR法) 目的 本ルーチンは対称行列のすべての固有値, および, 必要により固有ベクトルを Implicit QL または QR法により求める. dsytrd, dsptrd または dsbtrd を使用して三重対角形に変換した場合は, 元の対称行列または対称帯行列の固有ベクトルを求めることもできる. 引数 [in] compz = 'N': 固有値のみ求める. = 'V': 元の対称行列の固有値・固有ベクトルを求める. 元の行列から三重対角形への変換を行った直交行列を z[][] に入れておかなければならない. = 'I': 三重対角行列の固有値・固有ベクトルを求める. z[][] は単位行列に初期化される. [in] n 行列の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in,out] d[] 配列 d[ld] (ld >= n) [in] 三重対角行列の対角要素. [out] info = 0 であれば固有値が昇順に入る. [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' の場合, 三重対角形への変換を行った直交行列を z[][] に入れる. [out] info = 0 であれば, compz = 'V' の場合, 元の対称行列の正規直交化した固有ベクトルが z[][] に入る. compz = 'I' の場合, 対称三重対角行列の正規直交化した固有ベクトルが z[][] に入る. compz = 'N' の場合, z[][] は参照されない. [out] work[] 配列 work[lwork] (lwork >= max(1, 2*n - 2)) 作業領域. compz = 'N' の場合, work は参照されない. [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: 30*n 回以内の反復ですべての固有値を見つけることはできなかった. e[] の i 個の要素が 0 に収束しなかった. 終了時, d[] と e[] には元の行列に直交相似な対称三重対角行列の要素が入る. 出典 LAPACK