Numerical Scattering Amplitudes with pySecDec

We present a major update of the program pySecDec, a toolbox for the evaluation of dimensionally regulated parameter integrals. The new version enables the evaluation of multi-loop integrals as well as amplitudes in a highly distributed and flexible way, optionally on GPUs. The program has been optimised and runs up to an order of magnitude faster than the previous release. A new integration procedure that utilises construction-free median Quasi-Monte Carlo rules is implemented. The median lattice rules can outperform our previous component-by-component rules by a factor of 5 and remove the limitation on the maximum number of sampling points. The expansion by regions procedures have been extended to support Feynman integrals with numerators, and functions for automatically determining when and how analytic regulators should be introduced are now available. The new features and performance are illustrated with several examples.