From 925e475674e3b43758d6ebc5af3bcc25d2927a3f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 10 Apr 2011 10:04:52 +0000
Subject: [PATCH] bugfix: forgot a factor
[[Imported from SVN: r7125]]
---
dune/gfe/localgeodesicfefunction.hh | 2 ++
1 file changed, 2 insertions(+)
diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 2045b0de..83187d0d 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -317,6 +317,8 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
= TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[coefficient], q);
+ dvDqF = w[coefficient] * dvDqF;
+
// Put it all together
for (size_t i=0; i<result.size(); i++)
result[i] = dFdqPseudoInv * ( dvDqF[i] * dFdqPseudoInv * dFdw - dpDwF[i]) * B;
--
GitLab