Mercurial > octave-dspies
diff liboctave/dbleQR.cc @ 1921:ce0db0ca0729
[project @ 1996-02-11 02:31:58 by jwe]
author | jwe |
---|---|
date | Sun, 11 Feb 1996 02:31:58 +0000 |
parents | 1281a23a34dd |
children | e379c781ae7e |
line wrap: on
line diff
--- a/liboctave/dbleQR.cc Sun Feb 11 02:23:22 1996 +0000 +++ b/liboctave/dbleQR.cc Sun Feb 11 02:31:58 1996 +0000 @@ -47,6 +47,10 @@ QR::QR (const Matrix& a, QR::type qr_type) { + tau = 0; + work = 0; + tmp_data = 0; + int m = a.rows (); int n = a.cols (); @@ -56,11 +60,12 @@ return; } - double *tmp_data; int min_mn = m < n ? m : n; - double *tau = new double[min_mn]; + tau = new double[min_mn]; + int lwork = 32*n; - double *work = new double[lwork]; + work = new double[lwork]; + int info = 0; if (m > n) @@ -71,52 +76,72 @@ else tmp_data = dup (a.data (), a.length ()); - F77_FCN (dgeqrf, DGEQRF) (m, n, tmp_data, m, tau, work, lwork, info); + F77_XFCN (dgeqrf, DGEQRF, (m, n, tmp_data, m, tau, work, lwork, info)); delete [] work; + work = 0; - if (qr_type == QR::raw) - { - for (int j = 0; j < min_mn; j++) - { - int limit = j < min_mn - 1 ? j : min_mn - 1; - for (int i = limit + 1; i < m; i++) - tmp_data[m*j+i] *= tau[j]; - } - } + if (f77_exception_encountered) + (*current_liboctave_error_handler) ("unrecoverable error in dgeqrf"); else { - int n2; - if (qr_type == QR::economy && m > n) + if (qr_type == QR::raw) { - n2 = n; - r.resize (n, n, 0.0); + for (int j = 0; j < min_mn; j++) + { + int limit = j < min_mn - 1 ? j : min_mn - 1; + for (int i = limit + 1; i < m; i++) + tmp_data[m*j+i] *= tau[j]; + } } else { - n2 = m; - r.resize (m, n, 0.0); - } + volatile int n2; - for (int j = 0; j < n; j++) - { - int limit = j < min_mn-1 ? j : min_mn-1; - for (int i = 0; i <= limit; i++) - r.elem (i, j) = tmp_data[m*j+i]; - } + if (qr_type == QR::economy && m > n) + { + n2 = n; + r.resize (n, n, 0.0); + } + else + { + n2 = m; + r.resize (m, n, 0.0); + } + + for (int j = 0; j < n; j++) + { + int limit = j < min_mn-1 ? j : min_mn-1; + for (int i = 0; i <= limit; i++) + r.elem (i, j) = tmp_data[m*j+i]; + } + + lwork = 32*m; + work = new double[lwork]; - lwork = 32*m; - work = new double[lwork]; + F77_XFCN (dorgqr, 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); + delete [] work; + work = 0; - q = Matrix (tmp_data, m, m); - q.resize (m, n2); + if (f77_exception_encountered) + (*current_liboctave_error_handler) + ("unrecoverable error in dorgqr"); + else + { + q = Matrix (tmp_data, m, m); + tmp_data = 0; + q.resize (m, n2); + } + } + } - delete [] tau; - delete [] work; - } + delete [] tau; + tau = 0; + + delete [] tmp_data; + tmp_data = 0; } /*