8 void ldvmat(
double *a,
double *v,
int n)
10 double *p0, *q0, *p, *q, *qq;
16 for (i = 0, mm = n * n, q = v; i < mm; ++i)
22 p0 = a + n * n - n - n - 1;
23 for (i = n - 2, mm = 1; i > 0; --i, p0 -= n + 1, q0 -= n + 1, ++mm) {
24 if (*(p0 - 1) != 0.) {
25 for (j = 0, p = p0, h = 1.; j < mm; ++j, ++p)
29 for (j = 0, q = q0 + n, p = p0; j < mm; ++j, q += n)
31 for (k = i + 1, q = q0 + 1; k < n; ++k) {
32 for (j = 0, qq = q + n, p = p0, s = 0.; j < mm; ++j, qq += n)
35 for (j = 0, qq = q + n, p = p0; j < mm; ++j, qq += n)
42 for (j = 0, p = q0 + 1, q = q0 + n; j < mm; ++j, q += n)
void ldvmat(double *a, double *v, int n)