Public Member Functions | |
QGaussLog (const unsigned int n, const bool revert=false) | |
template<> | |
std::vector< double > | set_quadrature_points (const unsigned int) const |
template<> | |
std::vector< double > | set_quadrature_weights (const unsigned int) const |
template<> | |
QGaussLog (const unsigned int n, const bool revert) | |
Protected Member Functions | |
std::vector< double > | set_quadrature_points (const unsigned int n) const |
std::vector< double > | set_quadrature_weights (const unsigned int n) const |
ln|x|*f(x)
on the interval [0,1]
, where f is a smooth function without singularities. The collection of quadrature points and weights has been obtained using Numerical Recipes
.
Notice that only the function f(x)
should be provided, i.e., . Setting the
revert
flag to true at construction time switches the weight from ln|x|
to ln|1-x|
.
The weights and functions have been tabulated up to order 12.
Generate a formula with n
quadrature points
std::vector<double> QGaussLog< dim >::set_quadrature_points | ( | const unsigned int | n | ) | const [protected] |
Sets the points of the quadrature formula.
std::vector<double> QGaussLog< dim >::set_quadrature_weights | ( | const unsigned int | n | ) | const [protected] |
Sets the weights of the quadrature formula.
std::vector< double > QGaussLog< 1 >::set_quadrature_points | ( | const unsigned | int | ) | const [inline] |
std::vector< double > QGaussLog< 1 >::set_quadrature_weights | ( | const unsigned | int | ) | const [inline] |