Cartan-Khaneja-Glaser decomposition of $SU(2^n)$ via involutive automorphisms

We present a novel algorithm for performing the Cartan-Khaneja-Glaser decomposition of unitary matrices in $\mathrm{SU}\left(2^n\right)$, a critical task for efficient quantum circuit design. Building upon the approach introduced by Sá Earp and Pachos (2005), we overcome key limitations of their method, such as reliance on ill-defined matrix logarithms and the convergence issues of truncated Baker-Campbell-Hausdorff (BCH) series. Our reformulation leverages the algebraic structure of involutive automorphisms and symmetric Lie algebra decompositions to yield a stable and recursive factorization process. We provide a full Python implementation of the algorithm, available in an open-source repository, and validate its performance on matrices in $\mathrm{SU}(8)$ and $\mathrm{SU}(16)$ using random unitary benchmarks. The algorithm produces decompositions that are directly suited to practical quantum hardware, with factors that can be implemented near-optimally using standard gate sets.