🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: sci.math.num-analysis
Article #1133

Re: Numerical solution of ODEs with Bulirsch-Stoer extrapolation

#1133
From: Mike Hosea
Date: Thu, 28 Jan 1993 17:32
30 lines
2328 bytes
In article <C1H5uF.5HC@rice.edu> dharry@cyclone.rice.edu (Dennis Harry) writes:
>In article <1993Jan20.054301.455@seas.smu.edu>, mhosea@seas.smu.edu (Mike Hosea) writes:
> 
>|> The reason extrapolation methods get very accurate results with large step
>|> sizes is because they can use very high order approximations.  However,
>|> my research has shown that when compared to conventional Runge-Kutta
>|> and Adams methods of the same order, the extrapolation formulas are not
>|> as efficient in terms of function evaluations per unit step.
>                                                 ^^^^^^^^^^^^^
>Yes, but the point is that they may take much larger steps.  Ideally,
>their relative inefficiency per unit step is (more than) offset by having
>to take fewer steps?

No, it's not.  Not when you consider the competition.  It is true that
an extrapolation code may perform better than a fixed order Runge-Kutta
code when the tolerance is very stringent, but this is only because they
can go to higher order.  However, like the extrapolation formulas, the
Adams-Moulton formulas are variable order, and at each order (especially
high orders) the Adams-Moulton formulas are much more efficient.  If you
compare the performances of codes you may not see this since the Adams
code might not attept to use the higher orders--this would have nothing
to do with the methods themselves.  Anyway, I am not arguing from intuition.
Check out the mathematics in my report or any reputable numerical
comparison like

  L.F. Shampine and L.S. Baca  1986  "Fixed versus variable order
  Runge-Kutta", ACM TOMS, 12, no.  pp. 1-23.
-- 
Regards,
Mike Hosea (mhosea@seas.smu.edu)
Technical Details
Message-ID: <1993Jan28.173203.7930@seas.smu.edu>
Path: arc1-ca.pugleaf.net!sparky!uunet!europa.eng.gtefsd.com!emory!gatech!darwin.sura.net!sgiblab!nec-gw!netkeeper!vivaldi!seas.smu.edu!mhosea
References: <1993Jan20.024502.23189@leland.Stanford.EDU> <1993Jan20.054301.455@seas.smu.edu> <C1H5uF.5HC@rice.edu>

Show Full Headers
Complete Article Headers
Newsgroups: sci.math.num-analysis
Path: arc1-ca.pugleaf.net!sparky!uunet!europa.eng.gtefsd.com!emory!gatech!darwin.sura.net!sgiblab!nec-gw!netkeeper!vivaldi!seas.smu.edu!mhosea
From: mhosea@seas.smu.edu (Mike Hosea)
Subject: Re: Numerical solution of ODEs with Bulirsch-Stoer extrapolation
Message-ID: <1993Jan28.173203.7930@seas.smu.edu>
Sender: mhosea@seas.smu.edu
Nntp-Posting-Host: express.seas.smu.edu
Organization: SMU - School of Engineering & Applied Science - Dallas
References: <1993Jan20.024502.23189@leland.Stanford.EDU> <1993Jan20.054301.455@seas.smu.edu> <C1H5uF.5HC@rice.edu>
Date: Thu, 28 Jan 1993 17:32:03 GMT
Lines: 30