|
|
◆ _dsygv()
| void _dsygv |
( |
int |
itype, |
|
|
char |
jobz, |
|
|
char |
uplo, |
|
|
int |
n, |
|
|
int |
lda, |
|
|
double |
a[], |
|
|
int |
ldb, |
|
|
double |
b[], |
|
|
double |
w[], |
|
|
double |
work[], |
|
|
int |
lwork, |
|
|
int * |
info |
|
) |
| |
(シンプルドライバ) 一般化固有値問題 (対称行列)
- 目的
- 本ルーチンは実対称行列の一般化固有値問題
Ax = λBx, ABx = λx または BAx = λx
- 引数
-
| [in] | itype | 解くべき問題のタイプを指定.
= 1: Ax = λBx.
= 2: ABx = λx.
= 3: BAx = λx. |
| [in] | jobz | = 'N': 固有値のみ求める.
= 'V': 固有値と固有ベクトルを求める. |
| [in] | uplo | = 'U': A および B の上三角部分を格納する.
= 'L': A および B の下三角部分を格納する. |
| [in] | n | 行列 A および B の行および列数. (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 の場合, 次のように正規化された固有ベクトルからなる行列 Z を a[][] に返す.
itype = 1 または 2: Z^T*B*Z = I
itype = 3: Z^T*inv(B)*Z = I
jobz = 'N': a[][] の上三角部分(uplo = 'U' の場合)あるいは下三角部分(uplo = 'L' の場合)は, 対角部分を含め壊される. |
| [in] | ldb | 二次元配列 b[][] の整合寸法. (ldb >= max(1, n)) |
| [in,out] | b[][] | 配列 b[lb][ldb] (lb >= n)
[in] 正定値対称行列 B. uplo = 'U' の場合, b[][] の n x n 上三角部分に行列 B の上三角部分を格納する. uplo = 'L' の場合, b[][] の n x n 下三角部分に行列 B の下三角部分を格納する.
[out] info <= n の場合, b[][] の行列を格納している部分はコレスキー分解 B = U^T*U あるいは B = L*L^T の三角行列 U あるいは L により上書きされる. |
| [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[0] にその値を返す. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ itype の誤り (itype < 1 または itype > 3)
= -2: 入力パラメータ jobz の誤り (jobz != 'V' および 'N')
= -3: 入力パラメータ uplo の誤り (uplo != 'U' および 'L')
= -4: 入力パラメータ n の誤り (n < 0)
= -5: 入力パラメータ lda の誤り (lda < max(1, n))
= -7: 入力パラメータ ldb の誤り (ldb < max(1, n))
= -11: 入力パラメータ lwork の誤り (lworkが小さすぎる)
= i (0 < i <= n): dsyev が収束しなかった. 中間の3重対角形の非対角要素のうち i 個が 0 に収束しなかった.
= i (i > n): B の i-n 次小行列が正定値でない. B の分解が完了できず, 固有値・固有ベクトルは計算されなかった. |
- 出典
- LAPACK
|
1.9.6