|
|
◆ dpttrs()
| void dpttrs |
( |
int |
n, |
|
|
int |
nrhs, |
|
|
double |
d[], |
|
|
double |
e[], |
|
|
int |
ldb, |
|
|
double |
b[], |
|
|
int * |
info |
|
) |
| |
分解済の連立一次方程式 AX = B の解 (正定値対称3重対角行列)
- 目的
- 本ルーチンは正定値対称3重対角形の連立一次方程式 をdpttrfにより計算されたAのL*D*L^T分解を用いて解く. Dはd[]に格納された対角行列, Lは対角要素が1で副対角要素がe[]に格納された2重対角行列, また, XおよびBはn×nrhs行列である.
- 引数
-
| [in] | n | 行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) |
| [in] | nrhs | 右辺の数, すなわち, 行列Bの列数. (nrhs >= 0) (nrhs = 0 の場合, 処理を行わずに戻る) |
| [in] | d[] | 配列 d[ld] (ld >= n)
AのL*D*L^T分解結果の対角行列Dのn個の対角要素. |
| [in] | e[] | 配列 e[le] (le >= n - 1)
AのL*D*L^T分解の2重対角行列Lのn-1個の下副対角要素. e[]はまた, 分解 A = U^T*D*U の2重対角行列Uの上副対角要素とみなすこともできる. |
| [in] | ldb | 二次元配列b[][]の整合寸法. (ldb >= max(1, n)) |
| [in,out] | b[][] | 配列 b[lb][ldb] (lb >= nrhs)
[in] n×nrhs右辺行列 B.
[out] info = 0 の場合, n×nrhs解行列 X. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ n の誤り (n < 0)
= -2: 入力パラメータ nrhs の誤り (nrhs < 0)
= -5: 入力パラメータ ldb の誤り (ldb < max(1, n)) |
- 出典
- LAPACK
|