Different approximations have been proposed in the literature to reduce the time needed for SOC integrals evaluation . We propose here a new approach, which allows for significant speed ups of the calculations without loss of accuracy. This can achieved by the use of the Cholesky representation for the 2e-SOC integrals recently presented in the literature by some of us :