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);

Reply via email to