The minimal polynomial is
x^6 - 4750778730825177725463920948909726618214491718039471628856160047142000*x^5
+ 1247257156019654977752984724237035223986874495851369721538788760906203681342078134750000*x^4
- 3711837295929728841959711983585876317003061224025827146858434575359795370054804675000000000*x^3
+ 3160150517834696901784153329054103916683908447320628422609243420587305378116705388187500000000*x^2
- 700265800610377949731030701279743787907896705372015247301698864378963568331565476375000000000000*x
+ 327451677250026694198133278336402500637269636726169887842595798606445909367058576435277640625000000000000
It has degree $6$, not $3$, because $\sqrt{-652}$ generates
the quadratic order of discriminant $-4^2 163$ which has index $4$
in the full ring of integers, and this order has class number $6$.
The index-$2$ order is generated by $\sqrt{-163}$, and $j(\sqrt{-163})$
has cubic minimal polynomial
x^3 - 68925893036109279891085639286946000*x^2
+ 102561728837719322645921325412908000000*x
- 18095625621665522953693950872675200892692248000000000
gp's default precision is not nearly enough to find these polynomials.
200 digits suffice for the cubic, and 1200 for the sextic, e.g.
\p 1200
algdep(ellj(sqrt(-652)),6)
The cubic can also be obtained by plugging
$j\bigl(\frac12(1+\sqrt{-163})\bigr) = -640320^3$ into the modular polynomial
$\Phi_2(j,x)$; the sextic can be computed similarly with a bit more work.