diff --git a/tests/test_py_hipert_two_fluids.py b/tests/test_py_hipert_two_fluids.py
index 10690fc1..078a1007 100644
--- a/tests/test_py_hipert_two_fluids.py
+++ b/tests/test_py_hipert_two_fluids.py
@@ -112,3 +112,11 @@ def spec_params(Omegars=1.0e-5, w=1.0e-3, E0=1.0):
         Pk0 = s.eval(lnk, np.log(w))
         Pk1 = s_calib.eval(lnk, np.log(w))
         assert_allclose(Pk0, Pk1, rtol=1.0e-3)
+
+
+def test_evolve_array(two_fluids, cosmo_qgrw):
+    """Test NcHIPertTwoFluids evolve_array."""
+    two_fluids.props.reltol = 1.0e-9
+    m = two_fluids.evolve_array(cosmo=cosmo_qgrw, alphaf=-1.0e-1)
+
+    assert all(np.isfinite(m.dup_array()))