diff --git a/dumux/porousmediumflow/richards/implicit/newtoncontroller.hh b/dumux/porousmediumflow/richards/implicit/newtoncontroller.hh index 0c2db36640c4d9c9a33f495a7fcac94b1bdd2a73..e33313d4501e577d2205a09f50ccc613b72c0826 100644 --- a/dumux/porousmediumflow/richards/implicit/newtoncontroller.hh +++ b/dumux/porousmediumflow/richards/implicit/newtoncontroller.hh @@ -80,7 +80,8 @@ public: { ParentType::newtonUpdate(uCurrentIter, uLastIter, deltaU); - if ((!GET_PARAM_FROM_GROUP(TypeTag, bool, Newton, UseLineSearch))) + if ((!GET_PARAM_FROM_GROUP(TypeTag, bool, Newton, UseLineSearch)) + && GET_PARAM_FROM_GROUP(TypeTag, bool, Newton, EnableChop )) { // do not clamp anything after 5 iterations if (this->numSteps_ > 4) diff --git a/test/porousmediumflow/richards/implicit/test_boxrichards.input b/test/porousmediumflow/richards/implicit/test_boxrichards.input index 41619629976aa26cac47fa78c24208efecbe4abc..60d877dac8979f214dbf6243cc9de97635d9fd54 100644 --- a/test/porousmediumflow/richards/implicit/test_boxrichards.input +++ b/test/porousmediumflow/richards/implicit/test_boxrichards.input @@ -11,8 +11,9 @@ Name = richardslensbox EnableGravity = 1 # enable gravity [Implicit] -EnablePartialReassemble = 1 # enable partial reassembly of the Jacobian matrix -EnableJacobianRecycling = 1 # +EnablePartialReassemble = true # enable partial reassembly of the Jacobian matrix +EnableJacobianRecycling = true # [Newton] TargetSteps = 18 # set the "desireable" number of newton iterations of a time step +EnableChop = true # chop for better convergence diff --git a/test/porousmediumflow/richards/implicit/test_ccrichards.input b/test/porousmediumflow/richards/implicit/test_ccrichards.input index effd06edc0c35f4788f9859c01abbdb7b81e6881..5ed3375d5d264884f0d2326f019c2edec3262ab8 100644 --- a/test/porousmediumflow/richards/implicit/test_ccrichards.input +++ b/test/porousmediumflow/richards/implicit/test_ccrichards.input @@ -11,8 +11,9 @@ Name = richardslenscc EnableGravity = 1 # enable gravity [Implicit] -EnablePartialReassemble = 1 # enable partial reassembly of the Jacobian matrix -EnableJacobianRecycling = 1 # +EnablePartialReassemble = true # enable partial reassembly of the Jacobian matrix +EnableJacobianRecycling = true # [Newton] TargetSteps = 18 # set the "desireable" number of newton iterations of a time step +EnableChop = true # chop for better convergence