Implement degree 0 to 2 macro polynomials on a triangle (#687)

* Hard code degree 1 macro polynomials on a triangle * cell types for iso * correct basis functions * decimals not std::sqrt * Add degree 2 macro edge polyset on a triangle
FEniCS · Aug 16, 2023 · fd06036 · fd06036
1 parent 482741c
commit fd06036
Show file tree

Hide file tree

Showing 6 changed files with 538 additions and 7 deletions.
diff --git a/cpp/basix/e-lagrange.cpp b/cpp/basix/e-lagrange.cpp
@@ -1435,10 +1435,11 @@ FiniteElement<T> basix::element::create_iso(cell::type celltype, int degree,
                                             bool discontinuous)
 {
   if (celltype != cell::type::interval && celltype != cell::type::quadrilateral
-      && celltype != cell::type::hexahedron)
+      && celltype != cell::type::hexahedron && celltype != cell::type::triangle)
   {
-    throw std::runtime_error("Can currently only create iso elements on "
-                             "intervals, quadrilaterals, and hexahedra");
+    throw std::runtime_error(
+        "Can currently only create iso elements on "
+        "intervals, triangles, quadrilaterals, and hexahedra");
   }
 
   if (variant == lagrange_variant::unset)