The rod forces need to be flipped. That used to be in getResultantForce, but...

The rod forces need to be flipped.  That used to be in getResultantForce, but here is a better place

[[Imported from SVN: r5107]]

dirneucoupling.cc

@@ -364,6 +364,10 @@ int main (int argc, char *argv[]) try
         FieldVector<double,dim> resultantForce, resultantTorque;
         resultantForce  = rodAssembler.getResultantForce(couplingBoundary, rodX, resultantTorque);
 
+        // Flip orientation
+        resultantForce  *= -1;
+        resultantTorque *= -1;
+        
         std::cout << "resultant force: " << resultantForce << std::endl;
         std::cout << "resultant torque: " << resultantTorque << std::endl;