Here we describe how to generate all required SOC matrix elements from a minimal amount of calculations. The theory is formulated in terms of reduced one-particle density matrices, such that it is ansatz-agnostic and can be applied to any electronic structure method that can furnish required transition densities. We illustrate this approach by using the EOM-CC method with single and double excitations (EOM-CCSD) as an example.