7310595895

Committed 23 Dec 2023 09:32PM UTC coverage: 55.715% (+0.4%) from 55.278%

Build # 7310595895

Build Type

Pull #93

github

Committed by

web-flow

Commit Message

Merge 898c46308 into eb5e3f39f

Pull Request Pull Request #93: Galaxy WL reformulation

Run Details

7420 of 12717 branches covered (0.0%)

Branch coverage included in aggregate %.

1155 of 1436 new or added lines in 21 files covered. (80.43%)

39 existing lines in 4 files now uncovered.

49122 of 88767 relevant lines covered (55.34%)

2014774.86 hits per line

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

73.13

/numcosmo/data/nc_data_cluster_wll.c

/***************************************************************************
 *            nc_data_cluster_wll.c
 *
 *  Tue Jun 15 16:00:13 2023
 *  Copyright  2023  Caio Lima de Oliveira
 *  <caiolimadeoliveira@pm.me>
 ****************************************************************************/
/*
 * nc_data_cluster_wll.c
 * Copyright (C) 2023 Caio Lima de Oliveira <caiolimadeoliveira@pm.me>
 *
 * 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_wll
 * @title: NcDataClusterWLL
 * @short_description: Cluster weak lensing likelihood.
 * @stability: Unstable
 *
 * This class implements the weak lensing likelihood for galaxy clusters.
 *
 */

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

#include "data/nc_data_cluster_wll.h"

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

struct _NcDataClusterWLLPrivate
{
  NcmObjArray *galaxy_array;
  gdouble z_cluster;
  gboolean kde;
};

enum
{
  PROP_0,
  PROP_GALAXY_ARRAY,
  PROP_Z_CLUSTER,
  PROP_KDE,
  PROP_SIZE,
};

G_DEFINE_TYPE_WITH_PRIVATE (NcDataClusterWLL, nc_data_cluster_wll, NCM_TYPE_DATA);

static void
nc_data_cluster_wll_init (NcDataClusterWLL *dcwll)
{
  NcDataClusterWLLPrivate * const self = dcwll->priv = nc_data_cluster_wll_get_instance_private (dcwll);

  self->galaxy_array = ncm_obj_array_new ();
  self->z_cluster    = 0.0;
  self->kde          = FALSE;
}

static void
nc_data_cluster_wll_set_property (GObject *object, guint prop_id, const GValue *value, GParamSpec *pspec)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (object);
  NcDataClusterWLLPrivate * const self = dcwll->priv;

  g_return_if_fail (NC_IS_DATA_CLUSTER_WLL (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_Z_CLUSTER:
      self->z_cluster = g_value_get_double (value);
      break;
    case PROP_KDE:
      nc_data_cluster_wll_set_kde (dcwll, g_value_get_boolean (value));
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_cluster_wll_get_property (GObject *object, guint prop_id, GValue *value, GParamSpec *pspec)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (object);
  NcDataClusterWLLPrivate * const self = dcwll->priv;

  g_return_if_fail (NC_IS_DATA_CLUSTER_WLL (object));

  switch (prop_id)
  {
    case PROP_GALAXY_ARRAY:
      g_value_set_boxed (value, self->galaxy_array);
      break;
    case PROP_Z_CLUSTER:
      g_value_set_double (value, self->z_cluster);
      break;
    case PROP_KDE:
      g_value_set_boolean (value, self->kde);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
nc_data_cluster_wll_dispose (GObject *object)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (object);
  NcDataClusterWLLPrivate * const self = dcwll->priv;

  ncm_obj_array_clear (&self->galaxy_array);

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

static void
nc_data_cluster_wll_finalize (GObject *object)
{
  /*NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);*/

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

static void _nc_data_cluster_wll_m2lnL_val (NcmData *data, NcmMSet *mset, gdouble *m2lnL);
static guint _nc_data_cluster_wll_get_len (NcmData *data);
static void _nc_data_cluster_wll_prepare (NcmData *data, NcmMSet *mset);

static void
nc_data_cluster_wll_class_init (NcDataClusterWLLClass *klass)
{
  GObjectClass *object_class = G_OBJECT_CLASS (klass);
  NcmDataClass *data_class   = NCM_DATA_CLASS (klass);

  object_class->set_property = nc_data_cluster_wll_set_property;
  object_class->get_property = nc_data_cluster_wll_get_property;
  object_class->dispose      = nc_data_cluster_wll_dispose;
  object_class->finalize     = nc_data_cluster_wll_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_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_KDE,
                                   g_param_spec_boolean ("kde",
                                                         NULL,
                                                         "Whether to use KDE method",
                                                         TRUE,
                                                         G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));


  data_class->m2lnL_val  = &_nc_data_cluster_wll_m2lnL_val;
  data_class->get_length = &_nc_data_cluster_wll_get_len;
  data_class->prepare    = &_nc_data_cluster_wll_prepare;
}

static void
_nc_data_cluster_wll_m2lnL_val (NcmData *data, NcmMSet *mset, gdouble *m2lnL)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (data);
  NcDataClusterWLLPrivate * const self = dcwll->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++)
  {
    NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));
    NcmMatrix *gal_obs           = nc_galaxy_wl_likelihood_peek_obs (gwll_i);

    if (self->kde)
      m2lnL[0] += nc_galaxy_wl_likelihood_kde_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster, gal_obs, NULL);
    else
      m2lnL[0] += nc_galaxy_wl_likelihood_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster, gal_obs, NULL);
  }

  return;
}

