Mercurial > hg > octave-lojdl
changeset 1246:35950c102ac3
[project @ 1995-04-11 00:49:24 by jwe]
author | jwe |
---|---|
date | Tue, 11 Apr 1995 00:49:24 +0000 |
parents | 85d1899047e1 |
children | 7bddfd169d27 |
files | src/expm.cc |
diffstat | 1 files changed, 6 insertions(+), 8 deletions(-) [+] |
line wrap: on
line diff
--- a/src/expm.cc +++ b/src/expm.cc @@ -46,11 +46,11 @@ extern "C" { - double F77_FCN (dlange) (const char*, const int*, const int*, - const double*, const int*, double*); + double F77_FCN (dlange) (const char*, const int&, const int&, + const double*, const int&, double*); - double F77_FCN (zlange) (const char*, const int*, const int*, - const Complex*, const int*, double*); + double F77_FCN (zlange) (const char*, const int&, const int&, + const Complex*, const int&, double*); } DEFUN_DLD_BUILTIN ("expm", Fexpm, Sexpm, 2, 1, @@ -135,8 +135,7 @@ // Preconditioning step 3: scaling. ColumnVector work(nc); - inf_norm = F77_FCN (dlange) ("I", &nc, &nc, - m.fortran_vec (), &nc, + inf_norm = F77_FCN (dlange) ("I", nc, nc, m.fortran_vec (), nc, work.fortran_vec ()); sqpow = (int) (1.0 + log (inf_norm) / log (2.0)); @@ -228,8 +227,7 @@ // Preconditioning step 3: scaling. ColumnVector work (nc); - inf_norm = F77_FCN (zlange) ("I", &nc, &nc, m. - fortran_vec (), &nc, + inf_norm = F77_FCN (zlange) ("I", nc, nc, m.fortran_vec (), nc, work.fortran_vec ()); sqpow = (int) (1.0 + log (inf_norm) / log (2.0));