x = (
                np.sum(ss0 ** 2, 1) * (ss1[:, Y] - vv[:, Y])
                + np.sum(ss1 ** 2, 1) * (vv[:, Y] - ss0[:, Y])
                + np.sum(vv ** 2, 1) * (ss0[:, Y] - ss1[:, Y])