Zpbcon Zpbtrf Zpbtrs Zptcon Zpttrf Zpttrs

Zpbcon() Sub Zpbcon ( Uplo As  String, N As  Long, Kd As  Long, Ab() As  Complex, ANorm As  Double, RCond As  Double, Info As  Long  ) 行列の条件数 (正定値エルミート帯行列) 目的 本ルーチンはZpbtrfにより計算されたコレスキー分解 A = U^H*U または A = L*L^H を用いて正定値エルミート帯行列の(1-ノルムによる)条件数の推定を行う. norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する. RCond = 1 / (norm(A) * norm(inv(A))) 引数 [in] Uplo = 'U': 上三角行列UをAb()に格納. = 'L': 下三角行列LをAb()に格納. [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, RCond = 1 を返す) [in] Kd 行列Aの上帯幅(Uplo = "U"の場合)あるいは下帯幅(Uplo = "L"の場合). (Kd >= 0) [in] Ab() 配列 Ab(LAb1 - 1, LAb2 - 1) (LAb1 >= Kd + 1, LAb2 >= N) 正定値エルミート帯行列Aのコレスキー分解 A = U^H*U または A = L*L^H の三角行列 U または L (Kd+1×N対称帯行列形式). [in] ANorm 正定値エルミート帯行列Aの1ノルム(または無限大ノルム). [out] RCond 行列Aの条件数の逆数. 次のように計算する.   RCond = 1/(ANorm * Ainvnm) ここで, Ainvnmは本ルーチン内で求められたinv(A)の1-ノルムの推定値である. [out] Info = 0: 正常終了. = -1: パラメータ Uplo の誤り. (Uplo <> "U"および"L") = -2: パラメータ N の誤り. (N < 0) = -2: パラメータ Kd の誤り. (Kd < 0) = -4: パラメータ Ab() の誤り. = -5: パラメータ ANorm の誤り. (ANorm < 0) 出典 LAPACK 使用例 Zpbsvの使用例を参照せよ.