Mercurial > hg > octave-nkf
diff src/expm.cc @ 52:9eda1009cf95
[project @ 1993-08-11 21:05:59 by jwe]
| author | jwe |
|---|---|
| date | Wed, 11 Aug 1993 21:05:59 +0000 |
| parents | 6028dcac27ef |
| children | cd4df7ad58fa |
line wrap: on
line diff
--- a/src/expm.cc +++ b/src/expm.cc @@ -60,7 +60,7 @@ extern int empty_arg (tree_constant&); tree_constant -matrix_exp (const tree_constant& a) +matrix_exp (tree_constant& a) { tree_constant retval; tree_constant tmp = a.make_numeric (); @@ -132,12 +132,12 @@ fortran_vec (), &n_cols, work.fortran_vec ()); - sqpow = 1.0 + log (inf_norm) / log (2.0); + sqpow = (int) (1.0 + log (inf_norm) / log (2.0)); // Check whether we need to square at all. - if (sqpow <= 0.0) - sqpow = 0.0; + if (sqpow < 0) + sqpow = 0; else { for (inf_norm = 1.0, i = 0; i < sqpow; i++) @@ -235,12 +235,12 @@ m.fortran_vec (), &n_cols, work.fortran_vec ()); - sqpow = 1.0 + log (inf_norm) / log (2.0); + sqpow = (int) (1.0 + log (inf_norm) / log (2.0)); // Check whether we need to square at all. - if (sqpow <= 0.0) - sqpow = 0.0; + if (sqpow < 0) + sqpow = 0; else { for (inf_norm = 1.0, i = 0; i < sqpow; i++)