Abstract:Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid ( G2P2G ) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array ( AoSoA ) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 1024 3 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Accelerating AP3M-Based Computational Astrophysics Simulations with Reconfigurable Clusters

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

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

Efficient Particle-Grid Space Interpolation of an FPGA-Accelerated Particle-in-Cell Plasma Simulation

FP-AMR: A Reconfigurable Fabric Framework for Adaptive Mesh Refinement Applications

Multi-GPU Hybrid Programming Accelerated Three-Dimensional Phase-Field Model in Binary Alloy

FPGA Acceleration of TreePM N-body Simulations for Modified Newton Dynamics.

Design and implementation of reconfigurable acceleration for in-memory distributed big data computing.

Solving the Global Atmospheric Equations Through Heterogeneous Reconfigurable Platforms.

Accelerating Astrophysical Particle Simulations with Programmable Hardware (FPGA and GPU)

Acceleration for CFD Applications on Large GPU Clusters: an NPB Case Study

Efficient Ab-Initio Molecular Dynamic Simulations by Offloading Fast Fourier Transformations to FPGAs

Molecular Dynamics Simulations on High-Performance Reconfigurable Computing Systems

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

Design and Implementation of a Heterogeneous High-performance Computing Framework using Dynamic and Partial Reconfigurable FPGAs

SuperDragon: A Heterogeneous Parallel System for Accelerating 3D Reconstruction of Cryo-Electron Microscopy Images

A massively parallel and scalable multi-GPU material point method

An MPI+OpenACC-based PRM Scalar Advection Scheme in the GRAPES Model over a Cluster with Multiple CPUs and GPUs

Accelerating the 3D Euler Atmospheric Solver Through Heterogeneous CPU-GPU Platforms

Heterogeneous Reconfigurable Design for TreePM N-body Simulation

Performance Optimization of Non-Equilibrium Ionization Simulations from MapReduce and GPU Acceleration