The trick is that for the expression $ x^3 - 3Dx $, the substitution $ x = t + \frac{D}{t} $ changes it to $ t^3 + \frac{D^3}{t^3} $ hence completely eliminating the 'linear term'. With that substitution, the equation becomes a quadratic $ t^6 - 2Dt^3 + D^3 = 0 $.
Solve the quadratic to get $ t^3 = D \pm \sqrt{D^2 - D^3} $.
If $ a_+ = D + \sqrt{D^2 - D^3} $ and $ a_- = D - \sqrt{D^2-D^3} $ then observe that $ a_+ a_- = D^3 $. The roots of the equations $ t^3 = a_+ $ and $ t^3 = a_- $ are of the form $ \alpha, \omega \alpha, \omega^2 \alpha $ and $ \beta, \omega \beta, \omega^2 \beta $ respectively and we have $ (\alpha \beta)^3 = a_+a_- = D^3 $ therefore we may fix $ \alpha $ and $ \beta $ so that $ \alpha \beta = D $.
Then we see that the roots of the original equation are $ \alpha + \frac{D}{\alpha}, \omega \alpha + \frac{D}{\omega \alpha} $ and $ \omega^2 \alpha + \frac{D}{\omega^2 \alpha} $. (The list $ \beta + \frac{D}{\beta}, \omega \beta + \frac{D}{\omega \beta} $ and $ \omega^2 \beta + \frac{D}{\omega^2 \beta} $ is a permutation of the previous list, where the last two roots are flipped)
One has to be careful about which cube root of $ t^3 = a_+ $ you pick for the expression of $ x_m^* $ that you write.