We always use TNNMG

src/one-body-sample.cc

@@ -230,7 +230,7 @@ int main(int argc, char *argv[]) {
         std::cout.flush();
       }
 
-      if (parset.get<bool>("solver.tnnmg.use")) {
+      {
         assemble_neumann<GridType, GridView, SmallVector, P1Basis>(
             leafView, p1Basis, neumannNodes, b4,
             functions.get("neumannCondition"), h * run);

src/one-body-sample.parset

@@ -24,7 +24,6 @@ nu = 0.3 # Needs to be strictly below 0.5; the closer we get, the
 tolerance = 1e-14
 
 [solver.tnnmg]
-use = true
 maxiterations = 1000000
 fixed_point_iterations = 20
 fixed_point_tolerance = 1e-14