Fewer temporary variables

src/bisection-example-new.cc

@@ -41,13 +41,13 @@ class SampleFunctional {
 
        since A is symmetric.
      */
-    SmallVector tmp2;
-    A_.mv(descDir, tmp2);
-    double const rest_A = tmp2 * descDir;
+    SmallVector tmp;
 
-    SmallVector tmp3;
-    A_.mv(x, tmp3);
-    double const rest_b = (b_ - tmp3) * descDir;
+    A_.mv(descDir, tmp);
+    double const rest_A = tmp * descDir;
+
+    A_.mv(x, tmp);
+    double const rest_b = (b_ - tmp) * descDir;
 
     typedef MyNonlinearity<dimension, Function> MyNonlinearityType;
     MyNonlinearityType phi;