I will discuss a framework to compute two-loop multiscale QCD amplitudes that combines integrand reduction, integration-by-part identities, numerical evaluation over finite fields and analytic reconstruction techniques. Recent applications of this framework to compute two-loop 2 to 3 scattering processes with an offshell leg, as well as an analytic computation of two-loop gg->tt amplitude involving massive top loop will also be presented.