7311805470

Committed 24 Dec 2023 01:51AM UTC coverage: 55.27% (+0.004%) from 55.266%

Build # 7311805470

Build Type

Pull #85

github

Committed by

web-flow

Commit Message

Merge 88d822566 into 052525214

Pull Request Pull Request #85: Xcor cmp

Run Details

7284 of 12566 branches covered (0.0%)

Branch coverage included in aggregate %.

21 of 69 new or added lines in 2 files covered. (30.43%)

2 existing lines in 1 file now uncovered.

47975 of 87415 relevant lines covered (54.88%)

2251552.19 hits per line

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

28.0

/numcosmo/xcor/nc_xcor_limber_kernel_tSZ.c

/***************************************************************************
 *            nc_xcor_limber_kernel_tSZ.c
 *
 *  Tue January 10 12:00:00 2022
 *  Copyright  2022  Arthur de Souza Molina
 *  <arthur.souza.molina@uel.br>
 ****************************************************************************/
/*
 * numcosmo
 * Copyright (C) 2022 Arthur de Souza Molina <arthur.souza.molina@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_xcor_limber_kernel_tSZ
 * @title: NcXcorLimberKerneltSZ
 * @short_description: implementation of #NcXcorLimberKernel with the thermal Sunyaev Zel'dovich effect
 *
 * The thermal Sunyaev Zel'dovich (tSZ) effect is a modification in the observed temperature of the cosmic microwave 
 * background (CMB) due to the inverse Compton scattering of CMB photons with high-energy electrons along the line-of-sight. 
 * These electrons are present in the intracluster medium (ICM) of galaxy clusters, for example.   
 *   Compton-y parameter
 * The kernel is given by
 * \begin{equation}
 *    W(\\chi) = \\frac{\\sigma_T}{m_ec^2} \\frac{1}{1+z}.
 * \end{equation}
 *
 */

#ifdef HAVE_CONFIG_H
#include "config.h"
#endif /* HAVE_CONFIG_H */
#include "build_cfg.h"
#include "math/ncm_cfg.h"
#include "xcor/nc_xcor_limber_kernel_tSZ.h"
#include "xcor/nc_xcor.h"

#ifndef NUMCOSMO_GIR_SCAN

#endif /* NUMCOSMO_GIR_SCAN */

G_DEFINE_TYPE (NcXcorLimberKerneltSZ, nc_xcor_limber_kernel_tsz, NC_TYPE_XCOR_LIMBER_KERNEL);

#define VECTOR (NCM_MODEL (xclkl)->params)

enum
{
  PROP_0,
  PROP_ZMAX,
  PROP_SIZE,
};

static void
nc_xcor_limber_kernel_tsz_init (NcXcorLimberKerneltSZ *xclkl)
{
  xclkl->Zmax = 0;
}

static void
_nc_xcor_limber_kernel_tsz_set_property (GObject *object, guint prop_id, const GValue *value, GParamSpec *pspec)
{
  NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (object);
  
  g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (object));
  
  switch (prop_id)
  {
    case PROP_ZMAX:
      xclkl->Zmax = g_value_get_double(value);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
_nc_xcor_limber_kernel_tsz_get_property (GObject *object, guint prop_id, GValue *value, GParamSpec *pspec)
{
  NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (object);
  
  g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (object));
  
  switch (prop_id)
  {
    case PROP_ZMAX:
      g_value_set_double (value, xclkl->Zmax);
      break;
    default:
      G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);
      break;
  }
}

