From 7617b8f95f7a9d29474ea2fe6d9cc39b8c2a840d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 24 Apr 2012 14:51:18 +0000
Subject: [PATCH] I suspect that I swapped to indices when computing the
gradient of Neumann boundary forces. But I am not sure and I want to
postpone the issue. Therefore this patch adds a cpp warning.
[[Imported from SVN: r8625]]
---
dune/gfe/cosseratenergystiffness.hh | 1 +
1 file changed, 1 insertion(+)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 63e5a1eb..1efb76e0 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -875,6 +875,7 @@ assembleGradient(const Entity& element,
// Only translational dofs are affected by the Neumann force
for (size_t v_i=0; v_i<3; v_i++)
for (size_t j=0; j<3; j++)
+#warning Try whether the arguments of derOfValueWRTCoefficient are swapped
embeddedLocalGradient[i][v_i] += thickness_ * (neumannValue[j] * derOfValueWRTCoefficient[j][v_i]) * quad[pt].weight() * integrationElement;
}
--
GitLab