Mercurial > hg > octave-lyh
diff liboctave/CmplxQRP.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/CmplxQRP.cc +++ b/liboctave/CmplxQRP.cc @@ -34,11 +34,13 @@ extern "C" { - int F77_FCN (zgeqpf) (const int&, const int&, Complex*, const int&, - int*, Complex*, Complex*, double*, int&); + int F77_FCN (zgeqpf, ZGEQPF) (const int&, const int&, Complex*, + const int&, int*, Complex*, Complex*, + double*, int&); - int F77_FCN (zungqr) (const int&, const int&, const int&, Complex*, - const int&, Complex*, Complex*, const int&, int&); + int F77_FCN (zungqr, ZUNGQR) (const int&, const int&, const int&, + Complex*, const int&, Complex*, + Complex*, const int&, int&); } // It would be best to share some of this code with ComplexQR class... @@ -83,7 +85,8 @@ for (int i = 0; i < n; i++) jpvt[i] = 0; - F77_FCN (zgeqpf) (m, n, tmp_data, m, jpvt, tau, work, rwork, info); + F77_FCN (zgeqpf, ZGEQPF) (m, n, tmp_data, m, jpvt, tau, work, rwork, + info); // Form Permutation matrix (if economy is requested, return the // indices only!) @@ -126,7 +129,8 @@ lwork = 32*m; work = new Complex[lwork]; - F77_FCN (zungqr) (m, m, min_mn, tmp_data, m, tau, work, lwork, info); + F77_FCN (zungqr, ZUNGQR) (m, m, min_mn, tmp_data, m, tau, work, + lwork, info); q = ComplexMatrix (tmp_data, m, m); q.resize (m, n2);