◆ _dgecon()

行列の条件数 (一般行列)

目的: 本ルーチンは一般行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する. 計算にはdgetrf, dgesv または dgesvxにより求めた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]	lda	二次元配列a[][]の整合寸法. (lda >= max(1, n))
[in]	a[][]	配列 a[la][lda] (la >= n) dgetrf, dgesv または dgesvxにより計算された分解 A = PLU のLおよびU.
[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 >= 4*n) 作業領域.
[out]	iwork[]	配列 iwork[liwork] (liwork >= n) 作業領域.
[out]	info	= 0: 正常終了 = -1: 入力パラメータ norm の誤り (norm != '1', 'O'および'I') = -2: 入力パラメータ n の誤り (n < 0) = -3: 入力パラメータ lda の誤り (lda < max(1, n)) = -5: 入力パラメータ anorm の誤り (anorm < 0)