diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index a8062aaa2ecac781355fdcd3f0e6e50715f056cd..6a19bf65c141a3f22f17c137164692bb486fa5f1 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -479,7 +479,7 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
             for (size_t k=0; k<coefficients_.size(); k++) {
                 dqdwF_times_dvq(i, j, k) = 0;
                 for (int l=0; l<embeddedDim; l++)
-                    dqdwF_times_dvq(i, j, k) += dqdwF(l, j, k) * dvq[l][i];
+                    dqdwF_times_dvq(i, j, k) += dqdwF(l, j, k) * dvq[i][l];
             }
 
     Tensor3<double, embeddedDim,embeddedDim,dim> foo;