17 #ifndef __deal2__quadrature_lib_h
18 #define __deal2__quadrature_lib_h
21 #include <deal.II/base/config.h>
22 #include <deal.II/base/quadrature.h>
48 QGauss (
const unsigned int n);
102 std::vector<long double>
105 const int beta)
const;
116 std::vector<long double>
119 const int beta)
const;
131 long double JacobiP(
const long double x,
134 const unsigned int n)
const;
141 long double gamma(
const unsigned int n)
const;
227 const bool revert=
false);
310 const double alpha = 1,
311 const bool factor_out_singular_weight=
false);
387 const bool factor_out_singular_weight=
false);
425 const unsigned int vertex_index,
426 const bool factor_out_singular_weight=
false);
433 const unsigned int n);
452 JacobiP(
const long double,
const int,
const int,
const unsigned int)
const;
472 DEAL_II_NAMESPACE_CLOSE
QGaussLog(const unsigned int n, const bool revert=false)
QGaussOneOverR(const unsigned int n, const Point< dim > singularity, const bool factor_out_singular_weight=false)
long double gamma(const unsigned int n) const
long double JacobiP(const long double x, const int alpha, const int beta, const unsigned int n) const
std::vector< long double > compute_quadrature_points(const unsigned int q, const int alpha, const int beta) const
QGaussLobatto(const unsigned int n)
std::vector< long double > compute_quadrature_weights(const std::vector< long double > &x, const int alpha, const int beta) const
QGauss(const unsigned int n)
std::vector< double > set_quadrature_weights(const unsigned int n) const
QGaussLogR(const unsigned int n, const Point< dim > x0=Point< dim >(), const double alpha=1, const bool factor_out_singular_weight=false)
static unsigned int quad_size(const Point< dim > singularity, const unsigned int n)
std::vector< double > set_quadrature_points(const unsigned int n) const