k1 = np.flatnonzero(1 - np.isfinite(rhs))