From b25cc485b5a674413fddb11605460c3a13bef00c Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 30 Apr 2009 15:29:06 +0000
Subject: [PATCH] use multigrid
[[Imported from SVN: r4134]]
---
rod-eoc.cc | 11 ++++++-----
1 file changed, 6 insertions(+), 5 deletions(-)
diff --git a/rod-eoc.cc b/rod-eoc.cc
index d99a6dcc..4c5e1f36 100644
--- a/rod-eoc.cc
+++ b/rod-eoc.cc
@@ -30,7 +30,8 @@ using namespace Dune;
using std::string;
void solve (const GridType& grid,
- SolutionType& x, int numLevels,
+ SolutionType& x,
+ int numLevels,
const TargetSpace& dirichletValue,
ConfigParser& parameters)
{
@@ -84,7 +85,7 @@ void solve (const GridType& grid,
RodAssembler<GridType> rodAssembler(grid, &localStiffness);
RiemannianTrustRegionSolver<GridType,RigidBodyMotion<3> > rodSolver;
-#if 0
+#if 1
rodSolver.setup(grid,
&rodAssembler,
x,
@@ -94,9 +95,9 @@ void solve (const GridType& grid,
initialTrustRegionRadius,
multigridIterations,
innerTolerance,
- mu, nu1, nu2,
- baseIterations,
- baseTolerance,
+ 1, 3, 3,
+ 100, // iterations of the base solver
+ 1e-8, // base tolerance
false); // instrumentation
#else
rodSolver.setupTCG(grid,
--
GitLab