[petsc-users] TS and SNES

[Jed Brown]

On Tue, Nov 8, 2011 at 04:14, Konstantinos Kontzialis <ckontzialis at lycos.com
> wrote:

>    ierr = TSSetSolution(sys.ts, sys.gsv);
>    CHKERRQ(ierr);
>
>    ierr = TSSetIFunction(sys.ts, PETSC_NULL, base_residual_implicit, &sys);
>    CHKERRQ(ierr);
>

Provide a residual vector here (instead of PETSC_NULL). Since you provided
a state Vec in TSSetSolution(), it is possible to create a Vec for the
residual. I'll add that logic and improve the error message, but for now,
just pass in the residual Vec.


>
>    ierr = TSGetSNES(sys.ts, &sys.snes);
>    CHKERRQ(ierr);
>
>    ierr = MatCreateSNESMF(sys.snes, &sys.J);
>    CHKERRQ(ierr);
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111108/c2bfbdc7/attachment.htm>