Fix Neumann point conditions for nonlinear problems

Whoops, needed to reset the integrated values vector to zero before it is calculated.

Also, currently the integration is done each time solver_matrices_static/dynamic_assemble is called, which is every time the residual is calculate for a nonlinear problem. I don't think this is necessary? Is there a better place this could be done?