Dgecon Dgecov Dgels Dgesv Dlange Dlansy Dpocon Dposv Dsyev

Dgecon() Sub Dgecon ( Norm As  String, N As  Long, A() As  Double, ANorm As  Double, RCond As  Double, Info As  Long  ) 行列の条件数 (一般行列) 目的 本ルーチンは一般行列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] A() 配列 A(LA1 - 1, LA2 - 1) (LA1 >= N, LA2 >= N) Dgetrf, Dgesv または Dgesvxにより計算された分解 A = P*L*U の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] Info = 0: 正常終了. = -1: パラメータ Norm の誤り. (Norm <> "1", "O"および"I") = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ A() の誤り. = -4: パラメータ ANorm の誤り. (ANorm < 0) 出典 LAPACK 使用例 Dgesvの使用例を参照せよ.