[SciPy-user] Eigenvalues of the Floquet matrix

Rob Clewley rob.clewley@gmail....
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...
>
> Nils
>

... 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!

-Rob


More information about the SciPy-user mailing list