From f3a0092b42e16b850efd18c1a108f7985b12ad39 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Thu, 5 Jun 2014 22:40:20 +0000 Subject: [PATCH] Fixed a error introduced by Speth that caused the brine phase pop problems to fail. --- src/equil/vcs_phaseStability.cpp | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/equil/vcs_phaseStability.cpp b/src/equil/vcs_phaseStability.cpp index 98a29285b..33b5b0b45 100644 --- a/src/equil/vcs_phaseStability.cpp +++ b/src/equil/vcs_phaseStability.cpp @@ -601,7 +601,9 @@ double VCS_SOLVE::vcs_phaseStabilityTest(const size_t iph) } for (size_t k = 0; k < nsp; k++) { - fracDelta_new[k] = std::min(fracDelta_new[k], 1e-13); + if (fracDelta_new[k] < 1.0E-13) { + fracDelta_new[k] =1.0E-13; + } } bool converged = false; for (int its = 0; its < 200 && (!converged); its++) {