  We present the first numerical computation of two-loop amplitudes based on
the unitarity method. As a proof of principle, we compute the four-gluon
process. We discuss the new method, analyze its numerical properties and apply
it to reconstruct the analytic form of the amplitudes. The numerical method is
universal, and can be automated to provide multi-scale two-loop computations
for phenomenologically relevant signatures at hadron colliders.
