diff --git a/src/samplefunctional.hh b/src/samplefunctional.hh
index 4186838c799ba08bbe06af1f46f8d542af073dc9..7a2f63bac59770c3e948ff7ec8a0397358403ec3 100644
--- a/src/samplefunctional.hh
+++ b/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);