Cleanup

src/samplefunctional.hh

@@ -155,8 +155,9 @@ void minimise(Functional const J, typename Functional::SmallVector &x,
     J.A.mv(descDir, tmp);                // Av
     double const JRestA = tmp * descDir; // <Av,v>
 
-    J.A.mv(x, tmp);                              // Au
-    double const JRestb = (J.b - tmp) * descDir; // <b-Au,v>
+    tmp = J.b;
+    J.A.mmv(x, tmp);                     //  b-Au
+    double const JRestb = tmp * descDir; // <b-Au,v>
 
     typedef typename Functional::NonlinearityType MyNonlinearityType;
     MyNonlinearityType phi = J.phi;