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