changed where setLinearizationPoint is happening to avoid refactoring…

… solves during repeated hessian evaluation

soupy/modeling/controlCostFunctional.py

@@ -159,6 +159,8 @@ def computeComponents(self, z, order=0):
             self.model.evalGradientControl(self.x, self.grad_objective)
             self.has_adjoint_solve = True
 
+        self.model.setLinearizationPoint(self.x)
+
         self.has_forward_solve = True
 
 
@@ -182,6 +184,7 @@ def cost(self, z, order=0, **kwargs):
             penalization = self.penalization.cost(z)
         return objective + penalization 
 
+
     def grad(self, g):
         """
         Computes the gradient of the cost functional 
@@ -212,7 +215,7 @@ def hessian(self, zhat, Hzhat):
         
         .. note:: Assumes :code:`self.cost` has been called with :code:`order >= 2`
         """
-        self.model.setLinearizationPoint(self.x)
+        # self.model.setLinearizationPoint(self.x)
         self.model.applyCz(zhat, self.rhs_fwd)
         self.model.solveFwdIncremental(self.uhat, self.rhs_fwd)
         self.model.applyWuu(self.uhat, self.rhs_adj)