From 0f6de80cb86902e8331c058e686d95b72944e576 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 16 Sep 2011 23:45:11 +0200
Subject: [PATCH] Use phi from the functional, do not reconstruct it

---
 src/samplefunctional.hh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/samplefunctional.hh b/src/samplefunctional.hh
index 4186838c..7a2f63ba 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);
-- 
GitLab