[SciPy-user] Eigenvalues of the Floquet matrix

[Nils Wagner]

Hi all,

I am interested in the stability of time periodic ODE's of
the form

\dot{y} = A(t) y,  A(t)=A(t+T)          (1)

I have used scipy.integrate.ode to compute the eigenvalues
of the Floquet matrix. See attachment for details.

The eigenvalues \lambda of the Floquet matrix are called 
multipliers.
Multipliers of system (1) possess symmetry: If \lambda is 
the
multiplier, then 1/\lambda is also a multiplier.

The system is stable iff all eigenvalues are inside
the unit circle.

However, I cannot reproduce the symmetry of the
numerical multipliers computed by scipy, e.g.

>>> evals
array([ -2.55239771e-02+0.j        , 
 -2.39509330e-02+0.99921743j,
         -2.39509330e-02-0.99921743j,  -3.85602951e+01+0.j 
       ])
>>> 1./evals
array([ -3.91788472e+01-0.j        , 
 -2.39746890e-02-1.00020852j,
         -2.39746890e-02+1.00020852j,  -2.59334115e-02-0.j 
       ])

How can I improve the numerical results concerning the
symmetry ?

The computation of a stability chart is a time-consuming
task even for low dimensional systems.

How can one accelerate the process ?

Is it possible to parallelize the integration, I mean
each processor could solve (1) for a new set of
initial conditions ?

Any pointer would be appreciated.

BTW, the example is taken from a recent paper by 
Seyranian.

Thanks in advance.

                       Nils

  

-------------- next part --------------
A non-text attachment was scrubbed...
Name: seyranian.py
Type: text/x-python
Size: 1941 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20090314/6e78dd1b/attachment.py>