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

rename a few variables for easier legibility

[[Imported from SVN: r6759]]
parent feee12fc
No related branches found
No related tags found
No related merge requests found
......@@ -389,15 +389,15 @@ void RodContinuumSteklovPoincareStep<RodGridType,ContinuumGridType>::iterate(Rig
couplingBitfield[0] = true;
LeafBoundaryPatch<RodGridType> couplingBoundary(*complex_.rodGrids_["rod"], couplingBitfield);
Dune::FieldVector<double,dim> resultantForce, resultantTorque;
resultantForce = rodAssembler_->getResultantForce(couplingBoundary, rodX, resultantTorque);
Dune::FieldVector<double,dim> rodForce, rodTorque;
rodForce = rodAssembler_->getResultantForce(couplingBoundary, rodX, rodTorque);
// Flip orientation
resultantForce *= -1;
resultantTorque *= -1;
rodForce *= -1;
rodTorque *= -1;
std::cout << "resultant force: " << resultantForce << std::endl;
std::cout << "resultant torque: " << resultantTorque << std::endl;
std::cout << "resultant force: " << rodForce << ", "
<< "resultant torque: " << rodTorque << std::endl;
///////////////////////////////////////////////////////////////////
// Evaluate the Dirichlet-to-Neumann map for the continuum
......@@ -445,8 +445,8 @@ void RodContinuumSteklovPoincareStep<RodGridType,ContinuumGridType>::iterate(Rig
// Compute the overall Steklov-Poincare residual
///////////////////////////////////////////////////////////////
Dune::FieldVector<double,3> residualForce = resultantForce + continuumForce;
Dune::FieldVector<double,3> residualTorque = resultantTorque + continuumTorque;
Dune::FieldVector<double,3> residualForce = rodForce + continuumForce;
Dune::FieldVector<double,3> residualTorque = rodTorque + continuumTorque;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment