Mercurial > hg > octave-avbm
changeset 3375:1cdb3ab4180a
[project @ 1999-12-09 04:58:34 by jwe]
author | jwe |
---|---|
date | Thu, 09 Dec 1999 04:58:35 +0000 |
parents | 89405d9a9b0b |
children | 6c6cec6d4252 |
files | liboctave/ChangeLog liboctave/CmplxLU.cc liboctave/dbleLU.cc |
diffstat | 3 files changed, 12 insertions(+), 10 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/ChangeLog +++ b/liboctave/ChangeLog @@ -1,3 +1,9 @@ +1999-12-08 John W. Eaton <jwe@bevo.che.wisc.edu> + + * dbleLU.cc (LU::LU): Call DGETRF directly instead of calling DGESV. + * CmplxLU.cc (ComplexLU::ComplexLU): Likewise, call ZGETRF + directly, instead of calling ZGESV. + 1999-11-18 John W. Eaton <jwe@bevo.che.wisc.edu> * data-conv.cc (init_sized_type_lookup_table): New function.
--- a/liboctave/CmplxLU.cc +++ b/liboctave/CmplxLU.cc @@ -43,9 +43,8 @@ extern "C" { - int F77_FCN (zgesv, ZGESV) (const int&, const int&, Complex*, - const int&, int*, Complex*, const int&, - int&); + int F77_FCN (zgetrf, ZGETRF) (const int&, const int&, Complex*, + const int&, int*, int&); } ComplexLU::ComplexLU (const ComplexMatrix& a) @@ -68,9 +67,8 @@ Complex *tmp_data = a_fact.fortran_vec (); int info = 0; - Complex *dummy = 0; - F77_XFCN (zgesv, ZGESV, (n, 0, tmp_data, n, pipvt, dummy, n, info)); + F77_XFCN (zgetrf, ZGETRF, (n, n, tmp_data, n, pipvt, info)); if (f77_exception_encountered) (*current_liboctave_error_handler) ("unrecoverable error in zgesv");
--- a/liboctave/dbleLU.cc +++ b/liboctave/dbleLU.cc @@ -43,9 +43,8 @@ extern "C" { - int F77_FCN (dgesv, DGESV) (const int&, const int&, double*, - const int&, int*, double&, const int&, - int&); + int F77_FCN (dgetrf, DGETRF) (const int&, const int&, double*, + const int&, int*, int&); } LU::LU (const Matrix& a) @@ -68,9 +67,8 @@ double *tmp_data = a_fact.fortran_vec (); int info = 0; - double dummy = 0; - F77_XFCN (dgesv, DGESV, (n, 0, tmp_data, n, pipvt, dummy, n, info)); + F77_XFCN (dgetrf, DGETRF, (n, n, tmp_data, n, pipvt, info)); if (f77_exception_encountered) (*current_liboctave_error_handler) ("unrecoverable error in dgesv");