Extending the perturbed matrix method beyond the dipolar approximation: comparison of different levels of theory