Simultaneous decomplexification of a pair of complex matrices via a unitary similarity transformation

Let $(A,B)$ be a given pair of complex $n\times n$ matrices and let at least one of these matrices be unitarily irreducible. An algorithm for verifying whether $A$ and $B$ can be made real via the same unitary similarity transformation is proposed and justified.