7250689656

Committed 18 Dec 2023 04:07PM UTC coverage: 55.112% (+0.1%) from 54.995%

Build # 7250689656

Build Type

Pull #126

github

Committed by

web-flow

Commit Message

Merge d240d76d3 into bffca6a86

Pull Request Pull Request #126: More objects encapsulation

Run Details

7259 of 12530 branches covered (0.0%)

Branch coverage included in aggregate %.

379 of 1021 new or added lines in 90 files covered. (37.12%)

33 existing lines in 10 files now uncovered.

47434 of 86709 relevant lines covered (54.7%)

1624823.85 hits per line

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

23.81

/numcosmo/math/ncm_mset_trans_kern_flat.c

/***************************************************************************
 *            ncm_mset_trans_kern_flat.c
 *
 *  Thu October 02 13:37:11 2014
 *  Copyright  2014  Sandro Dias Pinto Vitenti
 *  <vitenti@uel.br>
 ****************************************************************************/
/*
 * ncm_mset_trans_kern_flat.c
 * Copyright (C) 2014 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:ncm_mset_trans_kern_flat
 * @title: NcmMSetTransKernFlat
 * @short_description: Multivariate flat sampler
 *
 * This object subclasses NcmMSetTransKern and implements a multivariate flat sampler.
 *
 * Implementation of a multivariate flat sampler, offering a straightforward method for
 * generating random parameter vectors with multivariate parameters. This sampler
 * generates vectors uniformly distributed in the hypercube defined by the specified
 * bounds for each parameter.
 *
 * **Key Functionality:**
 *
 * - Generates random parameter vectors with multivariate parameters.
 * - Utilizes a simple flat sampling method within the hypercube defined by parameter
 *   bounds.
 *
 * This implementation is particularly useful when a basic flat sampling approach is
 * needed for generating random parameter vectors with multivariate parameters,
 * especially in scenarios where limited information is available about the posterior
 * distribution.
 *
 */

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

#include "math/ncm_mset_trans_kern_flat.h"
#include "math/ncm_c.h"

#ifndef NUMCOSMO_GIR_SCAN
#include <gsl/gsl_linalg.h>
#include <gsl/gsl_randist.h>
#endif /* NUMCOSMO_GIR_SCAN */

enum
{
  PROP_0,
};

G_DEFINE_TYPE (NcmMSetTransKernFlat, ncm_mset_trans_kern_flat, NCM_TYPE_MSET_TRANS_KERN)

static void
ncm_mset_trans_kern_flat_init (NcmMSetTransKernFlat *tkernf)
{
  tkernf->parea = 0.0;
}

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

static void _ncm_mset_trans_kern_flat_set_mset (NcmMSetTransKern *tkern, NcmMSet *mset);
static void _ncm_mset_trans_kern_flat_generate (NcmMSetTransKern *tkern, NcmVector *theta, NcmVector *thetastar, NcmRNG *rng);
static gdouble _ncm_mset_trans_kern_flat_pdf (NcmMSetTransKern *tkern, NcmVector *theta, NcmVector *thetastar);
static const gchar *_ncm_mset_trans_kern_flat_get_name (NcmMSetTransKern *tkern);

static void
ncm_mset_trans_kern_flat_class_init (NcmMSetTransKernFlatClass *klass)
{
  GObjectClass *object_class         = G_OBJECT_CLASS (klass);
  NcmMSetTransKernClass *tkern_class = NCM_MSET_TRANS_KERN_CLASS (klass);

  object_class->finalize = &_ncm_mset_trans_kern_flat_finalize;

  tkern_class->set_mset = &_ncm_mset_trans_kern_flat_set_mset;
  tkern_class->generate = &_ncm_mset_trans_kern_flat_generate;
  tkern_class->pdf      = &_ncm_mset_trans_kern_flat_pdf;
  tkern_class->get_name = &_ncm_mset_trans_kern_flat_get_name;
}

static void
_ncm_mset_trans_kern_flat_set_mset (NcmMSetTransKern *tkern, NcmMSet *mset)
{
  NcmMSetTransKernFlat *tkernf = NCM_MSET_TRANS_KERN_FLAT (tkern);
  guint fparam_len             = ncm_mset_fparam_len (tkern->mset);
  guint i;

  tkernf->parea = 1.0;

  for (i = 0; i < fparam_len; i++)
  {
    const gdouble lb = ncm_mset_fparam_get_lower_bound (tkern->mset, i);
    const gdouble ub = ncm_mset_fparam_get_upper_bound (tkern->mset, i);

    tkernf->parea *= ub - lb;
    g_assert_cmpfloat (tkernf->parea, >, 0.0);
  }
}

