|
|
◆ _dsbevd()
| void _dsbevd |
( |
char |
jobz, |
|
|
char |
uplo, |
|
|
int |
n, |
|
|
int |
kd, |
|
|
int |
ldab, |
|
|
double |
ab[], |
|
|
double |
w[], |
|
|
int |
ldz, |
|
|
double |
z[], |
|
|
double |
work[], |
|
|
int |
lwork, |
|
|
int |
iwork[], |
|
|
int |
liwork, |
|
|
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] = Aij. ただし, max(0, j - kd - 1) <= i <= j <= n - 1.
uplo = 'L': ab[j][i - j] = Aij. ただし, 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]
作業領域.
info = 0 の場合, work[0] に lwork の最適値を返す. |
| [in] | lwork | 配列 work[]のサイズ (lwork >= 1 (n <= 1 の場合), 2*n (jobz = 'N' の場合), 2*n^2 + 5*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: 入力パラメータ uplo の誤り (uplo != 'U' および 'L')
= -3: 入力パラメータ n の誤り (n < 0)
= -4: 入力パラメータ kd の誤り (kd < 0)
= -5: 入力パラメータ ldab の誤り (ldab < kd + 1)
= -8: 入力パラメータ ldz の誤り (ldz が小さすぎる)
= -11: 入力パラメータ lwork の誤り. (lwork が小さすぎる)
= -13: 入力パラメータ liwork の誤り. (liwork が小さすぎる)
= i > 0: 収束しなかった. 中間結果の3重対角形の非対角要素のうち i 個が 0 にならなかった. |
- 出典
- LAPACK
|
1.9.6