WZgbsv WZgbsv2 WZgesv WZgesv2 WZgtsv WZgtsv2 WZsysv WZsysv2 WZtrtrs WZtrtrs2

WZgtsv2() Function WZgtsv2 ( N As  Long, Dl As  Variant, D As  Variant, Du As  Variant, B As  Variant, Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B の解 (複素3重対角行列) (実数部／虚数部分離形式) 目的 WZgtsv2は次の連立一次方程式を解く. A * X = B ここで, AはN×N 3重対角行列である. ピボットの部分選択を行うガウスの消去法を使用する. 引数DuとDlの順番を入れ替えることにより, 方程式 A^T*X = B を解くこともできる. 複素数を表現するために実数部と虚数部を隣り合ったセルに格納する(左が実数部, 右が虚数部). 得られた解も実数部と虚数部が隣り合った別々のセルに出力される. 戻り値 N+1 × 2Nrhs 列1 列2 ・・・ 列2Nrhs 行1〜N 解行列 X (実数部と虚数部を隣り合った列に格納(左が実数部, 右が虚数部)) 行N+1 1/条件数 リターンコード ・・・ 0 リターンコード. = 0: 正常終了. = i > 0: 行列のi番目のピボットがゼロになった. (行列 A は特異) 引数 [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] Dl (N-1 x 2) 方程式の係数行列 A の下副対角要素. [in] D (N x 2) 方程式の係数行列 A の対角要素. [in] Du (N-1 x 2) 方程式の係数行列 A の上副対角要素. [in] B (N × 2Nrhs) N×Nrhs右辺行列 B. [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, Aは3重対角行列で ( 0.57-0.91i 0.56+0.92i 0 ) A = ( -1.45-0.81i 1.74-0.93i -0.19-0.16i ) ( 0 0.28+0.09i 0.10+0.15i ) ( -1.0924+1.0032i ) B = ( 1.8436+1.5703i ) ( 0.2057+0.2156i ) とする.