On Thu, Nov 14, 2013 at 7:32 PM, Geoffrey Irving <[email protected]> wrote:
> 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); > No, its just idempotent in that case. Matt Geoffrey > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
