7062407362

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

Build # 7062407362

Build Type

Pull #120

github

Committed by

web-flow

Commit Message

Merge 182f70ced into 24e5f6a9b

Pull Request Pull Request #120: Moving to meson

Run Details

14317 of 71092 branches covered (0.0%)

Branch coverage included in aggregate %.

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

3942 existing lines in 211 files now uncovered.

57520 of 154555 relevant lines covered (37.22%)

2231793.7 hits per line

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

41.07

/numcosmo/model/nc_hiprim_sbpl.c

/***************************************************************************
 *            nc_hiprim_sbpl.c
 *
 *  Thu July 20 13:06:51 2017
 *  Copyright  2017  Sandro Dias Pinto Vitenti
 *  <vitenti@uel.br>
 ****************************************************************************/
/*
 * nc_hiprim_sbpl.c
 * Copyright (C) 2016 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_hiprim_sbpl
 * @title: NcHIPrimSBPL
 * @short_description: Smooth Broken power law modification of the power law primordial spectrum
 * 
 * FIXME
 * 
 */

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

#include "model/nc_hiprim_sbpl.h"

G_DEFINE_TYPE (NcHIPrimSBPL, nc_hiprim_sbpl, NC_TYPE_HIPRIM)

enum {
  PROP_0,
  PROP_SIZE,
};

static void
nc_hiprim_sbpl_init (NcHIPrimSBPL *nc_hiprim_sbpl)
{
}

static void
nc_hiprim_sbpl_finalize (GObject *object)
{

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

static gdouble _nc_hiprim_sbpl_lnSA_powespec_lnk (NcHIPrim *prim, const gdouble lnk);
static gdouble _nc_hiprim_sbpl_lnT_powespec_lnk (NcHIPrim *prim, const gdouble lnk);

static void
nc_hiprim_sbpl_class_init (NcHIPrimSBPLClass *klass)
{
  GObjectClass* object_class = G_OBJECT_CLASS (klass);
  NcHIPrimClass *prim_class  = NC_HIPRIM_CLASS (klass);
  NcmModelClass *model_class = NCM_MODEL_CLASS (klass);

  object_class->finalize = nc_hiprim_sbpl_finalize;

  ncm_model_class_set_name_nick (model_class, "Broken power law model for primordial spectra", "BPL");
  ncm_model_class_add_params (model_class, NC_HIPRIM_SBPL_SPARAM_LEN, 0, PROP_SIZE);

  /* Set ln10e10ASA param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LN10E10ASA, "\\log(10^{10}A_{SA})", "ln10e10ASA",
                              0.0, 5.0, 1.0e0,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LN10E10ASA,
                              NCM_PARAM_TYPE_FIXED);
  /* Set N_SA param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_N_SA, "n_{\\mathrm{SA}}", "n_SA",
                              0.5, 1.5, 1.0e-1,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_N_SA,
                              NCM_PARAM_TYPE_FIXED);
  /* Set lambdac param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_DELTA, "\\delta", "delta",
                              0.0, 8.0, 1.0,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_DELTA,
                              NCM_PARAM_TYPE_FIXED);
  /* Set lambdac param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_RA, "R_\\mathrm{A}", "RA",
                              0.0, 4.0, 1.0,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_RA,
                              NCM_PARAM_TYPE_FIXED);
  /* Set lnkc param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LNKB, "\\ln(k_\\mathrm{b})", "lnkb",
                              -12.0, -3.0, 1.0,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LNKB,
                              NCM_PARAM_TYPE_FIXED);
  /* Set lnkc param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LAMBDA, "\\lambda", "lambda",
                              1.0e-3, 30.0, 1.0,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LAMBDA,
                              NCM_PARAM_TYPE_FIXED);
  /* Set T_SA_ratio param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_T_SA_RATIO, "A_T/A_{\\mathrm{SA}}", "T_SA_ratio",
                              0.0, 10.0, 1.0e-1,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_T_SA_RATIO,
                              NCM_PARAM_TYPE_FIXED);
  /* Set N_T param info */
  ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_N_T, "n_{\\mathrm{T}}", "n_T",
                              -0.5, 0.5, 1.0e-2,
                              NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_N_T,
                              NCM_PARAM_TYPE_FIXED);

  /* Check for errors in parameters initialization */
  ncm_model_class_check_params_info (model_class);

  nc_hiprim_set_lnSA_powspec_lnk_impl (prim_class, &_nc_hiprim_sbpl_lnSA_powespec_lnk);
  nc_hiprim_set_lnT_powspec_lnk_impl  (prim_class, &_nc_hiprim_sbpl_lnT_powespec_lnk);
}

