About SNES solver of TS objects

[Hong Zhang]

Can you send us a simplified code on what you've tried?
What ts_type do you use?

Hong

On Thu, 11 Sep 2008, tsjb00 wrote:

>
> Hi!
>
> I use TS object to solve a non-linear problem du/dt=f(u,t). I would like to check the non-linear solver performance in TSStep. I tried to output the residuals by setting a monitor. So far my attempt didn't work. Please let me know if there are functions that calculate residuals of iterations in each TS time step. If not, how can I output the solver iteration/residuals correctly while using TS objects.
>
> Many thanks!
>
> P.S. The following are the code I tried:
>
>     PetscViewerASCIIOpen(PETSC_COMM_WORLD,"snes.log",&rviewer);
>     ierr = TSGetSNES(ts,&ts_snes);
>     ierr = SNESMonitorSet(ts_snes,SNESMonitorDefault,&rviewer,PETSC_NULL);
>
> I did get residual output in snes.log, but the format was very weird. Each line started with countless blank/space and at the very end was the SNES iteration info.  As a result, the file took a lot of memory space even with one line of outputs.
>
> I also tried:
>     PetscViewerASCIIMonitorCreate(PETSC_COMM_WORLD,"snes.log",0,&rviewer);
>     ierr = TSGetSNES(ts,&ts_snes);
>     ierr = SNESMonitorSet(ts_snes,SNESMonitorDefault,&rviewer,PETSC_NULL);
> which didn't work at all.
> _________________________________________________________________
> 一点即聊，MSN推出新功能“点我！”
> http://im.live.cn/click/
>
>