except:
            par_cov = np.full(hessian.shape, np.nan)