Currently using scipy.optimize.linprog is too slow to use with expanded confidence intervals. See if the Hessian has special structure.