home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / sci / math / numanal / 2553 < prev    next >
Encoding:
Internet Message Format  |  1992-08-26  |  2.1 KB
  1. Path: sparky!uunet!usc!sdd.hp.com!swrinde!mips!darwin.sura.net!convex!news.utdallas.edu!hermes.chpc.utexas.edu!jonathan
  2. From: jonathan@chpc.utexas.edu (Jonathan Thornburg)
  3. Newsgroups: sci.math.num-analysis
  4. Subject: Re: predictor corrector methods for ODE's
  5. Summary: PC methods are *not* obselete, many *excellent* codes use them
  6. Keywords: predictor corrector ODE integrate LSODE Hindmarsh DE/STEP Shampine Gordon
  7. Message-ID: <1992Aug26.184246.11311@chpc.utexas.edu>
  8. Date: 26 Aug 92 18:42:46 GMT
  9. References: <26AUG199218380252@robin.rp.csiro.au>
  10. Sender: jonathan@einstein.ph.utexas.edu
  11. Organization: U of Texas at Austin / Physics Dept / Center for Relativity
  12. Lines: 31
  13.  
  14. In article <26AUG199218380252@robin.rp.csiro.au>
  15. awhichel@robin.rp.csiro.au (Adrian Whichello) writes:
  16. >    I just read in Numerical Recipes that predictor-corrector methods for
  17. >solving initial value ODE's are dead, or at least on their way out. This
  18. >strikes me as rather strong, but I would like some other opinions. What do you
  19. >think? 
  20.  
  21. Numerical Recipes is (very) wrong here.  Predictor-Corrector (PC)
  22. methods for non-stiff initial value ODE integration are alive and
  23. well.  ODE experts (which I'm not) recommend them as being robust
  24. and efficient for a wide range of problems.  See, for example, the
  25. review and discussion in:
  26. \mybibitem{Gear, C. W.}{1981}
  27.         ``Numerical Solution of Ordinary Differential Equations:
  28.           Is There Anything Left to Do?''
  29.         SIAM Review {\bf 23}(1), 10--24.
  30.  
  31. Indeed, of the current high quality "black box" ODE codes, a large
  32. fraction use PC methods.  Notable examples include Hindmarsh's LSODE
  33. suite, and the Shampine & Gordon DE/STEP package.
  34.  
  35. PC methods are also "culturally similar" to the backwards
  36. differentiation multistep methods which seem to be the current
  37. favorites for stiff ODEs, which makes it easier to combine both in
  38. a single general purpose code or code suite.  (For example, LSODE
  39. does this.)
  40.  
  41. - Jonathan Thornburg
  42.   <jonathan@einstein.ph.utexas.edu> or <jonathan@hermes.chpc.utexas.edu>
  43.   University of Texas at Austin / Physics Dept / Center for Relativity
  44.   and (for a few more months) U of British Columbia / {Astronomy,Physics}
  45.