Mercurial > octave-antonio
changeset 8757:79576d40acb6
fix unscaling in expm
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Mon, 16 Feb 2009 12:45:15 +0100 |
parents | d0755c9db5ed |
children | 83c9d60c3c47 |
files | scripts/ChangeLog scripts/linear-algebra/expm.m |
diffstat | 2 files changed, 6 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/ChangeLog Mon Feb 16 12:41:55 2009 +0100 +++ b/scripts/ChangeLog Mon Feb 16 12:45:15 2009 +0100 @@ -1,3 +1,7 @@ +2009-02-16 Jaroslav Hajek <highegg@gmail.com> + + * linear-algebra/expm.m: Fix invalid unscaling. + 2009-02-16 John W. Eaton <jwe@octave.org> * help/help.m: Shorten default help text.
--- a/scripts/linear-algebra/expm.m Mon Feb 16 12:41:55 2009 +0100 +++ b/scripts/linear-algebra/expm.m Mon Feb 16 12:45:15 2009 +0100 @@ -149,8 +149,8 @@ endfor ## inverse balancing. - ds = diag (s); - r = ds * r / ds; + d = diag (d); + r = d * r / d; r = r(p, p); ## Inverse trace reduction. if (trshift >0)