WDgbsv WDgesv WDgtsv WDtrtrs

WDtrtrs() Function WDtrtrs ( Uplo As  String, N As  Long, A As  Variant, B As  Variant, Optional Trans As  String = "N", Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B または ATX = B の解 (三角行列) 目的 WDtrtrsは三角行列の連立一次方程式 A * X = B または A^T * X = B を解く. ここで, AはN×N三角行列, また, BはN×Nrhs行列である. Aが特異でないことを確認するためのチェックを行う. 戻り値 N+2 × Nrhs 列1 列2 ・・・ 列Nrhs 行1〜N 解行列 X 行N+1 1/条件数 0 ・・・ 0 行N+2 リターンコード 0 ・・・ 0 リターンコード. = 0: 正常終了. = i > 0: 行列のi番目のピボットがゼロになった. (行列 A は特異) 引数 [in] Uplo = "U": Aは上三角行列. = "L": Aは下三角行列. [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] A (N × N) N×N係数行列 A. (Uplo = 0のとき下三角要素のみ使用, Uplo = 1のとき上三角要素のみ使用) [in] B (N × Nrhs) N×Nrhs右辺行列 B. [in] Trans (省略可) 連立方程式の形を指定. (省略時 = "N") = "N": A * X = B. (転置なし) = "T": A^T * X = B. (転置あり) [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, Aは三角行列で ( -1.13 0 0 ) ( 0.0452 ) A = ( 0.26 -1.98 0 ), B = ( -0.4856 ) ( -0.96 0.30 -2.32 ) ( 1.2472 ) とする.