  Multi-loop scattering amplitudes are difficult to evaluate due to
singularities of the integrals involved, especially with increasing number of
loops, external legs, and mass scales. For the first time at two loops, we
enable their direct numerical integration by tackling infrared, ultraviolet and
threshold singularities simultaneously using local subtractions. We demonstrate
the feasibility of our approach by calculating previously unknown perturbative
corrections for processes of interest to the Large Hadron Collider.
