Check the number of iterations needed by newton

Closes: #11

src/CMakeLists.txt

@@ -31,7 +31,7 @@ dune_add_test(
   COMMAND sh
   CMD_ARGS -exc [[
 ./solution01 -output.filename solution-matrix-free \
-             -solver.preconditioner.type none -operator.type matrix-free
+             -solver.preconditioner.type none -newton.Operator.type matrix-free
 cmp solution-matrix-free.vtu ref.vtu
 ]])
 
@@ -64,7 +64,7 @@ dune_add_test(
   COMMAND sh
   CMD_ARGS -exc [[
 ./nitschesolution01 -output.filename nitschesolution-matrix-free \
-                    -solver.preconditioner.type none -operator.type matrix-free
+                    -solver.preconditioner.type none -newton.Operator.type matrix-free
 cmp nitschesolution-matrix-free.vtu nitscheref.vtu
 ]])
 

src/driver.hh

@@ -100,7 +100,7 @@ namespace PPS {
     newton.setMinLinearReduction(1e-4); // min. red. in lin. solve
     newton.setMaxIterations(25);        // limit number of its
     newton.setLineSearchMaxIterations(10); // limit line search
-    newton.setOperatorParams(ptree.sub("operator")); // select operator type
+    newton.setParameters(ptree.sub("newton"));
 
     // solve nonlinear problem
     newton.apply();

src/newton.hh

@@ -764,7 +764,7 @@ namespace PPS
     */
     void setParameters(const Dune::ParameterTree & param)
     {
-      typedef typename TstV::ElementType RFType;
+      typedef typename TstV::field_type RFType;
       if (param.hasKey("VerbosityLevel"))
         this->setVerbosityLevel(
           param.get<unsigned int>("VerbosityLevel"));

src/octagon.ini

@@ -15,7 +15,12 @@ stab = 100
 [problem]
 eta=100.0
 
-[operator]
+[newton]
+# The 3 is an empirical value, it is expected that the needed number of
+# iterations may vary to some degree
+MaxIterations = 3
+
+[newton.Operator]
 type = matrix-free
 
 [solver]

src/tutorial01.ini

@@ -15,6 +15,11 @@ stab = 100
 [problem]
 eta=100.0
 
+[newton]
+# The 3 is an empirical value, it is expected that the needed number of
+# iterations may vary to some degree
+MaxIterations = 3
+
 [solver]
 #preconditioner.type = none
 maxiter = 100