diff --git a/1badaptmeasure.cc b/1badaptmeasure.cc index a9d563a28238fda005e775461d5bb331c14feac5..c8fb17ffe979fe3a87af2371988937bf1e415ada 100644 --- a/1badaptmeasure.cc +++ b/1badaptmeasure.cc @@ -145,7 +145,7 @@ void compute (const GridType& grid, multigridStep.setMGType(1, 3, 3); // A V(3,3) cycle multigridStep.setIgnore(dirichletNodes); multigridStep.setBaseSolver(baseSolver); - multigridStep.setSmoother(&presmoother, &postsmoother); + multigridStep.setSmoother(presmoother, postsmoother); multigridStep.setHasObstacle(obsNodes); multigridStep.setObstacles(contactAssembler.getObstacles()); multigridStep.setVerbosity(NumProc::QUIET); diff --git a/1bcontact.cc b/1bcontact.cc index b7abc4c448e88104ca2c203284114b9345975c3a..4a1962eca7c6460bb582cf1c67485f32a1666a29 100644 --- a/1bcontact.cc +++ b/1bcontact.cc @@ -231,15 +231,14 @@ int main (int argc, char *argv[]) try QuadraticIPOptSolver<MatrixType,VectorType> baseSolver(baseTolerance, 100, NumProc::QUIET); // Make pre and postsmoothers - ProjectedBlockGSStep<MatrixType, VectorType> presmoother; - ProjectedBlockGSStep<MatrixType, VectorType> postsmoother; + ProjectedBlockGSStep<MatrixType, VectorType> presmoother, postsmoother; MonotoneMGStep<MatrixType, VectorType> multigridStep(stiffnessMatrix, x, rhs); multigridStep.setMGType(mu, nu1, nu2); multigridStep.setIgnore(dirichletNodes); multigridStep.setBaseSolver(baseSolver); - multigridStep.setSmoother(&presmoother, &postsmoother); + multigridStep.setSmoother(presmoother, postsmoother); multigridStep.setHasObstacles(obsNodes); multigridStep.setObstacles(contactAssembler.getObstacles()); multigridStep.setObstacleRestrictor(ContactObsRestriction<VectorType>());