diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 3cbdc37024806f8218d9f0f7cdb2aec30845ac68..b17e057bdd057beec152e45ff46674c023db58a8 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -262,7 +262,7 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local, const TargetSpace
       for (size_t j=0; j<RHS.M(); j++) {
         RHS[i][j] = 0;
         for (size_t k=0; k<dFdw.M(); k++)
-          RHS[i][j] += dFdw[i][k]*B[k][0][j];
+          RHS[i][j] += dFdw[i][k][0][0]*B[k][0][j];
       }
 
     // the actual system matrix