static void
_nc_xcor_limber_kernel_tsz_dispose (GObject *object)
{

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

static void
_nc_xcor_limber_kernel_tsz_finalize (GObject *object)
{
  /* Chain up : end */
  G_OBJECT_CLASS (nc_xcor_limber_kernel_tsz_parent_class)->finalize (object);
}

static gdouble _nc_xcor_limber_kernel_tsz_eval (NcXcorLimberKernel *xclk, NcHICosmo *cosmo, gdouble z, const NcXcorKinetic *xck, gint l);
static void _nc_xcor_limber_kernel_tsz_prepare (NcXcorLimberKernel *xclk, NcHICosmo *cosmo);
static void _nc_xcor_limber_kernel_tsz_add_noise (NcXcorLimberKernel *xclk, NcmVector *vp1, NcmVector *vp2, guint lmin);
static guint _nc_xcor_limber_kernel_tsz_obs_len (NcXcorLimberKernel *xclk);
static guint _nc_xcor_limber_kernel_tsz_obs_params_len (NcXcorLimberKernel *xclk);

static void
nc_xcor_limber_kernel_tsz_class_init (NcXcorLimberKerneltSZClass *klass)
{
  GObjectClass *object_class            = G_OBJECT_CLASS (klass);
  NcXcorLimberKernelClass *parent_class = NC_XCOR_LIMBER_KERNEL_CLASS (klass);
  NcmModelClass *model_class            = NCM_MODEL_CLASS (klass);
  
  object_class->finalize    = &_nc_xcor_limber_kernel_tsz_finalize;
  object_class->dispose     = &_nc_xcor_limber_kernel_tsz_dispose;
  model_class->set_property = &_nc_xcor_limber_kernel_tsz_set_property;
  model_class->get_property = &_nc_xcor_limber_kernel_tsz_get_property;
  
  ncm_model_class_set_name_nick (model_class, "Xcor with the thermal Sunyaev Zel'dovich Compton-y parameter", "Xcor-tSZ");
  ncm_model_class_add_params (model_class, 0, 0, PROP_SIZE);
  
  /**
   * NcXcorLimberKerneltSZ:Zmax:
   *
   * FIXME Set correct values (limits)
   */
  g_object_class_install_property (object_class,
                                   PROP_ZMAX,
                                                                   g_param_spec_double ("Zmax",
                                                        NULL,
                                                        "Maximum redshift up to which we define the kernel",
                                                        0,6.0, 6.0,
                                                        G_PARAM_READWRITE | G_PARAM_STATIC_NAME | G_PARAM_STATIC_BLURB));

  /* Check for errors in parameters initialization */
  ncm_model_class_check_params_info (model_class);
  
  parent_class->eval    = &_nc_xcor_limber_kernel_tsz_eval;
  parent_class->prepare = &_nc_xcor_limber_kernel_tsz_prepare;
  parent_class->add_noise = &_nc_xcor_limber_kernel_tsz_add_noise;
  
  parent_class->obs_len        = &_nc_xcor_limber_kernel_tsz_obs_len;
  parent_class->obs_params_len = &_nc_xcor_limber_kernel_tsz_obs_params_len;
  
  ncm_model_class_add_impl_flag (model_class, NC_XCOR_LIMBER_KERNEL_IMPL_ALL);
}

/**
 * nc_xcor_limber_kernel_tsz_new:
 * @Zmax: a #Gdouble
 *
 * FIXME
 *
 * Returns: FIXME
 *
 */
NcXcorLimberKerneltSZ *
nc_xcor_limber_kernel_tsz_new (gdouble Zmax)
{
  NcXcorLimberKerneltSZ *xclkl = g_object_new (NC_TYPE_XCOR_LIMBER_KERNEL_TSZ,
                                                      "Zmax", Zmax,
                                                      NULL);
  
  return xclkl;
}

/*NcmVector *ncm_vector_linear_space ( gdouble start, gdouble stop, gdouble value)
{
        NcmVector *LinearSpace = ncm_vector_new (stop - start);

        gdouble mean = value / (stop - start);

        for (int i=start;i<stop;i++){

                 ncm_vector_set (LinearSpace, i, mean * i );
        }

        return LinearSpace;
}*/

static gdouble
_nc_xcor_limber_kernel_tsz_eval (NcXcorLimberKernel *xclk, NcHICosmo *cosmo, gdouble z, const NcXcorKinetic *xck, gint l)
{

  NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);

  NCM_UNUSED (l);

  const gdouble nc_prefac = ncm_c_thomson_cs() / ( ncm_c_mass_e() * ncm_c_c() * ncm_c_c() );

  if (nc_prefac * (1.0 / (1+ xclkl->Zmax)) < 0)
          return - nc_prefac * (1.0 / (1+ xclkl->Zmax));
  else {
          return nc_prefac * (1.0 / (1+ xclkl->Zmax));
  }
}

