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.

A Parallel 1-D FFT Implementation Method for Multi-core Vector Processors

Efficient Large-Scale 1D FFT Vectorization on Multi-Core Vector Accelerator

Software and Hardware Cooperate for 1-D FFT Algorithm Optimization on Multicore Processors

A Parallel Radix-2k FFT Processor Using Single-Port Merged-Bank Memory

An Optimized Parallel FFT Algorithm on Multiprocessors with Cache Technology in Linux

The Implementation and Optimization of Parallel Linpack on Multi-Core Vector Accelerator

Memory Efficient Two-Pass 3D FFT Algorithm for Intel® Xeon Phi TM Coprocessor

Parallel FFT Implementation Based on Multi-Core DSPs

A multi-core mapping implementation of 3780-point FFT

Task Parallel Implementation of Matrix Multiplication on Multi-socket Multi-core Architectures.

An Efficient SIMD Parallel Memory Structure for Radix-2 FFT Computation

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

Adaptive Hybrid FFT: A Novel Pipeline and Memory-Based Architecture for Radix-$2^k$ FFT in Large Size Processing

MFFT: A GPU Accelerated Highly Efficient Mixed-precision Large-scale FFT Framework

Algorithm and Implementation of a Large-Point FFT under the Master-Slave Parallel Multi-Processor Architecture

Parallel computing method for two-dimensional matrix convolution

Efficient implementation of 3780-point FFT on a 16-core processor

Implementation of Prime-Point DFTs with Programable Processors

Some New Parallel Fast Fourier Transform Algorithms

FT-Matrix: A Coordination-Aware Architecture for Signal Processing