From 5ecb12beac39e2467c1a51bf6f4420f6b26ef366 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 22 Oct 2009 14:20:58 +0000
Subject: [PATCH] remove a mystery minus sign which should really be part of
the Dirichlet-Neumann implementations
[[Imported from SVN: r5105]]
---
src/rodassembler.cc | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index 6b5f3f9d..4687812c 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -243,10 +243,8 @@ getResultantForce(const BoundaryPatchBase<PatchGridView>& boundary,
// assert( std::abs(localStress[2]-canonicalStress*sol[0].q.director(2)) < 1e-6 );
// Multiply force times boundary normal to get the transmitted force
- /** \todo The minus sign comes from the coupling conditions. It
- should really be in the Dirichlet-Neumann code. */
- canonicalStress *= -it->unitOuterNormal(FieldVector<double,0>(0))[0];
- canonicalTorque *= -it->unitOuterNormal(FieldVector<double,0>(0))[0];
+ canonicalStress *= it->unitOuterNormal(FieldVector<double,0>(0))[0];
+ canonicalTorque *= it->unitOuterNormal(FieldVector<double,0>(0))[0];
}
--
GitLab