7176062912

Committed 12 Dec 2023 02:19AM UTC coverage: 55.458% (+0.5%) from 54.987%

Build # 7176062912

Build Type

Pull #93

github

Committed by

web-flow

Commit Message

Merge 8c3e4a862 into 04c61dd4e

Pull Request Pull Request #93: Galaxy WL reformulation

Run Details

7391 of 12715 branches covered (0.0%)

Branch coverage included in aggregate %.

1146 of 1428 new or added lines in 22 files covered. (80.25%)

11 existing lines in 5 files now uncovered.

48606 of 88256 relevant lines covered (55.07%)

1641662.95 hits per line

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

92.0

/numcosmo/galaxy/nc_galaxy_sd_position_lsst_srd.c

/***************************************************************************
 *            nc_galaxy_sd_position_lsst_srd.c
 *
 *  Tue June 22 14:59:27 2023
 *  Copyright  2023  Caio Lima de Oliveira
 *  <caiolimadeoliveira@pm.me>
 ****************************************************************************/
/*
 * nc_galaxy_sd_position_lsst_srd.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_galaxy_sd_position_lsst_srd
 * @title: NcGalaxySDPositionLSSTSRD
 * @short_description: Class describing galaxy sample position distributions as in LSST-SRD
 * @stability: Unstable
 *
 * Class defining a galaxy sample position distribution as described in the LSST
 * Science Roadmap Document, it includes the redshift distribution and a angular radius
 * distribution.
 *
 */

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

/* #include "nc_enum_types.h" */
#include "galaxy/nc_galaxy_sd_position_lsst_srd.h"
#include "galaxy/nc_galaxy_sd_position.h"
#include "math/ncm_stats_dist1d.h"
#include "math/ncm_stats_dist1d_spline.h"
#include "math/ncm_dtuple.h"
#include "math/ncm_spline.h"
#include "math/ncm_spline_func.h"
#include "math/ncm_spline_cubic.h"
#include "math/ncm_spline_cubic_notaknot.h"
#include <math.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_randist.h>

typedef struct _NcGalaxySDPositionLSSTSRDPrivate
{
  gdouble z_lb;
  gdouble z_ub;
  gdouble r_norm;
  gdouble r_lb;
  gdouble r_ub;
  gdouble r_lb2;
  gdouble r_ub2;
} NcGalaxySDPositionLSSTSRDPrivate;

struct _NcGalaxySDPositionLSSTSRD
{
  NcGalaxySDPosition parent_instance;
};

enum
{
  PROP_0,
  PROP_LEN,
};

G_DEFINE_TYPE_WITH_PRIVATE (NcGalaxySDPositionLSSTSRD, nc_galaxy_sd_position_lsst_srd, NC_TYPE_GALAXY_SD_POSITION);

static void
nc_galaxy_sd_position_lsst_srd_init (NcGalaxySDPositionLSSTSRD *gsdplsst)
{
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);

  self->z_lb   = 0.0;
  self->z_ub   = 0.0;
  self->r_norm = 0.0;
  self->r_lb   = 0.0;
  self->r_ub   = 0.0;
  self->r_lb2  = 0.0;
  self->r_ub2  = 0.0;
}

static void
_nc_galaxy_sd_position_lsst_srd_dispose (GObject *object)
{
  /* Chain up : end */
  G_OBJECT_CLASS (nc_galaxy_sd_position_lsst_srd_parent_class)->dispose (object);
}

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

static gdouble _nc_galaxy_sd_position_lsst_srd_gen_r (NcGalaxySDPosition *gsdp, NcmRNG *rng);
static gdouble _nc_galaxy_sd_position_lsst_srd_gen_z (NcGalaxySDPosition *gsdp, NcmRNG *rng);
static gdouble _nc_galaxy_sd_position_lsst_srd_integ (NcGalaxySDPosition *gsdp, const gdouble r, const gdouble z);
static void _nc_galaxy_sd_position_lsst_srd_set_z_lim (NcGalaxySDPosition *gsdp, gdouble z_min, gdouble z_max);
static void _nc_galaxy_sd_position_lsst_srd_get_z_lim (NcGalaxySDPosition *gsdp, gdouble *z_min, gdouble *z_max);
static void _nc_galaxy_sd_position_lsst_srd_set_r_lim (NcGalaxySDPosition *gsdp, gdouble r_min, gdouble r_max);
static void _nc_galaxy_sd_position_lsst_srd_get_r_lim (NcGalaxySDPosition *gsdp, gdouble *r_min, gdouble *r_max);

