[ODE] what integration method is used ?
nlin at nlin.net
Mon Nov 18 03:35:02 2002
Nguyen Binh <ngbinh at glassegg.com> wrote:
> On Fri, Nov 15, 2002 at 02:28:59PM +0100, Torsten Hans wrote:
> TH> what basic integration scheme does ODE use ? is it an
> TH> explicit or implicit integration method ? does it use
> TH> runge kutta or something else ?
> Refer to the 0.3 docs, you'll see that ODE curently use basic Euler
Not quite. Somewhere in the docs is mentioned that ODE uses a
_time_ _stepping_ Euler integration scheme. This is a bit more hairy
to extend to higher order schemes. For the gory details and some
links to research papers see the thread started in April 2002 at
> TH> is it possible to use a different integration scheme ?
> Sure, But maybe you'll have to do it yourself or waiting for future
> ODE release.
A discussion on this topic in March 2002:
nlin> The problem with all of this ... is the interplay of all of this
nlin> with the time-stepping scheme.
russ>actually i'm not sure what a higher order time stepping scheme
russ>will look like in this case. i'd have to sit down and read some
russ>papers and write some equations. it's not a problem i have
Don't let this stop you - I'd also be interested in seeing the performance
and accuracy of higher-order integration schemes with ODE - but be aware
that this isn't a completely trivial thing to do because of the time-stepping
(rather than instantaneous) view of the equations to be solved.