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

Re: NDSolve precision and velocity problem

• To: mathgroup at smc.vnet.net
• Subject: [mg91318] Re: NDSolve precision and velocity problem
• From: Joerg Schaber <schaber at molgen.mpg.de>
• Date: Fri, 15 Aug 2008 06:53:50 -0400 (EDT)
• References: <g81378$o$1@smc.vnet.net> <48A43754.1000105@gmail.com>

Thanks.

Actually, I tried Method->"Adams", but after 10 minutes waiting I gave up.
BDF seems to work fine though. I did't see it in the Documentation.

Jean-Marc Gulliet wrote:
> 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]
>
> Regards,
> - Jean-Marc
>
> [1] "LSODE", _Serial Fortran Solvers for ODE Initial Value Problems_,
> https://computation.llnl.gov/casc/odepack/odepack_home.html
>
> [2] _IDA Method for NDSolve_,
> http://reference.wolfram.com/mathematica/tutorial/NDSolveIDAMethod.html
>
> [3] "Adams methods", _NDSolve Method Plug-in Framework_,
> http://reference.wolfram.com/mathematica/tutorial/NDSolvePlugIns.html
>
> [4] _Monitoring and Selecting Algorithms_,
> http://reference.wolfram.com/mathematica/tutorial/MonitoringAndSelectingAlgorithms.html 
>

-- 
----------------------------------------------
Jörg Schaber
Max Planck Institute for Molecular Genetics
Ihnestrasse 63-73, 14195 Berlin, Germany
Phone: +49 30 804093 19, Fax: +49 30 804093 22
----------------------------------------------