static void
_nc_xcor_limber_kernel_tsz_prepare (NcXcorLimberKernel *xclk, NcHICosmo *cosmo)
{
  NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);
  g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (xclkl));

}

static void
_nc_xcor_limber_kernel_tsz_add_noise (NcXcorLimberKernel *xclk, NcmVector *vp1, NcmVector *vp2, guint lmin)
{
  NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);

  NCM_UNUSED (xclk);
  NCM_UNUSED (xclkl);
}

static guint
_nc_xcor_limber_kernel_tsz_obs_len (NcXcorLimberKernel *xclk)
{
  NCM_UNUSED (xclk);
  
  return 1;
}

static guint
_nc_xcor_limber_kernel_tsz_obs_params_len (NcXcorLimberKernel *xclk)
{
  NCM_UNUSED (xclk);
  
  return 0;
}

1	/***************************************************************************
2	* nc_xcor_limber_kernel_tSZ.c
3	*
4	* Tue January 10 12:00:00 2022
5	* Copyright 2022 Arthur de Souza Molina
6	* <arthur.souza.molina@uel.br>
7	****************************************************************************/
8	/*
9	* numcosmo
10	* Copyright (C) 2022 Arthur de Souza Molina <arthur.souza.molina@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_xcor_limber_kernel_tSZ
28	* @title: NcXcorLimberKerneltSZ
29	* @short_description: implementation of #NcXcorLimberKernel with the thermal Sunyaev Zel'dovich effect
30	*
31	* The thermal Sunyaev Zel'dovich (tSZ) effect is a modification in the observed temperature of the cosmic microwave
32	* background (CMB) due to the inverse Compton scattering of CMB photons with high-energy electrons along the line-of-sight.
33	* These electrons are present in the intracluster medium (ICM) of galaxy clusters, for example.
34	* Compton-y parameter
35	* The kernel is given by
36	* \begin{equation}
37	* W(\\chi) = \\frac{\\sigma_T}{m_ec^2} \\frac{1}{1+z}.
38	* \end{equation}
39	*
40	*/
41
42	#ifdef HAVE_CONFIG_H
43	#include "config.h"
44	#endif /* HAVE_CONFIG_H */
45	#include "build_cfg.h"
46	#include "math/ncm_cfg.h"
47	#include "xcor/nc_xcor_limber_kernel_tSZ.h"
48	#include "xcor/nc_xcor.h"
49
50	#ifndef NUMCOSMO_GIR_SCAN
51
52	#endif /* NUMCOSMO_GIR_SCAN */
53
54	G_DEFINE_TYPE (NcXcorLimberKerneltSZ, nc_xcor_limber_kernel_tsz, NC_TYPE_XCOR_LIMBER_KERNEL);	672✔
55
56	#define VECTOR (NCM_MODEL (xclkl)->params)
57
58	enum
59	{
60	PROP_0,
61	PROP_ZMAX,
62	PROP_SIZE,
63	};
64
65	static void
NEW 66	nc_xcor_limber_kernel_tsz_init (NcXcorLimberKerneltSZ *xclkl)	×
67	{
NEW 68	xclkl->Zmax = 0;	×
69	}
70
71	static void
NEW 72	_nc_xcor_limber_kernel_tsz_set_property (GObject object, guint prop_id, const GValue value, GParamSpec *pspec)	×
73	{
NEW 74	NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (object);	×
75
NEW 76	g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (object));	×
77
NEW 78	switch (prop_id)	×
79	{
NEW 80	case PROP_ZMAX:	×
NEW 81	xclkl->Zmax = g_value_get_double(value);	×
NEW 82	break;	×
NEW 83	default:	×
NEW 84	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);	×
NEW 85	break;	×
86	}
87	}
88
89	static void
NEW 90	_nc_xcor_limber_kernel_tsz_get_property (GObject object, guint prop_id, GValue value, GParamSpec *pspec)	×
91	{
NEW 92	NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (object);	×
93
NEW 94	g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (object));	×
95
NEW 96	switch (prop_id)	×
97	{
NEW 98	case PROP_ZMAX:	×
NEW 99	g_value_set_double (value, xclkl->Zmax);	×
NEW 100	break;	×
NEW 101	default:	×
NEW 102	G_OBJECT_WARN_INVALID_PROPERTY_ID (object, prop_id, pspec);	×
NEW 103	break;	×
104	}
105	}
106
107	static void
NEW 108	_nc_xcor_limber_kernel_tsz_dispose (GObject *object)	×
109	{
110
111	/* Chain up : end */
NEW 112	G_OBJECT_CLASS (nc_xcor_limber_kernel_tsz_parent_class)->dispose (object);	×
113	}
114
115	static void
NEW 116	_nc_xcor_limber_kernel_tsz_finalize (GObject *object)	×
117	{
118	/* Chain up : end */
NEW 119	G_OBJECT_CLASS (nc_xcor_limber_kernel_tsz_parent_class)->finalize (object);	×
120	}
121
122	static gdouble _nc_xcor_limber_kernel_tsz_eval (NcXcorLimberKernel xclk, NcHICosmo cosmo, gdouble z, const NcXcorKinetic *xck, gint l);
123	static void _nc_xcor_limber_kernel_tsz_prepare (NcXcorLimberKernel xclk, NcHICosmo cosmo);
124	static void _nc_xcor_limber_kernel_tsz_add_noise (NcXcorLimberKernel xclk, NcmVector vp1, NcmVector *vp2, guint lmin);
125	static guint _nc_xcor_limber_kernel_tsz_obs_len (NcXcorLimberKernel *xclk);
126	static guint _nc_xcor_limber_kernel_tsz_obs_params_len (NcXcorLimberKernel *xclk);
127
128	static void
129	nc_xcor_limber_kernel_tsz_class_init (NcXcorLimberKerneltSZClass *klass)	112✔
130	{
131	GObjectClass *object_class = G_OBJECT_CLASS (klass);	112✔
132	NcXcorLimberKernelClass *parent_class = NC_XCOR_LIMBER_KERNEL_CLASS (klass);	112✔
133	NcmModelClass *model_class = NCM_MODEL_CLASS (klass);	112✔
134
135	object_class->finalize = &_nc_xcor_limber_kernel_tsz_finalize;	112✔
136	object_class->dispose = &_nc_xcor_limber_kernel_tsz_dispose;	112✔
137	model_class->set_property = &_nc_xcor_limber_kernel_tsz_set_property;	112✔
138	model_class->get_property = &_nc_xcor_limber_kernel_tsz_get_property;	112✔
139
140	ncm_model_class_set_name_nick (model_class, "Xcor with the thermal Sunyaev Zel'dovich Compton-y parameter", "Xcor-tSZ");	112✔
141	ncm_model_class_add_params (model_class, 0, 0, PROP_SIZE);	112✔
142
143	/**
144	* NcXcorLimberKerneltSZ:Zmax:
145	*
146	* FIXME Set correct values (limits)
147	*/
148	g_object_class_install_property (object_class,	112✔
149	PROP_ZMAX,
150	g_param_spec_double ("Zmax",
151	NULL,
152	"Maximum redshift up to which we define the kernel",
153	0,6.0, 6.0,
154	G_PARAM_READWRITE \| G_PARAM_STATIC_NAME \| G_PARAM_STATIC_BLURB));
155
156	/* Check for errors in parameters initialization */
157	ncm_model_class_check_params_info (model_class);	112✔
158
159	parent_class->eval = &_nc_xcor_limber_kernel_tsz_eval;	112✔
160	parent_class->prepare = &_nc_xcor_limber_kernel_tsz_prepare;	112✔
161	parent_class->add_noise = &_nc_xcor_limber_kernel_tsz_add_noise;	112✔
162
163	parent_class->obs_len = &_nc_xcor_limber_kernel_tsz_obs_len;	112✔
164	parent_class->obs_params_len = &_nc_xcor_limber_kernel_tsz_obs_params_len;	112✔
165
166	ncm_model_class_add_impl_flag (model_class, NC_XCOR_LIMBER_KERNEL_IMPL_ALL);	112✔
167	}	112✔
168
169	/**
170	* nc_xcor_limber_kernel_tsz_new:
171	* @Zmax: a #Gdouble
172	*
173	* FIXME
174	*
175	* Returns: FIXME
176	*
177	*/
178	NcXcorLimberKerneltSZ *
NEW 179	nc_xcor_limber_kernel_tsz_new (gdouble Zmax)	×
180	{
NEW 181	NcXcorLimberKerneltSZ *xclkl = g_object_new (NC_TYPE_XCOR_LIMBER_KERNEL_TSZ,	×
182	"Zmax", Zmax,
183	NULL);
184
NEW 185	return xclkl;	×
186	}
187
188	/NcmVector ncm_vector_linear_space ( gdouble start, gdouble stop, gdouble value)
189	{
190	NcmVector *LinearSpace = ncm_vector_new (stop - start);
191
192	gdouble mean = value / (stop - start);
193
194	for (int i=start;i<stop;i++){
195
196	ncm_vector_set (LinearSpace, i, mean * i );
197	}
198
199	return LinearSpace;
200	}*/
201
202	static gdouble
NEW 203	_nc_xcor_limber_kernel_tsz_eval (NcXcorLimberKernel xclk, NcHICosmo cosmo, gdouble z, const NcXcorKinetic *xck, gint l)	×
204	{
205
NEW 206	NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);	×
207
NEW 208	NCM_UNUSED (l);	×
209
NEW 210	const gdouble nc_prefac = ncm_c_thomson_cs() / ( ncm_c_mass_e() * ncm_c_c() * ncm_c_c() );	×
211
NEW 212	if (nc_prefac * (1.0 / (1+ xclkl->Zmax)) < 0)	×
NEW 213	return - nc_prefac * (1.0 / (1+ xclkl->Zmax));	×
214	else {
215	return nc_prefac * (1.0 / (1+ xclkl->Zmax));
216	}
217	}
218
219	static void
NEW 220	_nc_xcor_limber_kernel_tsz_prepare (NcXcorLimberKernel xclk, NcHICosmo cosmo)	×
221	{
NEW 222	NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);	×
NEW 223	g_return_if_fail (NC_IS_XCOR_LIMBER_KERNEL_TSZ (xclkl));	×
224
225	}
226
227	static void
NEW 228	_nc_xcor_limber_kernel_tsz_add_noise (NcXcorLimberKernel xclk, NcmVector vp1, NcmVector *vp2, guint lmin)	×
229	{
NEW 230	NcXcorLimberKerneltSZ *xclkl = NC_XCOR_LIMBER_KERNEL_TSZ (xclk);	×
231
NEW 232	NCM_UNUSED (xclk);	×
NEW 233	NCM_UNUSED (xclkl);	×
234	}
235
236	static guint
NEW 237	_nc_xcor_limber_kernel_tsz_obs_len (NcXcorLimberKernel *xclk)	×
238	{
NEW 239	NCM_UNUSED (xclk);	×
240
NEW 241	return 1;	×
242	}
243
244	static guint
NEW 245	_nc_xcor_limber_kernel_tsz_obs_params_len (NcXcorLimberKernel *xclk)	×
246	{
NEW 247	NCM_UNUSED (xclk);	×
248
NEW 249	return 0;	×
250	}

NumCosmo / NumCosmo / 7311805470

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