An Algorithm Architecture for Radio Interferometric Data Processing

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Bhatnagar, S., Rau, U., Hsieh, M., Kern, J., Xue, R.
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866917435298283520
author Bhatnagar, S.
Rau, U.
Hsieh, M.
Kern, J.
Xue, R.
author_facet Bhatnagar, S.
Rau, U.
Hsieh, M.
Kern, J.
Xue, R.
contents We present a foundational, scalable algorithm architecture for processing data from aperture synthesis radio telescopes. The analysis leading to the architecture is rooted in the theory of aperture synthesis, signal processing and numerical optimization keeping it scalable for variations in computing load, algorithmic complexity, and accommodate the continuing evolution of algorithms. It also adheres to scientific software design principles and use of modern performance engineering techniques providing a stable foundation for long-term scalability, performance, and development cost. We first show that algorithms for both calibration and imaging algorithms share a common mathematical foundation and can be expressed as numerical optimization problems. We then decompose the resulting mathematical framework into fundamental conceptual architectural components, and assemble calibration and imaging algorithms from these foundational components. For a physical architectural view, we used a library of algorithms implemented in the LibRA software for the various architectural components, and used the Kokkos framework in the compute-intensive components for performance portable implementation. This was deployed on hardware ranging from desktop-class computers to multiple super-computer class high-performance computing (HPC) and high-throughput computing (HTC) platforms with a variety of CPU and GPU architectures, and job schedulers (HTCondor and Slurm). As a test, we imaged archival data from the NSF's Karl G. Jansky Very Large Array (VLA) telescope in the A-array configuration for the Hubble Ultra Deep Field. Using over 100 GPUs we achieve a processing rate of ~2 Terabyte per hour to make one of the deepest images in the 2 -- 4 GHz band with an RMS noise of ~1 microJy/beam.
format Preprint
id arxiv_https___arxiv_org_abs_2508_17187
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle An Algorithm Architecture for Radio Interferometric Data Processing
Bhatnagar, S.
Rau, U.
Hsieh, M.
Kern, J.
Xue, R.
Instrumentation and Methods for Astrophysics
We present a foundational, scalable algorithm architecture for processing data from aperture synthesis radio telescopes. The analysis leading to the architecture is rooted in the theory of aperture synthesis, signal processing and numerical optimization keeping it scalable for variations in computing load, algorithmic complexity, and accommodate the continuing evolution of algorithms. It also adheres to scientific software design principles and use of modern performance engineering techniques providing a stable foundation for long-term scalability, performance, and development cost. We first show that algorithms for both calibration and imaging algorithms share a common mathematical foundation and can be expressed as numerical optimization problems. We then decompose the resulting mathematical framework into fundamental conceptual architectural components, and assemble calibration and imaging algorithms from these foundational components. For a physical architectural view, we used a library of algorithms implemented in the LibRA software for the various architectural components, and used the Kokkos framework in the compute-intensive components for performance portable implementation. This was deployed on hardware ranging from desktop-class computers to multiple super-computer class high-performance computing (HPC) and high-throughput computing (HTC) platforms with a variety of CPU and GPU architectures, and job schedulers (HTCondor and Slurm). As a test, we imaged archival data from the NSF's Karl G. Jansky Very Large Array (VLA) telescope in the A-array configuration for the Hubble Ultra Deep Field. Using over 100 GPUs we achieve a processing rate of ~2 Terabyte per hour to make one of the deepest images in the 2 -- 4 GHz band with an RMS noise of ~1 microJy/beam.
title An Algorithm Architecture for Radio Interferometric Data Processing
topic Instrumentation and Methods for Astrophysics
url https://arxiv.org/abs/2508.17187