From 9503c13242c78c65119ce934dacc818c5a0a6d2f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 4 Oct 2015 14:11:13 +0200
Subject: [PATCH] Bugfix: 'hesse' is a symmetric matrix, so no all entries can
actually be accessed by operator()
---
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 8cff20de..10f79490 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -490,7 +490,7 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
for (size_t k=0; k<coefficients_.size(); k++) {
Dune::SymmetricMatrix<RT,embeddedDim> hesse = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[k], q);
for (int i=0; i<embeddedDim; i++)
- for (int j=0; j<embeddedDim; j++)
+ for (int j=0; j<i; j++)
dqdwF(i, j, k) = dqdwF(j, i, k) = hesse(i,j);
}
--
GitLab