From 6aaa804a0f6010f1bb878dfa927a9fd25d618160 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 3 Sep 2013 16:30:24 +0000
Subject: [PATCH] Do not hand the number of levels to the MonotoneMGStep
This number is nowadays infered from the number of transfer operators.
[[Imported from SVN: r9432]]
---
dune/gfe/riemanniantrsolver.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index b86fd19b..d0b3e37c 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -83,7 +83,7 @@ setup(const GridType& grid,
TrustRegionGSStep<MatrixType, CorrectionType>* presmoother = new TrustRegionGSStep<MatrixType, CorrectionType>;
TrustRegionGSStep<MatrixType, CorrectionType>* postsmoother = new TrustRegionGSStep<MatrixType, CorrectionType>;
- MonotoneMGStep<MatrixType, CorrectionType>* mmgStep = new MonotoneMGStep<MatrixType, CorrectionType>(numLevels);
+ MonotoneMGStep<MatrixType, CorrectionType>* mmgStep = new MonotoneMGStep<MatrixType, CorrectionType>;
mmgStep->setMGType(mu, nu1, nu2);
mmgStep->ignoreNodes_ = &dirichletNodes;
@@ -288,7 +288,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
}
- mgStep->setProblem(*hessianMatrix_, corr, rhs, grid_->maxLevel()+1);
+ mgStep->setProblem(*hessianMatrix_, corr, rhs);
trustRegionObstacles.back() = trustRegion.obstacles();
mgStep->obstacles_ = &trustRegionObstacles;
--
GitLab