/* 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