Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion FIAT/quadrature_schemes.py
Original file line number Diff line number Diff line change
Expand Up @@ -56,7 +56,7 @@ def create_quadrature(ref_el, degree, scheme="default", entity=None):
:kwarg scheme: The quadrature scheme, can be choosen from ["default", "canonical", "KMV"]
"default" -> hard-coded scheme for low degree and collapsed Gauss scheme for high degree,
"canonical" -> collapsed Gauss scheme,
"KMV" -> spectral lumped scheme for low degree (<=5 on triangles, <=3 on tetrahedra).
"KMV" -> spectral lumped scheme for low degree (<=6 on triangles, <=3 on tetrahedra).
:kwarg entity: A tuple of entity dimension and entity id specifying the
integration domain. If not provided, the domain is the entire cell.
"""
Expand Down
6 changes: 5 additions & 1 deletion finat/quadrature.py
Original file line number Diff line number Diff line change
Expand Up @@ -18,14 +18,18 @@ def make_quadrature(ref_el, degree, scheme="default"):
Generate quadrature rule for given reference element
that will integrate an polynomial of order 'degree' exactly.

For low-degree (<=6) polynomials on triangles and tetrahedra, this
For low-degree polynomials on triangles (<=50) and tetrahedra (<=15), this
uses hard-coded rules, otherwise it falls back to a collapsed
Gauss scheme on simplices. On tensor-product cells, it is a
tensor-product quadrature rule of the subcells.

:arg ref_el: The FIAT cell to create the quadrature for.
:arg degree: The degree of polynomial that the rule should
integrate exactly.
:kwarg scheme: The quadrature scheme, can be choosen from ["default", "canonical", "KMV"]
"default" -> hard-coded scheme for low degree and collapsed Gauss scheme for high degree,
"canonical" -> collapsed Gauss scheme,
"KMV" -> spectral lumped scheme for low degree (<=6 on triangles, <=3 on tetrahedra).
"""
if ref_el.get_shape() == TENSORPRODUCT:
try:
Expand Down