11 double z, *p, *q, *
r, *s, *
t;
15 for (j = 0, p = v; j < n; ++j, p += n + 1) {
16 for (q = v + j * n; q < p; ++q)
21 for (k = j + 1, q = p + n; k < n; ++k, q += n) {
22 for (r = v + j * n, s = v + k * n, z = 0.; r < p;)
29 for (j = 0, p = v; j < n; ++j, p += n + 1) {
31 for (q = v + j, t = v; q < p; t += n + 1, q += n) {
32 for (s = q, r = t, z = 0.; s < p; s += n)
37 for (j = 0, p = v; j < n; ++j, p += n + 1) {
38 for (q = v + j, t = p - j; q <= p; q += n) {
39 for (k = j, r = p, s = q, z = 0.; k < n; ++k)
void trnm(double *a, int n)
int psinv(double *v, int n)