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>
