diff --git a/src/bisection-example-new.cc b/src/bisection-example-new.cc
index c789cff35ceca310939eab18ce95889a1a1bcf33..59035dfae3aedf5bd9ae44c2bff813b3c2dd1c4b 100644
--- a/src/bisection-example-new.cc
+++ b/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;