fix compilation
testing neumann boundary condition, the solution is still divided by 2 and the non linear solver still doesn't converge