From 41a79985935a0c2053c988a579d38b1439c751e8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 4 Nov 2014 14:12:49 +0000
Subject: [PATCH] Bugfix: the Neumann vector field enters the energy with a
_negative_ sign
[[Imported from SVN: r9949]]
---
dune/gfe/cosseratenergystiffness.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 6d7f7414..755f749b 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -412,7 +412,7 @@ energy(const Entity& element,
// Only translational dofs are affected by the Neumann force
for (size_t i=0; i<neumannValue.size(); i++)
- energy += thickness_ * (neumannValue[i] * value.r[i]) * quad[pt].weight() * integrationElement;
+ energy -= thickness_ * (neumannValue[i] * value.r[i]) * quad[pt].weight() * integrationElement;
}
--
GitLab