static void
nc_galaxy_sd_position_lsst_srd_class_init (NcGalaxySDPositionLSSTSRDClass *klass)
{
  NcGalaxySDPositionClass *sd_position_class = NC_GALAXY_SD_POSITION_CLASS (klass);
  GObjectClass *object_class                 = G_OBJECT_CLASS (klass);
  NcmModelClass *model_class                 = NCM_MODEL_CLASS (klass);

  object_class->dispose  = &_nc_galaxy_sd_position_lsst_srd_dispose;
  object_class->finalize = &_nc_galaxy_sd_position_lsst_srd_finalize;

  ncm_model_class_set_name_nick (model_class, "LSST SRD Galaxy Distribution", "LSST_SRD");
  ncm_model_class_add_params (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_SPARAM_LEN, 0, PROP_LEN);

  /**
   * NcGalaxySDPositionLSSTSRD:alpha:
   *
   * The redshift exponential slope.
   *
   */
  ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_ALPHA, "\\alpha", "alpha",
                              1e-8,  1.0, 1.0e-2,
                              NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_ALPHA,
                              NCM_PARAM_TYPE_FIXED);

  /**
   * NcGalaxySDPositionLSSTSRD:beta:
   *
   * The redshift power law slope.
   *
   */
  ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_BETA, "\\beta", "beta",
                              1e-8,  5.0, 1.0e-1,
                              NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_BETA,
                              NCM_PARAM_TYPE_FIXED);

  /**
   * NcGalaxySDPositionLSSTSRD:z0:
   *
   * The redshift pivot.
   *
   */
  ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_Z0, "z_0", "z0",
                              1e-8,  10.0, 1.0e-2,
                              NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_Z0,
                              NCM_PARAM_TYPE_FIXED);

  ncm_model_class_check_params_info (model_class);

  sd_position_class->gen_r     = &_nc_galaxy_sd_position_lsst_srd_gen_r;
  sd_position_class->gen_z     = &_nc_galaxy_sd_position_lsst_srd_gen_z;
  sd_position_class->integ     = &_nc_galaxy_sd_position_lsst_srd_integ;
  sd_position_class->set_z_lim = &_nc_galaxy_sd_position_lsst_srd_set_z_lim;
  sd_position_class->get_z_lim = &_nc_galaxy_sd_position_lsst_srd_get_z_lim;
  sd_position_class->set_r_lim = &_nc_galaxy_sd_position_lsst_srd_set_r_lim;
  sd_position_class->get_r_lim = &_nc_galaxy_sd_position_lsst_srd_get_r_lim;
}

#define VECTOR (NCM_MODEL (gsdp))
#define ALPHA  (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_ALPHA))
#define BETA   (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_BETA))
#define Z0     (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_Z0))

static gdouble
_nc_galaxy_sd_position_lsst_srd_gen_r (NcGalaxySDPosition *gsdp, NcmRNG *rng)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);
  const gdouble cumul_gen                       = ncm_rng_uniform_gen (rng, 0.0, 1.0);

  return sqrt (cumul_gen * 2.0 / self->r_norm + self->r_lb2);
}

static gdouble
_nc_galaxy_sd_position_lsst_srd_gen_z (NcGalaxySDPosition *gsdp, NcmRNG *rng)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);
  const gdouble alpha                           = ALPHA;
  const gdouble beta                            = BETA;
  const gdouble z0                              = Z0;
  const gdouble y0                              = pow (z0, alpha);
  const gdouble gamma_a                         = (1.0 + beta) / alpha;

  gdouble z;

  do {
    const gdouble gen_y = ncm_rng_gamma_gen (rng, gamma_a, y0);

    z = pow (gen_y, 1.0 / alpha);
  } while (z < self->z_lb || z > self->z_ub);

  return z;
}

