i_fixed = np.flatnonzero(np.isnan(fit_dist.par_fix))[-1]