◆ zgbtrs()

LU分解済の連立一次方程式 AX = B, A^TX = B または A^HX = B の解 (複素帯行列)

目的: 本ルーチンは複素帯行列Aの連立一次方程式
A * X = B, A^T * X = B または A^H * X = B

をzgbtrfにより計算された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]	kl	Aの下帯幅. (kl >= 0)
[in]	ku	Aの上帯幅. (ku >= 0)
[in]	nrhs	右辺の数, すなわち, 行列Bの列数. (nrhs >= 0) (nrhs = 0 の場合, 処理を行わずに戻る)
[in]	ldab	二次元配列ab[][]の整合寸法. (ldab >= 2kl + ku + 1)
[in]	ab[][]	配列 ab[lab][ldab] (lab >= n) zgbtrfで計算された帯行列AのLU分解結果. Uは上帯幅kl+kuの上三角帯行列として第1〜kl+ku+1行に格納される. また, 分解中に使われた乗数が第kl+ku+2〜2kl+ku+1行に格納される.
[in]	ipiv[]	配列 ipiv[lipiv] (lipiv >= n) ピボットインデックス. 1 <= i <= nに対して, 第i行が第ipiv[i-1]行と交換されたことを表す.
[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: 入力パラメータ kl の誤り (kl < 0) = -4: 入力パラメータ ku の誤り (ku < 0) = -5: 入力パラメータ nrhs の誤り (nrhs < 0) = -6: 入力パラメータ ldab の誤り (ldab < 2kl+ku+1) = -9: 入力パラメータ ldb の誤り (ldb < max(1, n))