Comments

dune/tectonic/samplefunctional.hh

@@ -29,9 +29,9 @@ template <int dim> class SampleFunctional {
 
   double operator()(SmallVector const &v) const {
     SmallVector y;
-    A.mv(v, y);            // y = Av
-    y /= 2;                // y = 1/2 Av
-    y -= b;                // y = 1/2 Av - b
+    A.mv(v, y);            //      Av
+    y /= 2;                //  1/2 Av
+    y -= b;                //  1/2 Av - b
     return y * v + phi(v); // <1/2 Av - b,v> + H(|v|)
   }