i = numpy.where(abs(numpy.real(w) - 1.0) < 1e-8)[0]