7062620033

Committed 01 Dec 2023 04:31PM UTC coverage: 31.848% (-21.6%) from 53.447%

Build # 7062620033

Build Type

Pull #120

github

Committed by

web-flow

Commit Message

Merge 411ae5c9b into 24e5f6a9b

Pull Request Pull Request #120: Moving to meson

Run Details

14336 of 71092 branches covered (0.0%)

Branch coverage included in aggregate %.

1956 of 3488 new or added lines in 312 files covered. (56.08%)

3938 existing lines in 209 files now uncovered.

57528 of 154555 relevant lines covered (37.22%)

2504997.96 hits per line

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

12.84

/numcosmo/data/nc_data_cluster_wl.c

/***************************************************************************
 *            nc_data_cluster_wl.c
 *
 *  Mon Jul 27 16:10:25 2020
 *  Copyright  2020  Mariana Penna Lima
 *  <pennalima@gmail.com>
 ****************************************************************************/
/*
 * nc_data_cluster_wl.c
 * Copyright (C) 2020 Mariana Penna Lima <pennalima@gmail.com>
 *
 * 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_cluster_wl
 * @title: NcDataClusterWL
 * @short_description: Cluster weak lensing likelihood.
 * @stability: Unstable
 *
 * FIXME
 *
 */

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

#include "data/nc_data_cluster_wl.h"

#include "nc_hicosmo.h"
#include "lss/nc_halo_density_profile.h"
#include "lss/nc_wl_surface_mass_density.h"
#include "lss/nc_galaxy_wl.h"
#include "math/ncm_cfg.h"
#include "math/ncm_util.h"
#include "math/ncm_integrate.h"
#include "math/ncm_memory_pool.h"

struct _NcDataClusterWLPrivate
{
  NcmObjArray *galaxy_array;
  gdouble psf_size;
  gdouble z_cluster;
  gdouble ra_cluster;
  gdouble dec_cluster;
};

enum
{
  PROP_0,
  PROP_GALAXY_ARRAY,
  PROP_PSF_SIZE,
  PROP_Z_CLUSTER,
  PROP_RA_CLUSTER,
  PROP_DEC_CLUSTER,
  PROP_SIZE,
};

G_DEFINE_TYPE_WITH_PRIVATE (NcDataClusterWL, nc_data_cluster_wl, NCM_TYPE_DATA)

static void
nc_data_cluster_wl_init (NcDataClusterWL *dcwl)
{
  NcDataClusterWLPrivate * const self = dcwl->priv = nc_data_cluster_wl_get_instance_private (dcwl);

  self->galaxy_array = ncm_obj_array_new ();
  self->z_cluster    = 0.0;
  self->ra_cluster   = 0.0;
  self->dec_cluster  = 0.0;
  self->psf_size     = 0.0;
}

