Add regular sphere with quadrilateral cells (#28)

* use macros in tests

* add regular sphere meshed using quads

* correct faces of cube

* remove print

Author: mscroggs

ndelement/examples/lagrange_element.rs

@@ -22,7 +22,7 @@ fn main() {
     element.tabulate(&points, 0, &mut basis_values);
     println!(
         "The values of the basis functions at the point (1/3, 1/3) are: {:?}",
-        basis_values.data()
+        basis_values.data().unwrap()
     );
 
     // Set point to [1, 0]
@@ -32,6 +32,6 @@ fn main() {
     element.tabulate(&points, 0, &mut basis_values);
     println!(
         "The values of the basis functions at the point (1, 0) are: {:?}",
-        basis_values.data()
+        basis_values.data().unwrap()
     );
 }

ndfunctionspace/examples/test_mass_matrix.rs

@@ -16,7 +16,7 @@ use rlst::{DynArray, rlst_dynamic_array};
 
 /// Test values in Lagrange mass matrix
 fn test_lagrange_mass_matrix() {
-    let grid = regular_sphere(0);
+    let grid = regular_sphere(0, ReferenceCellType::Triangle);
 
     let family = LagrangeElementFamily::<f64>::new(1, Continuity::Standard);
     let space = FunctionSpaceImpl::new(&grid, &family);
@@ -84,7 +84,7 @@ fn test_lagrange_mass_matrix() {
 
 /// Test values in Raviart-Thomas mass matrix
 fn test_rt_mass_matrix() {
-    let grid = regular_sphere(0);
+    let grid = regular_sphere(0, ReferenceCellType::Triangle);
 
     let family = RaviartThomasElementFamily::<f64>::new(1, Continuity::Standard);
     let space = FunctionSpaceImpl::new(&grid, &family);

ndgrid/src/io/gmsh.rs

@@ -776,8 +776,11 @@ mod test {
 
     #[test]
     fn test_regular_sphere_gmsh_io() {
-        let g = regular_sphere::<f64>(2);
+        let g = regular_sphere::<f64>(2, ReferenceCellType::Triangle);
         g.export_as_gmsh("_test_io_sphere.msh");
+
+        let g = regular_sphere::<f64>(2, ReferenceCellType::Quadrilateral);
+        g.export_as_gmsh("_test_io_sphere_quads.msh");
     }
 
     #[test]

ndgrid/src/io/ron.rs

@@ -40,7 +40,7 @@ mod test {
 
     #[test]
     fn test_ron_export_and_import() {
-        let g = regular_sphere::<f64>(1);
+        let g = regular_sphere::<f64>(1, ReferenceCellType::Triangle);
         let n = g.entity_count(ReferenceCellType::Interval);
         g.export_as_ron("_test_export.ron");