|
|
◆ dpotrf()
| void dpotrf |
( |
char |
uplo, |
|
|
int |
n, |
|
|
int |
lda, |
|
|
double |
a[], |
|
|
int * |
info |
|
) |
| |
係数行列のコレスキー分解 (正定値対称行列)
- 目的
- 本ルーチンは正定値対称行列Aのコレスキー分解を計算する. 分解は次の形式である.
A = U^T * U (uplo = 'U'の場合)
A = L * L^T (uplo = 'L'の場合)
ここで, Uは上三角行列, Lは下三角行列である.
本ルーチンはLevel 3 BLASを呼び出すブロック版のアルゴリズムを使用する.
- 引数
-
| [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] n×n正定値対称行列 A. uploに従い上三角部分あるいは下三角部分が参照される.
[out] Info = 0の場合, コレスキー分解 A = U^T*U または A = L*L^T のUまたはL. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ uplo の誤り (uplo != 'U'および'L')
= -2: 入力パラメータ n の誤り (n < 0)
= -3: 入力パラメータ lda の誤り (lda < max(1, n))
= i > 0: i×i首座小行列が正定値でないため分解を完了できなかった. |
- 出典
- LAPACK
|