On Jun 23, 2010, at 12:14 PM, Jed Brown wrote:
> On Wed, 23 Jun 2010 10:22:19 -0400, Xuan YU <xxy113 at psu.edu> wrote: >> I got the correct results! Thanks! > > If you are feeling like experimenting, you could try -ts_type gl which > is a very new adaptive order adaptive step method for stiff > systems. It > seems to be okay on your problem, but the controller is currently > pretty > fragile (it's possible for step size to go to zero due to poor error > estimates, or for length/order to oscillate spuriously with several > step > period). > > Jed > I tried these method. Pseudo and Beuler work well. gl kept running after t=1.665007e+01. Actually, I will solve a global ODE system in form of y'=f(t,y;x), x is the forcing. The problem has been solved by cvode(sundials). I want to use petsc to improve efficiency.
