From 658e02b88496f8bae34ddecea713bbe9dda062a4 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 19 Jan 2011 19:32:29 +0000
Subject: [PATCH] bugfix: rod interface boundary wasn't written to the
RodContinuumComplex
[[Imported from SVN: r6803]]
---
dirneucoupling.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 1678e342..9bcd4fa1 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -213,7 +213,7 @@ int main (int argc, char *argv[]) try
BitSetVector<1> rodCouplingBitfield(rodX.size(),false);
// Using that index 0 is always the left boundary for a uniformly refined OneDGrid
rodCouplingBitfield[0] = true;
- LeafBoundaryPatch<RodGridType> rodCouplingBoundary(*complex.rodGrids_["rod"], rodCouplingBitfield);
+ complex.couplings_[interfaceName].rodInterfaceBoundary_.setup(*complex.rodGrids_["rod"], rodCouplingBitfield);
// then for the continuum
LevelBoundaryPatch<GridType> coarseInterfaceBoundary(*complex.continuumGrids_["continuum"], 0);
@@ -385,7 +385,7 @@ int main (int argc, char *argv[]) try
// ///////////////////////////////////////////////////////////
FieldVector<double,dim> resultantForce, resultantTorque;
- resultantForce = rodAssembler.getResultantForce(rodCouplingBoundary, rodX, resultantTorque);
+ resultantForce = rodAssembler.getResultantForce(complex.couplings_[interfaceName].rodInterfaceBoundary_, rodX, resultantTorque);
// Flip orientation
resultantForce *= -1;
--
GitLab