6749245209

Committed 03 Nov 2023 07:10PM UTC coverage: 52.661% (+1.4%) from 51.227%

Build # 6749245209

Build Type

Pull #116

github

Committed by

web-flow

Commit Message

Merge 7138feffd into 69afcab6b

Pull Request Pull Request #116: Improving documentation and encapsulating objects

Run Details

8575 of 20622 branches covered (0.0%)

Branch coverage included in aggregate %.

1538 of 2274 new or added lines in 59 files covered. (67.63%)

251 existing lines in 35 files now uncovered.

43676 of 78600 relevant lines covered (55.57%)

3707329.33 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

15.38

/numcosmo/data/nc_data_bao_dv.c

/***************************************************************************
 *            nc_data_bao_dv.c
 *
 *  Thu Apr 22 15:31:32 2010
 *  Copyright  2010  Sandro Dias Pinto Vitenti
 *  <vitenti@uel.br>
 ****************************************************************************/
/*
 * numcosmo
 * Copyright (C) 2012 Sandro Dias Pinto Vitenti <vitenti@uel.br>
 *
 * numcosmo is free software: you can redistribute it and/or modify it
 * under the terms of the GNU General Public License as published by the
 * Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * numcosmo is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 * See the GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License along
 * with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

/**
 * SECTION:nc_data_bao_dv
 * @title: NcDataBaoDV
 * @short_description: Baryon oscillation data -- volume mean $D_V$.
 *
 * See [Eisenstein et al. (2005)][XEisenstein2005].
 *
 * The data is stored in a #NcDataBaoDV object. The data is stored in a
 * #NcDataGaussDiag base class object, which is a subclass of #NcmData.
 * The data represents the mean values of the volume mean $D_V$ at the redshift
 * $z$.
 *
 */

#ifdef HAVE_CONFIG_H
#  include "config.h"
#endif /* HAVE_CONFIG_H */
#include "build_cfg.h"

#include "data/nc_data_bao_dv.h"
#include "nc_enum_types.h"
#include "math/ncm_cfg.h"

enum
{
  PROP_0,
  PROP_DIST,
  PROP_Z,
  PROP_SIZE,
};

G_DEFINE_TYPE (NcDataBaoDV, nc_data_bao_dv, NCM_TYPE_DATA_GAUSS_DIAG);

static void
nc_data_bao_dv_init (NcDataBaoDV *bao_dv)
{
  bao_dv->dist = nc_distance_new (2.0);
  bao_dv->x    = NULL;
}

static void
nc_data_bao_dv_set_property (GObject *object, guint prop_id, const GValue *value, GParamSpec *pspec)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);

  g_return_if_fail (NC_IS_DATA_BAO_DV (object));

  switch (prop_id)
  {
    case PROP_DIST:
      nc_data_bao_dv_set_dist (bao_dv, g_value_get_object (value));
      break;
    case PROP_Z:
      ncm_vector_substitute (&bao_dv->x, g_value_get_object (value), TRUE);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_bao_dv_get_property (GObject *object, guint prop_id, GValue *value, GParamSpec *pspec)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);

  g_return_if_fail (NC_IS_DATA_BAO_DV (object));

  switch (prop_id)
  {
    case PROP_DIST:
      g_value_set_object (value, bao_dv->dist);
      break;
    case PROP_Z:
      g_value_set_object (value, bao_dv->x);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_bao_dv_dispose (GObject *object)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);

  nc_distance_clear (&bao_dv->dist);
  ncm_vector_clear (&bao_dv->x);

  /* Chain up : end */
  G_OBJECT_CLASS (nc_data_bao_dv_parent_class)->dispose (object);
}

static void
nc_data_bao_dv_finalize (GObject *object)
{
  /* Chain up : start */
  G_OBJECT_CLASS (nc_data_bao_dv_parent_class)->finalize (object);
}

static void _nc_data_bao_dv_prepare (NcmData *data, NcmMSet *mset);
static void _nc_data_bao_dv_mean_func (NcmDataGaussDiag *diag, NcmMSet *mset, NcmVector *vp);
static void _nc_data_bao_dv_set_size (NcmDataGaussDiag *diag, guint np);

