From d3205dc95ad73c2afd3872b17408ee72023d6a2d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 27 Oct 2011 09:47:39 +0000
Subject: [PATCH] Fix index bug.  New the unit test in
 localgeodesicfefunctiontest.cc passes

[[Imported from SVN: r8059]]
---
 dune/gfe/localgeodesicfefunction.hh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index a8062aaa..6a19bf65 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;
-- 
GitLab