From e6804a84db5a3c58e4a6af4c861e53c5e945b89b Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 14 Jan 2011 10:32:27 +0000
Subject: [PATCH] use the RodFactory to compute starting values for the rod
Dirichlet problems
[[Imported from SVN: r6746]]
---
dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index faffd175..4b5c6ec2 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -368,10 +368,15 @@ void RodContinuumSteklovPoincareStep<RodGridType,ContinuumGridType>::iterate(Rig
// solve a Dirichlet problem for the rod
RodSolutionType rodX;
- /** \todo Using that index 0 is always the left boundary for a uniformly refined OneDGrid */
- rodX[0] = lambda;
+ /** \todo Using that the coupling boundary is the one with the lower coordinate */
#warning Dirichlet boundary not properly set
- rodX.back() = lambda;
+ RigidBodyMotion<3> rodDirichletValue;
+ rodDirichletValue.r = 0;
+ rodDirichletValue.q = Rotation<3,double>::identity();
+
+ RodFactory<typename RodGridType::LeafGridView> rodFactory(complex_.rodGrids_["rod"]->leafView());
+ rodFactory.create(rodX,lambda,rodDirichletValue);
+
rodSolver_->setInitialSolution(rodX);
rodSolver_->solve();
--
GitLab