static void
nc_data_bao_dv_class_init (NcDataBaoDVClass *klass)
{
  GObjectClass *object_class        = G_OBJECT_CLASS (klass);
  NcmDataClass *data_class          = NCM_DATA_CLASS (klass);
  NcmDataGaussDiagClass *diag_class = NCM_DATA_GAUSS_DIAG_CLASS (klass);

  object_class->set_property = &nc_data_bao_dv_set_property;
  object_class->get_property = &nc_data_bao_dv_get_property;
  object_class->dispose      = &nc_data_bao_dv_dispose;
  object_class->finalize     = &nc_data_bao_dv_finalize;

  g_object_class_install_property (object_class,
                                   PROP_DIST,
                                   g_param_spec_object ("dist",
                                                        NULL,
                                                        "Distance object",
                                                        NC_TYPE_DISTANCE,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));

  g_object_class_install_property (object_class,
                                   PROP_Z,
                                   g_param_spec_object ("z",
                                                        NULL,
                                                        "Data redshift",
                                                        NCM_TYPE_VECTOR,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));

  data_class->prepare   = &_nc_data_bao_dv_prepare;
  diag_class->mean_func = &_nc_data_bao_dv_mean_func;
  diag_class->set_size  = &_nc_data_bao_dv_set_size;
}

static void
_nc_data_bao_dv_prepare (NcmData *data, NcmMSet *mset)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (data);
  NcHICosmo *cosmo    = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));

  nc_distance_prepare_if_needed (bao_dv->dist, cosmo);
}

static void
_nc_data_bao_dv_mean_func (NcmDataGaussDiag *diag, NcmMSet *mset, NcmVector *vp)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (diag);
  NcHICosmo *cosmo    = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));
  const guint np      = ncm_data_gauss_diag_get_size (diag);
  guint i;

  for (i = 0; i < np; i++)
  {
    const gdouble z  = ncm_vector_get (bao_dv->x, i);
    const gdouble DV = nc_distance_dilation_scale (bao_dv->dist, cosmo, z) * nc_hicosmo_RH_Mpc (cosmo);

    ncm_vector_set (vp, i, DV);
  }
}

static void
_nc_data_bao_dv_set_size (NcmDataGaussDiag *diag, guint np)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (diag);
  const guint cnp     = ncm_data_gauss_diag_get_size (diag);

  if ((np == 0) || (np != cnp))
    ncm_vector_clear (&bao_dv->x);

  if ((np != 0) && (np != cnp))
    bao_dv->x = ncm_vector_new (np);

  /* Chain up : end */
  NCM_DATA_GAUSS_DIAG_CLASS (nc_data_bao_dv_parent_class)->set_size (diag, np);
}

/**
 * nc_data_bao_dv_new_from_file:
 * @filename: file containing a serialized #NcDataBaoDV.
 *
 * Creates a new #NcDataBaoDV from @filename.
 *
 * Returns: (transfer full): the newly created #NcDataBaoDV.
 */
NcDataBaoDV *
nc_data_bao_dv_new_from_file (const gchar *filename)
{
  NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (ncm_serialize_global_from_file (filename));

  g_assert (NC_IS_DATA_BAO_DV (bao_dv));

  return bao_dv;
}

/**
 * nc_data_bao_dv_new_from_id:
 * @dist: a #NcDistance
 * @id: a #NcDataBaoId
 *
 * Creates a new acustic scale data object #NcDataBaoDV from @id.
 * This object requires a #NcDistance object to be set.
 *
 * Returns: a #NcDataBaoDV
 */
NcDataBaoDV *
nc_data_bao_dv_new_from_id (NcDistance *dist, NcDataBaoId id)
{
  NcDataBaoDV *bao_dv;
  gchar *filename;

  switch (id)
  {
    case NC_DATA_BAO_DV_EISENSTEIN2005:
      filename = ncm_cfg_get_data_filename ("nc_data_bao_dv_eisenstein2005.obj", TRUE);
      break;
    default:
      g_error ("nc_data_bao_dv_new_from_id: id %d not recognized.", id);
      break;
  }

  bao_dv = nc_data_bao_dv_new_from_file (filename);
  nc_data_bao_dv_set_dist (bao_dv, dist);
  g_free (filename);

  return bao_dv;
}

/**
 * nc_data_bao_dv_set_dist:
 * @bao_dv: a #NcDataBaoDV
 * @dist: a #NcDistance
 *
 * Sets the distance object.
 *
 */
void
nc_data_bao_dv_set_dist (NcDataBaoDV *bao_dv, NcDistance *dist)
{
  nc_distance_clear (&bao_dv->dist);
  bao_dv->dist = nc_distance_ref (dist);
}


