Banfac Banslv Bint4 Bintk Bspldr Bsplev Bsplpp Bsplvd Bsplvn Bvalue Interv Ppvalu

Banfac() Sub Banfac ( N As  Long, Kl As  Long, Ku As  Long, Ab() As  Double, Info As  Long  ) 連立一次方程式の係数行列(帯行列)のLU分解 (Bint4およびBintkの補助ルーチン) 目的 本ルーチンはN次の(Kl+1+Ku)帯または対角行列AのLU分解を求める(ピボットの選択なし). LU分解された行列は連立一次方程式を解いてB-スプライン係数を求めるためにBanslvにより使用される. BanfacおよびBanslvはBint4およびBintkの補助ルーチンである. 引数 [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] Kl Aの下帯幅. (Kl >= 0) [in] Ku Aの上帯幅. (Ku >= 0) [in,out] Ab() 配列 Ab(LAb1 - 1, LAb2 - 1) (LAb1 >= Kl + Ku + 1, LAb2 >= N) (2D 配列) または Ab(LAb - 1) (LAb >= (Kl + Ku + 1)*N) (1D 配列) [in] 帯行列形式の行列A. [out] Aの対角要素が1の下三角行列Lおよび上三角行列U(両方とも帯行列)への分解結果. Aの各要素に上書き保存される. [out] Info = 0: 正常終了. = -1: パラメータ N の誤り. (N < 1) = -2: パラメータ Kl の誤り. (K < 1) = -3: パラメータ Ku の誤り. (Ku < 1) = -4: パラメータ Ab() の誤り. = 1: 分解に失敗した. これは, Aが特異であることを示している. 出典 SLATEC