![]() |
![]() |
![]() |
NumCosmo Reference Manual | ![]() |
---|---|---|---|---|
Top | Description |
Synopsis
#define NCM_INTEGRAL_ABS_ERROR #define NCM_INTEGRAL_ALG #define NCM_INTEGRAL_ERROR #define NCM_INTEGRAL_PARTITION struct NcmIntegralFixed; struct NcmIntegrand2dim; gdouble (*_NcmIntegrand2dimFunc) (gdouble x
,gdouble y
,gpointer userdata
); gint ncm_integral_cached_0_x (NcmFunctionCache *cache
,gsl_function *F
,gdouble x
,gdouble *result
,gdouble *error
); gint ncm_integral_cached_x_inf (NcmFunctionCache *cache
,gsl_function *F
,gdouble x
,gdouble *result
,gdouble *error
); void ncm_integral_fixed_calc_nodes (NcmIntegralFixed *intf
,gsl_function *F
); void ncm_integral_fixed_free (NcmIntegralFixed *intf
); gdouble ncm_integral_fixed_integ_mult (NcmIntegralFixed *intf
,gsl_function *F
); gdouble ncm_integral_fixed_integ_posdef_mult (NcmIntegralFixed *intf
,gsl_function *F
,gdouble max
,gdouble reltol
); NcmIntegralFixed * ncm_integral_fixed_new (gulong n_nodes
,gulong rule_n
,gdouble xl
,gdouble xu
); gdouble ncm_integral_fixed_nodes_eval (NcmIntegralFixed *intf
); gsl_integration_workspace ** ncm_integral_get_workspace (); gint ncm_integral_locked_a_b (gsl_function *F
,gdouble a
,gdouble b
,gdouble abstol
,gdouble reltol
,gdouble *result
,gdouble *error
); gint ncm_integral_locked_a_inf (gsl_function *F
,gdouble a
,gdouble abstol
,gdouble reltol
,gdouble *result
,gdouble *error
); gboolean ncm_integrate_2dim (NcmIntegrand2dim *integ
,gdouble xi
,gdouble yi
,gdouble xf
,gdouble yf
,gdouble epsrel
,gdouble epsabs
,gdouble *result
,gdouble *error
);
Details
_NcmIntegrand2dimFunc ()
gdouble (*_NcmIntegrand2dimFunc) (gdouble x
,gdouble y
,gpointer userdata
);
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.
|
a pointer to NcmFunctionCache. |
|
a gsl_function wich is the integrand. |
|
upper integration limit. |
|
a pointer to a gdouble in which the function stores the result. |
|
a pointer to a gdouble in which the function stores the estimated error. |
Returns : |
the error code returned by gsl_integration_qag. |
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 \f$ [x, x_{near}] \f$ interval. This result is summed to that obtained at \f$ [x_{near}, \infty] \f$ and then it is saved in the cache.
|
a pointer to NcmFunctionCache. |
|
a gsl_function wich is the integrand. |
|
lower integration limit. |
|
a pointer to a gdouble in which the function stores the result. |
|
a pointer to a gdouble in which the function stores the estimated error. |
Returns : |
the error code returned by gsl_integration_qagiu. |
ncm_integral_fixed_calc_nodes ()
void ncm_integral_fixed_calc_nodes (NcmIntegralFixed *intf
,gsl_function *F
);
This function FIXME
|
a pointer to NcmIntegralFixed. |
|
a pointer to a gsl_function. |
ncm_integral_fixed_free ()
void ncm_integral_fixed_free (NcmIntegralFixed *intf
);
This function frees the memory associated to NcmIntegralFixed.
|
a pointer to NcmIntegralFixed. |
ncm_integral_fixed_integ_mult ()
gdouble ncm_integral_fixed_integ_mult (NcmIntegralFixed *intf
,gsl_function *F
);
This function
|
a pointer to NcmIntegralFixed. |
|
a pointer to gsl_function. |
Returns : |
FIXME |
ncm_integral_fixed_integ_posdef_mult ()
gdouble ncm_integral_fixed_integ_posdef_mult (NcmIntegralFixed *intf
,gsl_function *F
,gdouble max
,gdouble reltol
);
This function
|
a pointer to NcmIntegralFixed. |
|
a pointer to gsl_function. |
|
FIXME |
|
FIXME |
Returns : |
FIXME |
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.
|
number of nodes in the full interval. |
|
order of the Gauss-Legendre integration rule to be applied in each interval. |
|
the interval lower limit. |
|
the interval upper limit. |
Returns : |
a pointer to the newly created NcmIntegralFixed structure. |
ncm_integral_fixed_nodes_eval ()
gdouble ncm_integral_fixed_nodes_eval (NcmIntegralFixed *intf
);
This function
|
a pointer to NcmIntegralFixed. |
Returns : |
FIXME |
ncm_integral_get_workspace ()
gsl_integration_workspace ** ncm_integral_get_workspace ();
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.
Returns : |
a pointer to gsl_integration_workspace structure. |
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.
|
a gsl_function wich is the integrand. |
|
lower integration limit. |
|
upper integration limit. |
|
absolute tolerance. |
|
relative tolerance. |
|
a pointer to a gdouble in which the function stores the result. |
|
a pointer to a gdouble in which the function stores the estimated error. |
Returns : |
the error code returned by gsl_integration_qag. |
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 \f$ [a, \infty] \f$ interval.
|
a gsl_function wich is the integrand. |
|
lower integration limit. |
|
absolute tolerance. |
|
relative tolerance. |
|
a pointer to a gdouble in which the function stores the result. |
|
a pointer to a gdouble in which the function stores the estimated error. |
Returns : |
the error code returned by gsl_integration_qagiu. |
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
|
a pointer to NcmIntegrand2dim. |
|
gbouble which is the lower integration limit of variable x. |
|
gbouble which is the lower integration limit of variable y. |
|
gbouble which is the upper integration limit of variable x. |
|
gbouble which is the upper integration limit of variable y. |
|
relative error |
|
absolute error |
|
a pointer to a gdouble in which the function stores the result. |
|
a pointer to a gdouble in which the function stores the estimated error. |
Returns : |
a gboolean |