Dpbcon Dpbtrf Dpbtrs Dptcon Dpttrf Dpttrs

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