|
|
◆ 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
|