Accurate yet feasible post-Hartree-Fock computation of magnetic interactions in large biradicals through a combined variational/perturbative approach: Setup and validation.