7100283118

Committed 05 Dec 2023 11:49AM UTC coverage: 54.434% (-0.3%) from 54.692%

Build # 7100283118

Build Type

Pull #93

github

Committed by

web-flow

Commit Message

Merge c9e8ceb33 into 03fd86611

Pull Request Pull Request #93: Galaxy WL reformulation

Run Details

7129 of 12482 branches covered (0.0%)

Branch coverage included in aggregate %.

295 of 983 new or added lines in 15 files covered. (30.01%)

7 existing lines in 3 files now uncovered.

47105 of 87150 relevant lines covered (54.05%)

1642370.35 hits per line

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

12.6

/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));

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

  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);	672✔
69
70	static void
NEW 71	nc_data_cluster_wll_init (NcDataClusterWLL *dcwll)	×
72	{
NEW 73	NcDataClusterWLLPrivate * const self = dcwll->priv = nc_data_cluster_wll_get_instance_private (dcwll);	×
74
NEW 75	self->galaxy_array = ncm_obj_array_new ();	×
NEW 76	self->z_cluster = 0.0;	×
NEW 77	self->kde = FALSE;	×
78	}
79
80	static void
NEW 81	nc_data_cluster_wll_set_property (GObject object, guint prop_id, const GValue value, GParamSpec *pspec)	×
82	{
NEW 83	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);	×
NEW 84	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
85
NEW 86	g_return_if_fail (NC_IS_DATA_CLUSTER_WLL (object));	×
87
NEW 88	switch (prop_id)	×
89	{
NEW 90	case PROP_GALAXY_ARRAY:	×
91	{
NEW 92	NcmObjArray *galaxy_array = g_value_get_boxed (value);	×
93
NEW 94	g_clear_pointer (&self->galaxy_array, ncm_obj_array_unref);	×
NEW 95	self->galaxy_array = ncm_obj_array_ref (galaxy_array);	×
96
NEW 97	break;	×
98	}
NEW 99	case PROP_Z_CLUSTER:	×
NEW 100	self->z_cluster = g_value_get_double (value);	×
NEW 101	break;	×
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
NEW 137	nc_data_cluster_wll_dispose (GObject *object)	×
138	{
NEW 139	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (object);	×
NEW 140	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
141
NEW 142	ncm_obj_array_clear (&self->galaxy_array);	×
143
144	/* Chain up : end */
NEW 145	G_OBJECT_CLASS (nc_data_cluster_wll_parent_class)->dispose (object);	×
146	}
147
148	static void
NEW 149	nc_data_cluster_wll_finalize (GObject *object)	×
150	{
151	/NcDataClusterWLL dcwll = NC_DATA_CLUSTER_WLL (object);*/
152
153	/* Chain up : end */
NEW 154	G_OBJECT_CLASS (nc_data_cluster_wll_parent_class)->finalize (object);	×
155	}
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)	112✔
163	{
164	GObjectClass *object_class = G_OBJECT_CLASS (klass);	112✔
165	NcmDataClass *data_class = NCM_DATA_CLASS (klass);	112✔
166
167	object_class->set_property = nc_data_cluster_wll_set_property;	112✔
168	object_class->get_property = nc_data_cluster_wll_get_property;	112✔
169	object_class->dispose = nc_data_cluster_wll_dispose;	112✔
170	object_class->finalize = nc_data_cluster_wll_finalize;	112✔
171
172	g_object_class_install_property (object_class,	112✔
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,	112✔
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,	112✔
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;	112✔
198	data_class->get_length = &_nc_data_cluster_wll_get_len;	112✔
199	data_class->prepare = &_nc_data_cluster_wll_prepare;	112✔
200	}	112✔
201
202	static void
NEW 203	_nc_data_cluster_wll_m2lnL_val (NcmData data, NcmMSet mset, gdouble *m2lnL)	×
204	{
NEW 205	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	×
NEW 206	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
NEW 207	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
NEW 208	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	×
NEW 209	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	×
NEW 210	const guint ngal = self->galaxy_array->len;	×
NEW 211	guint i;	×
212
NEW 213	m2lnL[0] = 0.0;	×
214
NEW 215	for (i = 0; i < ngal; i++)	×
216	{
NEW 217	NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));	×
218
NEW 219	if (self->kde)	×
NEW 220	m2lnL[0] += nc_galaxy_wl_likelihood_kde_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster);	×
221	else
NEW 222	m2lnL[0] += nc_galaxy_wl_likelihood_eval_m2lnP (gwll_i, cosmo, dp, smd, self->z_cluster);	×
223	}
224
NEW 225	return;	×
226	}
227
228	static guint
NEW 229	_nc_data_cluster_wll_get_len (NcmData *data)	×
230	{
NEW 231	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	×
NEW 232	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
233
NEW 234	if (self->galaxy_array != NULL)	×
235	{
NEW 236	const guint ngal = self->galaxy_array->len;	×
NEW 237	guint len = 0;	×
NEW 238	guint i;	×
239
NEW 240	for (i = 0; i < ngal; i++)	×
241	{
NEW 242	NcGalaxyWLLikelihood *gwll_i = NC_GALAXY_WL_LIKELIHOOD (ncm_obj_array_peek (self->galaxy_array, i));	×
243
NEW 244	len += nc_galaxy_wl_likelihood_len (gwll_i);	×
245	}
246
NEW 247	return len;	×
248	}
249	else
250	{
251	return 0;
252	}
253	}
254
255	static void
NEW 256	_nc_data_cluster_wll_prepare (NcmData data, NcmMSet mset)	×
257	{
NEW 258	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	×
NEW 259	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
NEW 260	NcHICosmo *cosmo = NC_HICOSMO (ncm_mset_peek (mset, nc_hicosmo_id ()));	×
NEW 261	NcWLSurfaceMassDensity *smd = NC_WL_SURFACE_MASS_DENSITY (ncm_mset_peek (mset, nc_wl_surface_mass_density_id ()));	×
NEW 262	NcHaloDensityProfile *dp = NC_HALO_DENSITY_PROFILE (ncm_mset_peek (mset, nc_halo_density_profile_id ()));	×
263
NEW 264	g_assert ((cosmo != NULL) && (smd != NULL) && (dp != NULL));	×
265
NEW 266	nc_wl_surface_mass_density_prepare_if_needed (smd, cosmo);	×
267	}
268
269	/**
270	* nc_data_cluster_wll_new:
271	*
272	* Creates a new #NcDataClusterWLL.
273	*
274	* Returns: the newly created #NcDataClusterWLL.
275	*/
276	NcDataClusterWLL *
NEW 277	nc_data_cluster_wll_new (void)	×
278	{
NEW 279	NcDataClusterWLL *dcwll = g_object_new (NC_TYPE_DATA_CLUSTER_WLL,	×
280	NULL);
281
NEW 282	return dcwll;	×
283	}
284
285	/**
286	* nc_data_cluster_wll_new_from_file:
287	* @filename: file containing a serialized #NcDataClusterWLL
288	*
289	* Creates a new #NcDataClusterWLL from @filename.
290	*
291	* Returns: (transfer full): the newly created #NcDataClusterWLL.
292	*/
293	NcDataClusterWLL *
NEW 294	nc_data_cluster_wll_new_from_file (const gchar *filename)	×
295	{
NEW 296	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (ncm_serialize_global_from_file (filename));	×
297
NEW 298	g_assert (NC_IS_DATA_CLUSTER_WLL (dcwll));	×
299
NEW 300	return dcwll;	×
301	}
302
303	/**
304	* nc_data_cluster_wll_ref:
305	* @dcwll: a #NcDataClusterWLL
306	*
307	* Increases the reference count of @dcwll by one.
308	*
309	* Returns: (transfer full): @dcwll
310	*/
311	NcDataClusterWLL *
NEW 312	nc_data_cluster_wll_ref (NcDataClusterWLL *dcwll)	×
313	{
NEW 314	return g_object_ref (dcwll);	×
315	}
316
317	/**
318	* nc_data_cluster_wll_free:
319	* @dcwll: a #NcDataClusterWLL
320	*
321	* Atomically decrements the reference count of @dcwll by one. If the reference count drops to 0,
322	* all memory allocated by @dcwll is released.
323	*
324	*/
325	void
NEW 326	nc_data_cluster_wll_free (NcDataClusterWLL *dcwll)	×
327	{
NEW 328	g_object_unref (dcwll);	×
329	}
330
331	/**
332	* nc_data_cluster_wll_clear:
333	* @dcwll: a #NcDataClusterWLL
334	*
335	* The reference count of @dcwll is decreased and the pointer is set to NULL.
336	*
337	*/
338	void
NEW 339	nc_data_cluster_wll_clear (NcDataClusterWLL **dcwll)	×
340	{
NEW 341	g_clear_object (dcwll);	×
342	}
343
344	/**
345	* nc_data_cluster_wll_set_kde:
346	* @dcwll: a #NcDataClusterWLL
347	* @kde: whether to use KDE method
348	*
349	* The reference count of @dcwll is decreased and the pointer is set to NULL.
350	*
351	*/
352	void
NEW 353	nc_data_cluster_wll_set_kde (NcDataClusterWLL *data, gboolean kde)	×
354	{
NEW 355	NcDataClusterWLL *dcwll = NC_DATA_CLUSTER_WLL (data);	×
NEW 356	NcDataClusterWLLPrivate * const self = dcwll->priv;	×
357
NEW 358	self->kde = kde;	×
359	}
360

NumCosmo / NumCosmo / 7100283118

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