From a62b39265d00ffb613e90c0038d67fd26e81d3d9 Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Wed, 8 Jul 2020 00:34:48 +0200
Subject: [PATCH] Reduce quadrature order, before 16 quadrature points were
used, an amount that's sufficient to integrate over a 3D-element, but as the
integral is on the boundary, now 9 quadrature points, which are enough for
2D, are used.
---
dune/gfe/surfacecosseratenergy.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/surfacecosseratenergy.hh b/dune/gfe/surfacecosseratenergy.hh
index 101c37b9..099f97d7 100644
--- a/dune/gfe/surfacecosseratenergy.hh
+++ b/dune/gfe/surfacecosseratenergy.hh
@@ -282,7 +282,7 @@ RT energy(const typename Basis::LocalView& localView,
auto id = idSet.subId(it.inside(), it.indexInInside(), 1);
auto boundaryGeometry = geometriesOnShellBoundary_.at(id);
auto quadOrder = (it.type().isSimplex()) ? localFiniteElement.localBasis().order()
- : localFiniteElement.localBasis().order() * gridDim;
+ : localFiniteElement.localBasis().order() * boundaryDim;
const auto& quad = Dune::QuadratureRules<DT, boundaryDim>::rule(it.type(), quadOrder);
for (size_t pt=0; pt<quad.size(); pt++) {
--
GitLab