expose bobyqa hessian

joezuntz · Dec 9, 2024 · 33a782c · 33a782c
1 parent 711efc1
commit 33a782c
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/cosmosis/samplers/maxlike/maxlike_sampler.py b/cosmosis/samplers/maxlike/maxlike_sampler.py
@@ -144,6 +144,10 @@ def run_optimizer(self, inputs):
             }
             optimizer_result = pybobyqa.solve(likefn, start_vector, **kw)
             opt_norm = optimizer_result.x
+            # bobyqa calls it .hessian but scipy calls it .hess, so copy it here
+            # if available
+            if optimizer_result.hessian is not None:
+                optimizer_result.hess = optimizer_result.hessian
         else:
             # Use scipy mainimizer instead
             optimizer_result = scipy.optimize.minimize(likefn, start_vector, method=self.method,