Fast Power Spectrum Estimation
Astrophysics
2009-11-07 v1
Abstract
We present a new algorithm to rapidly and optimally compute power spectra. This new algorithm is based on a generalization of iterative multigrid, and has computational cost O(N log N), compared to the standard brute force approach which costs O(N^3). The procedure retains this speed on the full sky and for ill-conditioned matrices. It is applicable to galaxy power spectra, CMB, polarization and weak lensing data. We present a mathematical convergence analysis, and performance results.
Cite
@article{arxiv.astro-ph/0304513,
title = {Fast Power Spectrum Estimation},
author = {Ue-Li Pen},
journal= {arXiv preprint arXiv:astro-ph/0304513},
year = {2009}
}
Comments
11 pages, 3 figures