Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive

MathGroup Archive 2008

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

Search the Archive

Re: NDSolve precision and velocity problem

  • To: mathgroup at
  • Subject: [mg91312] Re: NDSolve precision and velocity problem
  • From: Jean-Marc Gulliet <jeanmarc.gulliet at>
  • Date: Fri, 15 Aug 2008 06:52:44 -0400 (EDT)
  • Organization: The Open University, Milton Keynes, UK
  • References: <g81378$o$>

Joerg Schaber wrote:

> I have an ODE system that takes lots of time to solve with NDSolve 
> (several minutes),
> whereas the same system with the FORTRAN solver LSODA
> solves the problem in an second.
> When I decrease the precision, e.g. AccuracyGoal->6, NDSolve crashes. It 
> often crashes, because it finds complex values in a less-than-comparison.
> Maybe it's a method problem.
> How do I tell NDSolve, to use the same methods as the LSODEA solver?

If by the above you mean using either Adams or Backward Differentiation 
Formula (BDF) methods [1], for nonstiff or stiff cases, respectively, 
you can use the option *Method*.

The following might be worth reading (or at least browsing): [2, 3, 4]

- Jean-Marc

[1] "LSODE", _Serial Fortran Solvers for ODE Initial Value Problems_,

[2] _IDA Method for NDSolve_,

[3] "Adams methods", _NDSolve Method Plug-in Framework_,

[4] _Monitoring and Selecting Algorithms_,

  • Prev by Date: Re: ODE in mathematica. what is wrong?
  • Next by Date: Re: Another Out of Memory Problem
  • Previous by thread: Re: NDSolve precision and velocity problem
  • Next by thread: Re: NDSolve precision and velocity problem