[Getfem-commits] (no subject)

[Konstantinos Poulios]

branch: master
commit 3480c3bde8804e87ae15bf952cf9a7068ea35387
Author: Konstantinos Poulios <address@hidden>
Date:   Fri Sep 15 12:46:58 2017 +0200

    Account for pyramid elements in classic integration methods
---
 src/getfem_integration.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/getfem_integration.cc b/src/getfem_integration.cc
index 437f816..7131cc4 100644
--- a/src/getfem_integration.cc
+++ b/src/getfem_integration.cc
@@ -1263,6 +1263,9 @@ namespace getfem {
       }
       GMM_ASSERT1(false, "could not find an " << name.str()
                   << " of degree >= " << int(degree));
+    } else if (bgeot::basic_structure(cvs) == bgeot::pyramid_QK_structure(1)) {
+      GMM_ASSERT1(n == 3, "Wrong dimension");
+      name << "IM_PYRAMID(IM_GAUSS_PARALLELEPIPED(3," << degree << "))";
     } else if (cvs->is_product(&a,&b) ||
                (bgeot::basic_structure(cvs).get() &&
                 bgeot::basic_structure(cvs)->is_product(&a,&b))) {