/* this is to be tested, HHC, Mar 9, 1999 */
int main() {
int m = 4, n = 3;
array double a[4][3] = {-1,5, 6,
3,-6, 1,
1, 3, 4,
6, 8, 9}; /* mxn matrix, m > n */
int mn = min(m, n);
array double l[m][mn], u[mn][n]; /* ??? */
int status;
printf("warning: ludecomp_2 is to be tested\n");
exit(0); /* ??? */
status = ludecomp(a, l, u);
if(status == 0) {
printf("l =\n%f\n", l);
printf("u =\n%f\n", u);
printf("l*u =\n%f\n", l*u);
}
else
printf("error: numerical error in ludecomp()\n");
}
warning: ludecomp_2 is to be tested