#1569 Add a wrapper to set the divergence tolerance for SNES
Sometimes the penalization term (boundary condition, incompressility etc) can make the norm of the residual too high relative to the first iteration. This causes issues for adaptative loading strategies where the norm of the residual of the first iteration can be quite low : the default value (1e4) for the divergence tolerance makes the solver fail.