The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison