Return
An Algorithm Architecture for Radio Interferometric Data Processing
DOI:10.3847/1538-3881/adfe61.png)
Abstract
En 中文
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 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 supercomputer class high-performance computing 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 Very Large Array telescope in the A-array configuration for the Hubble Ultra Deep Field. Using over 100 GPUs, we achieve a processing rate of similar to 2 TB hr-1 to make one of the deepest images in the 2-4 GHz band with an rms noise of similar to 1 mu Jy beam-1.
Keywords:
DECONVOLUTION ALGORITHM
CALIBRATION
POLARIMETRY
Journal
IF:
5.1
Papers:
1.4W
Citations:
5.0W

