Ztbcon Ztbtrs Ztpcon Ztptri Ztptrs Ztrcon Ztrtri Ztrtrs

Ztbcon() Sub Ztbcon ( Norm As  String, Uplo As  String, Diag As  String, N As  Long, Kd As  Long, Ab() As  Complex, RCond As  Double, Info As  Long  ) 行列の条件数 (複素三角帯行列) 目的 本ルーチンは三角帯行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する. Aのノルムを計算し, norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する. RCond = 1 / (norm(A) * norm(inv(A))) 引数 [in] Norm 1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定. = "1"または"O": 1ノルム. = "I": 無限ノルム. [in] Uplo = "U": Aは上三角行列. = "L": Aは下三角行列. [in] Diag = "N": Aは単位三角行列ではない. = "U": Aは単位三角行列である. (Ab()の対角要素を参照せず, 1とみなす) [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, RCond = 1 を返す) [in] Kd 三角帯行列Aの上または下帯幅. (Kd >= 0) [in] Ab() 配列 Ab(LAb1 - 1, LAb2 - 1) (LAb1 >= Kd + 1, LAb2 >= N) Kd+1×N対称帯行列形式のN×N三角帯行列 A. Uploに従って上または下三角部分を格納. [out] RCond 行列Aの条件数の逆数. (RCond = 1/(norm(A) * norm(inv(A)))) [out] Info = 0: 正常終了. = -1: パラメータ Norm の誤り. (Norm <> "1", "O"および"I") = -2: パラメータ Uplo の誤り. (Uplo <> "U"および"L") = -3: パラメータ Diag の誤り. (Diag <> "N"および"U") = -4: パラメータ N の誤り. (N < 0) = -5: パラメータ Kd の誤り. (Kd < 0) = -6: パラメータ Ab() の誤り. 出典 LAPACK 使用例 Ztbtrsの使用例を参照せよ.