[SciPy-user] Eigenvalues of the Floquet matrix
Sat Mar 14 17:07:12 CDT 2009
> Here is the product of X^T Y
> [[ 0.99682762 0.00427521 -0.00417167 0.00414188]
> [ 0.00212284 0.99717958 0.00329442 -0.00330161]
> [ 0.00122548 -0.00120682 0.98516109 0.01544185]
> [-0.00108133 0.00108144 0.01406736 0.98537091]]
> Again, I cannot reproduce the theory. The off-diagonal
> elements should be zero...
... and the on-diagonal elements should be 1. I suspect you basically
have it set up right, it's just that there must be some numerical
ill-conditioning somewhere in the matrix solving (or inaccuracy in the
integration) that leads to these fairly large errors. I'll take a look
when I have a chance in the next couple of days!
More information about the SciPy-user