1	/***************************************************************************
2	* nc_data_bao_dv.c
3	*
4	* Thu Apr 22 15:31:32 2010
5	* Copyright 2010 Sandro Dias Pinto Vitenti
6	* <vitenti@uel.br>
7	****************************************************************************/
8	/*
9	* numcosmo
10	* Copyright (C) 2012 Sandro Dias Pinto Vitenti <vitenti@uel.br>
11	*
12	* numcosmo is free software: you can redistribute it and/or modify it
13	* under the terms of the GNU General Public License as published by the
14	* Free Software Foundation, either version 3 of the License, or
15	* (at your option) any later version.
16	*
17	* numcosmo is distributed in the hope that it will be useful, but
18	* WITHOUT ANY WARRANTY; without even the implied warranty of
19	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
20	* See the GNU General Public License for more details.
21	*
22	* You should have received a copy of the GNU General Public License along
23	* with this program. If not, see <http://www.gnu.org/licenses/>.
24	*/
25
26	/**
27	* SECTION:nc_data_bao_dv
28	* @title: NcDataBaoDV
29	* @short_description: Baryon oscillation data -- volume mean $D_V$.
30	*
31	* See [Eisenstein et al. (2005)][XEisenstein2005].
32	*
33	* The data is stored in a #NcDataBaoDV object. The data is stored in a
34	* #NcDataGaussDiag base class object, which is a subclass of #NcmData.
35	* The data represents the mean values of the volume mean $D_V$ at the redshift
36	* $z$.
37	*
38	*/
39
40	#ifdef HAVE_CONFIG_H
41	# include "config.h"
42	#endif /* HAVE_CONFIG_H */
43	#include "build_cfg.h"
44
45	#include "data/nc_data_bao_dv.h"
46	#include "nc_enum_types.h"
47	#include "math/ncm_cfg.h"
48
49	enum
50	{
51	PROP_0,
52	PROP_DIST,
53	PROP_Z,
54	PROP_SIZE,
55	};
56
57	G_DEFINE_TYPE (NcDataBaoDV, nc_data_bao_dv, NCM_TYPE_DATA_GAUSS_DIAG);	64!
58
59	static void
60	nc_data_bao_dv_init (NcDataBaoDV *bao_dv)	×
61	{
62	bao_dv->dist = nc_distance_new (2.0);	×
63	bao_dv->x = NULL;	×
64	}	×
65
66	static void
67	nc_data_bao_dv_set_property (GObject object, guint prop_id, const GValue value, GParamSpec *pspec)	×
68	{
69	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);	×
70
UNCOV 71	g_return_if_fail (NC_IS_DATA_BAO_DV (object));	×
72
73	switch (prop_id)	×
74	{
75	case PROP_DIST:	×
76	nc_data_bao_dv_set_dist (bao_dv, g_value_get_object (value));	×
77	break;	×
78	case PROP_Z:	×
79	ncm_vector_substitute (&bao_dv->x, g_value_get_object (value), TRUE);	×
80	break;	×
81	default:	×
82	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
83	break;	×
84	}
85	}
86
87	static void
88	nc_data_bao_dv_get_property (GObject object, guint prop_id, GValue value, GParamSpec *pspec)	×
89	{
90	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);	×
91
UNCOV 92	g_return_if_fail (NC_IS_DATA_BAO_DV (object));	×
93
94	switch (prop_id)	×
95	{
96	case PROP_DIST:	×
97	g_value_set_object (value, bao_dv->dist);	×
98	break;	×
99	case PROP_Z:	×
100	g_value_set_object (value, bao_dv->x);	×
101	break;	×
102	default:	×
103	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
104	break;	×
105	}
106	}
107
108	static void
109	nc_data_bao_dv_dispose (GObject *object)	×
110	{
111	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (object);	×
112
113	nc_distance_clear (&bao_dv->dist);	×
114	ncm_vector_clear (&bao_dv->x);	×
115
116	/* Chain up : end */
117	G_OBJECT_CLASS (nc_data_bao_dv_parent_class)->dispose (object);	×
118	}	×
119
120	static void
UNCOV 121	nc_data_bao_dv_finalize (GObject *object)	×
122	{
123	/* Chain up : start */
124	G_OBJECT_CLASS (nc_data_bao_dv_parent_class)->finalize (object);	×
125	}	×
126
127	static void _nc_data_bao_dv_prepare (NcmData data, NcmMSet mset);
128	static void _nc_data_bao_dv_mean_func (NcmDataGaussDiag diag, NcmMSet mset, NcmVector *vp);
129	static void _nc_data_bao_dv_set_size (NcmDataGaussDiag *diag, guint np);
130
131	static void
132	nc_data_bao_dv_class_init (NcDataBaoDVClass *klass)	64✔
133	{
134	GObjectClass *object_class = G_OBJECT_CLASS (klass);	64✔
135	NcmDataClass *data_class = NCM_DATA_CLASS (klass);	64✔
136	NcmDataGaussDiagClass *diag_class = NCM_DATA_GAUSS_DIAG_CLASS (klass);	64✔
137
138	object_class->set_property = &nc_data_bao_dv_set_property;	64✔
139	object_class->get_property = &nc_data_bao_dv_get_property;	64✔
140	object_class->dispose = &nc_data_bao_dv_dispose;	64✔
141	object_class->finalize = &nc_data_bao_dv_finalize;	64✔
142
143	g_object_class_install_property (object_class,	64✔
144	PROP_DIST,
145	g_param_spec_object ("dist",
146	NULL,
147	"Distance object",
148	NC_TYPE_DISTANCE,
149	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
150
151	g_object_class_install_property (object_class,	64✔
152	PROP_Z,
153	g_param_spec_object ("z",
154	NULL,
155	"Data redshift",
156	NCM_TYPE_VECTOR,
157	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
158
159	data_class->prepare = &_nc_data_bao_dv_prepare;	64✔
160	diag_class->mean_func = &_nc_data_bao_dv_mean_func;	64✔
161	diag_class->set_size = &_nc_data_bao_dv_set_size;	64✔
162	}	64✔
163
164	static void
165	_nc_data_bao_dv_prepare (NcmData data, NcmMSet mset)	×
166	{
167	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (data);	×
NEW 168	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
169
170	nc_distance_prepare_if_needed (bao_dv->dist, cosmo);	×
171	}	×
172
173	static void
174	_nc_data_bao_dv_mean_func (NcmDataGaussDiag diag, NcmMSet mset, NcmVector *vp)	×
175	{
176	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (diag);	×
177	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
178	const guint np = ncm_data_gauss_diag_get_size (diag);	×
179	guint i;
180
181	for (i = 0; i < np; i++)	×
182	{
NEW 183	const gdouble z = ncm_vector_get (bao_dv->x, i);	×
184	const gdouble DV = nc_distance_dilation_scale (bao_dv->dist, cosmo, z) * nc_hicosmo_RH_Mpc (cosmo);	×
185
UNCOV 186	ncm_vector_set (vp, i, DV);	×
187	}
188	}	×
189
190	static void
191	_nc_data_bao_dv_set_size (NcmDataGaussDiag *diag, guint np)	×
192	{
193	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (diag);	×
194	const guint cnp = ncm_data_gauss_diag_get_size (diag);	×
195
196	if ((np == 0) \|\| (np != cnp))	×
197	ncm_vector_clear (&bao_dv->x);	×
198
199	if ((np != 0) && (np != cnp))	×
200	bao_dv->x = ncm_vector_new (np);	×
201
202	/* Chain up : end */
203	NCM_DATA_GAUSS_DIAG_CLASS (nc_data_bao_dv_parent_class)->set_size (diag, np);	×
204	}	×
205
206	/**
207	* nc_data_bao_dv_new_from_file:
208	* @filename: file containing a serialized #NcDataBaoDV.
209	*
210	* Creates a new #NcDataBaoDV from @filename.
211	*
212	* Returns: (transfer full): the newly created #NcDataBaoDV.
213	*/
214	NcDataBaoDV *
215	nc_data_bao_dv_new_from_file (const gchar *filename)	×
216	{
217	NcDataBaoDV *bao_dv = NC_DATA_BAO_DV (ncm_serialize_global_from_file (filename));	×
218
UNCOV 219	g_assert (NC_IS_DATA_BAO_DV (bao_dv));	×
220
221	return bao_dv;	×
222	}
223
224	/**
225	* nc_data_bao_dv_new_from_id:
226	* @dist: a #NcDistance
227	* @id: a #NcDataBaoId
228	*
229	* Creates a new acustic scale data object #NcDataBaoDV from @id.
230	* This object requires a #NcDistance object to be set.
231	*
232	* Returns: a #NcDataBaoDV
233	*/
234	NcDataBaoDV *
235	nc_data_bao_dv_new_from_id (NcDistance *dist, NcDataBaoId id)	×
236	{
237	NcDataBaoDV *bao_dv;
238	gchar *filename;
239
240	switch (id)	×
241	{
242	case NC_DATA_BAO_DV_EISENSTEIN2005:	×
243	filename = ncm_cfg_get_data_filename ("nc_data_bao_dv_eisenstein2005.obj", TRUE);	×
244	break;	×
245	default:	×
246	g_error ("nc_data_bao_dv_new_from_id: id %d not recognized.", id);
247	break;
248	}
249
250	bao_dv = nc_data_bao_dv_new_from_file (filename);	×
251	nc_data_bao_dv_set_dist (bao_dv, dist);	×
252	g_free (filename);	×
253
254	return bao_dv;	×
255	}
256
257	/**
258	* nc_data_bao_dv_set_dist:
259	* @bao_dv: a #NcDataBaoDV
260	* @dist: a #NcDistance
261	*
262	* Sets the distance object.
263	*
264	*/
265	void
266	nc_data_bao_dv_set_dist (NcDataBaoDV bao_dv, NcDistance dist)	×
267	{
268	nc_distance_clear (&bao_dv->dist);	×
269	bao_dv->dist = nc_distance_ref (dist);	×
270	}	×
271

NumCosmo / NumCosmo / 6749245209

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous