MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Solving very large systems of ODEs with NDSolve

  • To: mathgroup at
  • Subject: [mg115331] Solving very large systems of ODEs with NDSolve
  • From: Gabriel Landi <gtlandi at>
  • Date: Fri, 7 Jan 2011 04:13:58 -0500 (EST)

Greetings all.

The problem I have in my hands is to solve a very large (500 to 10 000 or +)
system of non-linear dense ODEs (the Landau-Lifshitz equations in

In general, I think what I perhaps would like to ask is general advice on
implementing large systems in the NDSolve Framework
In particular, some of the problems that arose are:

 - In certain cases, I am not interested in storing the solution to all the
variables for all t. All I need is the value of each dependent variable at
the final integration time. I have been using combos of DependentVariables
and the EventLocator method but so far I haven't been able to formulate the
problem in a efficient way.

 - Each equation has a invariant and using the projection method to control
the error on these invariants turned out to be a very nice thing.
Unfortunately, if the system gets larger than 100 or so equations, the
kernel immediately crashes when I call this method (with no message at all).
I am not sure why since, when memory is the problem, a message is issued.
Anyone have any thoughts on this?

Thanks in advance for any suggestions,


Gabriel Landi

  • Prev by Date: Re: check if a square matrix is diagonal
  • Next by Date: Re: FillingStyle -> Striped ?
  • Previous by thread: Re: Extracting pixel values from an image under a polygon (region of interest)
  • Next by thread: Re: Solving very large systems of ODEs with NDSolve