Article View: sci.math.num-analysis
Article #1133Re: Numerical solution of ODEs with Bulirsch-Stoer extrapolation
From: Mike Hosea
Date: Thu, 28 Jan 1993 17:32
Date: Thu, 28 Jan 1993 17:32
30 lines
2328 bytes
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!mhoseaReferences:
<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