Fast Cosmological Parameter Estimation from Microwave Background Temperature and Polarization Power Spectra

We improve the algorithm of Kosowsky, Milosavljevic, and Jimenez (2002) for computing power spectra of the cosmic microwave background. The present algorithm computes not only the temperature power spectrum but also the E-mode polarization and the temperature-polarization cross power spectra, providing the accuracy required for current cosmological parameter estimation. We refine the optimum set of cosmological parameters for computing the power spectra as perturbations around a fiducial model, leading to an accuracy better than 0.5% for the temperature power spectrum throughout the region of parameter space within WMAP's first-year 3sigma confidence region. This accuracy is comparable to the difference between the widely-used CMBFAST code (Seljak and Zaldarriaga 1996) and Boltzmann codes. Our algorithm (CMBwarp) makes possible a full exploration of the likelihood region for eight cosmological parameters in about one hour on a laptop computer. We provide the code to compute power spectra as well as the Markov Chain Monte Carlo algorithm for cosmological parameters estimation at http://www.physics.upenn.edu/~raulj/CMBwarp