These are the Pauli $X$, $Y$ and $Z$ matrices respectively:
$$X=\begin{bmatrix}0&1\\1&0\end{bmatrix},\ Y=\begin{bmatrix}0&-i\\i&0\end{bmatrix} \text{ and } Z = \begin{bmatrix}1&0\\0&-1\end{bmatrix}.$$
I'm trying to find the $2\times 2$ unitary matrices $U_Y$ and $U_Z$ that satisfy
$$U_YXU_Y^\dagger = Y \text{ and } U_ZXU_Z^\dagger = Z.$$
What would be a quick algorithmic method to calculate $U_Y$ and $U_Z$? This is the context.
Another approach is to keep in mind that these are similarity relations.
For example, $U_Z X U_Z^\dagger = Z$ is the same as $X=U_Z^\dagger Z U_Z$, where $Z$ is diagonal. This means that this expression is essentially the eigendecomposition of $X$, and that the elements of $Z$ are the eigenvalues of $X$.
Thus, $U_Z$ must be the set of eigenvectors of $X$. More precisely, the columns of $U_Z^\dagger$ are the eigenvectors of $X$. Note that via this simple observation you immediately get that $U_Z$ must be Hadamard $H=(X+Z)/\sqrt2$, as expected.
You can follow similar ideas for the other case. $U_Y X U_Y^\dagger=Y$ is equivalent to $X=U_Y^\dagger Y U_Y$, but also $Y=VZV^\dagger$ where $$V=\frac{1}{\sqrt2}\begin{pmatrix}1 & 1 \\ i & -i\end{pmatrix},$$ which you know immediately if you know the eigenvalues and eigenvectors of $Y$. Then, $X=U_Y^\dagger VZV^\dagger U_Y$. But then again, this means that $U_Y^\dagger V=H$ by the first argument in this answer. You conclude that $U_Y^\dagger = HV^\dagger$.