  We describe a new method to perform NLO calculations, combining real and
virtual amplitudes at the integrand level, with a fully local compensation
between them in the IR, and between the virtual integrand and properly defined
counter-terms in the UV, in such a way that physical observables can be
computed in 4 dimensions. One of the advantages of the method is that all the
scattering amplitudes are integrated simultaneously, without the need for
tensor reduction, or projection onto sets of master integrals. As such, it
could offer great progress in the automation of NLO calculations, where the
actual bottle-necks are the complexity of the analytical calculations for
multi-leg processes, and the numerical stability of the result.
