FIAT

FIAT is also able of generating arbitrary order instances of Jacobi-type quadrature rules on the same element shapes.

Further, H(div) and H(curl) conforming finite element spaces such as the families of Raviart-Thomas, Brezzi-Douglas-Marini and Nedelec are supported on triangles and tetrahedra.

Upcoming versions will also support Hermite and nonconforming elements.

While FIAT may be used in a stand-alone mode by writers of finite element code, we suggest using it underneath a code generation system such as FFC or PETSc expression graphs.