[SciPy-user] ode on complex systems

Nils Wagner nwagner@iam.uni-stuttgart...
Fri Nov 16 09:25:53 CST 2007


On Fri, 16 Nov 2007 16:19:08 +0100
  Fabrice Silva <Fabrice.Silva@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@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



More information about the SciPy-user mailing list