Efficient wide-field radio interferometry response

Abstract

Radio interferometers do not measure the sky brightness distribution directly but rather a modified Fourier transform of it. Imaging algorithms, thus, need a computational representation of the linear measurement operator and its adjoint, irrespective of the specific chosen imaging algorithm. In this paper, we present a C++ implementation of the radio interferometric measurement operator for wide-field measurements which is based on "improved $w$-stacking". It can provide high accuracy (down to $\approx 10^{-12}$), is based on a new gridding kernel which allows smaller kernel support for given accuracy, dynamically chooses kernel, kernel support and oversampling factor for maximum performance, uses piece-wise polynomial approximation for cheap evaluations of the gridding kernel, treats the visibilities in cache-friendly order, uses explicit vectorisation if available and comes with a parallelisation scheme which scales well also in the adjoint direction (which is a problem for many previous implementations). The implementation has a small memory footprint in the sense that temporary internal data structures are much smaller than the respective input and output data, allowing in-memory processing of data sets which needed to be read from disk or distributed across several compute nodes before.

arXiv comment13 pages, 8 figures

References in corpus (4)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• Image Domain Gridding: a fast method for convolutional resampling of visibilities
• Comparison of classical and Bayesian imaging in radio interferometry
• Optimal gridding and degridding in radio interferometry imaging

Cited by in corpus (24)

• The MPIfR-MeerKAT Galactic Plane survey I -- System setup and early results
• Improved cosmic microwave background (de-)lensing using general spherical harmonic transforms
• Deep Images of the Galactic Center with GRAVITY
• High accuracy wide field imaging method in radio interferometry
• Forward-fitting STIX visibilities
• Bayesian radio interferometric imaging with direction-dependent calibration
• Data Processing Pipeline For Tianlai Experiment
• Ultra-Steep Spectrum Radio `Jellyfish' Uncovered in Abell 2877
• The merging galaxy cluster Abell 3266 at low radio frequencies
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• The prototypical major cluster merger Abell 754. I. Calibration of MeerKAT data and radio/X-ray spectral mapping of the cluster
• Transient RFI environment of LOFAR-LBA at 72-75 MHz: Impact on ultra-widefield AARTFAAC Cosmic Explorer observations of the redshifted 21-cm signal
• OmniUV: A Multi-Purpose Simulation Toolkit for VLBI Observation
• 1-arcsecond imaging of ELAIS-N1 field at 144MHz using the LoTSS survey with international LOFAR telescope
• Validation Solutions to the Full-Sky Radio Interferometry Measurement Equation for Diffuse Emission
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• First upper limits on the 21-cm signal power spectrum of neutral hydrogen at $z=9.16$ from the LOFAR 3C196 field
• A novel radio imaging method for physical spectral index modelling
• Revealing the intricacies of radio galaxies and filaments in the merging galaxy cluster Abell 2255. I. Insights from deep LOFAR-VLBI sub-arcsecond resolution images
• WS-Snapshot: An effective algorithm for wide-field and large-scale imaging
• Relighting the fire in Hickson Compact Group (HCG) 15: magnetised fossil plasma revealed by the SKA Pathfinders & Precursors
• Imaging swiFTly: streaming widefield Fourier Transforms for large-scale interferometry
• A practical preconditioner for wide-field continuum imaging of radio interferometric data
• Cosmic-ray electron propagation in NGC 3044 from radio continuum observations