With a key improvement, the auxiliary mass flow method is now able to compute many badly-needed Feynman integrals encountered in cutting-edge collider processes. We have successfully applied it to two-loop electroweak correction to $e^+e^-to HZ$, two-loop QCD corrections to $3j$, $W/Z/H+2j$, $tbar{t}H$ and $4j$ production at hadron colliders, and three-loop QCD correction to $tbar{t}$ production at hadron colliders, all of which are crucial for precision test in the following decade. Our results are important building blocks and benchmarks for future study of these processes.