Mercurial > hg > octave-lyh
diff liboctave/dbleQRP.cc @ 1253:bb67a902760b
[project @ 1995-04-11 16:35:23 by jwe]
| author | jwe |
|---|---|
| date | Tue, 11 Apr 1995 16:35:23 +0000 |
| parents | 97eac19837dc |
| children | f93b7fa5e113 |
line wrap: on
line diff
--- a/liboctave/dbleQRP.cc +++ b/liboctave/dbleQRP.cc @@ -34,11 +34,13 @@ extern "C" { - int F77_FCN (dgeqpf) (const int&, const int&, double*, const int&, - int*, double*, double*, int&); + int F77_FCN (dgeqpf, DGEQPF) (const int&, const int&, double*, + const int&, int*, double*, double*, + int&); - int F77_FCN (dorgqr) (const int&, const int&, const int&, double*, - const int&, double*, double*, const int&, int&); + int F77_FCN (dorgqr, DORGQR) (const int&, const int&, const int&, + double*, const int&, double*, double*, + const int&, int&); } // It would be best to share some of this code with QR class... @@ -79,7 +81,7 @@ for (int i = 0; i < n; i++) jpvt[i] = 0; - F77_FCN (dgeqpf) (m, n, tmp_data, m, jpvt, tau, work, info); + F77_FCN (dgeqpf, DGEQPF) (m, n, tmp_data, m, jpvt, tau, work, info); // Form Permutation matrix (if economy is requested, return the // indices only!) @@ -122,7 +124,8 @@ lwork = 32*m; work = new double[lwork]; - F77_FCN (dorgqr) (m, m, min_mn, tmp_data, m, tau, work, lwork, info); + F77_FCN (dorgqr, DORGQR) (m, m, min_mn, tmp_data, m, tau, work, + lwork, info); q = Matrix (tmp_data, m, m); q.resize (m, n2);