English

Automated Translation and Accelerated Solving of Differential Equations on Multiple GPU Platforms

Distributed, Parallel, and Cluster Computing 2023-11-20 v3 Mathematical Software Numerical Analysis Numerical Analysis

Abstract

We demonstrate a high-performance vendor-agnostic method for massively parallel solving of ensembles of ordinary differential equations (ODEs) and stochastic differential equations (SDEs) on GPUs. The method is integrated with a widely used differential equation solver library in a high-level language (Julia's DifferentialEquations.jl) and enables GPU acceleration without requiring code changes by the user. Our approach achieves state-of-the-art performance compared to hand-optimized CUDA-C++ kernels while performing 20--100×\times faster than the vectorizing map (vmap) approach implemented in JAX and PyTorch. Performance evaluation on NVIDIA, AMD, Intel, and Apple GPUs demonstrates performance portability and vendor-agnosticism. We show composability with MPI to enable distributed multi-GPU workflows. The implemented solvers are fully featured -- supporting event handling, automatic differentiation, and incorporation of datasets via the GPU's texture memory -- allowing scientists to take advantage of GPU acceleration on all major current architectures without changing their model code and without loss of performance. We distribute the software as an open-source library https://github.com/SciML/DiffEqGPU.jl

Keywords

Cite

@article{arxiv.2304.06835,
  title  = {Automated Translation and Accelerated Solving of Differential Equations on Multiple GPU Platforms},
  author = {Utkarsh Utkarsh and Valentin Churavy and Yingbo Ma and Tim Besard and Prakitr Srisuma and Tim Gymnich and Adam R. Gerlach and Alan Edelman and George Barbastathis and Richard D. Braatz and Christopher Rackauckas},
  journal= {arXiv preprint arXiv:2304.06835},
  year   = {2023}
}

Comments

14 figures

R2 v1 2026-06-28T10:05:28.898Z