From a9889474cff171220afb0fe7a23ec635790a8275 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 19 Feb 2010 11:43:55 +0000
Subject: [PATCH] use the AverageDistanceAssembler to assemble the Hessian of
the average distance functional
[[Imported from SVN: r5596]]
---
src/localgeodesicfefunction.hh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/localgeodesicfefunction.hh b/src/localgeodesicfefunction.hh
index 6b63c877..f465e3fa 100644
--- a/src/localgeodesicfefunction.hh
+++ b/src/localgeodesicfefunction.hh
@@ -130,10 +130,10 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local)
// the actual system matrix
Dune::FieldVector<ctype,dim+1> w = barycentricCoordinates(local);
-
+ AverageDistanceAssembler<TargetSpace> assembler(coefficients_, w);
+
Dune::FieldMatrix<ctype,dim+1,dim+1> dFdq(0);
- for (int i=0; i<dim+1; i++)
- dFdq.axpy(w[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], q));
+ assembler.assembleHessian(q,dFdq);
// ////////////////////////////////////
// solve the system
--
GitLab