From 6bbf585ef80487b9c424c8e82fc80c5bd6631c0f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 19 Nov 2010 10:16:30 +0000
Subject: [PATCH] bugfix
[[Imported from SVN: r6546]]
---
dune/gfe/localgeodesicfefunction.hh | 14 +++++---------
1 file changed, 5 insertions(+), 9 deletions(-)
diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index f79f2c3b..bb2098ef 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -313,17 +313,13 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
// best way, though.
Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> dFdqPseudoInv = pseudoInverse(dFdq);
- // Put it all together
- for (size_t i=0; i<result.size(); i++) {
-
- //
- Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
- = TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[i], q);
-
+ //
+ Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
+ = TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[coefficient], q);
+ // Put it all together
+ for (size_t i=0; i<result.size(); i++)
result[i] = dFdqPseudoInv * ( dvDqF[i] * dFdqPseudoInv * dFdw - dpDwF[i]) * B;
-
- }
}
#endif
--
GitLab