static gdouble
_nc_galaxy_sd_position_lsst_srd_integ (NcGalaxySDPosition *gsdp, const gdouble r, const gdouble z)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);
  const gdouble alpha                           = ALPHA;
  const gdouble beta                            = BETA;
  const gdouble z0                              = Z0;
  const gdouble y0                              = pow (z0, alpha);
  const gdouble gamma_a                         = (1.0 + beta) / alpha;
  const gdouble y                               = pow (z, alpha);

  return gsl_ran_gamma_pdf (y, gamma_a, y0) * alpha * y / z * r * self->r_norm;
}

static void
_nc_galaxy_sd_position_lsst_srd_set_z_lim (NcGalaxySDPosition *gsdp, gdouble z_min, gdouble z_max)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);

  g_assert_cmpfloat (z_min, <, z_max);

  self->z_lb = z_min;
  self->z_ub = z_max;
}

static void
_nc_galaxy_sd_position_lsst_srd_get_z_lim (NcGalaxySDPosition *gsdp, gdouble *z_min, gdouble *z_max)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);

  g_assert_nonnull (z_min);
  g_assert_nonnull (z_max);

  *z_min = self->z_lb;
  *z_max = self->z_ub;
}

static void
_nc_galaxy_sd_position_lsst_srd_set_r_lim (NcGalaxySDPosition *gsdp, gdouble r_min, gdouble r_max)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);

  g_assert_cmpfloat (r_min, <, r_max);

  self->r_lb = r_min;
  self->r_ub = r_max;

  self->r_lb2  = self->r_lb * self->r_lb;
  self->r_ub2  = self->r_ub * self->r_ub;
  self->r_norm = 2.0 / (self->r_ub2 - self->r_lb2);
}

static void
_nc_galaxy_sd_position_lsst_srd_get_r_lim (NcGalaxySDPosition *gsdp, gdouble *r_min, gdouble *r_max)
{
  NcGalaxySDPositionLSSTSRD *gsdplsst           = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);
  NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);

  g_assert_nonnull (r_min);
  g_assert_nonnull (r_max);

  *r_min = self->r_lb;
  *r_max = self->r_ub;
}

/**
 * nc_galaxy_sd_position_lsst_srd_new:
 * @z_min: the minimum redshift
 * @z_max: the maximum redshift
 * @r_min: the minimum angular radius
 * @r_max: the maximum angular radius
 *
 * Creates a new #NcGalaxySDPositionLSSTSRD, the parameter values correspond to the
 * LSST SRD year 1.
 *
 * Returns: (transfer full): a new NcGalaxySDPositionLSSTSRD.
 */
NcGalaxySDPositionLSSTSRD *
nc_galaxy_sd_position_lsst_srd_new (const gdouble z_min, const gdouble z_max, const gdouble r_min, const gdouble r_max)
{
  NcmDTuple2 z_lim                    = NCM_DTUPLE2_STATIC_INIT (z_min, z_max);
  NcmDTuple2 r_lim                    = NCM_DTUPLE2_STATIC_INIT (r_min, r_max);
  NcGalaxySDPositionLSSTSRD *gsdplsst = g_object_new (NC_TYPE_GALAXY_SD_POSITION_LSST_SRD,
                                                      "z-lim", &z_lim,
                                                      "r-lim", &r_lim,
                                                      NULL);

  return gsdplsst;
}

/**
 * nc_galaxy_sd_position_lsst_srd_new_y10:
 * @z_min: the minimum redshift
 * @z_max: the maximum redshift
 * @r_min: the minimum angular radius
 * @r_max: the maximum angular radius
 *
 * Creates a new #NcGalaxySDPositionLSSTSRD, the parameter values correspond to the
 * LSST SRD year 10.
 *
 * Returns: (transfer full): a new NcGalaxySDPositionLSSTSRD.
 */
