Add Gauss quadrature for pyramids up to degree 20

src/analysis/CMakeLists.txt

@@ -3,6 +3,7 @@
 add_library(
   PugsAnalysis
   CubeGaussQuadrature.cpp
+  PyramidGaussQuadrature.cpp
   SquareGaussQuadrature.cpp
   TensorialGaussLegendreQuadrature.cpp
   TensorialGaussLobattoQuadrature.cpp

src/analysis/PyramidGaussQuadrature.hpp

+#ifndef PYRAMID_GAUSS_QUADRATURE_HPP
+#define PYRAMID_GAUSS_QUADRATURE_HPP
+
+#include <analysis/QuadratureFormula.hpp>
+
+/**
+ * Defines Gauss quadrature on the reference pyramid element
+ *
+ * \note formulae are provided by
+ *
+ * 'High-order symmetric cubature rules for tetrahedra and pyramids'
+ * Jan Jaśkowiec & N. Sukumar (2020).
+ */
+class PyramidGaussQuadrature final : public QuadratureForumla<3>
+{
+ private:
+  void _buildPointAndWeightLists();
+
+ public:
+  PyramidGaussQuadrature(PyramidGaussQuadrature&&)      = default;
+  PyramidGaussQuadrature(const PyramidGaussQuadrature&) = default;
+
+  explicit PyramidGaussQuadrature(const size_t degree) : QuadratureForumla<3>(degree)
+  {
+    this->_buildPointAndWeightLists();
+  }
+
+  PyramidGaussQuadrature()  = delete;
+  ~PyramidGaussQuadrature() = default;
+};
+
+#endif   // PYRAMID_GAUSS_QUADRATURE_HPP

tests/CMakeLists.txt

@@ -98,6 +98,7 @@ add_executable (unit_tests
   test_PugsAssert.cpp
   test_PugsFunctionAdapter.cpp
   test_PugsUtils.cpp
+  test_PyramidGaussQuadrature.cpp
   test_Q1Transformation.cpp
   test_RevisionInfo.cpp
   test_SmallArray.cpp