Fix Neumann point conditions for nonlinear problems

Need to reset the integrated values vector to zero before it is
calculated.
Add load incremented Neumann point boundary conditions
Fix conflict