From 403d695edf7d07e6b711c60a48b62c9ce5f4e4a4 Mon Sep 17 00:00:00 2001 From: Oliver Sander <sander@igpm.rwth-aachen.de> Date: Sun, 12 Jun 2011 10:05:55 +0000 Subject: [PATCH] use only first-order quadrature rule [[Imported from SVN: r7416]] --- 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 f24a7676..ab730dbb 100644 --- a/dune/gfe/harmonicenergystiffness.hh +++ b/dune/gfe/harmonicenergystiffness.hh @@ -57,7 +57,7 @@ energy(const Entity& element, LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution); - int quadOrder = gridDim; + int quadOrder = 1;//gridDim; const Dune::QuadratureRule<double, gridDim>& quad = Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder); @@ -107,7 +107,7 @@ assembleEmbeddedGradient(const Entity& element, LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution); // I am not sure about the correct quadrature order - int quadOrder = gridDim; + int quadOrder = 1;//gridDim; // numerical quadrature loop const Dune::QuadratureRule<double, gridDim>& quad -- GitLab