Massively Parallelized Interpolated Factored Green Function Method

TL;DR

This paper introduces a highly scalable parallel implementation of the Interpolated Factored Green Function (IFGF) method, enabling efficient large-scale wave scattering computations without FFT reliance, demonstrating excellent parallel performance up to thousands of cores.

Contribution

It presents the first parallel, MPI-OpenMP implementation of the IFGF method, maintaining O(N log N) complexity and superior scalability for large wave scattering problems.

Findings

Achieves excellent weak and strong parallel scaling.

Handles problems up to 4,096 wavelengths in size.

Scales efficiently from 1 to 1,680 cores.

Abstract

This paper presents the first parallel implementation of the novel "Interpolated Factored Green Function" (IFGF) method introduced recently for the accelerated evaluation of discrete integral operators arising in wave scattering and other areas (Bauinger and Bruno, Jour. Computat. Phys., 2021). On the basis of the hierarchical IFGF interpolation strategy, the proposed (hybrid MPI-OpenMP) parallel implementation results in highly efficient data communication, and it exhibits in practice excellent parallel scaling up to large numbers of cores - without any hard limitations on the number of cores concurrently employed with high efficiency. Moreover, on any given number of cores, the proposed parallel approach preserves the O(N log N) computing cost inherent in the sequential version of the IFGF algorithm. Unlike other approaches, the IFGF method does not utilize the Fast Fourier Transform (FFT), and it is thus better suited than other methods for efficient parallelization in distributed-memory computer systems. In particular, the IFGF method relies on a "peer-to-peer" strategy wherein, at every level, field propagation is directly enacted via "exchanges" between "peer" polynomials of low and constant degree, without data accumulation in large-scale "telephone-central" mathematical constructs such as those in the Fast Multipole Method (FMM) or pure FFT-based approaches. A variety of numerical results presented in this paper illustrate the character of the proposed parallel algorithm, including excellent weak and strong parallel scaling properties in all cases considered - for problems of up to 4,096 wavelengths in acoustic size, and scaling tests spanning from 1 compute core to all 1,680 cores available in the High Performance Computing cluster used.