Periodic Poisson compatibility ($\int f=0$)

On a periodic domain, integrating $\Delta \Phi =f$ gives $\int f=0$ because boundary fluxes cancel. On a periodic grid the row sum of the discrete Laplacian is zero, so ${\sum }_j{f}_j=0$ is necessary. The constant gravitational potential remains an arbitrary gauge.