_dsbev _dsbevd _dsbevx _dspev _dspevd _dspevx _dstev _dstevd _dstevr _dstevx _dsyev _dsyevd _dsyevr _dsyevx

_dstevd() void _dstevd ( char  jobz, int  n, double  d[], double  e[], int  ldz, double  z[], double  work[], int  lwork, int  iwork[], int  liwork, int *  info  ) (分割統治法ドライバ) 固有値・固有ベクトル (対称3重対角行列) 目的 本ルーチンは実対称3重対角行列のすべての固有値, および, 必要により固有ベクトルを求める. 固有ベクトルを求める場合には分割統治法を使用する. 引数 [in] jobz = 'N': 固有値のみ求める. = 'V': 固有値と固有ベクトルを求める. [in] n 行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in,out] d[] 配列 d[ld] (ld >= n) [in] 対称3重対角行列 A の n 個の対角要素. [out] info = の場合, 求められた固有値(昇順). [in,out] e[] 配列 e[le] (le >= n - 1) [in] 対称3重対角行列 A の n - 1 個の副対角要素(e[0] 〜 e[n-2]に格納). [out] e[] の内容は壊される. [in] ldz 配列z[][]の整合寸法. (ldz >= 1 (jobz='N'), ldz >= max(1,n) (jobz='V')) [out] z[][] 配列 z[lz][ldz] (lz >= n) jobz = 'V': info = 0 の場合, z[][] に正規直交固有ベクトルが入る. d[i] に関連する固有ベクトルが i 行に入る. jobz = 'N': z[][] は参照されない. [out] work[] 配列 work[lwork] 作業領域. info = 0の場合, work[0]にlworkの最適値を返す. [in] lwork 配列 work[]のサイズ (lwork >= 1 (n <= 1 の場合), 1 (jobz = 'N' の場合), n^2 + 4*n + 1 (jobz = 'V' の場合)) lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[]とiwork[]の最適サイズを求める計算だけを行い, work[0]とiwork[0]にそれらの値を返す. [out] iwork[] 配列 iwork[liwork] 整数作業領域. info = 0の場合, iwork[0]にliworkの最適値を返す. [in] liwork 配列 iwork[]のサイズ (liwork >= 1 (n <= 1 の場合), 1 (jobz = 'N' の場合), 5*n + 3 (jobz = 'V' の場合)) lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[]とiwork[]の最適サイズを求める計算だけを行い, work[0]とiwork[0]にそれらの値を返す. [out] info = 0: 正常終了 = -1: 入力パラメータ jobz の誤り (jobz != 'V'および'N') = -2: 入力パラメータ n の誤り (n < 0) = -5: 入力パラメータ ldz の誤り (ldzが小さすぎる) = -8: 入力パラメータ lwork の誤り. (lworkが小さすぎる) = -10: 入力パラメータ liwork の誤り. (liworkが小さすぎる) = i > 0: 収束しなかった. 非対角要素 e[] のうち i 個が 0 にならなかった. 出典 LAPACK