static void
nc_data_cluster_wl_set_property (GObject *object, guint prop_id, const GValue *value, GParamSpec *pspec)
{
  NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (object);
  NcDataClusterWLPrivate * const self = dcwl->priv;

  g_return_if_fail (NC_IS_DATA_CLUSTER_WL (object));

  switch (prop_id)
  {
    case PROP_GALAXY_ARRAY:
    {
      NcmObjArray *galaxy_array = g_value_get_boxed (value);

      g_clear_pointer (&self->galaxy_array, ncm_obj_array_unref);
      self->galaxy_array = ncm_obj_array_ref (galaxy_array);

      break;
    }
    case PROP_PSF_SIZE:
      self->psf_size = g_value_get_double (value);
      break;
    case PROP_Z_CLUSTER:
      self->z_cluster = g_value_get_double (value);
      break;
    case PROP_RA_CLUSTER:
      self->ra_cluster = g_value_get_double (value);
      break;
    case PROP_DEC_CLUSTER:
      self->dec_cluster = g_value_get_double (value);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_cluster_wl_get_property (GObject *object, guint prop_id, GValue *value, GParamSpec *pspec)
{
  NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (object);
  NcDataClusterWLPrivate * const self = dcwl->priv;

  g_return_if_fail (NC_IS_DATA_CLUSTER_WL (object));

  switch (prop_id)
  {
    case PROP_GALAXY_ARRAY:
      g_value_set_boxed (value, self->galaxy_array);
      break;
    case PROP_PSF_SIZE:
      g_value_set_double (value, self->psf_size);
      break;
    case PROP_Z_CLUSTER:
      g_value_set_double (value, self->z_cluster);
      break;
    case PROP_RA_CLUSTER:
      g_value_set_double (value, self->ra_cluster);
      break;
    case PROP_DEC_CLUSTER:
      g_value_set_double (value, self->dec_cluster);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_cluster_wl_dispose (GObject *object)
{
  NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (object);
  NcDataClusterWLPrivate * const self = dcwl->priv;

  ncm_obj_array_clear (&self->galaxy_array);

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

static void
nc_data_cluster_wl_finalize (GObject *object)
{
  /*NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (object);*/

  /* Chain up : end */
  G_OBJECT_CLASS (nc_data_cluster_wl_parent_class)->finalize (object);
}

static void _nc_data_cluster_wl_m2lnL_val (NcmData *data, NcmMSet *mset, gdouble *m2lnL);
static guint _nc_data_cluster_wl_get_len (NcmData *data);
static void _nc_data_cluster_wl_prepare (NcmData *data, NcmMSet *mset);

static void
nc_data_cluster_wl_class_init (NcDataClusterWLClass *klass)
{
  GObjectClass *object_class = G_OBJECT_CLASS (klass);
  NcmDataClass *data_class   = NCM_DATA_CLASS (klass);

  object_class->set_property = nc_data_cluster_wl_set_property;
  object_class->get_property = nc_data_cluster_wl_get_property;
  object_class->dispose      = nc_data_cluster_wl_dispose;
  object_class->finalize     = nc_data_cluster_wl_finalize;

  g_object_class_install_property (object_class,
                                   PROP_GALAXY_ARRAY,
                                   g_param_spec_boxed ("galaxy-array",
                                                       NULL,
                                                       "Array of galaxy weak lensing objects",
                                                       NCM_TYPE_OBJ_ARRAY,
                                                       G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));
  g_object_class_install_property (object_class,
                                   PROP_PSF_SIZE,
                                   g_param_spec_double ("psf-size",
                                                        NULL,
                                                        "PSF size",
                                                        0.0, G_MAXDOUBLE, 0.0,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));

  g_object_class_install_property (object_class,
                                   PROP_Z_CLUSTER,
                                   g_param_spec_double ("z-cluster",
                                                        NULL,
                                                        "Cluster (halo) redshift",
                                                        0.0, G_MAXDOUBLE, 0.0,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));
  g_object_class_install_property (object_class,
                                   PROP_RA_CLUSTER,
                                   g_param_spec_double ("ra-cluster",
                                                        NULL,
                                                        "Cluster (halo) RA",
                                                        -360.0, 360.0, 0.0,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));
  g_object_class_install_property (object_class,
                                   PROP_DEC_CLUSTER,
                                   g_param_spec_double ("dec-cluster",
                                                        NULL,
                                                        "Cluster (halo) DEC",
                                                        -360.0, 360.0, 0.0,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));

  data_class->m2lnL_val  = &_nc_data_cluster_wl_m2lnL_val;
  data_class->get_length = &_nc_data_cluster_wl_get_len;
  data_class->prepare    = &_nc_data_cluster_wl_prepare;
}

static void
_nc_data_cluster_wl_m2lnL_val (NcmData *data, NcmMSet *mset, gdouble *m2lnL)
{
  NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (data);
  NcDataClusterWLPrivate * const self = dcwl->priv;
  NcHICosmo *cosmo                    = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));
  NcWLSurfaceMassDensity *smd         = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));
  NcHaloDensityProfile *dp            = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));
  const guint ngal                    = self->galaxy_array->len;
  guint i;

  m2lnL[0] = 0.0;

  for (i = 0; i < ngal; i++)
  {
    NcGalaxyWL *gwl_i = NC_GALAXY_WL (ncm_obj_array_peek (self->galaxy_array, i));

    m2lnL[0] += nc_galaxy_wl_eval_m2lnP (gwl_i, cosmo, dp, smd, self->z_cluster);
  }

  return;
}

