basix::quadrature Namespace Reference

Quadrature rules. More...

Enumerations
enum	type {   Default = 0, gauss_jacobi = 1, gll = 2, xiao_gimbutas = 3,   zienkiewicz_taylor = 20, keast = 21, strang_fix = 22 }
	Quadrature type.

Functions
std::array< std::vector< double >, 2 >	make_quadrature (const quadrature::type rule, cell::type celltype, int m)
std::array< std::vector< double >, 2 >	make_quadrature (cell::type celltype, int m)
quadrature::type	get_default_rule (cell::type celltype, int m)
std::vector< double >	get_gll_points (int m)
std::vector< double >	get_gl_points (int m)

Detailed Description

Quadrature rules.

Function Documentation

get_default_rule()

quadrature::type basix::quadrature::get_default_rule	(	cell::type	celltype,
		int	m
	)

Get the default quadrature type for the given cell and order

Parameters

[in]	celltype	The cell type
[in]	m	Maximum degree of polynomial that this quadrature rule will integrate exactly

Returns

The quadrature type that will be used by default

get_gl_points()

std::vector< double > basix::quadrature::get_gl_points

(

int

m

)

Get Gauss-Legendre (GL) points on the interval [0, 1].

Parameters

[in]

m

The number of points

Returns

An array of GL points. Shape is (num points, gdim)

get_gll_points()

std::vector< double > basix::quadrature::get_gll_points

(

int

m

)

Get Gauss-Lobatto-Legendre (GLL) points on the interval [0, 1].

Parameters

[in]

m

The number of points

Returns

An array of GLL points. Shape is (num points, gdim)

make_quadrature() [1/2]

std::array< std::vector< double >, 2 > basix::quadrature::make_quadrature	(	cell::type	celltype,
		int	m
	)

Make a default quadrature rule on reference cell

Parameters

[in]	celltype	The cell type
[in]	m	Maximum degree of polynomial that this quadrature rule will integrate exactly

Returns

List of points and list of weights. The number of points arrays has shape (num points, gdim)

make_quadrature() [2/2]

std::array< std::vector< double >, 2 > basix::quadrature::make_quadrature	(	const quadrature::type	rule,
		cell::type	celltype,
		int	m
	)

Make a quadrature rule on a reference cell

Parameters

[in]	rule	Type of quadrature rule (or use quadrature::Default)
[in]	celltype	The cell type
[in]	m	Maximum degree of polynomial that this quadrature rule will integrate exactly

Returns

List of points and list of weights. The number of points arrays has shape (num points, gdim)