From a6dc994439d00c038e1d5d511f67285c46606deb Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Fri, 12 Oct 2012 20:34:59 +0000
Subject: [PATCH] Fixed a problem with the VCS equilibrium solver at constant
 HP or UP

A sign error when limiting the temperature step caused the solver to fail if
the initial temperature was higher than sqrt(Tlow*Thigh).
---
 src/equil/vcs_MultiPhaseEquil.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/equil/vcs_MultiPhaseEquil.cpp b/src/equil/vcs_MultiPhaseEquil.cpp
index aa9b7bc9e..72fa6a48b 100644
--- a/src/equil/vcs_MultiPhaseEquil.cpp
+++ b/src/equil/vcs_MultiPhaseEquil.cpp
@@ -294,7 +294,7 @@ int vcs_MultiPhaseEquil::equilibrate_HP(doublereal Htarget,
                 Tnew = sqrt(Tlow*Thigh);
                 dT = Tnew - Tnow;
                 if (dT < -200.) {
-                    dT = 200;
+                    dT = -200;
                 }
                 if (dT > 200.) {
                     dT = 200.;