diff --git a/src/harmonicenergystiffness.hh b/src/harmonicenergystiffness.hh
index c2b2780c454310ea0c186716185f9829f3945b10..c653f8c56354e57e5db26a13efc9f12f5ac33c13 100644
--- a/src/harmonicenergystiffness.hh
+++ b/src/harmonicenergystiffness.hh
@@ -49,8 +49,10 @@ energy(const Entity& element,
        const std::vector<TargetSpace>& localSolution) const
 {
     RT energy = 0;
+
+    assert(element.type().isSimplex());
     
-    LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(element.type(), localSolution);
+    LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
 
     int quadOrder = 1;//gridDim;