This repository was archived by the owner on Feb 21, 2022. It is now read-only.
This repository was archived by the owner on Feb 21, 2022. It is now read-only.
Default degree 3 quadrature on triangle inefficient? #27
Description
Original report by Miklós Homolya (Bitbucket: miklos1, GitHub: miklos1).
FIAT has generic rules for generating arbitrary order collapsed quadrature on triangles and tetrahedra ("canonical"). It also has a few hand-written rules for low-order quadrature on triangles and tetrahedra ("default"), which are presumably more efficient. The default rules indeed have fewer quadrature points than the canonical ones in all cases except one: degree 3 quadrature on the triangle.
The default scheme uses 6-point quadrature in this case, while the canonical scheme only uses 2 x 2 = 4 points.