JFNKgamma_lin = JFNKgamma_lin_max
JFNKgamma_lin = SEAICE_JFNKphi
JFNKgamma_lin = min(JFNKgamma_lin_max, JFNKgamma_lin)
JFNKgamma_lin = max(JFNKgamma_lin_min, JFNKgamma_lin)
FGMRESeps = JFNKgamma_lin * JFNKresidual
& JFNKgamma_lin,JFNKresidual
_RL JFNKgamma_lin