Add the 1-point 3D Q1 quadrature (for testing)

pism · Oct 14, 2020 · 46b5831 · 46b5831
1 parent ee471bd
commit 46b5831
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/src/util/fem/Quadrature.cc b/src/util/fem/Quadrature.cc
@@ -162,5 +162,10 @@ Q13DQuadrature8::Q13DQuadrature8() {
   }
 }
 
+Q13DQuadrature1::Q13DQuadrature1() {
+  m_points = {{0.0, 0.0, 0.0}};
+  m_weights = {8.0};
+}
+
 } // end of namespace fem
 } // end of namespace pism
diff --git a/src/util/fem/Quadrature.hh b/src/util/fem/Quadrature.hh
@@ -127,6 +127,11 @@ public:
   Q13DQuadrature8();
 };
 
+class Q13DQuadrature1 : public Quadrature {
+public:
+  Q13DQuadrature1();
+};
+
 } // end of namespace fem
 } // end of namespace pism