# HG changeset patch # User John W. Eaton # Date 1294427728 18000 # Node ID f617ce2472f82866e69a2c901f77e54c02df4f76 # Parent 93a039fe681e97ad32dfa183a35c6cbba60dec99 fix logm for identity matrices diff --git a/scripts/ChangeLog b/scripts/ChangeLog --- a/scripts/ChangeLog +++ b/scripts/ChangeLog @@ -1,3 +1,8 @@ +2011-01-07 John W. Eaton + + * linear-algebra/logm.m: Don't call logm_pade_pf unless m > 0. + Bug #32075. + 2011-01-07 John W. Eaton * linear-algebra/logm.m: Style fixes. diff --git a/scripts/linear-algebra/logm.m b/scripts/linear-algebra/logm.m --- a/scripts/linear-algebra/logm.m +++ b/scripts/linear-algebra/logm.m @@ -84,7 +84,11 @@ warning ("logm: maximum number of square roots exceeded; results may still be accurate"); endif - s = logm_pade_pf (s - eye (size (s)), m); + s -= eye (size (s)); + + if (m > 1) + s = logm_pade_pf (s, m); + endif s = 2^k * u * s * u';