[petsc-users] help about my first petsc program

[Xuan YU]

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.