use neumann boundary conditions in poission solver