diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index fc8d0abd4bde31e95b62c0efa6257507db2263d3..fd4d1f0bef13bc39f8d787d5cb6cd67b5d5629bc 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -559,12 +559,12 @@ evaluateFDDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>
         for (int l=0; l<dim; l++) {
                     
             for (int k=0; k<embeddedDim; k++)
-                foo[k] = result[k][j][l];
+                foo[k] = result[j][k][l];
 
             foo = q.projectOntoTangentSpace(foo);
 
             for (int k=0; k<embeddedDim; k++)
-                result[k][j][l] = foo[k];
+                result[j][k][l] = foo[k];
                     
         }