  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.
