diff --git a/dune/gfe/localgeodesicfestiffness.hh b/dune/gfe/localgeodesicfestiffness.hh
index 129c4faac5d22137290ac0495dd46aa022b6a552..a83533789739c9b035115054f4b477a5bade07e5 100644
--- a/dune/gfe/localgeodesicfestiffness.hh
+++ b/dune/gfe/localgeodesicfestiffness.hh
@@ -372,7 +372,7 @@ public:
         std::vector<TargetSpace> backwardSolution = localSolution;
         
         
-        for (int j=0; j<blocksize; j++) {
+        for (int j=0; j<embeddedBlocksize; j++) {
             
             // The return value does not have unit norm.  But assigning it to a UnitVector object
             // will normalize it.  This amounts to an extension of the energy functional