Zpbcon Zpbtrf Zpbtrs Zptcon Zpttrf Zpttrs

Zpttrs() Sub Zpttrs ( Uplo As  String, N As  Long, D() As  Double, E() As  Complex, B() As  Complex, Info As  Long, Optional Nrhs As  Long = 1  ) 分解済の連立一次方程式 AX = B の解 (正定値エルミート3重対角行列) 目的 本ルーチンは3重対角形の連立一次方程式 A * X = B をZpttrfにより計算された分解 A = U^H*D*U または A = L*D*L^H を用いて解く. Dは配列D()に格納された対角行列, U(またはL)は対角要素が1で副対角要素が配列E()に格納された2重対角行列, また, XおよびBはn×nrhs行列である. 引数 [in] Uplo 分解の形式, および, E()が上2重対角行列Uの上副対角要素であるか下2重対角行列Lの下副対角要素であるかを指定する. = "U": A = U^H*D*U, E()はUの上副対角要素. = "L": A = L*D*L^H, E()はLの下副対角要素. [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in] D() 配列 D(LD - 1) (LD >= N) 分解 A = U^H*D*U または A = L*D*L^H の対角行列DのN個の対角要素. [in] E() 配列 E(LE - 1) (LE >= N - 1) 分解 A = U^H*D*U または A = L*D*L^H の2重対角行列UまたはLのN-1個の上または下副対角要素. [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: パラメータ D() の誤り. = -4: パラメータ E() の誤り. = -5: パラメータ B() の誤り. = -7: パラメータ Nrhs の誤り. (Nrhs < 0) [in] Nrhs (省略可) 右辺の数, すなわち, 行列Bの列数. (Nrhs >= 0) (Nrhs = 0 の場合, 処理を行わずに戻る) (省略時 = 1) 出典 LAPACK 使用例 Zpttrfの使用例を参照せよ.