Use phi from the functional, do not reconstruct it

src/samplefunctional.hh

@@ -144,7 +144,7 @@ void minimise(const Functional J, const typename Functional::SmallVector x,
   double const JRestb = (J.b - tmp) * descDir; // <b-Au,v>
 
   typedef typename Functional::NonlinearityType MyNonlinearityType;
-  MyNonlinearityType phi;
+  MyNonlinearityType phi = J.phi;
   typedef DirectionalConvexFunction<MyNonlinearityType>
   MyDirectionalConvexFunctionType;
   MyDirectionalConvexFunctionType JRest(JRestA, JRestb, phi, x, descDir);