# HG changeset patch # User jwe # Date 797561364 0 # Node ID 35950c102ac3ca382a9215a5168df987be896b43 # Parent 85d1899047e1a7f27f026b17e46171443ced3c90 [project @ 1995-04-11 00:49:24 by jwe] diff -r 85d1899047e1 -r 35950c102ac3 src/expm.cc --- a/src/expm.cc Tue Apr 11 00:45:35 1995 +0000 +++ b/src/expm.cc Tue Apr 11 00:49:24 1995 +0000 @@ -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));