if somefixed:
            allfixed = np.all(np.isfinite(self.par_fix))
            if not allfixed:
                pcov = -pinv(hessian[self.i_notfixed, :][..., self.i_notfixed])
                for row, i in enumerate(list(self.i_notfixed)):