From 35f5bcb5d8a824971bcb5d6c1c668c07c7df232f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 14 Oct 2011 16:52:14 +0000
Subject: [PATCH] a more intelligent way to choose the quadrature order
[[Imported from SVN: r7923]]
---
dune/gfe/cosseratenergystiffness.hh | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 067120cb..18d9ea09 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -207,7 +207,8 @@ energy(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, LocalFiniteElement, TargetSpace> localGeodesicFEFunction(localFiniteElement,
localSolution);
- int quadOrder = 1;//gridDim;
+ int quadOrder = (element.type().isSimplex()) ? localFiniteElement.localBasis().order()
+ : localFiniteElement.localBasis().order() * gridDim;
const Dune::QuadratureRule<double, gridDim>& quad
= Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder);
--
GitLab