[petsc-users] sundials and beuler

[Xuan YU]

On Jul 28, 2010, at 11:15 AM, Hong Zhang wrote:

> Run ex4 with '-ts_type sundials -nox -ts_sundials_monitor_steps'
>

it does output all the  time steps. But they are different in time
beuler
Timestep 0: time = 0, 2-norm error = 1.01507e-15, max norm error =  
3.10862e-15
Timestep 1: time = 0.000143637, 2-norm error = 0.00112483, max norm  
error = 0.00162124
Timestep 2: time = 0.000287274, 2-norm error = 0.00213968, max norm  
error = 0.00308653
Timestep 3: time = 0.000430911, 2-norm error = 0.00305264, max norm  
error = 0.0044073
Timestep 4: time = 0.000574548, 2-norm error = 0.00387125, max norm  
error = 0.00559423
Timestep 5: time = 0.000718184, 2-norm error = 0.0046026, max norm  
error = 0.00665731
Timestep 6: time = 0.000861821, 2-norm error = 0.00525327, max norm  
error = 0.00760586
Timestep 7: time = 0.00100546, 2-norm error = 0.00582943, max norm  
error = 0.00844859
Timestep 8: time = 0.0011491, 2-norm error = 0.00633684, max norm  
error = 0.00919364
Timestep 9: time = 0.00129273, 2-norm error = 0.00678087, max norm  
error = 0.00984859

sundials
Timestep 0: time = 0, 2-norm error = 1.01507e-15, max norm error =  
3.10862e-15
Timestep 1: time = 4.23858e-06, 2-norm error = 9.80672e-06, max norm  
error = 1.41163e-05
Timestep 2: time = 8.47716e-06, 2-norm error = 1.95841e-05, max norm  
error = 2.8191e-05
Timestep 3: time = 1.94633e-05, 2-norm error = 4.45115e-05, max norm  
error = 6.40773e-05
Timestep 4: time = 3.04494e-05, 2-norm error = 6.80772e-05, max norm  
error = 9.80074e-05
Timestep 5: time = 4.14354e-05, 2-norm error = 9.10699e-05, max norm  
error = 0.000131116
Timestep 6: time = 5.24215e-05, 2-norm error = 0.000113753, max norm  
error = 0.000163784
Timestep 7: time = 7.6572e-05, 2-norm error = 0.000162885, max norm  
error = 0.000234557
Timestep 8: time = 0.000135787, 2-norm error = 0.000278648, max norm  
error = 0.000401418
Timestep 9: time = 0.000195002, 2-norm error = 0.000388956, max norm  
error = 0.000560541

Can I custom the output time of sundials?
> Hong
>
> On Wed, Jul 28, 2010 at 9:17 AM, Xuan YU <xxy113 at psu.edu> wrote:
>> Hi
>>
>> I am running the example: src/ts/examples/tutorials/ex4.c
>>
>> For this example,  why sundials(./ex4 -ts_type sundials) only  
>> output the
>> first time step and the last time step, but beuler can output every  
>> time
>> step?
>>
>> How can I make sundials output every time step?
>>
>> Thanks!
>>
>

Xuan YU (俞烜)
xxy113 at psu.edu




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100728/c75bdc1f/attachment.htm>