diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 1ac7a33509123e4144e6955c609f68ece42b19c6..59f503adae592a10eb9310e7aa401a9183046a6c 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -473,6 +473,7 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
 
     AverageDistanceAssembler<TargetSpace> assembler(coefficients_, w);
 
+    /** \todo Use a symmetric matrix here */
     Dune::FieldMatrix<RT,embeddedDim,embeddedDim> dFdq(0);
     assembler.assembleEmbeddedHessian(q,dFdq);