since I can easily manufacture solution to homogeneous Neumann on square, I add that geometry
because filename is a pointer, sizeof(filename)==8, so this fixes a weird bug
I think this successfully solves homogeneous Neumann Poisson equation