XLPack 7.0
XLPack 数値計算ライブラリ (C API) リファレンスマニュアル
読み取り中…
検索中…
一致する文字列を見つけられません

◆ dpttrs()

void dpttrs ( int  n,
int  nrhs,
double  d[],
double  e[],
int  ldb,
double  b[],
int *  info 
)

分解済の連立一次方程式 AX = B の解 (正定値対称3重対角行列)

目的
本ルーチンは正定値対称3重対角形の連立一次方程式
A * X = B
を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