i = np.flatnonzero(np.isfinite(pvec))