Constness

dune/tectonic/samplefunctional.hh

@@ -160,7 +160,7 @@ void minimise(Functional const J, typename Functional::SmallVector &x,
       J.A.mv(descDir, tmp);                //  Av
       double const JRestA = tmp * descDir; // <Av,v>
 
-      MyDirectionalConvexFunction<LocalNonlinearityType> JRest(
+      MyDirectionalConvexFunction<LocalNonlinearityType> const JRest(
           JRestA, JRestb, *J.phi, x, descDir);
       // }}}
 
@@ -197,8 +197,8 @@ void minimise(Functional const J, typename Functional::SmallVector &x,
       Bisection slowBisection(bisection);
       slowBisection.setFastQuadratic(false);
 
-      CircularConvexFunction<LocalNonlinearityType> JRest(J.A, J.b, *J.phi, x,
-                                                          descDir);
+      CircularConvexFunction<LocalNonlinearityType> const JRest(
+          J.A, J.b, *J.phi, x, descDir);
 
       int count;
       double const stepsize = slowBisection.minimize(JRest, 0.0, 1.0, count);