Another error is needed

usnistgov · Sep 7, 2023 · 92515dc · 92515dc
1 parent d6cacec
commit 92515dc
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/include/teqp/models/saft/polar_terms.hpp b/include/teqp/models/saft/polar_terms.hpp
@@ -934,6 +934,9 @@ class MultipolarContributionGrayGubbins {
             // The packing fraction is defined by eta = pi/6*rho^*, so use the (temperature-dependent) eta to obtain rho^*
             rhostar = forceeval(packing_fraction/(static_cast<double>(EIGEN_PI)/6.0));
         }
+        else{
+            throw teqp::InvalidArgument("The method used to determine rho^* is invalid");
+        }
         return rhostar;
     }