[petsc-dev] Picard with TS

[Mark Adams]

I tried using Picard inside of TS, and SNES diverged. Is this supposed to
be possible?  I tried to do do this by wrapping my operator in a SNES
version and simply adding it like this.

Mark

  ierr = TSSetIFunction(ts,NULL,LandIFunction,ctx);CHKERRQ(ierr);
  ierr = TSSetIJacobian(ts,J,J,LandIJacobian,ctx);CHKERRQ(ierr);
    SNES snes;
    Vec vec;
    ierr = VecDuplicate(uu,&vec);CHKERRQ(ierr);
    ierr = TSGetSNES(ts,&snes);CHKERRQ(ierr);
    ierr = SNESSetPicard(snes, vec, LandSNESFunction, J, J,
LandSNESJacobian, ctx);CHKERRQ(ierr);
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170101/a2a7bb3c/attachment.html>