Matthew Knepley <[email protected]> writes:

> When you and Peter changed the way SNES used to work, I asked you how I
> should
> change Plex. You told me to do this:
>
>     ierr = DMSNESSetFunctionLocal(dm,  (PetscErrorCode
> (*)(DM,Vec,Vec,void*)) DMPlexComputeResidualFEM, &user);CHKERRQ(ierr);
>     ierr = DMSNESSetJacobianLocal(dm,  (PetscErrorCode
> (*)(DM,Vec,Mat,Mat,MatStructure*,void*)) DMPlexComputeJacobianFEM,
> &user);CHKERRQ(ierr);

Can you find the email?  I don't think I would have recommended this.

> which I did. Moreover, this is how it always worked, namely specifying the
> local function.
> What exactly are you proposing to replace this?

DMPlexSNESSetFunctionFEM(), similar to DMDASNESSetFunctionLocal in
src/snes/examples/tutorials/ex5.c and ex19.c.

Attachment: pgpl8NOlk1Kqs.pgp
Description: PGP signature

Reply via email to