16 #ifndef dealii_quadrature_lib_h
17 #define dealii_quadrature_lib_h
47 QGauss(
const unsigned int n);
247 QGaussLog(
const unsigned int n,
const bool revert =
false);
253 static std::vector<double>
259 static std::vector<double>
314 const double alpha = 1,
315 const bool factor_out_singular_weight =
false);
329 const
double fraction;
394 const bool factor_out_singular_weight =
false);
430 const unsigned int vertex_index,
431 const bool factor_out_singular_weight =
false);
471 compare_weights(
const unsigned int a,
const unsigned int b)
const;
602 const unsigned int n,
708 template <
int spacedim = dim>
710 compute_affine_transformation(
725 template <
int spacedim = dim>
728 const std::vector<std::array<
Point<spacedim>, dim + 1>> &simplices)
const;
822 const double beta = 1.0);
831 QDuffy(
const unsigned int n,
const double beta);
957 const bool use_odd_order =
true);
972 const unsigned int n_copies);
987 explicit QGaussWedge(
const unsigned int n_points_1D);
1044 const unsigned int index,
QGaussLobatto(const unsigned int n)
QGaussLogR(const unsigned int n, const Point< dim > &x0=Point< dim >(), const double alpha=1, const bool factor_out_singular_weight=false)
QGaussLogR(QGaussLogR< dim > &&) noexcept=default
QGaussLog(const unsigned int n, const bool revert=false)
static std::vector< double > get_quadrature_points(const unsigned int n)
static std::vector< double > get_quadrature_weights(const unsigned int n)
QGaussOneOverR(const unsigned int n, const unsigned int vertex_index, const bool factor_out_singular_weight=false)
QGaussOneOverR(const unsigned int n, const Point< dim > &singularity, const bool factor_out_singular_weight=false)
static unsigned int quad_size(const Point< dim > &singularity, const unsigned int n)
QGaussRadauChebyshev(QGaussRadauChebyshev< dim > &&) noexcept=default
QGaussRadau(QGaussRadau< dim > &&) noexcept=default
QGaussRadau(const unsigned int n, const EndPoint end_point=QGaussRadau::EndPoint::left)
QGauss(const unsigned int n)
QTelles(const Quadrature< 1 > &base_quad, const Point< dim > &singularity)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
void split_point(const Point< dim1+dim2 > &source, Point< dim1 > &p1, Point< dim2 > &p2)