Hi, maybe some unpopular opinion here. Did you try to use Numba in the minimizer functions F, dF and H? I get speedups of a factor \times 15 in all of them. This is just an easy to implement/intermediate solution, because it seems that only Julia will be much faster, but in the meantime it may help a bit and you can still maintain your Python code.