Two polynomials f, g ∈ double-struck F[x1...,xn] are called shift-equivalent if there exists a vector (a1..., a n ∈ double-struck Fn such that the polynomial identity f(x1+a1,...,xn+an) ≡g(x1,...,xn) holds. Our main result is a new randomized algorithm that tests whether two given polynomials are shift equivalent. Our algorithm runs in time polynomial in the circuit size of the polynomials, to which it is given black box access. This complements a previous work of Grigoriev [Gri97 who gave a deterministic algorithm running in time nO(d) for degree d polynomials. Our algorithm uses randomness only to solve instances of the Polynomial Identity Testing (PIT) problem. Hence, if one could de-randomize PIT (a long-standing open problem in complexity) a de-randomization of our algorithm would follow. This establishes an equivalence between de-randomizing shift-equivalence testing and de-randomizing PIT (both in the black-box and the white-box setting). For certain restricted models, such as Read Once Branching Programs, we already obtain a deterministic algorithm using existing PIT results.