ix = np.flatnonzero(s > tol)