Allow up to 30 tries at the saturation solver

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-11-29 18:38:55 -05:00
parent e306a67f1e
commit ae6b5f85df

View File

@@ -1079,8 +1079,8 @@ void SaturationSolvers::saturation_T_pure_Maxwell(HelmholtzEOSMixtureBackend &HE
}
iter++;
if (iter > 20){
throw SolutionError(format("Maxwell solver did not converge after 20 iterations; rhoL: %0.16Lg rhoV: %0.16Lg error: %Lg dvL/vL: %Lg dvV/vV: %Lg pL: %Lg pV: %Lg\n", rhoL, rhoV, error, DeltavL/vL, DeltavV/vV, pL, pV));
if (iter > 30){
throw SolutionError(format("Maxwell solver did not converge after 30 iterations; rhoL: %0.16Lg rhoV: %0.16Lg error: %Lg dvL/vL: %Lg dvV/vV: %Lg pL: %Lg pV: %Lg\n", rhoL, rhoV, error, DeltavL/vL, DeltavV/vV, pL, pV));
}
}
while ((SatL->p() < 0) || (error > 1e-5 && small_step_count < 4));