Comments

src/bisection-example-new.cc

@@ -43,11 +43,11 @@ class SampleFunctional {
      */
     SmallVector tmp;
 
-    A_.mv(descDir, tmp);
-    double const rest_A = tmp * descDir;
+    A_.mv(descDir, tmp);                 // Av
+    double const rest_A = tmp * descDir; // <Av,v>
 
-    A_.mv(x, tmp);
-    double const rest_b = (b_ - tmp) * descDir;
+    A_.mv(x, tmp);                              // Au
+    double const rest_b = (b_ - tmp) * descDir; // <b-Au,v>
 
     typedef MyNonlinearity<dimension, Function> MyNonlinearityType;
     MyNonlinearityType phi;
@@ -55,7 +55,8 @@ class SampleFunctional {
     MyDirectionalConvexFunctionType;
     MyDirectionalConvexFunctionType rest(rest_A, rest_b, phi, x, descDir);
 
-    Bisection bisection; // FIXME: default values
+    // FIXME: default values are used
+    Bisection bisection;
     int count;
     // FIXME: does x_old = 1 make any sense?!
     double const m = bisection.minimize(rest, 0.0, 1.0, count);