9 void ldumat(
double *a,
double *u,
int m,
int n)
11 double *p0, *q0, *p, *q, *w;
17 w = (
double *)calloc(m,
sizeof(
double));
18 for (i = 0, mm = m * m, q = u; i < mm; ++i)
24 for (j = 0; j < mm; ++j, q0 -= m + 1)
33 for (; i >= 0; --i, ++mm, p0 -= n + 1, q0 -= m + 1) {
35 for (j = 0, p = p0 + n, h = 1.; j < mm; p += n)
39 for (j = 0, q = q0 + m; j < mm; q += m)
41 for (k = i + 1, q = q0 + 1; k < m; ++k) {
42 for (j = 0, p = q + m, s = 0.; j < mm; p += m)
45 for (j = 0, p = q + m; j < mm; p += m)
52 for (j = 0, p = q0 + 1, q = q0 + m; j < mm; ++j, q += m)
void ldumat(double *a, double *u, int m, int n)