On Thu, Nov 14, 2013 at 7:38 PM, Geoffrey Irving <[email protected]> wrote:
> On Thu, Nov 14, 2013 at 5:35 PM, Matthew Knepley <[email protected]> > wrote: > > 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. > > That's what I meant by fail to test: if there was a bug related to > null spaces, it wouldn't show up in this example. > I do not understand. If you have Neumann conditions, nothing happens, so it would check that case. Matt > Thanks, > 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
