From b3216db29bdfa8a8aa968a03e94c471e324bf7ab Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 23 Jan 2011 17:14:34 +0000
Subject: [PATCH] minor cleanup
[[Imported from SVN: r6850]]
---
.../gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 13 +------------
1 file changed, 1 insertion(+), 12 deletions(-)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index 9173dd8a..03afe5ba 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -494,12 +494,6 @@ linearizedRodNeumannToDirichletMap(const RodConfigurationType& currentIterate,
{
std::pair<std::string,std::string> interfaceName = std::make_pair("rod","continuum");
- Dune::FieldVector<double,3> force, torque;
- for (int i=0; i<3; i++) {
- force[i] = forceTorque[i];
- torque[i] = forceTorque[i+3];
- }
-
////////////////////////////////////////////////////
// Assemble the linearized rod problem
////////////////////////////////////////////////////
@@ -527,12 +521,7 @@ linearizedRodNeumannToDirichletMap(const RodConfigurationType& currentIterate,
/////////////////////////////////////////////////////////////////////
// The weak form of the Neumann data
- rhs[0][0] += force[0];
- rhs[0][1] += force[1];
- rhs[0][2] += force[2];
- rhs[0][3] += torque[0];
- rhs[0][4] += torque[1];
- rhs[0][5] += torque[2];
+ rhs[0] += forceTorque;
///////////////////////////////////////////////////////////
// Modify matrix and rhs to contain one Dirichlet node
--
GitLab