Wiener filter reloaded: fast signal reconstruction without preconditioning

Abstract

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

References in corpus (9)

• Planck 2018 results. VII. Isotropy and Statistics of the CMB
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Reconstructed Density and Velocity Fields from the 2MASS Redshift Survey
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• Estimation of Polarized Power Spectra by Gibbs sampling
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• Efficient Wiener filtering without preconditioning
• Pure E and B polarization maps via Wiener filtering
• Fast Wiener filtering of CMB maps

Cited by in corpus (15)

• Review: Deep Learning in Electron Microscopy
• Bayesian optimisation for likelihood-free cosmological inference
• Bayesian delensing of CMB temperature and polarization
• Efficient Optimal Reconstruction of Linear Fields and Band-powers from Cosmological Data
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Pure E and B polarization maps via Wiener filtering
• Optimal and fast E/B separation with a dual messenger field
• Cosmic Microwave Background Mapmaking with a Messenger Field
• Solving linear equations with messenger-field and conjugate gradients techniques - an application to CMB data analysis
• B-mode Power Spectrum of CMB via Polarized Compton Scattering
• Preconditioner-free Wiener filtering with a dense noise matrix
• Enhancing CMB map reconstruction and power spectrum estimation with convolutional neural networks
• Pixel domain multi-resolution minimum variance painting of CMB maps
• Cooling Improves Cosmic Microwave Background Map-Making When Low-Frequency Noise is Large