static guint
_nc_data_cluster_wll_get_len (NcmData *data)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (data);
  NcDataClusterWLLPrivate * const self = dcwll->priv;

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

    for (i = 0; i < ngal; i++)
    {
      NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));

      len += nc_galaxy_wl_likelihood_len (gwll_i);
    }

    return len;
  }
  else
  {
    return 0;
  }
}

static void
_nc_data_cluster_wll_prepare (NcmData *data, NcmMSet *mset)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (data);
  NcDataClusterWLLPrivate * const self = dcwll->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_wll_new:
 *
 * Creates a new #NcDataClusterWLL.
 *
 * Returns: the newly created #NcDataClusterWLL.
 */
NcDataClusterWLL *
nc_data_cluster_wll_new (void)
{
  NcDataClusterWLL *dcwll = g_object_new (NC_TYPE_DATA_CLUSTER_WLL,
                                          NULL);

  return dcwll;
}

/**
 * nc_data_cluster_wll_new_from_file:
 * @filename: file containing a serialized #NcDataClusterWLL
 *
 * Creates a new #NcDataClusterWLL from @filename.
 *
 * Returns: (transfer full): the newly created #NcDataClusterWLL.
 */
NcDataClusterWLL *
nc_data_cluster_wll_new_from_file (const gchar *filename)
{
  NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (ncm_serialize_global_from_file (filename));

  g_assert (NC_IS_DATA_CLUSTER_WLL (dcwll));

  return dcwll;
}

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

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

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

/**
 * nc_data_cluster_wll_set_kde:
 * @dcwll: a #NcDataClusterWLL
 * @kde: whether to use KDE method
 *
 * The reference count of @dcwll is decreased and the pointer is set to NULL.
 *
 */
void
nc_data_cluster_wll_set_kde (NcDataClusterWLL *data, gboolean kde)
{
  NcDataClusterWLL *dcwll              = NC_DATA_CLUSTER_WLL (data);
  NcDataClusterWLLPrivate * const self = dcwll->priv;

  self->kde = kde;
}


