How to reuse Jacobian and Inverse Hessian with SciPy minimize

Question

I am currently using the minimize(BFGS) function from SciPy.optimize to calibrate my model. Once calibrated the parameters will be just slightly perturbed from time to time