Suppose I have a skew-symmetric, real $4\times4$ matrix $M$. This implies that all its eigenvalues are imaginary and come in pairs. Let us call them $\pm i \lambda_1$ and $\pm i \lambda_2$ with $\lambda_1, \lambda_2 \in \mathbb{R}$ and let us call the corresponding eigenvectors $v^\pm_1$ and $v^\pm_2$ such that $M v^-_2= -i\lambda_2 v^-_2$, and so on.
Using these four eigenvectors one can construct the matrix $S = (v^+_1, v^-_1, v^+_2, v^-_2)$ such that $J = S^{-1}MS = \text{diag}( i\lambda_1, -i\lambda_1, i\lambda_2, -i\lambda_2 )$ is the diagonalized matrix. So far so good.
I seem to recall that under some additional circumstances (i.e. $\text{rank}M = 4$?) one can find a special orthogonal matrix $R \in SO(4)$ such that
$D = R^T M R = \begin{pmatrix} 0 & i\lambda_1 & 0 & 0 \\ -i\lambda_1 & 0 & 0 & 0 \\ 0 & 0 & 0 & i\lambda_2 \\ 0 & 0 & -i\lambda_2 & 0 \end{pmatrix}$.
//edit: The matrix should be brought to the form
$D = R^T M R = \begin{pmatrix} 0 & \lambda_1 & 0 & 0 \\ \lambda_1 & 0 & 0 & 0 \\ 0 & 0 & 0 & \lambda_2 \\ 0 & 0 & \lambda_2 & 0 \end{pmatrix}$.
My question: how do I construct this matrix $R$ from the four eigenvectors $v^\pm_{1,2}$ that I already know? This basis is sometimes called a Darboux basis I think.