1	/***************************************************************************
2	* nc_data_cluster_wll.c
3	*
4	* Tue Jun 15 16:00:13 2023
5	* Copyright 2023 Caio Lima de Oliveira
6	* <caiolimadeoliveira@pm.me>
7	****************************************************************************/
8	/*
9	* nc_data_cluster_wll.c
10	* Copyright (C) 2023 Caio Lima de Oliveira <caiolimadeoliveira@pm.me>
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_wll
28	* @title: NcDataClusterWLL
29	* @short_description: Cluster weak lensing likelihood.
30	* @stability: Unstable
31	*
32	* This class implements the weak lensing likelihood for galaxy clusters.
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_wll.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 "galaxy/nc_galaxy_wl_likelihood.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 _NcDataClusterWLLPrivate
53	{
54	NcmObjArray *galaxy_array;
55	gdouble z_cluster;
56	gboolean kde;
57	};
58
59	enum
60	{
61	PROP_0,
62	PROP_GALAXY_ARRAY,
63	PROP_Z_CLUSTER,
64	PROP_KDE,
65	PROP_SIZE,
66	};
67
68	G_DEFINE_TYPE_WITH_PRIVATE (NcDataClusterWLL, nc_data_cluster_wll, NCM_TYPE_DATA);	698✔
69
70	static void
71	nc_data_cluster_wll_init (NcDataClusterWLL *dcwll)	2✔
72	{
73	NcDataClusterWLLPrivate * const self = dcwll->priv = nc_data_cluster_wll_get_instance_private (dcwll);	2✔
74
75	self->galaxy_array = ncm_obj_array_new ();	2✔
76	self->z_cluster = 0.0;	2✔
77	self->kde = FALSE;	2✔
78	}	2✔
79
80	static void
81	nc_data_cluster_wll_set_property (GObject object, guint prop_id, const GValue value, GParamSpec *pspec)	4✔
82	{
83	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);	4✔
84	NcDataClusterWLLPrivate * const self = dcwll->priv;	4✔
85
86	g_return_if_fail (NC_IS_DATA_CLUSTER_WLL (object));	4✔
87
88	switch (prop_id)	4✔
89	{
90	case PROP_GALAXY_ARRAY:	2✔
91	{
92	NcmObjArray *galaxy_array = g_value_get_boxed (value);	2✔
93
94	g_clear_pointer (&self->galaxy_array, ncm_obj_array_unref);	2✔
95	self->galaxy_array = ncm_obj_array_ref (galaxy_array);	2✔
96
97	break;	2✔
98	}
99	case PROP_Z_CLUSTER:	2✔
100	self->z_cluster = g_value_get_double (value);	2✔
101	break;	2✔
NEW 102	case PROP_KDE:	×
NEW 103	nc_data_cluster_wll_set_kde (dcwll, g_value_get_boolean (value));	×
NEW 104	break;	×
NEW 105	default:	×
NEW 106	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);	×
NEW 107	break;	×
108	}
109	}
110
111	static void
NEW 112	nc_data_cluster_wll_get_property (GObject object, guint prop_id, GValue value, GParamSpec *pspec)	×
113	{
NEW 114	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);	×
NEW 115	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
116
NEW 117	g_return_if_fail (NC_IS_DATA_CLUSTER_WLL (object));	×
118
NEW 119	switch (prop_id)	×
120	{
NEW 121	case PROP_GALAXY_ARRAY:	×
NEW 122	g_value_set_boxed (value, self->galaxy_array);	×
NEW 123	break;	×
NEW 124	case PROP_Z_CLUSTER:	×
NEW 125	g_value_set_double (value, self->z_cluster);	×
NEW 126	break;	×
NEW 127	case PROP_KDE:	×
NEW 128	g_value_set_boolean (value, self->kde);	×
NEW 129	break;	×
NEW 130	default:	×
NEW 131	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);	×
NEW 132	break;	×
133	}
134	}
135
136	static void
137	nc_data_cluster_wll_dispose (GObject *object)	2✔
138	{
139	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);	2✔
140	NcDataClusterWLLPrivate * const self = dcwll->priv;	2✔
141
142	ncm_obj_array_clear (&self->galaxy_array);	2✔
143
144	/* Chain up : end */
145	G_OBJECT_CLASS (nc_data_cluster_wll_parent_class)->dispose (object);	2✔
146	}	2✔
147
148	static void
149	nc_data_cluster_wll_finalize (GObject *object)	2✔
150	{
151	/NcDataClusterWLL dcwll = NC_DATA_CLUSTER_WLL (object);*/
152
153	/* Chain up : end */
154	G_OBJECT_CLASS (nc_data_cluster_wll_parent_class)->finalize (object);	2✔
155	}	2✔
156
157	static void _nc_data_cluster_wll_m2lnL_val (NcmData data, NcmMSet mset, gdouble *m2lnL);
158	static guint _nc_data_cluster_wll_get_len (NcmData *data);
159	static void _nc_data_cluster_wll_prepare (NcmData data, NcmMSet mset);
160
161	static void
162	nc_data_cluster_wll_class_init (NcDataClusterWLLClass *klass)	114✔
163	{
164	GObjectClass *object_class = G_OBJECT_CLASS (klass);	114✔
165	NcmDataClass *data_class = NCM_DATA_CLASS (klass);	114✔
166
167	object_class->set_property = nc_data_cluster_wll_set_property;	114✔
168	object_class->get_property = nc_data_cluster_wll_get_property;	114✔
169	object_class->dispose = nc_data_cluster_wll_dispose;	114✔
170	object_class->finalize = nc_data_cluster_wll_finalize;	114✔
171
172	g_object_class_install_property (object_class,	114✔
173	PROP_GALAXY_ARRAY,
174	g_param_spec_boxed ("galaxy-array",
175	NULL,
176	"Array of galaxy weak lensing objects",
177	NCM_TYPE_OBJ_ARRAY,
178	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
179
180	g_object_class_install_property (object_class,	114✔
181	PROP_Z_CLUSTER,
182	g_param_spec_double ("z-cluster",
183	NULL,
184	"Cluster (halo) redshift",
185	0.0, G_MAXDOUBLE, 0.0,
186	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
187
188	g_object_class_install_property (object_class,	114✔
189	PROP_KDE,
190	g_param_spec_boolean ("kde",
191	NULL,
192	"Whether to use KDE method",
193	TRUE,
194	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
195
196
197	data_class->m2lnL_val = &_nc_data_cluster_wll_m2lnL_val;	114✔
198	data_class->get_length = &_nc_data_cluster_wll_get_len;	114✔
199	data_class->prepare = &_nc_data_cluster_wll_prepare;	114✔
200	}	114✔
201
202	static void
203	_nc_data_cluster_wll_m2lnL_val (NcmData data, NcmMSet mset, gdouble *m2lnL)	343✔
204	{
205	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	343✔
206	NcDataClusterWLLPrivate * const self = dcwll->priv;	343✔
207	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	343✔
208	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	343✔
209	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	343✔
210	const guint ngal = self->galaxy_array->len;	343✔
211	guint i;	343✔
212
213	m2lnL[0] = 0.0;	343✔
214
215	for (i = 0; i < ngal; i++)	686✔
216	{
217	NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));	343✔
218	NcmMatrix *gal_obs = nc_galaxy_wl_likelihood_peek_obs (gwll_i);	343✔
219
220	if (self->kde)	343✔
NEW 221	m2lnL[0] += nc_galaxy_wl_likelihood_kde_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster, gal_obs, NULL);	×
222	else
223	m2lnL[0] += nc_galaxy_wl_likelihood_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster, gal_obs, NULL);	343✔
224	}
225
226	return;	343✔
227	}
228
229	static guint
230	_nc_data_cluster_wll_get_len (NcmData *data)	12✔
231	{
232	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	12✔
233	NcDataClusterWLLPrivate * const self = dcwll->priv;	12✔
234
235	if (self->galaxy_array != NULL)	12✔
236	{
237	const guint ngal = self->galaxy_array->len;	12✔
238	guint len = 0;	12✔
239	guint i;	12✔
240
241	for (i = 0; i < ngal; i++)	24✔
242	{
243	NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));	12✔
244
245	len += nc_galaxy_wl_likelihood_len (gwll_i);	12✔
246	}
247
248	return len;	12✔
249	}
250	else
251	{
252	return 0;
253	}
254	}
255
256	static void
257	_nc_data_cluster_wll_prepare (NcmData data, NcmMSet mset)	343✔
258	{
259	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	343✔
260	NcDataClusterWLLPrivate * const self = dcwll->priv;	343✔
261	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	343✔
262	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	343✔
263	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	343✔
264
265	g_assert ((cosmo != NULL) && (smd != NULL) && (dp != NULL));	343✔
266
267	nc_wl_surface_mass_density_prepare_if_needed (smd, cosmo);	343✔
268	}	343✔
269
270	/**
271	* nc_data_cluster_wll_new:
272	*
273	* Creates a new #NcDataClusterWLL.
274	*
275	* Returns: the newly created #NcDataClusterWLL.
276	*/
277	NcDataClusterWLL *
278	nc_data_cluster_wll_new (void)	2✔
279	{
280	NcDataClusterWLL *dcwll = g_object_new (NC_TYPE_DATA_CLUSTER_WLL,	2✔
281	NULL);
282
283	return dcwll;	2✔
284	}
285
286	/**
287	* nc_data_cluster_wll_new_from_file:
288	* @filename: file containing a serialized #NcDataClusterWLL
289	*
290	* Creates a new #NcDataClusterWLL from @filename.
291	*
292	* Returns: (transfer full): the newly created #NcDataClusterWLL.
293	*/
294	NcDataClusterWLL *
NEW 295	nc_data_cluster_wll_new_from_file (const gchar *filename)	×
296	{
NEW 297	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (ncm_serialize_global_from_file (filename));	×
298
NEW 299	g_assert (NC_IS_DATA_CLUSTER_WLL (dcwll));	×
300
NEW 301	return dcwll;	×
302	}
303
304	/**
305	* nc_data_cluster_wll_ref:
306	* @dcwll: a #NcDataClusterWLL
307	*
308	* Increases the reference count of @dcwll by one.
309	*
310	* Returns: (transfer full): @dcwll
311	*/
312	NcDataClusterWLL *
NEW 313	nc_data_cluster_wll_ref (NcDataClusterWLL *dcwll)	×
314	{
NEW 315	return g_object_ref (dcwll);	×
316	}
317
318	/**
319	* nc_data_cluster_wll_free:
320	* @dcwll: a #NcDataClusterWLL
321	*
322	* Atomically decrements the reference count of @dcwll by one. If the reference count drops to 0,
323	* all memory allocated by @dcwll is released.
324	*
325	*/
326	void
327	nc_data_cluster_wll_free (NcDataClusterWLL *dcwll)	2✔
328	{
329	g_object_unref (dcwll);	2✔
330	}	2✔
331
332	/**
333	* nc_data_cluster_wll_clear:
334	* @dcwll: a #NcDataClusterWLL
335	*
336	* The reference count of @dcwll is decreased and the pointer is set to NULL.
337	*
338	*/
339	void
NEW 340	nc_data_cluster_wll_clear (NcDataClusterWLL **dcwll)	×
341	{
NEW 342	g_clear_object (dcwll);	×
343	}
344
345	/**
346	* nc_data_cluster_wll_set_kde:
347	* @dcwll: a #NcDataClusterWLL
348	* @kde: whether to use KDE method
349	*
350	* The reference count of @dcwll is decreased and the pointer is set to NULL.
351	*
352	*/
353	void
354	nc_data_cluster_wll_set_kde (NcDataClusterWLL *data, gboolean kde)	2✔
355	{
356	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	2✔
357	NcDataClusterWLLPrivate * const self = dcwll->priv;	2✔
358
359	self->kde = kde;	2✔
360	}	2✔
361

NumCosmo / NumCosmo / 7310595895

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