diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index 8c63941ad5d916ac780d6d5530971ba54a05d132..f8814cdbcba164afc5dc875bb20d9da2aab249e8 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -544,7 +544,6 @@ continuumDirichletToNeumannMap(const std::string& continuumName,
     dynamic_cast<LinearIterationStep<MatrixType,VectorType>* >(continuum(continuumName).solver_->iterationStep_)->setProblem(*continuum(continuumName).stiffnessMatrix_, x3d, rhs3d);
 
     continuum(continuumName).solver_->preprocess();
-    dynamic_cast<IterationStep<VectorType>* >(continuum(continuumName).solver_->iterationStep_)->preprocess();
         
     continuum(continuumName).solver_->solve();
         
@@ -795,7 +794,6 @@ linearizedContinuumNeumannToDirichletMap(const std::string& continuumName,
     dynamic_cast<LinearIterationStep<MatrixType,VectorType>* >(continuum(continuumName).solver_->iterationStep_)->setProblem(stiffnessMatrix, x, rhs);
 
     //solver.preprocess();
-    continuum(continuumName).solver_->iterationStep_->preprocess();
         
     continuum(continuumName).solver_->solve();
         
@@ -1125,7 +1123,6 @@ iterateWithContact(std::map<std::pair<std::string,std::string>, RigidBodyMotion<
     multigridStep->ignoreNodes_ = &totalDirichletNodes;
     
     contactSolver->preprocess();
-    multigridStep->preprocess();
             
     contactSolver->solve();
             
@@ -1299,7 +1296,6 @@ iterateWithContact(std::map<std::pair<std::string,std::string>, RigidBodyMotion<
         multigridStep->ignoreNodes_ = &totalDirichletNodes;
             
         contactSolver->preprocess();
-        multigridStep->preprocess();
             
         contactSolver->solve();