Wiener filter reloaded: fast signal reconstruction without preconditioning
| Autor: | Doogesh Kodi Ramanah, Benjamin D. Wandelt, Guilhem Lavaux |
|---|---|
| Přispěvatelé: | Institut d'Astrophysique de Paris (IAP), Institut national des sciences de l'Univers (INSU - CNRS)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS), Institut Lagrange de Paris, Sorbonne Universités, Institut d'Astrophysique de Paris ( IAP ), Université Pierre et Marie Curie - Paris 6 ( UPMC ) -Institut national des sciences de l'Univers ( INSU - CNRS ) -Centre National de la Recherche Scientifique ( CNRS ), Sorbonne Université (SU) |
| Jazyk: | angličtina |
| Rok vydání: | 2017 |
| Předmět: |
Cosmology and Nongalactic Astrophysics (astro-ph.CO)
[ PHYS.ASTR ] Physics [physics]/Astrophysics [astro-ph] FOS: Physical sciences cosmic background radiation 01 natural sciences symbols.namesake Conjugate gradient method 0103 physical sciences Convergence (routing) 010303 astronomy & astrophysics Physics methods: statistical Series (mathematics) 010308 nuclear & particles physics Noise (signal processing) Signal reconstruction Wiener filter Spectral density Astronomy and Astrophysics methods: data analysis Data set Space and Planetary Science cosmology: observations symbols [PHYS.ASTR]Physics [physics]/Astrophysics [astro-ph] Algorithm Astrophysics - Cosmology and Nongalactic Astrophysics |
| Zdroj: | Mon.Not.Roy.Astron.Soc. Mon.Not.Roy.Astron.Soc., 2017, 468 (2), pp.1782-1793. ⟨10.1093/mnras/stx527⟩ Mon.Not.Roy.Astron.Soc., 2017, 468 (2), pp.1782-1793. 〈10.1093/mnras/stx527〉 Monthly Notices of the Royal Astronomical Society Monthly Notices of the Royal Astronomical Society, Oxford University Press (OUP): Policy P-Oxford Open Option A, 2017, 468 (2), pp.1782-1793. ⟨10.1093/mnras/stx527⟩ |
| ISSN: | 0035-8711 1365-2966 |
| DOI: | 10.1093/mnras/stx527⟩ |
| Popis: | We present a high performance solution to the Wiener filtering problem via a formulation that is dual to the recently developed messenger technique. This new dual messenger algorithm, like its predecessor, efficiently calculates the Wiener filter solution of large and complex data sets without preconditioning and can account for inhomogeneous noise distributions and arbitrary mask geometries. We demonstrate the capabilities of this scheme in signal reconstruction by applying it on a simulated cosmic microwave background (CMB) temperature data set. The performance of this new method is compared to that of the standard messenger algorithm and the preconditioned conjugate gradient (PCG) approach, using a series of well-known convergence diagnostics and their processing times, for the particular problem under consideration. This variant of the messenger algorithm matches the performance of the PCG method in terms of the effectiveness of reconstruction of the input angular power spectrum and converges smoothly to the final solution. The dual messenger algorithm outperforms the standard messenger and PCG methods in terms of execution time, as it runs to completion around 2 and 3-4 times faster than the respective methods, for the specific problem considered. 13 pages, 10 figures. Accepted for publication in MNRAS main journal |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|