construct higher degree elements in example run, not test

examples/test_high_degree.rs

@@ -0,0 +1,56 @@
+//! Finite element definitions
+
+use ndelement::ciarlet::{lagrange, nedelec, raviart_thomas};
+use ndelement::types::{Continuity, ReferenceCellType};
+use paste::paste;
+
+fn main() {
+    macro_rules! construct_lagrange {
+        ($cell:ident, $max_degree:expr) => {
+            paste! {
+                for d in 1..[<$max_degree>] {
+                    println!("Constructing Lagrange(degree={d}, cell={:?})", ReferenceCellType::[<$cell>]);
+                    let _e = lagrange::create::<f64>(ReferenceCellType::[<$cell>], d, Continuity::Standard);
+                }
+            }
+        };
+    }
+
+    macro_rules! construct_raviart_thomas {
+        ($cell:ident, $max_degree:expr) => {
+            paste! {
+                for d in 1..[<$max_degree>] {
+                    println!("Constructing RaviartThomas(degree={d}, cell={:?})", ReferenceCellType::[<$cell>]);
+                    let _e = raviart_thomas::create::<f64>(ReferenceCellType::[<$cell>], d, Continuity::Standard);
+                }
+            }
+        };
+    }
+
+    macro_rules! construct_nedelec {
+        ($cell:ident, $max_degree:expr) => {
+            paste! {
+                for d in 1..[<$max_degree>] {
+                    println!("Constructing Nedelec(degree={d}, cell={:?})", ReferenceCellType::[<$cell>]);
+                    let _e = nedelec::create::<f64>(ReferenceCellType::[<$cell>], d, Continuity::Standard);
+                }
+            }
+        };
+    }
+
+    construct_lagrange!(Interval, 14);
+    construct_lagrange!(Triangle, 8);
+    construct_lagrange!(Quadrilateral, 8);
+    construct_lagrange!(Tetrahedron, 6);
+    construct_lagrange!(Hexahedron, 6);
+
+    construct_raviart_thomas!(Triangle, 8);
+    construct_raviart_thomas!(Quadrilateral, 8);
+    construct_raviart_thomas!(Tetrahedron, 6);
+    construct_raviart_thomas!(Hexahedron, 6);
+
+    construct_nedelec!(Triangle, 8);
+    construct_nedelec!(Quadrilateral, 8);
+    construct_nedelec!(Tetrahedron, 6);
+    construct_nedelec!(Hexahedron, 6);
+}

python/ndelement/ciarlet.py

@@ -228,38 +228,6 @@ def continuity(self) -> Continuity:
 
     def element(self, cell: ReferenceCellType) -> CiarletElement:
         """Create an element."""
-        # TODO: remove these error once https://github.com/linalg-rs/rlst/issues/98 is fixed
-        msg = "Cannot create element due to bug in RLST"
-        if self.family == Family.Lagrange:
-            if cell == ReferenceCellType.Interval and self.degree >= 99:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Triangle and self.degree >= 13:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Quadrilateral and self.degree >= 10:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Tetrahedron and self.degree >= 7:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Hexahedron and self.degree >= 5:
-                raise RuntimeError(msg)
-        if self.family == Family.RaviartThomas:
-            if cell == ReferenceCellType.Triangle and self.degree >= 10:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Quadrilateral and self.degree >= 7:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Tetrahedron and self.degree >= 5:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Hexahedron and self.degree >= 3:
-                raise RuntimeError(msg)
-        if self.family == Family.NedelecFirstKind:
-            if cell == ReferenceCellType.Triangle and self.degree >= 10:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Quadrilateral and self.degree >= 7:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Tetrahedron and self.degree >= 5:
-                raise RuntimeError(msg)
-            if cell == ReferenceCellType.Hexahedron and self.degree >= 3:
-                raise RuntimeError(msg)
-
         return CiarletElement(_lib.element_family_create_element(self._rs_family, cell.value))
 
 

src/ciarlet.rs

@@ -1275,5 +1275,4 @@ mod test {
     test_entity_closure_dofs_lagrange!(Tetrahedron, 5);
     test_entity_closure_dofs_lagrange!(Hexahedron, 2);
     test_entity_closure_dofs_lagrange!(Hexahedron, 3);
-    test_entity_closure_dofs_lagrange!(Hexahedron, 4);
 }