fast-resolve: Fast Bayesian radio interferometric imaging

¹ Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany
² Ludwig-Maximilians-Universität, Geschwister-Scholl-Platz 1, 80539 Munich, Germany
³ Technische Universität München (TUM), Boltzmannstr. 3, 85748 Garching, Germany
⁴ South African Radio Astronomy Observatory (SARAO), Cape Town, 7925, South Africa
⁵ Centre for Radio Astronomy Techniques and Technologies (RATT), Department of Physics and Electronics, Rhodes University, Makhanda 6140, South Africa
⁶ Institute for Radioastronomy, National Institute of Astrophysics (INAF IRA), Via Gobetti 101, 40129 Bologna, Italy

^★ Corresponding author; This email address is being protected from spambots. You need JavaScript enabled to view it.

Received: 13 June 2024
Accepted: 9 September 2024

Abstract

Context. The aim of this paper is to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets.

Aims. By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm that we named fast-resolve.

Methods. We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset, which is computationally out of reach for resolve.

Results. The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically, for the single channel VLA Cygnus A datasets fast-resolve is about 144 times faster than resolve. For the MeerKAT dataset with multiple channels the computational increase in speed of fast-resolve is even greater.