Saved in:
Bibliographic Details
Main Authors: Tonarelli, Melanie, Riva, Simone, Benedusi, Pietro, Ferrandi, Fabrizio, Krause, Rolf
Format: Preprint
Published: 2025
Subjects:
Online Access:https://arxiv.org/abs/2511.01573
Tags: Add Tag
No Tags, Be the first to tag this record!
Table of Contents:
  • We introduce a distributed adaptive quadrature method that formulates multidimensional integration as a hierarchical domain decomposition problem on multi-GPU architectures. The integration domain is recursively partitioned into subdomains whose refinement is guided by local error estimators. Each subdomain evolves independently on a GPU, which exposes a significant load imbalance as the adaptive process progresses. To address this challenge, we introduce a decentralised load redistribution schemes based on a cyclic round-robin policy. This strategy dynamically rebalance subdomains across devices through non-blocking, CUDA-aware MPI communication that overlaps with computation. The proposed strategy has two main advantages compared to a state-of-the-art GPU-tailored package: higher efficiency in high dimensions; and improved robustness w.r.t the integrand regularity and the target accuracy.