/**
 * nc_hiprim_sbpl_new: (constructor)
 *
 * FIXME
 *
 * Returns: (transfer full): FIXME
 */
NcHIPrimSBPL *
nc_hiprim_sbpl_new (void)
{
  NcHIPrimSBPL *prim_bpl = g_object_new (NC_TYPE_HIPRIM_SBPL,
                                          NULL);
  return prim_bpl;
}

#define VECTOR     (NCM_MODEL (prim)->params)
#define LN10E10ASA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LN10E10ASA))
#define N_SA       (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_N_SA))
#define DELTA      (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_DELTA))
#define RA         (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_RA))
#define LNKB       (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LNKB))
#define LAMBDA     (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LAMBDA))
#define T_SA_RATIO (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_T_SA_RATIO))
#define N_T        (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_N_T))

/****************************************************************************
 * Power spectrum
 ****************************************************************************/

static gdouble
_nc_hiprim_sbpl_lnSA_powespec_lnk (NcHIPrim *prim, const gdouble lnk)
{
  const gdouble ln_ka   = lnk - prim->lnk_pivot;
  const gdouble delta   = DELTA;
  const gdouble ra      = RA;
  const gdouble lnk_b   = LNKB;
  const gdouble lambda  = LAMBDA;

  const gdouble lnA_eff = LN10E10ASA - 10.0 * M_LN10;

  const gdouble lnPplaw = (N_SA - 1.0) * ln_ka + lnA_eff;

  const gdouble lnX     = (lnk - lnk_b);

  gdouble one_1px_plambda, one_1px_mlambda;
    
  if (lnX > 0.0)
  {
    const gdouble x_mlambda = exp (-lambda * lnX);

    one_1px_plambda = x_mlambda / (1.0 + x_mlambda);
    one_1px_mlambda = 1.0 / (1.0 + x_mlambda);
  }
  else
  {
    const gdouble x_plambda = exp (+lambda * lnX);

    one_1px_plambda = 1.0 / (1.0 + x_plambda);
    one_1px_mlambda = x_plambda / (1.0 + x_plambda);
  }

  {
    const gdouble lnxi = log (ra * exp (delta * ln_ka) * one_1px_plambda + one_1px_mlambda);
    return lnxi + lnPplaw;
  }
}

static gdouble
_nc_hiprim_sbpl_lnT_powespec_lnk (NcHIPrim *prim, const gdouble lnk)
{
  const gdouble ln_ka = lnk - prim->lnk_pivot;
  return N_T * ln_ka + LN10E10ASA - 10.0 * M_LN10 + log (T_SA_RATIO);
}


