From f85c1df4c035611c3f0f48853cec9c4449441320 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 24 Jan 2011 08:46:46 +0000
Subject: [PATCH] bugfix
[[Imported from SVN: r6857]]
---
dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index acd3a0d6..45f66d5b 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -588,7 +588,6 @@ linearizedRodNeumannToDirichletMap(const std::string& rodName,
if (couplingName.first != rodName)
continue;
- // Use 'forceTorque' as a Neumann value for the rod
const LeafBoundaryPatch<RodGridType>& interfaceBoundary = complex_.coupling(couplingName).rodInterfaceBoundary_;
/** \todo Use the BoundaryPatch iterator, which will be a lot faster
@@ -597,7 +596,7 @@ linearizedRodNeumannToDirichletMap(const std::string& rodName,
*/
for (size_t i=0; i<rhs.size(); i++)
if (interfaceBoundary.containsVertex(i))
- interfaceCorrection[couplingName] = rhs[i];
+ interfaceCorrection[couplingName] = x[i];
}
return interfaceCorrection;
--
GitLab