[petsc-users] snes ex12 question about DMPlexProjectFunctionLocal

[Geoffrey Irving]

If this line in ex12 is part of Dirichlet boundary conditions, why is
it done unconditionally even for Neumann?  Is this a mistake that
fails to test correct Neumann conditions?

    ierr = DMPlexProjectFunctionLocal(dm, user.fe, user.exactFuncs,
INSERT_BC_VALUES, userJ.u);CHKERRQ(ierr);

Geoffrey