Try to create a 3D Q1 element

pism · Oct 14, 2020 · 02a2859 · 02a2859
1 parent 377504a
commit 02a2859
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/stressbalance/blatter/Poisson3.cc b/src/stressbalance/blatter/Poisson3.cc
@@ -20,7 +20,7 @@
 #include <cassert>
 
 #include "Poisson3.hh"
-
+#include "pism/util/fem/FEM.hh"
 #include "pism/util/error_handling.hh"
 
 namespace pism {
@@ -58,6 +58,8 @@ void Poisson3::compute_residual(DMDALocalInfo *info,
     ym  = info->ym,
     zm  = info->zm;
 
+  fem::Q1Element3 E(*info, m_grid->dx(), m_grid->dy(), fem::Q13DQuadrature8());
+
   // loop over elements
   for (int k = gzs; k < zs + zm; k++) {
     for (int j = gys; j < ys + ym; j++) {