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

_dsbev() void _dsbev ( char  jobz, char  uplo, int  n, int  kd, int  ldab, double  ab[], double  w[], int  ldz, double  z[], double  work[], int *  info  ) (シンプルドライバ) 固有値・固有ベクトル (対称帯行列) 目的 本ルーチンは実対称帯行列 A のすべての固有値, および, 必要により固有ベクトルを求める. 固有値および固有ベクトルはQL法またはQR法により求める. 引数 [in] jobz = 'N': 固有値のみ求める. = 'V': 固有値と固有ベクトルを求める. [in] uplo = 'U': A の上三角部分を格納する. = 'L': A の下三角部分を格納する. [in] n 行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] kd 上帯幅(uplo = 'U' の場合)あるいは下帯幅(uplo = 'L' の場合). (kd >= 0) [in] ldab 二次元配列 ab[][] の整合寸法. (ldab >= kd + 1) [in,out] ab[][] 配列 ab[lab][ldab] (lab >= n) [in] 対称帯行列 A の上または下三角部分を配列の最初の kd+1 列に格納する. A のj列が配列 ab のj行に次のように格納される. uplo = 'U': ab[j][kd + i - j] = A(i, j). ただし, max(0, j - kd - 1) <= i <= j <= n - 1. uplo = 'L': ab[j][i - j] = A(i, j). ただし, 0 <= j <= i <= min(n - 1, j + kd - 1). [out] ab[][] は三重対角形への変換中に生成される値で上書きされる. uplo = 'U' の場合, 三重対角行列 T の上副対角要素および対角要素を ab[][] の kd および kd+1 番目の列に返す. uplo = 'L' の場合, T の対角要素および下副対角要素を ab[][] の最初の2列に返す. [out] w[] 配列 w[lw] (lw >= n) info = 0 の場合, 求められた固有値(昇順). [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[][] に正規直交固有ベクトルが入る. w[i] に関連する固有ベクトルが i 行に入る. jobz = 'N': z[][] は参照されない. [out] work[] 配列 work[lwork] (lwork >= max(1, 3*n - 2)) 作業領域. [out] info = 0: 正常終了 = -1: 入力パラメータ jobz の誤り (jobz != 'V' および 'N') = -2: 入力パラメータ uplo の誤り (uplo != 'U' および 'L') = -3: 入力パラメータ n の誤り (n < 0) = -4: 入力パラメータ kd の誤り (kd < 0) = -5: 入力パラメータ ldab の誤り (ldab < kd + 1) = -8: 入力パラメータ ldz の誤り (ldz が小さすぎる) = i > 0: 収束しなかった. 中間結果の3重対角形の非対角要素のうち i 個が 0 にならなかった. 出典 LAPACK