[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Expm giving NaN
From: |
Marco Caliari |
Subject: |
Expm giving NaN |
Date: |
Thu, 27 Sep 2007 09:51:26 +0200 (CEST) |
Dear maintainers,
the patch against CMatrix.cc in 2.9.14 fixing the "expm giving NaN" bug.
Here the result:
octave:1> version
ans = 2.9.14
octave:2> expm(1000*i-10)
ans = 2.5532e-05 + 3.7540e-05i
octave:3> expm([1000*i-10 0;0 1000*i-10])
ans =
2.5532e-05 + 3.7540e-05i 0.0000e+00 + 0.0000e+00i
0.0000e+00 + 0.0000e+00i 2.5532e-05 + 3.7540e-05i
I also rewrote the computation of the rational approximation avoiding a
matrix-scalar product and including a for loop. I'm not sure if it is
better.
Best regards,
Marco
CMatrix.patch
Description: Text document
- Expm giving NaN,
Marco Caliari <=