From 5b506c71f775c929600e1abe11dee1b1b94b3c72 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 10 Apr 2011 13:40:17 +0000
Subject: [PATCH] use a higher-order quadrature rule
[[Imported from SVN: r7131]]
---
dune/gfe/harmonicenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index 8a8ea8f7..7a85534b 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -55,7 +55,7 @@ energy(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
- int quadOrder = 1;//gridDim;
+ int quadOrder = gridDim;
const Dune::QuadratureRule<double, gridDim>& quad
= Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder);
@@ -121,7 +121,7 @@ assembleEmbeddedGradient(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
// I am not sure about the correct quadrature order
- int quadOrder = 1;//gridDim;
+ int quadOrder = gridDim;
// numerical quadrature loop
const Dune::QuadratureRule<double, gridDim>& quad
--
GitLab