From ed20a1b6f4c62dba37eeb72248a9729493256364 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 23 Jan 2011 16:56:49 +0000
Subject: [PATCH] hack removal: computeAveragePressure takes a single combined
force/torque vector now
[[Imported from SVN: r6847]]
---
dirneucoupling.cc | 17 ++++-------------
1 file changed, 4 insertions(+), 13 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index da870fcc..7684803d 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -388,24 +388,15 @@ int main (int argc, char *argv[]) try
RigidBodyMotion<3>::TangentVector resultantForceTorque
= rodAssembler.getResultantForce(complex.couplings_[interfaceName].rodInterfaceBoundary_, rodX);
- // separate into translational and rotational part
- FieldVector<double,dim> resultantForce, resultantTorque;
- for (int j=0; j<dim; j++) {
- resultantForce[j] = resultantForceTorque[j];
- resultantTorque[j] = resultantForceTorque[dim+j];
- }
-
// Flip orientation
- resultantForce *= -1;
- resultantTorque *= -1;
-
- std::cout << "resultant force: " << resultantForce << std::endl;
- std::cout << "resultant torque: " << resultantTorque << std::endl;
+ resultantForceTorque *= -1;
+
+ std::cout << "resultant force and torque: " << resultantForceTorque << std::endl;
VectorType neumannValues(rhs3d.size());
// Using that index 0 is always the left boundary for a uniformly refined OneDGrid
- computeAveragePressure<GridType::LeafGridView>(resultantForce, resultantTorque,
+ computeAveragePressure<GridType::LeafGridView>(resultantForceTorque,
complex.couplings_[interfaceName].continuumInterfaceBoundary_,
rodX[0].r,
neumannValues);
--
GitLab