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