From f5c69b34eeb3ca7513182c04699ba79f1995a92a Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Thu, 21 Dec 2017 11:07:28 +0100
Subject: [PATCH] Lower the quad order for cube elements
The new value is what you need for linear target spaces.
I am not convinced that this is sufficient, but it matches
what we do for simplices.
The whole quadrature business will need some theoretical
investigations eventually.
---
dune/gfe/harmonicenergystiffness.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index b39d882f..645b519a 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -38,7 +38,7 @@ energy(const typename Basis::LocalView& localView,
LocalInterpolationRule localInterpolationRule(localFiniteElement,localSolution);
int quadOrder = (localFiniteElement.type().isSimplex()) ? (localFiniteElement.localBasis().order()-1) * 2
- : localFiniteElement.localBasis().order() * 2 * gridDim;
+ : (localFiniteElement.localBasis().order() * gridDim - 1) * 2;
const auto element = localView.element();
--
GitLab