WZhesv WZhesv2 WZpbsv WZpbsv2 WZposv WZposv2 WZptsv WZptsv2

WZptsv2() Function WZptsv2 ( N As  Long, D As  Variant, E As  Variant, B As  Variant, Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B の解 (正定値エルミート3重対角行列) 目的 WZptsv2は次の複素連立一次方程式を解く. A * X = B ここで, AはN×N正定値エルミート3重対角行列, また, XおよびBはN×Nrhs行列である. Aを A = L*D*L^H と分解し, 次に分解されたAを用いて連立一次方程式の解を求める. 複素数を表現するために実数部と虚数部を隣り合ったセルに格納する(左が実数部, 右が虚数部). 得られた解も実数部と虚数部が隣り合った別々のセルに出力される. 戻り値 N+1 × 2Nrhs 列1 列2 ・・・ 列2Nrhs 行1〜N 解行列 X (実数部と虚数部を隣り合った列に格納(左が実数部, 右が虚数部)) 行N+1 1/条件数 リターンコード ・・・ 0 リターンコード. = 0: 正常終了. = i > 0: i×i首座小行列が正定値でないため分解を完了できなかった. 解は計算されなかった. i = N でない限り分解は完了していない. 引数 [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] D (N) N×N係数行列 A の対角要素. [in] E (N-1 × 2) N×N係数行列 A の副対角要素. [in] B (N × 2Nrhs) N×Nrhs右辺行列 B. [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, Aは正定値エルミート3重対角行列で ( 2.88 0.29-0.44i 0 ) A = ( 0.29+0.44i 0.62 -0.01-0.02i ) ( 0 -0.01+0.02i 0.46 ) ( 1.6236-0.7300i ) B = ( 0.1581+0.1537i ) ( 0.1132-0.2290i ) とする.