From 022dabcbf2650e69232ca9a7507ce3170a490775 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Fri, 27 Apr 2018 19:40:11 +0200
Subject: [PATCH] Remove Hessian method taking a FieldMatrix
These Hessians are always symmetric, and we should really
exploit this.
Additionally, it seems that the method was never used anyway.
---
dune/gfe/averagedistanceassembler.hh | 16 +---------------
1 file changed, 1 insertion(+), 15 deletions(-)
diff --git a/dune/gfe/averagedistanceassembler.hh b/dune/gfe/averagedistanceassembler.hh
index b4def6ba..3c23d232 100644
--- a/dune/gfe/averagedistanceassembler.hh
+++ b/dune/gfe/averagedistanceassembler.hh
@@ -69,21 +69,6 @@ public:
orthonormalFrame.mv(embeddedGradient,gradient);
}
- void assembleEmbeddedHessian(const TargetSpace& x,
- Dune::FieldMatrix<ctype,embeddedSize,embeddedSize>& matrix) const
- {
- matrix = 0;
- for (size_t i=0; i<coefficients_.size(); i++) {
- Dune::SymmetricMatrix<ctype,embeddedSize> tmp = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], x);
- for (int j=0; j<embeddedSize; j++)
- for (int k=0; k<=j; k++) {
- matrix[j][k] += weights_[i] * tmp(j,k);
- if (j!=k)
- matrix[k][j] += weights_[i] * tmp(j,k);
- }
- }
- }
-
void assembleEmbeddedHessian(const TargetSpace& x,
Dune::SymmetricMatrix<ctype,embeddedSize>& matrix) const
{
@@ -92,6 +77,7 @@ public:
matrix.axpy(weights_[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], x));
}
+ // TODO Use a Symmetric matrix for the result
void assembleHessian(const TargetSpace& x,
Dune::FieldMatrix<ctype,size,size>& matrix) const
{
--
GitLab