static void
_ncm_mset_trans_kern_flat_generate (NcmMSetTransKern *tkern, NcmVector *theta, NcmVector *thetastar, NcmRNG *rng)
{
  guint fparam_len = ncm_mset_fparam_len (tkern->mset);
  guint i;

  for (i = 0; i < fparam_len; i++)
  {
    const gdouble lb  = ncm_mset_fparam_get_lower_bound (tkern->mset, i);
    const gdouble ub  = ncm_mset_fparam_get_upper_bound (tkern->mset, i);
    const gdouble val = lb + (ub - lb) * ncm_rng_uniform01_pos_gen (rng);

    ncm_vector_set (thetastar, i, val);
  }
}

static gdouble
_ncm_mset_trans_kern_flat_pdf (NcmMSetTransKern *tkern, NcmVector *theta, NcmVector *thetastar)
{
  return 1.0 / NCM_MSET_TRANS_KERN_FLAT (tkern)->parea;
}

static const gchar *
_ncm_mset_trans_kern_flat_get_name (NcmMSetTransKern *tkern)
{
  return "Multivariate Flat Sampler";
}

/**
 * ncm_mset_trans_kern_flat_new:
 *
 * New NcmMSetTransKern flat.
 *
 * Returns: (transfer full): a new #NcmMSetTransKernFlat.
 *
 */
NcmMSetTransKernFlat *
ncm_mset_trans_kern_flat_new (void)
{
  NcmMSetTransKernFlat *tkernf = g_object_new (NCM_TYPE_MSET_TRANS_KERN_FLAT,
                                               NULL);

  return tkernf;
}


