◆ dgecon()

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

目的: dgeconは一般行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する. 計算にはdgetrf, dgesv または dgesvxにより求められたLU分解が使用される.

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

戻り値: (rcond, info)

rcond (float):
行列Aの条件数の逆数. (rcond = 1/(norm(A) * norm(inv(A))))

info (int):
= 0: 正常終了.
= -1: 入力パラメータ norm の誤り. (norm != '1', 'O'および'I')
= -2: 入力パラメータ n の誤り. (n < 0)
= -3: 入力パラメータ a の誤り.
= -4: 入力パラメータ anorm の誤り. (anorm < 0)

引数

[in]	norm	1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定. = '1'または'O': 1ノルム = 'I': 無限ノルム
[in]	n	行列Aの行および列数. (n >= 0) (n = 0 の場合, rcond = 1 を返す)
[in]	a	Numpy ndarray (2次元配列, float, n x n) dgetrf, dgesv または dgesvxにより計算された分解 A = PLU のLおよびU.
[in]	anorm	norm = '1'または'O': 分解前の行列Aの1ノルム (anorm >= 0) norm = 'I': 分解前の行列Aの無限ノルム (anorm >= 0)