<div dir="ltr"><div class="gmail_default" style="font-size:small">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!</div></div>