Edit: well, with boundary conditions, G will be more complicated. But for example, with Neumann boundary conditions, it is clear to see that the remaining surface integral should not be negative (by comparing with the wikipedia page for Helmholtz decomposition). It should have the opposite sign. The other part of the surface integral should also be the negative of what is written.