  Some of the tools and techniques that have recently been used to compute
Higgs boson pair production at NLO in QCD are discussed. The calculation relies
on the use of integral reduction, to reduce the number of integrals which must
be computed, and expressing the amplitude in terms of a quasi-finite basis,
which simplifies their numeric evaluation. Emphasis is placed on sector
decomposition and Quasi-Monte Carlo (QMC) integration which are used to
numerically compute the master integrals.
