Mercurial > hg > octave-nkf
diff liboctave/dbleHESS.cc @ 7482:29980c6b8604
don't check f77_exception_encountered
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 14 Feb 2008 21:57:50 -0500 |
parents | a1dbe9d80eee |
children | eb63fbe60fab |
line wrap: on
line diff
--- a/liboctave/dbleHESS.cc +++ b/liboctave/dbleHESS.cc @@ -87,58 +87,36 @@ n, h, n, ilo, ihi, pscale, info F77_CHAR_ARG_LEN (1))); - if (f77_exception_encountered) - (*current_liboctave_error_handler) ("unrecoverable error in dgebal"); - else - { - Array<double> tau (n-1); - double *ptau = tau.fortran_vec (); + Array<double> tau (n-1); + double *ptau = tau.fortran_vec (); - Array<double> work (lwork); - double *pwork = work.fortran_vec (); + Array<double> work (lwork); + double *pwork = work.fortran_vec (); - F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork, - lwork, info)); + F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork, + lwork, info)); - if (f77_exception_encountered) - (*current_liboctave_error_handler) ("unrecoverable error in dgehrd"); - else - { - unitary_hess_mat = hess_mat; - double *z = unitary_hess_mat.fortran_vec (); + unitary_hess_mat = hess_mat; + double *z = unitary_hess_mat.fortran_vec (); - F77_XFCN (dorghr, DORGHR, (n, ilo, ihi, z, n, ptau, pwork, - lwork, info)); + F77_XFCN (dorghr, DORGHR, (n, ilo, ihi, z, n, ptau, pwork, + lwork, info)); - if (f77_exception_encountered) - (*current_liboctave_error_handler) - ("unrecoverable error in dorghr"); - else - { - F77_XFCN (dgebak, DGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1), - F77_CONST_CHAR_ARG2 (&side, 1), - n, ilo, ihi, pscale, n, z, - n, info - F77_CHAR_ARG_LEN (1) - F77_CHAR_ARG_LEN (1))); + F77_XFCN (dgebak, DGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1), + F77_CONST_CHAR_ARG2 (&side, 1), + n, ilo, ihi, pscale, n, z, + n, info + F77_CHAR_ARG_LEN (1) + F77_CHAR_ARG_LEN (1))); - if (f77_exception_encountered) - (*current_liboctave_error_handler) - ("unrecoverable error in dgebak"); - else - { - // If someone thinks of a more graceful way of doing - // this (or faster for that matter :-)), please let - // me know! + // If someone thinks of a more graceful way of doing + // this (or faster for that matter :-)), please let + // me know! - if (n > 2) - for (octave_idx_type j = 0; j < a_nc; j++) - for (octave_idx_type i = j+2; i < a_nr; i++) - hess_mat.elem (i, j) = 0; - } - } - } - } + if (n > 2) + for (octave_idx_type j = 0; j < a_nc; j++) + for (octave_idx_type i = j+2; i < a_nr; i++) + hess_mat.elem (i, j) = 0; return info; }