Skip to content
Snippets Groups Projects
Commit 5ecb12be authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

remove a mystery minus sign which should really be part of the Dirichlet-Neumann implementations 

[[Imported from SVN: r5105]]
parent 2f214482
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 4 deletions
src/rodassembler.cc
+ 2
− 4
View file @ 5ecb12be
...@@ -243,10 +243,8 @@ getResultantForce(const BoundaryPatchBase<PatchGridView>& boundary, ...@@ -243,10 +243,8 @@ getResultantForce(const BoundaryPatchBase<PatchGridView>& boundary,
// assert( std::abs(localStress[2]-canonicalStress*sol[0].q.director(2)) < 1e-6 ); // assert( std::abs(localStress[2]-canonicalStress*sol[0].q.director(2)) < 1e-6 );
// Multiply force times boundary normal to get the transmitted force // Multiply force times boundary normal to get the transmitted force
/** \todo The minus sign comes from the coupling conditions. It canonicalStress *= it->unitOuterNormal(FieldVector<double,0>(0))[0];
should really be in the Dirichlet-Neumann code. */ canonicalTorque *= it->unitOuterNormal(FieldVector<double,0>(0))[0];
canonicalStress *= -it->unitOuterNormal(FieldVector<double,0>(0))[0];
canonicalTorque *= -it->unitOuterNormal(FieldVector<double,0>(0))[0];
} }
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment