Skip to content
Snippets Groups Projects

WIP: Add test computation for bending isometries using reducedcubichermitetrianglebasis

Open Klaus Böhnlein requested to merge feature/bendingIsometries into master
4 files
+ 60
55
Compare changes
  • Side-by-side
  • Inline
Files
4
@@ -46,7 +46,7 @@ namespace Dune::GFE
: (localFiniteElement.localBasis().order() * gridDim - 1) * 2;
const auto element = localView.element();
localFunction_.bind(element, localSolution);
localFunction_.bind(element);
const auto& quad = QuadratureRules<double, gridDim>::rule(localFiniteElement.type(), quadOrder);
@@ -59,7 +59,7 @@ namespace Dune::GFE
// Compute the Frobenius norm squared of the derivative -- i.e., the harmonic energy
// TODO: Replace this by the actual bending energy
energy += weight * jacobian.two_norm2();
energy += weight * jacobian.frobenius_norm2();
}
return 0.5 * energy;
Loading