static guint
_nc_data_cluster_wl_get_len (NcmData *data)
{
  NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (data);
  NcDataClusterWLPrivate * const self = dcwl->priv;

  if (self->galaxy_array != NULL)
  {
    const guint ngal = self->galaxy_array->len;
    guint len        = 0;
    guint i;

    for (i = 0; i < ngal; i++)
    {
      NcGalaxyWL *gwl_i = NC_GALAXY_WL (ncm_obj_array_peek (self->galaxy_array, i));

      len += nc_galaxy_wl_len (gwl_i);
    }

    return len;
  }
  else
  {
    return 0;
  }
}

static void
_nc_data_cluster_wl_prepare (NcmData *data, NcmMSet *mset)
{
  /*NcDataClusterWL *dcwl               = NC_DATA_CLUSTER_WL (data);*/
  /*NcDataClusterWLPrivate * const self = dcwl->priv;*/
  NcHICosmo *cosmo            = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));
  NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));
  NcHaloDensityProfile *dp    = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));

  g_assert ((cosmo != NULL) && (smd != NULL) && (dp != NULL));

  nc_wl_surface_mass_density_prepare_if_needed (smd, cosmo);
}

/**
 * nc_data_cluster_wl_new:
 *
 * Creates a new #NcDataClusterWL.
 *
 * Returns: the newly created #NcDataClusterWL.
 */
NcDataClusterWL *
nc_data_cluster_wl_new (void)
{
  NcDataClusterWL *dcwl = g_object_new (NC_TYPE_DATA_CLUSTER_WL,
                                        NULL);

  return dcwl;
}

/**
 * nc_data_cluster_wl_new_from_file:
 * @filename: file containing a serialized #NcDataClusterWL
 *
 * Creates a new #NcDataClusterWL from @filename.
 *
 * Returns: (transfer full): the newly created #NcDataClusterWL.
 */
NcDataClusterWL *
nc_data_cluster_wl_new_from_file (const gchar *filename)
{
  NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (ncm_serialize_global_from_file (filename));

  g_assert (NC_IS_DATA_CLUSTER_WL (dcwl));

  return dcwl;
}

/**
 * nc_data_cluster_wl_ref:
 * @dcwl: a #NcDataClusterWL
 *
 * Increases the reference count of @dcwl by one.
 *
 * Returns: (transfer full): @dcwl
 */
NcDataClusterWL *
nc_data_cluster_wl_ref (NcDataClusterWL *dcwl)
{
  return g_object_ref (dcwl);
}

/**
 * nc_data_cluster_wl_free:
 * @dcwl: a #NcDataClusterWL
 *
 * Atomically decrements the reference count of @dcwl by one. If the reference count drops to 0,
 * all memory allocated by @dcwl is released.
 *
 */
void
nc_data_cluster_wl_free (NcDataClusterWL *dcwl)
{
  g_object_unref (dcwl);
}

/**
 * nc_data_cluster_wl_clear:
 * @dcwl: a #NcDataClusterWL
 *
 * The reference count of @dcwl is decreased and the pointer is set to NULL.
 *
 */
void
nc_data_cluster_wl_clear (NcDataClusterWL **dcwl)
{
  g_clear_object (dcwl);
}


1	/***************************************************************************
2	* nc_data_cluster_wl.c
3	*
4	* Mon Jul 27 16:10:25 2020
5	* Copyright 2020 Mariana Penna Lima
6	* <pennalima@gmail.com>
7	****************************************************************************/
8	/*
9	* nc_data_cluster_wl.c
10	* Copyright (C) 2020 Mariana Penna Lima <pennalima@gmail.com>
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_cluster_wl
28	* @title: NcDataClusterWL
29	* @short_description: Cluster weak lensing likelihood.
30	* @stability: Unstable
31	*
32	* FIXME
33	*
34	*/
35
36	#ifdef HAVE_CONFIG_H
37	# include "config.h"
38	#endif /* HAVE_CONFIG_H */
39	#include "build_cfg.h"
40
41	#include "data/nc_data_cluster_wl.h"
42
43	#include "nc_hicosmo.h"
44	#include "lss/nc_halo_density_profile.h"
45	#include "lss/nc_wl_surface_mass_density.h"
46	#include "lss/nc_galaxy_wl.h"
47	#include "math/ncm_cfg.h"
48	#include "math/ncm_util.h"
49	#include "math/ncm_integrate.h"
50	#include "math/ncm_memory_pool.h"
51
52	struct _NcDataClusterWLPrivate
53	{
54	NcmObjArray *galaxy_array;
55	gdouble psf_size;
56	gdouble z_cluster;
57	gdouble ra_cluster;
58	gdouble dec_cluster;
59	};
60
61	enum
62	{
63	PROP_0,
64	PROP_GALAXY_ARRAY,
65	PROP_PSF_SIZE,
66	PROP_Z_CLUSTER,
67	PROP_RA_CLUSTER,
68	PROP_DEC_CLUSTER,
69	PROP_SIZE,
70	};
71
72	G_DEFINE_TYPE_WITH_PRIVATE (NcDataClusterWL, nc_data_cluster_wl, NCM_TYPE_DATA)	333!
73
74	static void
75	nc_data_cluster_wl_init (NcDataClusterWL *dcwl)	×
76	{
77	NcDataClusterWLPrivate * const self = dcwl->priv = nc_data_cluster_wl_get_instance_private (dcwl);	×
78
79	self->galaxy_array = ncm_obj_array_new ();	×
80	self->z_cluster = 0.0;	×
81	self->ra_cluster = 0.0;	×
82	self->dec_cluster = 0.0;	×
83	self->psf_size = 0.0;	×
84	}	×
85
86	static void
87	nc_data_cluster_wl_set_property (GObject object, guint prop_id, const GValue value, GParamSpec *pspec)	×
88	{
89	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (object);	×
90	NcDataClusterWLPrivate * const self = dcwl->priv;	×
91
92	g_return_if_fail (NC_IS_DATA_CLUSTER_WL (object));	×
93
94	switch (prop_id)	×
95	{
96	case PROP_GALAXY_ARRAY:	×
97	{
98	NcmObjArray *galaxy_array = g_value_get_boxed (value);	×
99
100	g_clear_pointer (&self->galaxy_array, ncm_obj_array_unref);	×
101	self->galaxy_array = ncm_obj_array_ref (galaxy_array);	×
102
103	break;	×
104	}
105	case PROP_PSF_SIZE:	×
106	self->psf_size = g_value_get_double (value);	×
107	break;	×
108	case PROP_Z_CLUSTER:	×
109	self->z_cluster = g_value_get_double (value);	×
110	break;	×
111	case PROP_RA_CLUSTER:	×
112	self->ra_cluster = g_value_get_double (value);	×
113	break;	×
114	case PROP_DEC_CLUSTER:	×
115	self->dec_cluster = g_value_get_double (value);	×
116	break;	×
117	default:	×
118	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
119	break;	×
120	}
121	}
122
123	static void
124	nc_data_cluster_wl_get_property (GObject object, guint prop_id, GValue value, GParamSpec *pspec)	×
125	{
126	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (object);	×
127	NcDataClusterWLPrivate * const self = dcwl->priv;	×
128
129	g_return_if_fail (NC_IS_DATA_CLUSTER_WL (object));	×
130
131	switch (prop_id)	×
132	{
133	case PROP_GALAXY_ARRAY:	×
134	g_value_set_boxed (value, self->galaxy_array);	×
135	break;	×
136	case PROP_PSF_SIZE:	×
137	g_value_set_double (value, self->psf_size);	×
138	break;	×
139	case PROP_Z_CLUSTER:	×
140	g_value_set_double (value, self->z_cluster);	×
141	break;	×
142	case PROP_RA_CLUSTER:	×
143	g_value_set_double (value, self->ra_cluster);	×
144	break;	×
145	case PROP_DEC_CLUSTER:	×
146	g_value_set_double (value, self->dec_cluster);	×
147	break;	×
148	default:	×
149	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
150	break;	×
151	}
152	}
153
154	static void
155	nc_data_cluster_wl_dispose (GObject *object)	×
156	{
157	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (object);	×
158	NcDataClusterWLPrivate * const self = dcwl->priv;	×
159
160	ncm_obj_array_clear (&self->galaxy_array);	×
161
162	/* Chain up : end */
163	G_OBJECT_CLASS (nc_data_cluster_wl_parent_class)->dispose (object);	×
164	}	×
165
166	static void
167	nc_data_cluster_wl_finalize (GObject *object)	×
168	{
169	/NcDataClusterWL dcwl = NC_DATA_CLUSTER_WL (object);*/
170
171	/* Chain up : end */
172	G_OBJECT_CLASS (nc_data_cluster_wl_parent_class)->finalize (object);	×
173	}	×
174
175	static void _nc_data_cluster_wl_m2lnL_val (NcmData data, NcmMSet mset, gdouble *m2lnL);
176	static guint _nc_data_cluster_wl_get_len (NcmData *data);
177	static void _nc_data_cluster_wl_prepare (NcmData data, NcmMSet mset);
178
179	static void
180	nc_data_cluster_wl_class_init (NcDataClusterWLClass *klass)	111✔
181	{
182	GObjectClass *object_class = G_OBJECT_CLASS (klass);	111✔
183	NcmDataClass *data_class = NCM_DATA_CLASS (klass);	111✔
184
185	object_class->set_property = nc_data_cluster_wl_set_property;	111✔
186	object_class->get_property = nc_data_cluster_wl_get_property;	111✔
187	object_class->dispose = nc_data_cluster_wl_dispose;	111✔
188	object_class->finalize = nc_data_cluster_wl_finalize;	111✔
189
190	g_object_class_install_property (object_class,	111✔
191	PROP_GALAXY_ARRAY,
192	g_param_spec_boxed ("galaxy-array",
193	NULL,
194	"Array of galaxy weak lensing objects",
195	NCM_TYPE_OBJ_ARRAY,
196	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
197	g_object_class_install_property (object_class,	111✔
198	PROP_PSF_SIZE,
199	g_param_spec_double ("psf-size",
200	NULL,
201	"PSF size",
202	0.0, G_MAXDOUBLE, 0.0,
203	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
204
205	g_object_class_install_property (object_class,	111✔
206	PROP_Z_CLUSTER,
207	g_param_spec_double ("z-cluster",
208	NULL,
209	"Cluster (halo) redshift",
210	0.0, G_MAXDOUBLE, 0.0,
211	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
212	g_object_class_install_property (object_class,	111✔
213	PROP_RA_CLUSTER,
214	g_param_spec_double ("ra-cluster",
215	NULL,
216	"Cluster (halo) RA",
217	-360.0, 360.0, 0.0,
218	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
219	g_object_class_install_property (object_class,	111✔
220	PROP_DEC_CLUSTER,
221	g_param_spec_double ("dec-cluster",
222	NULL,
223	"Cluster (halo) DEC",
224	-360.0, 360.0, 0.0,
225	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
226
227	data_class->m2lnL_val = &_nc_data_cluster_wl_m2lnL_val;	111✔
228	data_class->get_length = &_nc_data_cluster_wl_get_len;	111✔
229	data_class->prepare = &_nc_data_cluster_wl_prepare;	111✔
230	}	111✔
231
232	static void
233	_nc_data_cluster_wl_m2lnL_val (NcmData data, NcmMSet mset, gdouble *m2lnL)	×
234	{
235	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (data);	×
236	NcDataClusterWLPrivate * const self = dcwl->priv;	×
237	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
238	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	×
239	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	×
240	const guint ngal = self->galaxy_array->len;	×
NEW 241	guint i;	×
242
UNCOV 243	m2lnL[0] = 0.0;	×
244
245	for (i = 0; i < ngal; i++)	×
246	{
247	NcGalaxyWL *gwl_i = NC_GALAXY_WL (ncm_obj_array_peek (self->galaxy_array, i));	×
248
249	m2lnL[0] += nc_galaxy_wl_eval_m2lnP (gwl_i, cosmo, dp, smd, self->z_cluster);	×
250	}
251
252	return;	×
253	}
254
255	static guint
256	_nc_data_cluster_wl_get_len (NcmData *data)	×
257	{
258	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (data);	×
259	NcDataClusterWLPrivate * const self = dcwl->priv;	×
260
261	if (self->galaxy_array != NULL)	×
262	{
263	const guint ngal = self->galaxy_array->len;	×
264	guint len = 0;	×
NEW 265	guint i;	×
266
UNCOV 267	for (i = 0; i < ngal; i++)	×
268	{
269	NcGalaxyWL *gwl_i = NC_GALAXY_WL (ncm_obj_array_peek (self->galaxy_array, i));	×
270
271	len += nc_galaxy_wl_len (gwl_i);	×
272	}
273
274	return len;	×
275	}
276	else
277	{
278	return 0;
279	}
280	}
281
282	static void
283	_nc_data_cluster_wl_prepare (NcmData data, NcmMSet mset)	×
284	{
285	/NcDataClusterWL dcwl = NC_DATA_CLUSTER_WL (data);*/
286	/NcDataClusterWLPrivate const self = dcwl->priv;*/
287	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
288	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	×
289	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	×
290
291	g_assert ((cosmo != NULL) && (smd != NULL) && (dp != NULL));	×
292
293	nc_wl_surface_mass_density_prepare_if_needed (smd, cosmo);	×
294	}	×
295
296	/**
297	* nc_data_cluster_wl_new:
298	*
299	* Creates a new #NcDataClusterWL.
300	*
301	* Returns: the newly created #NcDataClusterWL.
302	*/
303	NcDataClusterWL *
304	nc_data_cluster_wl_new (void)	×
305	{
306	NcDataClusterWL *dcwl = g_object_new (NC_TYPE_DATA_CLUSTER_WL,	×
307	NULL);
308
309	return dcwl;	×
310	}
311
312	/**
313	* nc_data_cluster_wl_new_from_file:
314	* @filename: file containing a serialized #NcDataClusterWL
315	*
316	* Creates a new #NcDataClusterWL from @filename.
317	*
318	* Returns: (transfer full): the newly created #NcDataClusterWL.
319	*/
320	NcDataClusterWL *
321	nc_data_cluster_wl_new_from_file (const gchar *filename)	×
322	{
323	NcDataClusterWL *dcwl = NC_DATA_CLUSTER_WL (ncm_serialize_global_from_file (filename));	×
324
325	g_assert (NC_IS_DATA_CLUSTER_WL (dcwl));	×
326
327	return dcwl;	×
328	}
329
330	/**
331	* nc_data_cluster_wl_ref:
332	* @dcwl: a #NcDataClusterWL
333	*
334	* Increases the reference count of @dcwl by one.
335	*
336	* Returns: (transfer full): @dcwl
337	*/
338	NcDataClusterWL *
339	nc_data_cluster_wl_ref (NcDataClusterWL *dcwl)	×
340	{
341	return g_object_ref (dcwl);	×
342	}
343
344	/**
345	* nc_data_cluster_wl_free:
346	* @dcwl: a #NcDataClusterWL
347	*
348	* Atomically decrements the reference count of @dcwl by one. If the reference count drops to 0,
349	* all memory allocated by @dcwl is released.
350	*
351	*/
352	void
353	nc_data_cluster_wl_free (NcDataClusterWL *dcwl)	×
354	{
355	g_object_unref (dcwl);	×
356	}	×
357
358	/**
359	* nc_data_cluster_wl_clear:
360	* @dcwl: a #NcDataClusterWL
361	*
362	* The reference count of @dcwl is decreased and the pointer is set to NULL.
363	*
364	*/
365	void
366	nc_data_cluster_wl_clear (NcDataClusterWL **dcwl)	×
367	{
368	g_clear_object (dcwl);	×
369	}	×
370

NumCosmo / NumCosmo / 7062620033

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