Astronomy and Astrophysics – Astrophysics
Scientific paper
1999-02-11
Astronomy and Astrophysics
Astrophysics
9 pages LaTeX including 2 PostScript figures. Additional discussion of conjugate gradient chi-squared algorithm. Matches accep
Scientific paper
10.1086/312151
The derivation of cosmological parameters from astrophysical data sets routinely involves operations counts which scale as O(N^3) where N is the number of data points. Currently planned missions, including MAP and Planck, will generate sky maps with N_d = 10^6 or more pixels. Simple ``brute force'' analysis, applied to such mega-pixel data, would require years of computing even on the fastest computers. We describe an algorithm which allows estimation of the likelihood function in the direct pixel basis. The algorithm uses a conjugate gradient approach to evaluate chi-squared and a geometric approximation to evaluate the determinant. Monte Carlo simulations provide a correction to the determinant, yielding an unbiased estimate of the likelihood surface in an arbitrary region surrounding the likelihood peak. The algorithm requires O(N_d^{3/2}) operations and O(N_d) storage for each likelihood evaluation, and allows for significant parallel computation.
No associations
LandOfFree
Likelihood Analysis for Mega-Pixel Maps does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Likelihood Analysis for Mega-Pixel Maps, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Likelihood Analysis for Mega-Pixel Maps will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-461283