We present a specific prescription for the calculation of cosmological power
spectra, exploited here at two-loop order in perturbation theory (PT), based on
the multi-point propagator expansion. In this approach power spectra are
constructed from the regularized expressions of the propagators that reproduce
both the resummed behavior in the high-k limit and the standard PT results at
low-k. With the help of N-body simulations, we show that such a construction
gives robust and accurate predictions for both the density power spectrum and
the correlation function at percent-level in the weakly non-linear regime. We
then present an algorithm that allows accelerated evaluations of all the
required diagrams by reducing the computational tasks to one-dimensional
integrals. This is achieved by means of pre-computed kernel sets defined for
appropriately chosen fiducial models. The computational time for two-loop
results is then reduced from a few minutes, with the direct method, to a few
seconds with the fast one. The robustness and applicability of this method are
tested against the power spectrum cosmic emulator from which a wide variety of
cosmological models can be explored. The fortran program with which direct and
fast calculations of power spectra can be done, RegPT, is publicly released as
part of this paper.Comment: 28 pages, 15 figure