diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh index faffd1756c95c5b58ab21cd423016ab41430095c..4b5c6ec22605ce0deb5ed8b5baf23f9d30f33b3e 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();