Zpbcon Zpbtrf Zpbtrs Zptcon Zpttrf Zpttrs

Zpbtrs() Sub Zpbtrs ( Uplo As  String, N As  Long, Kd As  Long, Ab() As  Complex, B() As  Complex, Info As  Long, Optional Nrhs As  Long = 1  ) 分解済の連立一次方程式 AX = B の解 (正定値エルミート帯行列) 目的 本ルーチンは正定値エルミート帯行列Aの連立一次方程式 A * X = B を, Zpbtrfにより求められたコレスキー分解 A = U^H*U あるいは A = L*L^H を用いて解く. 引数 [in] Uplo = "U": Aの上三角部分を格納. = "L": Aの下三角部分を格納. [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [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,out] B() 配列 B(LB1 - 1, LB2 - 1) (LB1 >= max(1, N), LB2 >= Nrhs) (2次元配列) または B(LB - 1) (LB >= max(1, N), Nrhs = 1) (1次元配列) [in] N×Nrhs右辺行列 B. [out] Info = 0 の場合, N×Nrhs解行列 X. [out] Info = 0: 正常終了. = -1: パラメータ Uplo の誤り. (Uplo <> "U"および"L") = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ Kd の誤り. (Kd < 0) = -4: パラメータ Ab() の誤り. = -5: パラメータ B() の誤り. = -7: パラメータ Nrhs の誤り. (nrhs < 0) [in] Nrhs (省略可) 右辺の数, すなわち, 行列Bの列数. (Nrhs >= 0) (Nrhs = 0 の場合, 処理を行わずに戻る) (省略時 = 1) 出典 LAPACK 使用例 Zpbtrfの使用例を参照せよ.