Testing interface implementation for Vexp.

numcosmo/model/nc_hicosmo_qgw.c

@@ -72,36 +72,6 @@ nc_hicosmo_qgw_init (NcHICosmoQGW *cosmo_qgw)
   self->place_holder = 0;
 }
 
-static void
-_nc_hicosmo_qgw_set_property (GObject *object, guint prop_id, const GValue *value, GParamSpec *pspec)
-{
-  NcHICosmoQGW *qgw = NC_HICOSMO_QGW (object);
-
-  g_return_if_fail (NC_IS_HICOSMO_QGW (object));
-
-  switch (prop_id)
-  {
-    default:
-      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
-      break;
-  }
-}
-
-static void
-_nc_hicosmo_qgw_get_property (GObject *object, guint prop_id, GValue *value, GParamSpec *pspec)
-{
-  NcHICosmoQGW *qgw = NC_HICOSMO_QGW (object);
-
-  g_return_if_fail (NC_IS_HICOSMO_QGW (object));
-
-  switch (prop_id)
-  {
-    default:
-      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
-      break;
-  }
-}
-
 static void
 nc_hicosmo_qgw_finalize (GObject *object)
 {
@@ -127,9 +97,7 @@ nc_hicosmo_qgw_class_init (NcHICosmoQGWClass *klass)
   NcHICosmoClass *parent_class = NC_HICOSMO_CLASS (klass);
   NcmModelClass *model_class   = NCM_MODEL_CLASS (klass);
 
-  model_class->set_property = &_nc_hicosmo_qgw_set_property;
-  model_class->get_property = &_nc_hicosmo_qgw_get_property;
-  object_class->finalize    = nc_hicosmo_qgw_finalize;
+  object_class->finalize = nc_hicosmo_qgw_finalize;
 
   ncm_model_class_set_name_nick (model_class, "QGW", "QGW");
   ncm_model_class_add_params (model_class, NC_HICOSMO_QGW_SPARAM_LEN, 0, PROP_SIZE);

tests/test_py_hicosmo_Vexp.py

@@ -24,6 +24,7 @@
 
 """Tests on NcHICosmoVexp cosmological model."""
 
+from itertools import product
 from numpy.testing import assert_allclose
 import numpy as np
 
@@ -165,6 +166,78 @@ def test_serialize():
     assert vexp2.props.xb == vexp.props.xb  # pylint: disable=no-member
 
 
+def test_eval_at():
+    """Evaluate NcHICosmoVexp at a given time."""
+    vexp = Nc.HICosmoVexp()
+    assert vexp is not None
+
+    tau_min = vexp.tau_min()
+    tau_max = vexp.tau_max()
+
+    tau_a = np.linspace(tau_min, tau_max, 1000)
+
+    for tau in tau_a:
+        assert np.isfinite(vexp.x_tau(tau))
+        assert np.isfinite(vexp.alpha(tau))
+        assert np.isfinite(vexp.phi(tau))
+        assert np.isfinite(vexp.Ricci_scale(tau))
+        assert np.all(np.isfinite(vexp.x_y(tau)))
+
+
+def test_iadiab_eval_at():
+    """Evaluate NcHICosmoVexp implementation of NcHIPertAdiab at a given time."""
+    vexp = Nc.HICosmoVexp()
+    assert vexp is not None
+
+    tau_min = vexp.tau_min()
+    tau_max = vexp.tau_max()
+
+    tau_a = np.linspace(tau_min, tau_max, 1000)
+    k_a = np.geomspace(1.0e-3, 1.0e3, 5)
+
+    for tau, k in product(tau_a, k_a):
+        assert np.isfinite(Nc.HIPertIAdiab.eval_F1(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIAdiab.eval_m(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIAdiab.eval_nu(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIAdiab.eval_xi(vexp, tau, k))
+
+
+def test_igw_eval_at():
+    """Evaluate NcHICosmoVexp implementation of NcHIPertIGW at a given time."""
+    vexp = Nc.HICosmoVexp()
+    assert vexp is not None
+
+    tau_min = vexp.tau_min()
+    tau_max = vexp.tau_max()
+
+    tau_a = np.linspace(tau_min, tau_max, 1000)
+    k_a = np.geomspace(1.0e-3, 1.0e3, 5)
+
+    for tau, k in product(tau_a, k_a):
+        assert np.isfinite(Nc.HIPertIGW.eval_F1(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIGW.eval_m(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIGW.eval_nu(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIGW.eval_xi(vexp, tau, k))
+
+
+def test_iem_eval_at():
+    """Evaluate NcHICosmoVexp implementation of NcHIPertIEM at a given time."""
+    vexp = Nc.HICosmoVexp()
+    assert vexp is not None
+
+    tau_min = vexp.tau_min()
+    tau_max = vexp.tau_max()
+
+    tau_a = np.linspace(tau_min, tau_max, 1000)
+    k_a = np.geomspace(1.0e-3, 1.0e3, 5)
+
+    for tau, k in product(tau_a, k_a):
+        assert np.isfinite(Nc.HIPertIEM.eval_F1(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIEM.eval_m(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIEM.eval_nu(vexp, tau, k))
+        assert np.isfinite(Nc.HIPertIEM.eval_xi(vexp, tau, k))
+
+
 if __name__ == "__main__":
     test_init()
     test_hubble_negative_dphi()