Sense of parameters for sumColumns() reversed, to make it more intuitive.
PSSurface subclasses can now choose Neumann boundary condition.
Code compiles, has not yet been tested.  This is very experimental.
At least, it shouldnt break anything that doesnt use it.