◆ zgbcon()

行列の条件数 (複素帯行列)

目的: 本ルーチンは複素帯行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する. 計算にはzgbtrfにより求められたLU分解が使用される.

norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する.
rcond = 1 / (norm(A) * norm(inv(A)))

引数

[in]	norm	1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定 = '1'または'O': 1ノルム. = 'I': 無限ノルム.
[in]	n	行列Aの行および列数. (n >= 0) (n = 0 の場合, rcond = 1 を返す)
[in]	kl	行列Aの下帯幅. (kl >= 0)
[in]	ku	行列Aの上帯幅. (ku >= 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]	anorm	norm = '1'または'O': 分解前の行列Aの1ノルム. (anorm >= 0) norm = 'I': 分解前の行列Aの無限ノルム. (anorm >= 0)
[out]	rcond	行列Aの条件数の逆数. (rcond = 1/(norm(A) * norm(inv(A))))
[out]	work[]	配列 work[lwork] (lwork >= 2*n) 作業領域.
[out]	rwork[]	配列 rwork[lrwork] (lrwork >= n) 作業領域.
[out]	info	= 0: 正常終了 = -1: 入力パラメータ norm の誤り (norm != '1', 'O'および'I') = -2: 入力パラメータ n の誤り (n < 0) = -3: 入力パラメータ kl の誤り (kl < 0) = -4: 入力パラメータ ku の誤り (ku < 0) = -5: 入力パラメータ ldab の誤り (ldab < 2kl+ku+1) = -8: 入力パラメータ anorm の誤り (anorm < 0)