XLPack 6.1
C/C++ API リファレンスマニュアル
読み取り中…
検索中…
一致する文字列を見つけられません

◆ _dsyev()

void _dsyev ( char  jobz,
char  uplo,
int  n,
int  lda,
double  a[],
double  w[],
double  work[],
int  lwork,
int *  info 
)

(シンプルドライバ) 固有値・固有ベクトル (対称行列)

目的
本ルーチンは実対称行列 A のすべての固有値, および, 必要により固有ベクトルを求める. 固有値および固有ベクトルはQL法またはQR法により求める.
引数
[in]jobz= 'N': 固有値のみ求める.
= 'V': 固有値と固有ベクトルを求める.
[in]uplo= 'U': A の上三角部分を格納する.
= 'L': A の下三角部分を格納する.
[in]n行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]lda二次元配列 a[][] の整合寸法. (lda >= max(1, n))
[in,out]a[][]配列 a[la][lda] (la >= n)
[in] 対称行列 A. uplo = 'U' の場合, a[][] の n x n 上三角部分に行列 A の上三角部分を格納する. uplo = 'L' の場合, a[][] の n x n 下三角部分に行列 A の下三角部分を格納する.
[out] jobz = 'V': info = 0の場合, 行列 A の正規直交固有ベクトルを a[][] に返す.
  jobz = 'N': a[][] の下三角部分(uplo = 'L'の場合)あるいは上三角部分(uplo = 'U'の場合)は対角部分を含め壊される.
[out]w[]配列 w[lw] (lw >= n)
info = 0 の場合, 求められた固有値(昇順).
[out]work[]配列 work[lwork]
作業領域.
info = 0 の場合, work[0] に lwork の最適値を返す.
[in]lwork配列 work[] のサイズ (lwork >= max(1, 3*n-1))
最適効率のためには lwork >= (nb + 2)*n とせよ. ただし, nb は最適ブロックサイズである.
lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[] の最適サイズを求める計算だけを行い, work[0] にその値を返す.
[out]info= 0: 正常終了.
= -1: 入力パラメータ jobz の誤り. (jobz != 'V' および 'N')
= -2: 入力パラメータ uplo の誤り. (uplo != 'U 'および 'L')
= -3: 入力パラメータ n の誤り. (n < 0)
= -4: 入力パラメータ lda の誤り. (lda < max(1, n))
= -8: 入力パラメータ lwork の誤り. (lwork が小さすぎる)
= i > 0: 収束しなかった. 中間結果の3重対角形の非対角要素のうち i 個が 0 に収束しなかった.
出典
LAPACK