62 template<u
int dimension>
65 ASSERT_LT(order, quads.size()).error(
"Quadrature of given order is not implemented.");
66 auto &point_list = quads[order];
70 for (
int i=0; i<point_list->npoints; i++)
#define ASSERT_LT(a, b)
Definition of comparative assert macro (Less Than) only for debug mode.
#define ASSERT_EQ(a, b)
Definition of comparative assert macro (EQual) only for debug mode.
void append(const ArmaMat< Type, nr, nc > &item)
QGauss(unsigned int dim, const unsigned int order)
Create a formula of given order.
Base class for quadrature rules on simplices in arbitrary dimensions.
unsigned int size() const
Returns number of quadrature points.
std::vector< double > weights
List of weights to the quadrature points.
Armor::Array< double > quadrature_points
List of quadrature points.
Global macros to enhance readability and debugging, general constants.
typename arma::Col< Type >::template fixed< nr > ArmaVec
double __unit_cell_volume[]
std::vector< DimQuadList > __gauss_quadratures
std::vector< QUAD * > DimQuadList
Definitions of particular quadrature rules on simplices.