1	/***************************************************************************
2	* ncm_mset_trans_kern_flat.c
3	*
4	* Thu October 02 13:37:11 2014
5	* Copyright 2014 Sandro Dias Pinto Vitenti
6	* <vitenti@uel.br>
7	****************************************************************************/
8	/*
9	* ncm_mset_trans_kern_flat.c
10	* Copyright (C) 2014 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:ncm_mset_trans_kern_flat
28	* @title: NcmMSetTransKernFlat
29	* @short_description: Multivariate flat sampler
30	*
31	* This object subclasses NcmMSetTransKern and implements a multivariate flat sampler.
32	*
33	* Implementation of a multivariate flat sampler, offering a straightforward method for
34	* generating random parameter vectors with multivariate parameters. This sampler
35	* generates vectors uniformly distributed in the hypercube defined by the specified
36	* bounds for each parameter.
37	*
38	* Key Functionality:
39	*
40	* - Generates random parameter vectors with multivariate parameters.
41	* - Utilizes a simple flat sampling method within the hypercube defined by parameter
42	* bounds.
43	*
44	* This implementation is particularly useful when a basic flat sampling approach is
45	* needed for generating random parameter vectors with multivariate parameters,
46	* especially in scenarios where limited information is available about the posterior
47	* distribution.
48	*
49	*/
50
51	#ifdef HAVE_CONFIG_H
52	# include "config.h"
53	#endif /* HAVE_CONFIG_H */
54	#include "build_cfg.h"
55
56	#include "math/ncm_mset_trans_kern_flat.h"
57	#include "math/ncm_c.h"
58
59	#ifndef NUMCOSMO_GIR_SCAN
60	#include <gsl/gsl_linalg.h>
61	#include <gsl/gsl_randist.h>
62	#endif /* NUMCOSMO_GIR_SCAN */
63
64	enum
65	{
66	PROP_0,
67	};
68
69	G_DEFINE_TYPE (NcmMSetTransKernFlat, ncm_mset_trans_kern_flat, NCM_TYPE_MSET_TRANS_KERN)	12✔
70
71	static void
72	ncm_mset_trans_kern_flat_init (NcmMSetTransKernFlat *tkernf)	×
73	{
74	tkernf->parea = 0.0;	×
75	}
76
77	static void
78	_ncm_mset_trans_kern_flat_finalize (GObject *object)	×
79	{
80	/* Chain up : end */
81	G_OBJECT_CLASS (ncm_mset_trans_kern_flat_parent_class)->finalize (object);	×
82	}
83
84	static void _ncm_mset_trans_kern_flat_set_mset (NcmMSetTransKern tkern, NcmMSet mset);
85	static void _ncm_mset_trans_kern_flat_generate (NcmMSetTransKern tkern, NcmVector theta, NcmVector thetastar, NcmRNG rng);
86	static gdouble _ncm_mset_trans_kern_flat_pdf (NcmMSetTransKern tkern, NcmVector theta, NcmVector *thetastar);
87	static const gchar _ncm_mset_trans_kern_flat_get_name (NcmMSetTransKern tkern);
88
89	static void
90	ncm_mset_trans_kern_flat_class_init (NcmMSetTransKernFlatClass *klass)	2✔
91	{
92	GObjectClass *object_class = G_OBJECT_CLASS (klass);	2✔
93	NcmMSetTransKernClass *tkern_class = NCM_MSET_TRANS_KERN_CLASS (klass);	2✔
94
95	object_class->finalize = &_ncm_mset_trans_kern_flat_finalize;	2✔
96
97	tkern_class->set_mset = &_ncm_mset_trans_kern_flat_set_mset;	2✔
98	tkern_class->generate = &_ncm_mset_trans_kern_flat_generate;	2✔
99	tkern_class->pdf = &_ncm_mset_trans_kern_flat_pdf;	2✔
100	tkern_class->get_name = &_ncm_mset_trans_kern_flat_get_name;	2✔
101	}
102
103	static void
104	_ncm_mset_trans_kern_flat_set_mset (NcmMSetTransKern tkern, NcmMSet mset)	×
105	{
106	NcmMSetTransKernFlat *tkernf = NCM_MSET_TRANS_KERN_FLAT (tkern);	×
107	guint fparam_len = ncm_mset_fparam_len (tkern->mset);	×
108	guint i;	×
109
110	tkernf->parea = 1.0;	×
111
112	for (i = 0; i < fparam_len; i++)	×
113	{
114	const gdouble lb = ncm_mset_fparam_get_lower_bound (tkern->mset, i);	×
115	const gdouble ub = ncm_mset_fparam_get_upper_bound (tkern->mset, i);	×
116
117	tkernf->parea *= ub - lb;	×
118	g_assert_cmpfloat (tkernf->parea, >, 0.0);	×
119	}
120	}
121
122	static void
123	_ncm_mset_trans_kern_flat_generate (NcmMSetTransKern tkern, NcmVector theta, NcmVector thetastar, NcmRNG rng)	×
124	{
125	guint fparam_len = ncm_mset_fparam_len (tkern->mset);	×
126	guint i;	×
127
128	for (i = 0; i < fparam_len; i++)	×
129	{
130	const gdouble lb = ncm_mset_fparam_get_lower_bound (tkern->mset, i);	×
131	const gdouble ub = ncm_mset_fparam_get_upper_bound (tkern->mset, i);	×
NEW 132	const gdouble val = lb + (ub - lb) * ncm_rng_uniform01_pos_gen (rng);	×
133
134	ncm_vector_set (thetastar, i, val);	×
135	}
136	}
137
138	static gdouble
139	_ncm_mset_trans_kern_flat_pdf (NcmMSetTransKern tkern, NcmVector theta, NcmVector *thetastar)	×
140	{
141	return 1.0 / NCM_MSET_TRANS_KERN_FLAT (tkern)->parea;	×
142	}
143
144	static const gchar *
145	_ncm_mset_trans_kern_flat_get_name (NcmMSetTransKern *tkern)	×
146	{
147	return "Multivariate Flat Sampler";	×
148	}
149
150	/**
151	* ncm_mset_trans_kern_flat_new:
152	*
153	* New NcmMSetTransKern flat.
154	*
155	* Returns: (transfer full): a new #NcmMSetTransKernFlat.
156	*
157	*/
158	NcmMSetTransKernFlat *
159	ncm_mset_trans_kern_flat_new (void)	×
160	{
161	NcmMSetTransKernFlat *tkernf = g_object_new (NCM_TYPE_MSET_TRANS_KERN_FLAT,	×
162	NULL);
163
164	return tkernf;	×
165	}
166

NumCosmo / NumCosmo / 7250689656

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