Aperture synthesis by interferometry in radio astronomy allows observation of the sky by antenna arrays with otherwise inaccessible angular resolutions and sensitivities, providing a whole wealth of information for astrophysics and cosmology. At the target resolution and dynamic range of interest to upcoming telescopes, image cubes will reach close to Petabyte sizes, with data volumes orders of magnitude larger, possibly verging on the Exabyte scale. In this context, convex optimisation theory offers modern algorithmic structures that can handle extreme data volumes by distributing the computation across large computer systems, while explicitly incorporating complex image prior models to regularise the problem. Motivated by these advantages, we introduce a C++ library and set of applications, dubbed Puri-Psi, as a production implementation of the recently proposed RI imaging approach HyperSARA, previously available as a proof-of-concept MATLAB implementation. Puri-Psi enables scaling up to very large data volumes by parallelising across a large number processes and compute nodes. We validate Puri-Psi against the MATLAB implementation, and evaluate its performance in a representative high performance computing setting.
|Title of host publication||Proceedings of ISC HPC 2021|
|Number of pages||10|
|Publication status||Submitted - Dec 2020|