HCGrid: A Convolution-based Gridding Framework for RadioAstronomy in Hybrid Computing Environments
arXiv:2012.13173 · doi:10.1093/mnras/staa3800
Abstract
Gridding operation, which is to map non-uniform data samples onto a uniformly distributedgrid, is one of the key steps in radio astronomical data reduction process. One of the mainbottlenecks of gridding is the poor computing performance, and a typical solution for suchperformance issue is the implementation of multi-core CPU platforms. Although such amethod could usually achieve good results, in many cases, the performance of gridding is stillrestricted to an extent due to the limitations of CPU, since the main workload of gridding isa combination of a large number of single instruction, multi-data-stream operations, which ismore suitable for GPU, rather than CPU implementations. To meet the challenge of massivedata gridding for the modern large single-dish radio telescopes, e.g., the Five-hundred-meterAperture Spherical radio Telescope (FAST), inspired by existing multi-core CPU griddingalgorithms such as Cygrid, here we present an easy-to-install, high-performance, and open-source convolutional gridding framework, HCGrid,in CPU-GPU heterogeneous platforms. Itoptimises data search by employing multi-threading on CPU, and accelerates the convolutionprocess by utilising massive parallelisation of GPU. In order to make HCGrid a more adaptivesolution, we also propose the strategies of thread organisation and coarsening, as well as optimalparameter settings under various GPU architectures. A thorough analysis of computing timeand performance gain with several GPU parallel optimisation strategies show that it can leadto excellent performance in hybrid computing environments.
arXiv commentThis is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society, 11 pages, 13 figures
References in corpus (14)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Astropy: A Community Python Package for Astronomy
- The NumPy array: a structure for efficient numerical computation
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI: Final data release of the combined LDS and IAR surveys with improved stray-radiation corrections
- The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project
- The Arecibo Legacy Fast ALFA Survey: I. Science Goals, Survey Design and Strategy
- Representations of celestial coordinates in FITS
- Considerations for a Multi-beam Multi-purpose Survey with FAST
- GASS: The Parkes Galactic All-Sky Survey. II. Stray-Radiation Correction and Second Data Release
- The Effelsberg-Bonn HI Survey: Milky Way gas. First data release
- Deep residual detection of radio frequency interference for FAST
- Cygrid: A fast Cython-powered convolution-based gridding module for Python
- Faster GPU-based convolutional gridding via thread coarsening