Create bisection object just once

src/myblockproblem.hh

@@ -32,8 +32,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
   class IterateObject;
 
   MyBlockProblem(MyConvexProblemType& problem) : problem(problem) {
-    // TODO: Is it clever to create a bisection here?
-    bisection = Bisection(0.0, 1.0, 1e-15, true, 1e-14);
+    bisection = Bisection(0.0, 1.0, 1e-12, true, 0);
   };
 
   /** \brief Constructs and returns an iterate object */
@@ -119,7 +118,7 @@ class MyBlockProblem<MyConvexProblemTypeTEMPLATE>::IterateObject {
       Dune::SampleFunctional<block_size> localJ(*localA, localb, phi);
 
       LocalVectorType correction;
-      Dune::minimise(localJ, ui, 10); // FIXME: hardcoded value
+      Dune::minimise(localJ, ui, 10, bisection); // FIXME: hardcoded value
     }
   }