Sublinear scaling method for indirect spin–spin coupling constant calculations at Hartree–Fock and density functional theory level