From f81186f1c0e5b69d2e17924b56c422a4852a0e78 Mon Sep 17 00:00:00 2001
From: Ian Bell <ian.bell@nist.gov>
Date: Tue, 12 Oct 2021 19:06:35 -0400
Subject: [PATCH] Update tests for cubics

---
 src/tests/catch_test_cubics.cxx | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/src/tests/catch_test_cubics.cxx b/src/tests/catch_test_cubics.cxx
index 52e66b7..ff5f439 100644
--- a/src/tests/catch_test_cubics.cxx
+++ b/src/tests/catch_test_cubics.cxx
@@ -12,10 +12,13 @@ TEST_CASE("Test construction of cubic", "[cubic]")
     auto modelSRK = canonical_SRK(Tc_K, pc_Pa, acentric);
     auto modelPR = canonical_PR(Tc_K, pc_Pa, acentric);
 
-    double T = 300, rho = 300;
-    Eigen::ArrayXd molefrac(2); molefrac = 0.5;
+    double T = 800, rho = 5000;
+    auto molefrac = (Eigen::ArrayXd(3) << 0.5, 0.3, 0.2).finished();
     
     auto Ar02SRK = TDXDerivatives<decltype(modelSRK)>::get_Ar02(modelSRK, T, rho, molefrac);
+    auto Ar01PR = TDXDerivatives<decltype(modelPR)>::get_Ar01(modelPR, T, rho, molefrac);
     auto Ar02PR = TDXDerivatives<decltype(modelPR)>::get_Ar02(modelPR, T, rho, molefrac);
+    auto Ar03PR = TDXDerivatives<decltype(modelPR)>::get_Ar0n<3>(modelPR, T, rho, molefrac)[3];
+    auto Ar04PR = TDXDerivatives<decltype(modelPR)>::get_Ar0n<4>(modelPR, T, rho, molefrac)[4];
     int rr = 0;
 }
\ No newline at end of file
-- 
GitLab