Skip to content
Snippets Groups Projects
Commit 5b506c71 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

use a higher-order quadrature rule 

[[Imported from SVN: r7131]]
parent e57583ac
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/gfe/harmonicenergystiffness.hh
+ 2
− 2
View file @ 5b506c71
...@@ -55,7 +55,7 @@ energy(const Entity& element, ...@@ -55,7 +55,7 @@ energy(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution); LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
int quadOrder = 1;//gridDim; int quadOrder = gridDim;
const Dune::QuadratureRule<double, gridDim>& quad const Dune::QuadratureRule<double, gridDim>& quad
= Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder); = Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder);
...@@ -121,7 +121,7 @@ assembleEmbeddedGradient(const Entity& element, ...@@ -121,7 +121,7 @@ assembleEmbeddedGradient(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution); LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
// I am not sure about the correct quadrature order // I am not sure about the correct quadrature order
int quadOrder = 1;//gridDim; int quadOrder = gridDim;
// numerical quadrature loop // numerical quadrature loop
const Dune::QuadratureRule<double, gridDim>& quad const Dune::QuadratureRule<double, gridDim>& quad
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment