10 int sv2u1v(
double *d,
double *a,
int m,
double *v,
int n)
12 double *p, *p1, *q, *pp, *w, *e;
14 double s,
t, h,
r, sv;
16 int i, j, k, mm, nm, ms;
20 w = (
double *)calloc(m + n,
sizeof(
double));
22 for (i = 0, mm = m, p = a; i < n; ++i, --mm, p += n + 1) {
25 for (j = 0, q = p, s = 0.; j < mm; ++j, q += n) {
36 sv = 1. + fabs(*p / h);
37 for (k = 1, ms = n - i; k < ms; ++k) {
38 for (j = 0, q = p + k,
r = 0.; j < mm; q += n)
41 for (j = 0, q = p + k; j < mm; q += n)
44 for (j = 1, q = p; j < mm;)
45 *(q += n) = w[j++] *
t;
53 for (i = 0, q = v, p = a; i < n; ++i) {
54 for (j = 0; j < n; ++j, ++q, ++p) {
64 for (i = 0, mm = n, nm = n - 1, p = v; i < n; ++i, --mm, --nm, p += n + 1) {
67 for (j = 0, q = p, s = 0.; j < mm; ++j, q += n) {
78 sv = 1. + fabs(*p / h);
79 for (k = 1, ms = n - i; k < ms; ++k) {
80 for (j = 0, q = p + k,
r = 0.; j < mm; q += n)
82 for (j = 0, q = p + k,
r *= s; j < mm; q += n)
85 for (k = 0, p1 = a + i; k < m; ++k, p1 += n) {
86 for (j = 0, q = p1,
r = 0.; j < mm;)
88 for (j = 0, q = p1,
r *= s; j < mm;)
100 for (j = 0, q = p1, s = 0.; j < nm; ++j, ++q)
106 sv = 1. + fabs(*p1 / h);
110 for (k = n, ms = n * (n - i); k < ms; k += n) {
111 for (j = 0, q = p1, pp = p1 + k,
r = 0.; j < nm; ++j)
113 for (j = 0, q = p1, pp = p1 + k,
r *= s; j < nm; ++j)
116 for (j = 1, q = p1 + 1; j < nm; ++j)
127 for (i = 0; i < n; ++i) {
130 for (j = 0, p = v + i; j < n; ++j, p += n)
void atou1(double *a, int m, int n)
void atovm(double *v, int n)
int qrbdu1(double *d, double *e, double *u, int m, double *v, int n)
int sv2u1v(double *d, double *a, int m, double *v, int n)