and the full SO matrix is a product of Clebsh-Gordan's coefficients to these reduced matrix elements. The prefactors in the expressions above arise because the spin part is taken in spin-tensor form whereas the orbital momentum is taken in L + /L -, L z form, which is not a spin-tensor form. The overall workflow of the algorithm is shown in Fig. .