NcGalaxySDPositionLSSTSRD *
nc_galaxy_sd_position_lsst_srd_new_y10 (const gdouble z_min, const gdouble z_max, const gdouble r_min, const gdouble r_max)
{
  NcmDTuple2 z_lim                    = NCM_DTUPLE2_STATIC_INIT (z_min, z_max);
  NcmDTuple2 r_lim                    = NCM_DTUPLE2_STATIC_INIT (r_min, r_max);
  NcGalaxySDPositionLSSTSRD *gsdplsst = g_object_new (NC_TYPE_GALAXY_SD_POSITION_LSST_SRD,
                                                      "alpha", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_ALPHA,
                                                      "beta", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_BETA,
                                                      "z0", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_Z0,
                                                      "z-lim", &z_lim,
                                                      "r-lim", &r_lim,
                                                      NULL);

  return gsdplsst;
}

/**
 * nc_galaxy_sd_position_lsst_srd_ref:
 * @gsdplsst: a #NcGalaxySDPositionLSSTSRD
 *
 * Increase the reference of @gsdplsst by one.
 *
 * Returns: (transfer full): @gsdplsst.
 */
NcGalaxySDPositionLSSTSRD *
nc_galaxy_sd_position_lsst_srd_ref (NcGalaxySDPositionLSSTSRD *gsdplsst)
{
  return g_object_ref (gsdplsst);
}

/**
 * nc_galaxy_sd_position_lsst_srd_free:
 * @gsdplsst: a #NcGalaxySDPositionLSSTSRD
 *
 * Decrease the reference count of @gsdplsst by one.
 *
 */
void
nc_galaxy_sd_position_lsst_srd_free (NcGalaxySDPositionLSSTSRD *gsdplsst)
{
  g_object_unref (gsdplsst);
}

/**
 * nc_galaxy_sd_position_lsst_srd_clear:
 * @gsdplsst: a #NcGalaxySDPositionLSSTSRD
 *
 * Decrease the reference count of @gsdplsst by one, and sets the pointer *@gsdplsst to
 * NULL.
 *
 */
void
nc_galaxy_sd_position_lsst_srd_clear (NcGalaxySDPositionLSSTSRD **gsdplsst)
{
  g_clear_object (gsdplsst);
}