1	/***************************************************************************
2	* nc_hiprim_sbpl.c
3	*
4	* Thu July 20 13:06:51 2017
5	* Copyright 2017 Sandro Dias Pinto Vitenti
6	* <vitenti@uel.br>
7	****************************************************************************/
8	/*
9	* nc_hiprim_sbpl.c
10	* Copyright (C) 2016 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_hiprim_sbpl
28	* @title: NcHIPrimSBPL
29	* @short_description: Smooth Broken power law modification of the power law primordial spectrum
30	*
31	* FIXME
32	*
33	*/
34
35	#ifdef HAVE_CONFIG_H
36	# include "config.h"
37	#endif /* HAVE_CONFIG_H */
38	#include "build_cfg.h"
39
40	#include "model/nc_hiprim_sbpl.h"
41
42	G_DEFINE_TYPE (NcHIPrimSBPL, nc_hiprim_sbpl, NC_TYPE_HIPRIM)	333!
43
44	enum {
45	PROP_0,
46	PROP_SIZE,
47	};
48
49	static void
50	nc_hiprim_sbpl_init (NcHIPrimSBPL *nc_hiprim_sbpl)	×
51	{
52	}	×
53
54	static void
55	nc_hiprim_sbpl_finalize (GObject *object)	×
56	{
57
58	/* Chain up : end */
59	G_OBJECT_CLASS (nc_hiprim_sbpl_parent_class)->finalize (object);	×
60	}	×
61
62	static gdouble _nc_hiprim_sbpl_lnSA_powespec_lnk (NcHIPrim *prim, const gdouble lnk);
63	static gdouble _nc_hiprim_sbpl_lnT_powespec_lnk (NcHIPrim *prim, const gdouble lnk);
64
65	static void
66	nc_hiprim_sbpl_class_init (NcHIPrimSBPLClass *klass)	111✔
67	{
68	GObjectClass* object_class = G_OBJECT_CLASS (klass);	111✔
69	NcHIPrimClass *prim_class = NC_HIPRIM_CLASS (klass);	111✔
70	NcmModelClass *model_class = NCM_MODEL_CLASS (klass);	111✔
71
72	object_class->finalize = nc_hiprim_sbpl_finalize;	111✔
73
74	ncm_model_class_set_name_nick (model_class, "Broken power law model for primordial spectra", "BPL");	111✔
75	ncm_model_class_add_params (model_class, NC_HIPRIM_SBPL_SPARAM_LEN, 0, PROP_SIZE);	111✔
76
77	/* Set ln10e10ASA param info */
78	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LN10E10ASA, "\\log(10^{10}A_{SA})", "ln10e10ASA",	111✔
79	0.0, 5.0, 1.0e0,
80	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LN10E10ASA,
81	NCM_PARAM_TYPE_FIXED);
82	/* Set N_SA param info */
83	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_N_SA, "n_{\\mathrm{SA}}", "n_SA",	111✔
84	0.5, 1.5, 1.0e-1,
85	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_N_SA,
86	NCM_PARAM_TYPE_FIXED);
87	/* Set lambdac param info */
88	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_DELTA, "\\delta", "delta",	111✔
89	0.0, 8.0, 1.0,
90	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_DELTA,
91	NCM_PARAM_TYPE_FIXED);
92	/* Set lambdac param info */
93	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_RA, "R_\\mathrm{A}", "RA",	111✔
94	0.0, 4.0, 1.0,
95	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_RA,
96	NCM_PARAM_TYPE_FIXED);
97	/* Set lnkc param info */
98	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LNKB, "\\ln(k_\\mathrm{b})", "lnkb",	111✔
99	-12.0, -3.0, 1.0,
100	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LNKB,
101	NCM_PARAM_TYPE_FIXED);
102	/* Set lnkc param info */
103	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_LAMBDA, "\\lambda", "lambda",	111✔
104	1.0e-3, 30.0, 1.0,
105	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_LAMBDA,
106	NCM_PARAM_TYPE_FIXED);
107	/* Set T_SA_ratio param info */
108	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_T_SA_RATIO, "A_T/A_{\\mathrm{SA}}", "T_SA_ratio",	111✔
109	0.0, 10.0, 1.0e-1,
110	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_T_SA_RATIO,
111	NCM_PARAM_TYPE_FIXED);
112	/* Set N_T param info */
113	ncm_model_class_set_sparam (model_class, NC_HIPRIM_SBPL_N_T, "n_{\\mathrm{T}}", "n_T",	111✔
114	-0.5, 0.5, 1.0e-2,
115	NC_HIPRIM_DEFAULT_PARAMS_ABSTOL, NC_HIPRIM_SBPL_DEFAULT_N_T,
116	NCM_PARAM_TYPE_FIXED);
117
118	/* Check for errors in parameters initialization */
119	ncm_model_class_check_params_info (model_class);	111✔
120
121	nc_hiprim_set_lnSA_powspec_lnk_impl (prim_class, &_nc_hiprim_sbpl_lnSA_powespec_lnk);	111✔
122	nc_hiprim_set_lnT_powspec_lnk_impl (prim_class, &_nc_hiprim_sbpl_lnT_powespec_lnk);	111✔
123	}	111✔
124
125	/**
126	* nc_hiprim_sbpl_new: (constructor)
127	*
128	* FIXME
129	*
130	* Returns: (transfer full): FIXME
131	*/
132	NcHIPrimSBPL *
133	nc_hiprim_sbpl_new (void)	×
134	{
135	NcHIPrimSBPL *prim_bpl = g_object_new (NC_TYPE_HIPRIM_SBPL,	×
136	NULL);
137	return prim_bpl;	×
138	}
139
140	#define VECTOR (NCM_MODEL (prim)->params)
141	#define LN10E10ASA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LN10E10ASA))
142	#define N_SA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_N_SA))
143	#define DELTA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_DELTA))
144	#define RA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_RA))
145	#define LNKB (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LNKB))
146	#define LAMBDA (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_LAMBDA))
147	#define T_SA_RATIO (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_T_SA_RATIO))
148	#define N_T (ncm_vector_get (VECTOR, NC_HIPRIM_SBPL_N_T))
149
150	/****************************************************************************
151	* Power spectrum
152	****************************************************************************/
153
154	static gdouble
155	_nc_hiprim_sbpl_lnSA_powespec_lnk (NcHIPrim *prim, const gdouble lnk)	×
156	{
157	const gdouble ln_ka = lnk - prim->lnk_pivot;	×
158	const gdouble delta = DELTA;	×
159	const gdouble ra = RA;	×
160	const gdouble lnk_b = LNKB;	×
161	const gdouble lambda = LAMBDA;	×
162
163	const gdouble lnA_eff = LN10E10ASA - 10.0 * M_LN10;	×
164
165	const gdouble lnPplaw = (N_SA - 1.0) * ln_ka + lnA_eff;	×
166
167	const gdouble lnX = (lnk - lnk_b);	×
168
UNCOV 169	gdouble one_1px_plambda, one_1px_mlambda;	×
170
171	if (lnX > 0.0)	×
172	{
173	const gdouble x_mlambda = exp (-lambda * lnX);	×
174
175	one_1px_plambda = x_mlambda / (1.0 + x_mlambda);	×
176	one_1px_mlambda = 1.0 / (1.0 + x_mlambda);	×
177	}
178	else
179	{
180	const gdouble x_plambda = exp (+lambda * lnX);	×
181
182	one_1px_plambda = 1.0 / (1.0 + x_plambda);	×
183	one_1px_mlambda = x_plambda / (1.0 + x_plambda);	×
184	}
185
186	{
187	const gdouble lnxi = log (ra * exp (delta * ln_ka) * one_1px_plambda + one_1px_mlambda);	×
188	return lnxi + lnPplaw;	×
189	}
190	}
191
192	static gdouble
193	_nc_hiprim_sbpl_lnT_powespec_lnk (NcHIPrim *prim, const gdouble lnk)	×
194	{
195	const gdouble ln_ka = lnk - prim->lnk_pivot;	×
196	return N_T * ln_ka + LN10E10ASA - 10.0 * M_LN10 + log (T_SA_RATIO);	×
197	}
198

NumCosmo / NumCosmo / 7062407362

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