Minor improvements

thermo/eos.py

@@ -1303,7 +1303,7 @@ def set_from_PT(self, Vs, only_l=False, only_g=False):
             self.phase = "l/g"
         else:
             # Even in the case of three real roots, it is still the min/max that make sense
-            print([self.T, self.P, b, self.delta, self.epsilon, self.a_alpha, "coordinates of failure"])
+            # print([self.T, self.P, b, self.delta, self.epsilon, self.a_alpha, "coordinates of failure"])
             if self.multicomponent:
                 extra = f", zs is {self.zs}"
             else:

thermo/eos_volume.py

@@ -860,6 +860,7 @@ def volume_solutions_halley(T, P, b, delta, epsilon, a_alpha):
     # the value of the second term plus P is equal to P.
     if a_alpha/(b*(b + delta) + epsilon) + P == P:
         return (b + R*T/P, 0.0, 0.0)
+
     # Run this first, before the low P criteria
     if a_alpha > 1e4:
         V_possible = high_alpha_one_root(T, P, b, delta, epsilon, a_alpha)