Mercurial > hg > octave-nkf
diff liboctave/dbleLU.cc @ 5275:23b37da9fd5b
[project @ 2005-04-08 16:07:35 by jwe]
author | jwe |
---|---|
date | Fri, 08 Apr 2005 16:07:37 +0000 |
parents | e35b034d3523 |
children | 4c8a2e4e0717 |
line wrap: on
line diff
--- a/liboctave/dbleLU.cc +++ b/liboctave/dbleLU.cc @@ -40,30 +40,30 @@ extern "C" { F77_RET_T - F77_FUNC (dgetrf, DGETRF) (const int&, const int&, double*, - const int&, int*, int&); + F77_FUNC (dgetrf, DGETRF) (const octave_idx_type&, const octave_idx_type&, double*, + const octave_idx_type&, octave_idx_type*, octave_idx_type&); } LU::LU (const Matrix& a) { - int a_nr = a.rows (); - int a_nc = a.cols (); - int mn = (a_nr < a_nc ? a_nr : a_nc); + octave_idx_type a_nr = a.rows (); + octave_idx_type a_nc = a.cols (); + octave_idx_type mn = (a_nr < a_nc ? a_nr : a_nc); ipvt.resize (mn); - int *pipvt = ipvt.fortran_vec (); + octave_idx_type *pipvt = ipvt.fortran_vec (); a_fact = a; double *tmp_data = a_fact.fortran_vec (); - int info = 0; + octave_idx_type info = 0; F77_XFCN (dgetrf, DGETRF, (a_nr, a_nc, tmp_data, a_nr, pipvt, info)); if (f77_exception_encountered) (*current_liboctave_error_handler) ("unrecoverable error in dgetrf"); else - ipvt -= 1; + ipvt -= static_cast<octave_idx_type> (1); } /*