diff --git a/CMakeLists.txt b/CMakeLists.txt
index 2fce861ea504807e5e5be436d10944ec5f5ffd39..2b423a41f00e99475c64694c47288c2fa00534e7 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -390,9 +390,11 @@ add_test(test_general_box    bin/runTest.sh references/generallens_box-reference
 add_test(test_general_dec    bin/runTest.sh references/generallens_decoupled-reference.vtu generallens_decoupled-00003.vtu test/common/generalproblem/test_generalproblem_2p --decoupled 1e2 2e1)
 add_test(test_spline         test/common/spline/test_spline)
 
-add_test(test_fluidsystems test/material/fluidsystems/test_fluidsystems)
-add_test(test_ncpflash     test/material/ncpflash/test_ncpflash)
-add_test(test_tabulation   test/material/tabulation/test_tabulation)
+add_test(test_fluidsystems    test/material/fluidsystems/test_fluidsystems)
+add_test(test_immiscibleflash test/material/immiscibleflash/test_immiscibleflash)
+add_test(test_ncpflash        test/material/ncpflash/test_ncpflash)
+add_test(test_pengrobinson    test/material/pengrobinson/test_pengrobinson)
+add_test(test_tabulation      test/material/tabulation/test_tabulation)
 
 add_test(test_1p       bin/runTest.sh references/1ptest-reference.vtu        1ptest-00002.vtu        test/boxmodels/1p/test_1p             --parameter-file=parameters/test_1p.input       -gridFile test/boxmodels/1p/grids/test_1p_2d.dgf               -tEnd 1    -dtInitial 1)
 add_test(test_1p2c     bin/runTest.sh references/outflow-reference.vtu       outflow-00010.vtu       test/boxmodels/1p2c/test_1p2c         --parameter-file=parameters/test_1p2c.input     -gridFile test/boxmodels/1p2c/grids/test_1p2c.dgf              -tEnd 100  -dtInitial 1)
diff --git a/test/material/pengrobinson/test_pengrobinson.cc b/test/material/pengrobinson/test_pengrobinson.cc
index 82c135aab1fe53ba6bebf7c8d236fc941ac70fad..297b88decaccec0dc0634af5ed8841740b3b17a3 100644
--- a/test/material/pengrobinson/test_pengrobinson.cc
+++ b/test/material/pengrobinson/test_pengrobinson.cc
@@ -217,7 +217,7 @@ int main(int argc, char** argv)
     Scalar rho_oRef = surfaceFluidState.density(oPhaseIdx);
 
     std::cout << "alpha[-] p[Pa] S_g[-] rho_o[kg/m^3] rho_g[kg/m^3] <M_o>[kg/mol] <M_g>[kg/mol] R_s[m^3/m^3] B_g[-] B_o[-]\n";
-    int n = 1000;
+    int n = 50;
     for (int i = 0; i < n; ++i) {
         Scalar minAlpha = 0.98;
         Scalar maxAlpha = 2.0;