42 std::shared_ptr<dealii::Quadrature<dim>> quadrature;
43 if(element_type == ElementType::Hypercube)
45 quadrature = std::make_shared<dealii::QGauss<dim>>(n_q_points_1d);
47 else if(element_type == ElementType::Simplex)
49 quadrature = std::make_shared<dealii::QGaussSimplex<dim>>(n_q_points_1d);
53 AssertThrow(
false, ExcNotImplemented());
std::shared_ptr< dealii::Quadrature< dim > > create_quadrature(ElementType const &element_type, unsigned int const n_q_points_1d)
Definition quadrature.h:40