In order to exploit the full potential of the LHC and future colliders, high-precision calculations of a very wide range of observables are crucial.
Automated tools for next-to-next-to-leading order calculations of perturbative scattering amplitudes are therefore a highly desirable goal.
So far, we have developed several major ingredients of such a tool in the OpenLoops framework.

In our approach we split the calculation of scattering amplitudes into three components: loop momentum tensor integrals, the corresponding process-dependent tensor coefficients, and the interplay of $(D-4)$-dimensional parts of the integrand with divergences of the integrals.
In these proceedings, we present a new recursive algorithm to reduce arbitrary two-loop tensor integrals to scalar integrals, which has been implemented into an efficient numerical tool.
We also implemented a first version of the subsequent reduction to master integrals, which allows for a full validation of the algorithm.
We present a first successful validation computation and discuss its dependence on the precision of the externally computed master integrals.