prevent negatives

BioSTEAMDevelopmentGroup · Sep 17, 2024 · 9409981 · 9409981
1 parent db50bbb
commit 9409981
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/thermosteam/equilibrium/vle.py b/thermosteam/equilibrium/vle.py
@@ -29,6 +29,8 @@
 
 @njit(cache=True)
 def xy(x, Ks):
+    x[x < 0] = 1e-16
+    x /= x.sum()
     y = x * Ks
     y /= y.sum()
     return x, y