Skip to content
Snippets Groups Projects
Commit 9503c132 authored by Oliver Sander's avatar Oliver Sander
Browse files

Bugfix: 'hesse' is a symmetric matrix, so no all entries can actually be accessed by operator()

parent da4c3bcd
No related branches found
No related tags found
No related merge requests found
......@@ -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);
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment