i_fixed = np.flatnonzero(np.isfinite(par_fix))