[petsc-users] DMDA and boundary conditions

[zakaryah .]

For finite difference methods, using a DMDA with DMDA_BOUNDARY_GHOSTED is
great for Dirichlet boundary conditions, because I can just set the
Dirichlet values on the ghost points, correct?  It seems to me that the
analogous method for Neumann conditions would be to use a reflecting
boundary condition, and if the Neumann values are nonzero, these could be
added to the constant vector.  Is this correct?  Is DMDA_BOUNDARY_MIRROR
supposed to do this?  Is it implemented yet, or are there plans to
implement it soon?  Otherwise, is there any way to efficiently implement
Neumann conditions besides branching with if statements?  Thanks in advance!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171015/aa220090/attachment.html>