Mercurial > octave-nkf

changeset 8757:79576d40acb6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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)