Abstract:Gridding is the performance-critical step in the data reduction pipeline for radio astronomy research, allowing astronomers to create the correct sky images for further analysis. Like the 2D stencil computation, gridding iteratively updates the output cells by convolution, where the value at each output cell in the space is computed as a weighted sum of neighboring point values. Existing state-of-the-art works have achieved performance improvement of gridding by using multi-core CPUs and GPUs in real-world applications, and their study proved that gridding is a type of scientific computation with high-density computing characteristics. However, low computational performance or high power consumption becomes the main limitation for their processing of large-scale astronomical data. The high-density computing feature of gridding provides opportunities to accelerate it on the multi-core vector processor with vector-SIMD architectures. However, existing works’ (such as those implemented on CPUs or GPUs) task parallelization and data transfer strategies are inefficient to perform gridding directly on the vector processor without any dedicated mapping algorithm. M-DSP is a multi-core vector processor with vector-SIMD architectures designed for the next-generation exascale supercomputer, delivering high performance with ultra-low power consumption. In this paper, we present, for the first time, a novel method to achieve efficient gridding on the M-DSP. Specifically, we propose a gridding workflow designed for the vector-SIMD architectures and present a vectorized version of the gridding convolution algorithm to fully exploit the computational power of the M-DSP. In addition, centering on the processor architectures, we propose task-based parallelization strategies for block and line computing as well as different data loading strategies to achieve high parallel performance and high data transfer efficiency. Experimental results show that our work on M-DSP exhibits very competitive performance compared to other methods running on CPUs or GPUs. This demonstrates the efficiency of our method and the fact that the vector-SIMD architecture is beneficial for scientific computing with ”high density” characteristics, which can exploit its wide vector core and achieve higher performance than its competitors.

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

HCGrid: a Convolution-Based Gridding Framework for Radio Astronomy in Hybrid Computing Environments

HEGrid: A High Efficient Multi-Channel Radio Astronomical Data Gridding Framework in Heterogeneous Computing Environments

Efficient Radio Interferometric Imaging on the GPU

A Method for Efficient Radio Astronomical Data Gridding on Multi-Core Vector Processor

High-Performance Gridding For Radio Interferometric Image Synthesis

CPU/GPU Computing for Long-Wave Radiation Physics on Large GPU Clusters.

Heterogeneous CPU-GPU Accelerated Parallel Subgridding FDTD Algorithm

Parallel Multi-threaded Gridrec Algorithm for Computer Tomography on GPU for Edge Computing

Mapping the Sbr and Tw-Ildcs to Heterogeneous Cpu-Gpu Architecture for Fast Computation of Electromagnetic Scattering

GPU-Accelerated Radial Image Reconstruction with an Improved Parallel Gridding Method

FROM CPU TO GPU: GPU-BASED ELECTROMAGNETIC COMPUTING (GPUECO)

Design and Optimization of Remote Sensing Image Fusion Parallel Algorithms Based on CPU-GPU Heterogeneous Platforms

Numerical Study of Geometric Multigrid Methods on CPU--GPU Heterogeneous Computers

Exploring High Efficiency Hardware Accelerator for the Key Algorithm of Square Kilometer Array Telescope Data Processing

A Hierarchical Grid Algorithm for Accelerating High-Performance Conjugate Gradient Benchmark on Sunway Many-Core Processor

Accelerating Radio Astronomy Cross-Correlation with Graphics Processing Units

Accelerating Spectral Calculation Through Hybrid GPU-Based Computing

Parallel Implementation of W-Projection Wide-Field Imaging

GPU-RRTMG_SW: Accelerating a Shortwave Radiative Transfer Scheme on GPU

Accelerating Faceting Wide-Field Imaging Algorithm with FPGA for SKA Radio Telescope As a Vast Sensor Array.