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 (x - xp) ** 
    2 + (y - yp) ** 2 + (z - zp) ** 2 <= dist ** 2 + sys.float_info.epsilon)
понимание python с несколькими предложениями «для» и одним «если»