[SciPy-user] ode on complex systems

[Nils Wagner]

On Fri, 16 Nov 2007 16:19:08 +0100
  Fabrice Silva <Fabrice.Silva at crans.org> wrote:
> 
> Le vendredi 16 novembre 2007 à 15:53 +0100, lorenzo 
>bolla a écrit :
>> Hi all,
>> can I use any of scipy.integrate.ode/odeint etc. for 
>>complex systems? How?
> 
> You may need to separate real and imaginary parts of 
>your complex
> unknowns :
> you need to transform
> y(t)=[y0(t), y1(t)] with y0 and y1 complex signals
> into:
> y(t)=[real(y0(t)), imag(y0(t)), real (y1(t)), 
>imag(y1(t))]
> 
> and modify your matrix M as well.
> 
> Can anyone confirm?
> -- 
>Fabricio
> 
> 
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user

AFAIK, you can also apply the patch which is available
at

http://projects.scipy.org/scipy/scipy/ticket/334

  
So you can integrate complex ODE's directly without
doubling the size. Please let me know how it works
for you.

Nils