From a3938e126f021ab5665d3d9b68275494bc17725d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 19 Feb 2010 11:45:06 +0000
Subject: [PATCH] implement the actual Hessian
[[Imported from SVN: r5597]]
---
src/averagedistanceassembler.hh | 8 ++++++++
1 file changed, 8 insertions(+)
diff --git a/src/averagedistanceassembler.hh b/src/averagedistanceassembler.hh
index 0dfd06d2..18de1a9e 100644
--- a/src/averagedistanceassembler.hh
+++ b/src/averagedistanceassembler.hh
@@ -53,6 +53,14 @@ public:
matrix[i][j] = (i==j);
}
+ void assembleHessian(const TargetSpace& x,
+ Dune::FieldMatrix<double,size,size>& matrix) const
+ {
+ matrix = 0;
+ for (int i=0; i<dim+1; i++)
+ matrix.axpy(w[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], q));
+ }
+
const std::vector<TargetSpace> coefficients_;
const std::vector<double> weights_;
--
GitLab