From e42c84089ecda47277563e517bd77c196d1d9f04 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 15 Mar 2012 07:33:55 +0000
Subject: [PATCH] enable all three parts of the gradient
[[Imported from SVN: r8570]]
---
dune/gfe/cosseratenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 764d6c19..029689c5 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -816,7 +816,7 @@ assembleGradient(const Entity& element,
if (gridDim==2) {
typename TargetSpace::EmbeddedTangentVector tmp(0);
longQuadraticMembraneEnergyGradient(tmp,R,dR_dv,derivative,derOfGradientWRTCoefficient,U);
- //embeddedLocalGradient[i].axpy(weight * thickness_, tmp);
+ embeddedLocalGradient[i].axpy(weight * thickness_, tmp);
tmp = 0;
curvatureEnergyGradient(tmp,R,DR,dDR_dv);
@@ -824,7 +824,7 @@ assembleGradient(const Entity& element,
tmp = 0;
bendingEnergyGradient(tmp,R,dR_dv,DR,dDR_dv);
- //embeddedLocalGradient[i].axpy(weight * std::pow(thickness_,3) / 12.0, tmp);
+ embeddedLocalGradient[i].axpy(weight * std::pow(thickness_,3) / 12.0, tmp);
} else if (gridDim==3) {
assert(gridDim==2); // 3d not implemented yet
// energy += weight * quadraticMembraneEnergyGradient(U);
--
GitLab