From 5dddbaa2e71d27fa21c5571713cb8cc96174d65f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 30 Sep 2013 15:56:39 +0000
Subject: [PATCH] [bugfix] Properly implement the curvature correction
The Weingarten map is block-diagonal for a product space.
[[Imported from SVN: r9517]]
---
dune/gfe/localgeodesicfeadolcstiffness.hh | 32 +++++++++++++++--------
1 file changed, 21 insertions(+), 11 deletions(-)
diff --git a/dune/gfe/localgeodesicfeadolcstiffness.hh b/dune/gfe/localgeodesicfeadolcstiffness.hh
index 1dc7c18b..95470608 100644
--- a/dune/gfe/localgeodesicfeadolcstiffness.hh
+++ b/dune/gfe/localgeodesicfeadolcstiffness.hh
@@ -253,22 +253,32 @@ assembleHessian(const Entity& element,
this->A_[row][col][subRow][subCol] = tmp2[subRow];
}
- // + A_x (z, P_x^\orth \partial f)
- for (size_t row=0; row<nDofs; row++) {
+ }
- // Get the Euclidean gradient projected onto the normal space
- EmbeddedTangentVector porthGrad = localSolution[row].projectOntoNormalSpace(localEmbeddedGradient[row]);
+ }
- EmbeddedTangentVector tmp3 = localSolution[row].weingarten(z, porthGrad);
+ //////////////////////////////////////////////////////////////////////////////////////
+ // Further correction due to non-planar configuration space
+ // + \mathfrak{A}_x(z,P^\orth_x \partial f)
+ //////////////////////////////////////////////////////////////////////////////////////
- TangentVector tmp4;
- orthonormalFrame[row].mv(tmp3,tmp4);
+ // Project embedded gradient onto normal space
+ std::vector<typename TargetSpace::EmbeddedTangentVector> projectedGradient(localSolution.size());
+ for (size_t i=0; i<localSolution.size(); i++)
+ projectedGradient[i] = localSolution[i].projectOntoNormalSpace(localEmbeddedGradient[i]);
- for (int subRow=0; subRow<blocksize; subRow++)
- this->A_[row][col][subRow][subCol] += tmp4[subRow];
- }
+ for (size_t row=0; row<nDofs; row++) {
- }
+ for (size_t subRow=0; subRow<blocksize; subRow++) {
+
+ EmbeddedTangentVector z = orthonormalFrame[row][subRow];
+ EmbeddedTangentVector tmp1 = localSolution[row].weingarten(z,projectedGradient[row]);
+
+ TangentVector tmp2;
+ orthonormalFrame[row].mv(tmp1,tmp2);
+
+ this->A_[row][row][subRow] += tmp2;
+ }
}
--
GitLab