Skip to content
Snippets Groups Projects
Commit fbc97657 authored by Sander, Oliver's avatar Sander, Oliver
Browse files

Merge branch 'fix/use-correct-quadrature-order' into 'master'

Reduce quadrature order, before 16 quadrature points were used, an amount...

See merge request !46
parents 2008bdd7 a62b3926
No related branches found
No related tags found
1 merge request!46Reduce quadrature order, before 16 quadrature points were used, an amount...
Pipeline #4384 failed
......@@ -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++) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment