Abstract:SIAM Journal on Scientific Computing, Volume 43, Issue 1 , Page C31-C60, January 2021. A Fourier-based library of unbounded Poisson solvers (FLUPS) for 2D and 3D homogeneous distributed grids is presented. It is designed to handle every possible combination of periodic, symmetric, semi-unbounded, and fully unbounded boundary conditions for the Poisson equation on rectangular domains with uniform resolution. FLUPS leverages a dedicated implementation of 3D Fourier transforms to solve the Poisson equation using Green's functions in a fast and memory-efficient way. Several Green's functions are available, optionally with explicit regularization, spectral truncation, or using lattice Green's functions, and provide verified convergence orders from 2 to spectral-like. The algorithm depends on the FFTW library to perform 1D transforms, while message passing interface (MPI) communications enable the required remapping of data in memory. For the latter operation, a first available implementation resorts to the standard all-to-all routines. A second implementation, featuring non-blocking and persistent point-to-point communications, is however shown to be more efficient in a majority of cases and especially while taking advantage of the shared memory parallelism with OpenMP. The scalability of the algorithm, aimed at massively parallel architectures, is demonstrated up to 73720 cores. The results obtained with three different supercomputers show that the weak efficiency remains above and the strong efficiency above when the number of cores is multiplied by 16, for typical problems. These figures are slightly better than those expected from a third party 3D fast Fourier transform (FFT) tool, with which a longer execution time was also measured on average. From the outside, the solving procedure is fully automated so that the user benefits from the optimal performances while not having to handle the complexity associated with memory management, data mapping, and Fourier transform computation. The parallel code is available under Apache license 2.0 at github.com/vortexlab-uclouvain/flups.

FLUPS -- a flexible and performant massively parallel Fourier transform library

FLUPS: A Fourier-Based Library of Unbounded Poisson Solvers

FluidFFT: Common API (C++ and Python) for Fast Fourier Transform HPC Libraries

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

Flexible Multi-Dimensional FFTs for Plane Wave Density Functional Theory Codes

AccFFT: A library for distributed-memory FFT on CPU and GPU architectures

A FFT-based finite-difference solver for massively-parallel direct numerical simulations of turbulent flows

cuFINUFFT: a load-balanced GPU library for general-purpose nonuniform FFTs

Accelerating the computation of FLAPW methods on heterogeneous architectures

Hybrid CPU-GPU generation of the Hamiltonian and Overlap matrices in FLAPW methods

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

Scalable Multi-node Fast Fourier Transform on GPUs

VkFFT-A Performant, Cross-Platform and Open-Source GPU FFT Library

FFT for the APE Parallel Computer

SailFFish: A Lightweight, Parallelised Fast Poisson Solver Library

hyper.deal: An efficient, matrix-free finite-element library for high-dimensional partial differential equations

MPFFT：An Auto-Tuning FFT Library for OpenCL GPUs

The uniform sparse FFT with application to PDEs with random coefficients

Distributed-memory parallelization of the aggregated unfitted finite element method

A Massively Parallel Performance Portable Free-space Spectral Poisson Solver

A Flexible Framework for Parallel Multi-Dimensional DFTs