1	/***************************************************************************
2	* nc_galaxy_sd_position_lsst_srd.c
3	*
4	* Tue June 22 14:59:27 2023
5	* Copyright 2023 Caio Lima de Oliveira
6	* <caiolimadeoliveira@pm.me>
7	****************************************************************************/
8	/*
9	* nc_galaxy_sd_position_lsst_srd.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_galaxy_sd_position_lsst_srd
28	* @title: NcGalaxySDPositionLSSTSRD
29	* @short_description: Class describing galaxy sample position distributions as in LSST-SRD
30	* @stability: Unstable
31	*
32	* Class defining a galaxy sample position distribution as described in the LSST
33	* Science Roadmap Document, it includes the redshift distribution and a angular radius
34	* distribution.
35	*
36	*/
37
38	#ifdef HAVE_CONFIG_H
39	# include "config.h"
40	#endif /* HAVE_CONFIG_H */
41	#include "build_cfg.h"
42
43	/* #include "nc_enum_types.h" */
44	#include "galaxy/nc_galaxy_sd_position_lsst_srd.h"
45	#include "galaxy/nc_galaxy_sd_position.h"
46	#include "math/ncm_stats_dist1d.h"
47	#include "math/ncm_stats_dist1d_spline.h"
48	#include "math/ncm_dtuple.h"
49	#include "math/ncm_spline.h"
50	#include "math/ncm_spline_func.h"
51	#include "math/ncm_spline_cubic.h"
52	#include "math/ncm_spline_cubic_notaknot.h"
53	#include <math.h>
54	#include <gsl/gsl_math.h>
55	#include <gsl/gsl_randist.h>
56
57	typedef struct _NcGalaxySDPositionLSSTSRDPrivate
58	{
59	gdouble z_lb;
60	gdouble z_ub;
61	gdouble r_norm;
62	gdouble r_lb;
63	gdouble r_ub;
64	gdouble r_lb2;
65	gdouble r_ub2;
66	} NcGalaxySDPositionLSSTSRDPrivate;
67
68	struct _NcGalaxySDPositionLSSTSRD
69	{
70	NcGalaxySDPosition parent_instance;
71	};
72
73	enum
74	{
75	PROP_0,
76	PROP_LEN,
77	};
78
79	G_DEFINE_TYPE_WITH_PRIVATE (NcGalaxySDPositionLSSTSRD, nc_galaxy_sd_position_lsst_srd, NC_TYPE_GALAXY_SD_POSITION);	9,726,811✔
80
81	static void
82	nc_galaxy_sd_position_lsst_srd_init (NcGalaxySDPositionLSSTSRD *gsdplsst)	5✔
83	{
84	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	5✔
85
86	self->z_lb = 0.0;	5✔
87	self->z_ub = 0.0;	5✔
88	self->r_norm = 0.0;	5✔
89	self->r_lb = 0.0;	5✔
90	self->r_ub = 0.0;	5✔
91	self->r_lb2 = 0.0;	5✔
92	self->r_ub2 = 0.0;	5✔
93	}	5✔
94
95	static void
96	_nc_galaxy_sd_position_lsst_srd_dispose (GObject *object)	5✔
97	{
98	/* Chain up : end */
99	G_OBJECT_CLASS (nc_galaxy_sd_position_lsst_srd_parent_class)->dispose (object);	5✔
100	}	5✔
101
102	static void
103	_nc_galaxy_sd_position_lsst_srd_finalize (GObject *object)	5✔
104	{
105	/* Chain up : end */
106	G_OBJECT_CLASS (nc_galaxy_sd_position_lsst_srd_parent_class)->finalize (object);	5✔
107	}	5✔
108
109	static gdouble _nc_galaxy_sd_position_lsst_srd_gen_r (NcGalaxySDPosition gsdp, NcmRNG rng);
110	static gdouble _nc_galaxy_sd_position_lsst_srd_gen_z (NcGalaxySDPosition gsdp, NcmRNG rng);
111	static gdouble _nc_galaxy_sd_position_lsst_srd_integ (NcGalaxySDPosition *gsdp, const gdouble r, const gdouble z);
112	static void _nc_galaxy_sd_position_lsst_srd_set_z_lim (NcGalaxySDPosition *gsdp, gdouble z_min, gdouble z_max);
113	static void _nc_galaxy_sd_position_lsst_srd_get_z_lim (NcGalaxySDPosition gsdp, gdouble z_min, gdouble *z_max);
114	static void _nc_galaxy_sd_position_lsst_srd_set_r_lim (NcGalaxySDPosition *gsdp, gdouble r_min, gdouble r_max);
115	static void _nc_galaxy_sd_position_lsst_srd_get_r_lim (NcGalaxySDPosition gsdp, gdouble r_min, gdouble *r_max);
116
117	static void
118	nc_galaxy_sd_position_lsst_srd_class_init (NcGalaxySDPositionLSSTSRDClass *klass)	114✔
119	{
120	NcGalaxySDPositionClass *sd_position_class = NC_GALAXY_SD_POSITION_CLASS (klass);	114✔
121	GObjectClass *object_class = G_OBJECT_CLASS (klass);	114✔
122	NcmModelClass *model_class = NCM_MODEL_CLASS (klass);	114✔
123
124	object_class->dispose = &_nc_galaxy_sd_position_lsst_srd_dispose;	114✔
125	object_class->finalize = &_nc_galaxy_sd_position_lsst_srd_finalize;	114✔
126
127	ncm_model_class_set_name_nick (model_class, "LSST SRD Galaxy Distribution", "LSST_SRD");	114✔
128	ncm_model_class_add_params (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_SPARAM_LEN, 0, PROP_LEN);	114✔
129
130	/**
131	* NcGalaxySDPositionLSSTSRD:alpha:
132	*
133	* The redshift exponential slope.
134	*
135	*/
136	ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_ALPHA, "\\alpha", "alpha",	114✔
137	1e-8, 1.0, 1.0e-2,
138	NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_ALPHA,
139	NCM_PARAM_TYPE_FIXED);
140
141	/**
142	* NcGalaxySDPositionLSSTSRD:beta:
143	*
144	* The redshift power law slope.
145	*
146	*/
147	ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_BETA, "\\beta", "beta",	114✔
148	1e-8, 5.0, 1.0e-1,
149	NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_BETA,
150	NCM_PARAM_TYPE_FIXED);
151
152	/**
153	* NcGalaxySDPositionLSSTSRD:z0:
154	*
155	* The redshift pivot.
156	*
157	*/
158	ncm_model_class_set_sparam (model_class, NC_GALAXY_SD_POSITION_LSST_SRD_Z0, "z_0", "z0",	114✔
159	1e-8, 10.0, 1.0e-2,
160	NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_PARAMS_ABSTOL, NC_GALAXY_SD_POSITION_LSST_SRD_DEFAULT_Z0,
161	NCM_PARAM_TYPE_FIXED);
162
163	ncm_model_class_check_params_info (model_class);	114✔
164
165	sd_position_class->gen_r = &_nc_galaxy_sd_position_lsst_srd_gen_r;	114✔
166	sd_position_class->gen_z = &_nc_galaxy_sd_position_lsst_srd_gen_z;	114✔
167	sd_position_class->integ = &_nc_galaxy_sd_position_lsst_srd_integ;	114✔
168	sd_position_class->set_z_lim = &_nc_galaxy_sd_position_lsst_srd_set_z_lim;	114✔
169	sd_position_class->get_z_lim = &_nc_galaxy_sd_position_lsst_srd_get_z_lim;	114✔
170	sd_position_class->set_r_lim = &_nc_galaxy_sd_position_lsst_srd_set_r_lim;	114✔
171	sd_position_class->get_r_lim = &_nc_galaxy_sd_position_lsst_srd_get_r_lim;	114✔
172	}	114✔
173
174	#define VECTOR (NCM_MODEL (gsdp))
175	#define ALPHA (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_ALPHA))
176	#define BETA (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_BETA))
177	#define Z0 (ncm_model_orig_param_get (VECTOR, NC_GALAXY_SD_POSITION_LSST_SRD_Z0))
178
179	static gdouble
180	_nc_galaxy_sd_position_lsst_srd_gen_r (NcGalaxySDPosition gsdp, NcmRNG rng)	1,102,200✔
181	{
182	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	1,102,200✔
183	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	1,102,200✔
184	const gdouble cumul_gen = ncm_rng_uniform_gen (rng, 0.0, 1.0);	1,102,200✔
185
186	return sqrt (cumul_gen * 2.0 / self->r_norm + self->r_lb2);	1,102,200✔
187	}
188
189	static gdouble
190	_nc_galaxy_sd_position_lsst_srd_gen_z (NcGalaxySDPosition gsdp, NcmRNG rng)	1,104,264✔
191	{
192	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	1,104,264✔
193	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	1,104,264✔
194	const gdouble alpha = ALPHA;	1,104,264✔
195	const gdouble beta = BETA;	1,104,264✔
196	const gdouble z0 = Z0;	1,104,264✔
197	const gdouble y0 = pow (z0, alpha);	1,104,264✔
198	const gdouble gamma_a = (1.0 + beta) / alpha;	1,104,264✔
199
200	gdouble z;	1,357,236✔
201
202	do {	1,357,236✔
203	const gdouble gen_y = ncm_rng_gamma_gen (rng, gamma_a, y0);	1,357,236✔
204
205	z = pow (gen_y, 1.0 / alpha);	1,357,236✔
206	} while (z < self->z_lb \|\| z > self->z_ub);	1,357,236✔
207
208	return z;	1,104,264✔
209	}
210
211	static gdouble
212	_nc_galaxy_sd_position_lsst_srd_integ (NcGalaxySDPosition *gsdp, const gdouble r, const gdouble z)	7,475,627✔
213	{
214	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	7,475,627✔
215	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	7,475,627✔
216	const gdouble alpha = ALPHA;	7,475,627✔
217	const gdouble beta = BETA;	7,475,627✔
218	const gdouble z0 = Z0;	7,475,627✔
219	const gdouble y0 = pow (z0, alpha);	7,475,627✔
220	const gdouble gamma_a = (1.0 + beta) / alpha;	7,475,627✔
221	const gdouble y = pow (z, alpha);	7,475,627✔
222
223	return gsl_ran_gamma_pdf (y, gamma_a, y0) * alpha * y / z * r * self->r_norm;	7,475,627✔
224	}
225
226	static void
227	_nc_galaxy_sd_position_lsst_srd_set_z_lim (NcGalaxySDPosition *gsdp, gdouble z_min, gdouble z_max)	11,006✔
228	{
229	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	11,006✔
230	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	11,006✔
231
232	g_assert_cmpfloat (z_min, <, z_max);	11,006✔
233
234	self->z_lb = z_min;	11,006✔
235	self->z_ub = z_max;	11,006✔
236	}	11,006✔
237
238	static void
239	_nc_galaxy_sd_position_lsst_srd_get_z_lim (NcGalaxySDPosition gsdp, gdouble z_min, gdouble *z_max)	11,003✔
240	{
241	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	11,003✔
242	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	11,003✔
243
244	g_assert_nonnull (z_min);	11,003✔
245	g_assert_nonnull (z_max);	11,003✔
246
247	*z_min = self->z_lb;	11,003✔
248	*z_max = self->z_ub;	11,003✔
249	}	11,003✔
250
251	static void
252	_nc_galaxy_sd_position_lsst_srd_set_r_lim (NcGalaxySDPosition *gsdp, gdouble r_min, gdouble r_max)	11,006✔
253	{
254	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	11,006✔
255	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	11,006✔
256
257	g_assert_cmpfloat (r_min, <, r_max);	11,006✔
258
259	self->r_lb = r_min;	11,006✔
260	self->r_ub = r_max;	11,006✔
261
262	self->r_lb2 = self->r_lb * self->r_lb;	11,006✔
263	self->r_ub2 = self->r_ub * self->r_ub;	11,006✔
264	self->r_norm = 2.0 / (self->r_ub2 - self->r_lb2);	11,006✔
265	}	11,006✔
266
267	static void
268	_nc_galaxy_sd_position_lsst_srd_get_r_lim (NcGalaxySDPosition gsdp, gdouble r_min, gdouble *r_max)	11,000✔
269	{
270	NcGalaxySDPositionLSSTSRD *gsdplsst = NC_GALAXY_SD_POSITION_LSST_SRD (gsdp);	11,000✔
271	NcGalaxySDPositionLSSTSRDPrivate * const self = nc_galaxy_sd_position_lsst_srd_get_instance_private (gsdplsst);	11,000✔
272
273	g_assert_nonnull (r_min);	11,000✔
274	g_assert_nonnull (r_max);	11,000✔
275
276	*r_min = self->r_lb;	11,000✔
277	*r_max = self->r_ub;	11,000✔
278	}	11,000✔
279
280	/**
281	* nc_galaxy_sd_position_lsst_srd_new:
282	* @z_min: the minimum redshift
283	* @z_max: the maximum redshift
284	* @r_min: the minimum angular radius
285	* @r_max: the maximum angular radius
286	*
287	* Creates a new #NcGalaxySDPositionLSSTSRD, the parameter values correspond to the
288	* LSST SRD year 1.
289	*
290	* Returns: (transfer full): a new NcGalaxySDPositionLSSTSRD.
291	*/
292	NcGalaxySDPositionLSSTSRD *
293	nc_galaxy_sd_position_lsst_srd_new (const gdouble z_min, const gdouble z_max, const gdouble r_min, const gdouble r_max)	5✔
294	{
295	NcmDTuple2 z_lim = NCM_DTUPLE2_STATIC_INIT (z_min, z_max);	5✔
296	NcmDTuple2 r_lim = NCM_DTUPLE2_STATIC_INIT (r_min, r_max);	5✔
297	NcGalaxySDPositionLSSTSRD *gsdplsst = g_object_new (NC_TYPE_GALAXY_SD_POSITION_LSST_SRD,	5✔
298	"z-lim", &z_lim,
299	"r-lim", &r_lim,
300	NULL);
301
302	return gsdplsst;	5✔
303	}
304
305	/**
306	* nc_galaxy_sd_position_lsst_srd_new_y10:
307	* @z_min: the minimum redshift
308	* @z_max: the maximum redshift
309	* @r_min: the minimum angular radius
310	* @r_max: the maximum angular radius
311	*
312	* Creates a new #NcGalaxySDPositionLSSTSRD, the parameter values correspond to the
313	* LSST SRD year 10.
314	*
315	* Returns: (transfer full): a new NcGalaxySDPositionLSSTSRD.
316	*/
317	NcGalaxySDPositionLSSTSRD *
NEW 318	nc_galaxy_sd_position_lsst_srd_new_y10 (const gdouble z_min, const gdouble z_max, const gdouble r_min, const gdouble r_max)	×
319	{
NEW 320	NcmDTuple2 z_lim = NCM_DTUPLE2_STATIC_INIT (z_min, z_max);	×
NEW 321	NcmDTuple2 r_lim = NCM_DTUPLE2_STATIC_INIT (r_min, r_max);	×
NEW 322	NcGalaxySDPositionLSSTSRD *gsdplsst = g_object_new (NC_TYPE_GALAXY_SD_POSITION_LSST_SRD,	×
323	"alpha", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_ALPHA,
324	"beta", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_BETA,
325	"z0", NC_GALAXY_SD_POSITION_LSST_SRD_Y10_Z0,
326	"z-lim", &z_lim,
327	"r-lim", &r_lim,
328	NULL);
329
NEW 330	return gsdplsst;	×
331	}
332
333	/**
334	* nc_galaxy_sd_position_lsst_srd_ref:
335	* @gsdplsst: a #NcGalaxySDPositionLSSTSRD
336	*
337	* Increase the reference of @gsdplsst by one.
338	*
339	* Returns: (transfer full): @gsdplsst.
340	*/
341	NcGalaxySDPositionLSSTSRD *
NEW 342	nc_galaxy_sd_position_lsst_srd_ref (NcGalaxySDPositionLSSTSRD *gsdplsst)	×
343	{
NEW 344	return g_object_ref (gsdplsst);	×
345	}
346
347	/**
348	* nc_galaxy_sd_position_lsst_srd_free:
349	* @gsdplsst: a #NcGalaxySDPositionLSSTSRD
350	*
351	* Decrease the reference count of @gsdplsst by one.
352	*
353	*/
354	void
355	nc_galaxy_sd_position_lsst_srd_free (NcGalaxySDPositionLSSTSRD *gsdplsst)	2✔
356	{
357	g_object_unref (gsdplsst);	2✔
358	}	2✔
359
360	/**
361	* nc_galaxy_sd_position_lsst_srd_clear:
362	* @gsdplsst: a #NcGalaxySDPositionLSSTSRD
363	*
364	* Decrease the reference count of @gsdplsst by one, and sets the pointer *@gsdplsst to
365	* NULL.
366	*
367	*/
368	void
NEW 369	nc_galaxy_sd_position_lsst_srd_clear (NcGalaxySDPositionLSSTSRD **gsdplsst)	×
370	{
NEW 371	g_clear_object (gsdplsst);	×
372	}
373

NumCosmo / NumCosmo / 7176062912

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