Class inte_gauss56_cern (o2scl)

O2scl : Class List

template<class func_t = funct, class fp_t = double, class weights_t = inte_gauss56_coeffs_double>
class o2scl::inte_gauss56_cern : public o2scl::inte<funct, double>

5,6-point Gaussian quadrature (CERNLIB)

If \( I_5 \) is the 5-point approximation, and \( I_6 \) is the 6-point approximation to the integral, then integ_err() returns the result \( \frac{1}{2}(I_5+I_6) \) with uncertainty \( |I_5-I_6| \).

This class is based on the CERNLIB routines RGS56P and DGS56P were originally documented at http://wwwasdoc.web.cern.ch/wwwasdoc/shortwrupsdir/d106/top.html . (3/10/2020: The CERNLIB links are apparently dead and haven’t been supported since 2003.)

Public Functions

inte_gauss56_cern()
int integ_err(func_t &func, fp_t a, fp_t b, fp_t &res, fp_t &err)

Integrate function func from a to b giving result res and error err.

This function always returns success.

Protected Attributes

const fp_t *w5
const fp_t *x5
const fp_t *w6
const fp_t *x6
weights_t wgts