In $\mathbb R^4$ a reflect in a hyperplane through the origin $O$ is, $\forall q \in \mathbb R^4$ $$q \mapsto -u\overline{q}u$$ , where $u$ is a unit quaternion.
In $\mathbb R^3$ a rotation fixing the origin $O$ is $$p \mapsto x^{-1}px$$ , where $p$ is a pure quaternion, and $x$ is a unit quaternion.
Since a rotation is equivalent to 2 reflects chained up, given unit quaternion $u$ for a rotation $p \mapsto u^{-1}pu$, I'm trying to find unit quaternions $x$, $y$ so that the rotation can be decomposed to two reflects, first via $x$, then via $y$, i.e.
$$p \mapsto -x\overline{p}x \mapsto -y \overline{-x\overline{p}x} y = yx^{-1}px^{-1}y$$
Let this equals to $u^{-1}pu$ we get
$$x^{-1}y=u, \qquad yx^{-1}=u^{-1}$$
Please enlighten me how to find such a pair of $x$ and $y$ ?