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

dstein() void dstein ( int  n, double  d[], double  e[], int  m, double  w[], int  iblock[], int  isplit[], int  ldz, double  z[], double  work[], int  iwork[], int  ifail[], int *  info  ) 対称三重対角行列の固有ベクトル (逆反復法) 目的 本ルーチンは実対称三重対角行列 T の指定された固有値に対応する固有ベクトルを逆反復法により求める. 各固有ベクトルの計算に許される最大反復回数は内部パラメータ maxits (現在 5 に設定されている) で指定される. 引数 [in] n 三重対角行列の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] d[] 配列 d[ld] (ld >= n) 三重対角行列 T の n 個の対角要素. [in] e[] 配列 e[le] (le >= n - 1) 三重対角行列 T の n - 1 個の副対角要素. [in] m 求める固有ベクトルの数 (0 <= m <= n) [in] w[] 配列 w[lw] (lw >= n) w[] の先頭 m 個にその固有ベクトルを求める固有値を入れる. 固有値は分割されたブロックごとにグループ分けされ, ブロック内で最小から最大の順に並べられていること. (order = 'B' とした dstebz の出力である配列 w[] が想定されている.) [in] iblock[] 配列 iblock[liblock] (liblock >= n) w[] 内の固有値に対応する小行列のインデックス. 固有値 w[i] が先頭の小行列に属するとき iblock[i] = 1, w[i] が2つ目の小行列に属するときは 2, ... である. (dstebz の出力である配列 iblock[] が想定されている.) [in] isplit[] 配列 isplit[lisplit] (lisplit >= n) T が小行列に分割する分割点. 最初の小行列は行／列 1 〜 isplit[0], ２番目は行／列 isplit[0] + 1 〜 isplit[1], ... よりなる. (dstebz の出力である配列 isplit[] が想定されている.) [in] ldz 二次元配列 z[][] の整合寸法. (ldz >= max(1, n)) [out] z[][] 配列 z[lz][ldz] (lz >= max(1, m)) 求められた固有ベクトル. 固有値 w[i] に関連する固有ベクトルが z[] の i 行に入る. 収束に失敗したときは maxits 回の反復後時点のベクトルが設定される. [out] work[] 配列 work[lwork] (lwork >= 5*n) 作業領域. [out] iwork[] 配列 iwork[liwork] (liwork >= n) 整数作業領域. [out] ifail[] 配列 ifail[lifail] (lifail >= m) 正常終了時, ifail[] のすべての要素はゼロである. 1つあるいは複数の固有ベクトルが maxits 回の反復で収束しなかったならばそれらのインデックスが配列 ifail[] に格納される. [out] info = 0: 正常終了. = -1: 入力パラメータ n の誤り. (n < 0) = -4: 入力パラメータ m の誤り. (m < 0 または m > n) = -5: 入力パラメータ w の誤り. (要素が昇順でない) = -6: 入力パラメータ iblock の誤り. (要素が昇順でない) = -8: 入力パラメータ ldz の誤り. (ldz が小さすぎる) = i > 0: i 個の固有ベクトルが maxits 回の反復で収束しなかった. そのインデックスが ifail[] に格納される. 出典 LAPACK