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

◆ zgttrs()

void zgttrs ( char  trans,
int  n,
int  nrhs,
doublecomplex  dl[],
doublecomplex  d[],
doublecomplex  du[],
doublecomplex  du2[],
int  ipiv[],
int  ldb,
doublecomplex  b[],
int *  info 
)

LU分解済の連立一次方程式 AX = B, ATX = B または AHX = B の解 (複素3重対角行列)

目的
本ルーチンは3重対角行列Aの連立一次方程式
A * X = B, A^T * X = B または A^H * X = B
をzgttrfにより求められたLU分解を用いて解く.
引数
[in]trans連立方程式の形を指定
= 'N': A * X = B. (転置なし)
= 'T': A^T * X = B. (転置あり)
= 'C': A^H * X = B. (共役転置あり)
[in]n行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]nrhs右辺の数, すなわち, 行列Bの列数. (nrhs >= 0) (nrhs = 0 の場合, 処理を行わずに戻る)
[in]dl()配列 dl[ldl] (ldl >= n - 1)
AのLU分解結果の行列Lを定義するn-1個の乗数.
[in]d[]配列 d[ld] (ld >= n)
AのLU分解結果の上三角行列Uの対角要素.
[in]du[]配列 du[ldu] (ldu >= n - 1)
Uのn-1個の第1上副対角要素.
[in]du2[]配列 du2[ldu2] (ldu2 >= n - 2)
Uのn-2個の第2上副対角要素.
[out]ipiv[]配列 ipiv[lipiv] (lipiv >= n)
ピボットインデックス. 1 <= i <= n について, 行列の第i行は第ipiv[i-1]行と交換されたことを表す. ipiv[i-1]は常にiまたはi+1である. ipiv[i-1] = i は行の交換が不要であったことを示す.
[in]ldb二次元配列b[][]の整合寸法. (ldb >= max(1, n))
[in,out]b[][]配列 b[lb][ldb] (lb >= nrhs)
[in] 右辺行列 B.
[out] 解行列 X.
[out]info= 0: 正常終了
= -1: 入力パラメータ trans の誤り (trans != 'N', 'T'および'C')
= -2: 入力パラメータ n の誤り (n < 0)
= -3: 入力パラメータ nrhs の誤り (nrhs < 0)
= -9: 入力パラメータ ldb の誤り (ldb < max(1, n))
出典
LAPACK