it_0 = ((x + xp, y + yp, z + zp) for xp in range(-dist, dist + 1) for yp in
    range(-dist, dist + 1) for zp in range(-dist, dist + 1) if xp ** 2 + yp **
    2 + zp ** 2 <= dist ** 2 + sys.float_info.epsilon)
понимание python с несколькими предложениями «для» и одним «если»