Top | ![]() |
![]() |
![]() |
![]() |
Functions
Types and Values
struct | NcmIntegrand2dim |
struct | NcmIntegrand3dim |
struct | NcmIntegralFixed |
#define | NCM_INTEGRAL_PARTITION |
#define | NCM_INTEGRAL_ALG |
#define | NCM_INTEGRAL_ERROR |
#define | NCM_INTEGRAL_ABS_ERROR |
Functions
ncm_integral_get_workspace ()
gsl_integration_workspace **
ncm_integral_get_workspace (void
);
This function provides a workspace to be used by numerical integration functions of GSL. It keeps a internal pool of workspaces and allocate a new one if the function is called and the pool is empty. It is designed to be used in a multithread enviroment. The workspace must be unlocked in order to return to the pool. This must be done using the ncm_memory_pool_return.
[skip]
_NcmIntegrand2dimFunc ()
gdouble (*_NcmIntegrand2dimFunc) (gdouble x
,gdouble y
,gpointer userdata
);
NcmIntegralPeakfinder ()
void (*NcmIntegralPeakfinder) (const gint *ndim
,const gdouble b[]
,gint *n
,gdouble x[]
,void *userdata
);
_NcmIntegrand3dimFunc ()
gdouble (*_NcmIntegrand3dimFunc) (gdouble x
,gdouble y
,gdouble z
,gpointer userdata
);
ncm_integral_locked_a_b ()
gint ncm_integral_locked_a_b (gsl_function *F
,gdouble a
,gdouble b
,gdouble abstol
,gdouble reltol
,gdouble *result
,gdouble *error
);
This function uses a workspace from the pool and gsl_integration_qag function to perform the numerical integration in the [a, b] interval.
[skip]
Parameters
F |
a gsl_function wich is the integrand. |
|
a |
lower integration limit. |
|
b |
upper integration limit. |
|
abstol |
absolute tolerance. |
|
reltol |
relative tolerance. |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integral_locked_a_inf ()
gint ncm_integral_locked_a_inf (gsl_function *F
,gdouble a
,gdouble abstol
,gdouble reltol
,gdouble *result
,gdouble *error
);
This function uses a workspace from the pool and gsl_integration_qagiu function to perform the numerical integration in the $[a, \infty]$ interval.
[skip]
ncm_integral_cached_0_x ()
gint ncm_integral_cached_0_x (NcmFunctionCache *cache
,gsl_function *F
,gdouble x
,gdouble *result
,gdouble *error
);
This function searchs for the nearest x_near value previously chosed as the upper integration limit and perform the integration at [x_near, x] interval. This result is summed to that obtained at [0, x_near] and then it is saved in the cache.
[skip]
Parameters
cache |
a pointer to NcmFunctionCache. |
|
F |
a gsl_function wich is the integrand. |
|
x |
upper integration limit. |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integral_cached_x_inf ()
gint ncm_integral_cached_x_inf (NcmFunctionCache *cache
,gsl_function *F
,gdouble x
,gdouble *result
,gdouble *error
);
This function searchs for the nearest x_near value previously chosed as the lower integration limit and perform the integration at $[x, x_{near}]$ interval. This result is summed to that obtained at $[x_{near}, \infty]$ and then it is saved in the cache.
[skip]
Parameters
cache |
a pointer to NcmFunctionCache. |
|
F |
a gsl_function wich is the integrand. |
|
x |
lower integration limit. |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_2dim ()
gboolean ncm_integrate_2dim (NcmIntegrand2dim *integ
,gdouble xi
,gdouble yi
,gdouble xf
,gdouble yf
,gdouble epsrel
,gdouble epsabs
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand2dim. |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_2dim_divonne ()
gboolean ncm_integrate_2dim_divonne (NcmIntegrand2dim *integ
,gdouble xi
,gdouble yi
,gdouble xf
,gdouble yf
,gdouble epsrel
,gdouble epsabs
,const gint ngiven
,const gint ldxgiven
,gdouble xgiven[]
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand2dim |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
ngiven |
number of peaks |
|
ldxgiven |
the leading dimension of xgiven, i.e. the offset between one point and the next in memory (ref. libcuba documentation) |
|
xgiven |
list of points where the integrand might have peaks (ref. libcuba documentation) |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_2dim_divonne_peakfinder ()
gboolean ncm_integrate_2dim_divonne_peakfinder (NcmIntegrand2dim *integ
,gdouble xi
,gdouble yi
,gdouble xf
,gdouble yf
,gdouble epsrel
,gdouble epsabs
,const gint ngiven
,const gint ldxgiven
,gdouble xgiven[]
,const gint nextra
,NcmIntegralPeakfinder peakfinder
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand2dim |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
ngiven |
number of peaks |
|
ldxgiven |
the leading dimension of xgiven, i.e. the offset between one point and the next in memory (ref. libcuba documentation) |
|
xgiven |
list of points where the integrand might have peaks (ref. libcuba documentation) |
|
nextra |
the maximum number of extra points the peakfinder will return. |
|
peakfinder |
[scope call] | |
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_2dim_vegas ()
gboolean ncm_integrate_2dim_vegas (NcmIntegrand2dim *integ
,gdouble xi
,gdouble yi
,gdouble xf
,gdouble yf
,gdouble epsrel
,gdouble epsabs
,const gint nstart
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand2dim |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
nstart |
number of samples to start the first round of integration with. |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_3dim ()
gboolean ncm_integrate_3dim (NcmIntegrand3dim *integ
,gdouble xi
,gdouble yi
,gdouble zi
,gdouble xf
,gdouble yf
,gdouble zf
,gdouble epsrel
,gdouble epsabs
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand3dim. |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
zi |
gbouble which is the lower integration limit of variable z. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
zf |
gbouble which is the upper integration limit of variable z. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_3dim_divonne ()
gboolean ncm_integrate_3dim_divonne (NcmIntegrand3dim *integ
,gdouble xi
,gdouble yi
,gdouble zi
,gdouble xf
,gdouble yf
,gdouble zf
,gdouble epsrel
,gdouble epsabs
,const gint ngiven
,const gint ldxgiven
,gdouble xgiven[]
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand3dim |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
zi |
gbouble which is the lower integration limit of variable z. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
zf |
gbouble which is the upper integration limit of variable z. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
ngiven |
number of peaks |
|
ldxgiven |
the leading dimension of xgiven, i.e. the offset between one point and the next in memory (ref. libcuba documentation) |
|
xgiven |
list of points where the integrand might have peaks (ref. libcuba documentation) |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integrate_3dim_vegas ()
gboolean ncm_integrate_3dim_vegas (NcmIntegrand3dim *integ
,gdouble xi
,gdouble yi
,gdouble zi
,gdouble xf
,gdouble yf
,gdouble zf
,gdouble epsrel
,gdouble epsabs
,const gint nstart
,gdouble *result
,gdouble *error
);
This function FIXME
Parameters
integ |
a pointer to NcmIntegrand3dim |
|
xi |
gbouble which is the lower integration limit of variable x. |
|
yi |
gbouble which is the lower integration limit of variable y. |
|
zi |
gbouble which is the lower integration limit of variable z. |
|
xf |
gbouble which is the upper integration limit of variable x. |
|
yf |
gbouble which is the upper integration limit of variable y. |
|
zf |
gbouble which is the upper integration limit of variable z. |
|
epsrel |
relative error |
|
epsabs |
absolute error |
|
nstart |
number of samples to start the first round of integration with. |
|
result |
a pointer to a gdouble in which the function stores the result. |
|
error |
a pointer to a gdouble in which the function stores the estimated error. |
ncm_integral_fixed_new ()
NcmIntegralFixed * ncm_integral_fixed_new (gulong n_nodes
,gulong rule_n
,gdouble xl
,gdouble xu
);
This function prepares the NcmIntegralFixed with a grid with n_nodes - 1 intervals beteween xl and xu. In each interval it uses a fixed order (rule_n) Gauss-Legendre integration rule to determine the interval inner points. This results in a grid with (n_nodes - 1) * rule_n points.
[skip]
ncm_integral_fixed_free ()
void
ncm_integral_fixed_free (NcmIntegralFixed *intf
);
This function frees the memory associated to NcmIntegralFixed.
ncm_integral_fixed_calc_nodes ()
void ncm_integral_fixed_calc_nodes (NcmIntegralFixed *intf
,gsl_function *F
);
This function FIXME
[skip]
ncm_integral_fixed_nodes_eval ()
gdouble
ncm_integral_fixed_nodes_eval (NcmIntegralFixed *intf
);
This function
ncm_integral_fixed_integ_mult ()
gdouble ncm_integral_fixed_integ_mult (NcmIntegralFixed *intf
,gsl_function *F
);
This function
[skip]
ncm_integral_fixed_integ_posdef_mult ()
gdouble ncm_integral_fixed_integ_posdef_mult (NcmIntegralFixed *intf
,gsl_function *F
,gdouble max
,gdouble reltol
);
This function
[skip]