ludec2
n×n 行列を LU 分解するための倍精度サブルーチン．部分配列の使用を前提とすることにより，サブルーチンludecよりも引数を減らしてある．
サブルーチン名： ludec2
呼び出し方：call ludec2(adim,n,ip,epsln)
引 数：
引数
意 味
入出力の区分
備 考
adim(n,*)
LU分解される係数行列 [a
ij
] を格納しておく配列の配列名
入力＆出力
出力時には LU 分解された結果が格納される．
n
行列 [a
ij
] の寸法
入力
2≦nを満足するように選ぶ．
epsln
特異性判定子
入力
1.0×10
-14
程度．マシン・イプシロンでも良い．
ip(*)
インデックス・ベクトルを格納する１次元配列
出力
配列寸法は n 以上