All roots of $8x^3-6x+1$ are real. (*)
The discriminant of $8x^3-6x+1$ is $5184=72^2$ and so the splitting field of $8x^3-6x+1$ has degree $3$.
Therefore, all three roots can be expressed as polynomials in any one given root.
Indeed, if $a$ is a root, then the others are $2a^2-1$ and $-2a^2-a+1$.
This can be easily checked. But how can we find these expressions in the first place?
I've tried this: let $b,c$ be the other roots. Then from Vieta's formulas we get $$b+c=-a, \qquad bc=-\dfrac{1}{8a}$$
The corresponding quadratic has discriminant $3-3a^2$ but it is not at all clear that this is the square of a polynomial in $a$. (It is $(4 a^2+a-2 )^2=(b-c)^2$, as it should be.) I'm stuck here.
More generally, is there an algorithm that given a cubic with integer coefficients and having real splitting field of degree $3$, expresses all three roots as polynomials in any one given root?
(*) From the triple-angle formula $\sin (3\theta) = - 4\sin^3\theta + 3\sin\theta$ when $\sin (3\theta) = 1/2$, these roots are $\sin(10^\circ)$, $\sin(50^\circ)$, $\sin(-70^\circ)$, but perhaps this is immaterial here.