Document Type
Conference Paper
Publication Date
2022
Pages
1-19
Conference Name
International Conference on High Performance Computing (ISC High Performance 2022), May 29-June 2, 2022, Hamburg, Germany
Abstract
The task of multi-dimensional numerical integration is frequently encountered in physics and other scientific fields, e.g., in modeling the effects of systematic uncertainties in physical systems and in Bayesian parameter estimation. Multi-dimensional integration is often time-prohibitive on CPUs. Efficient implementation on many-core architectures is challenging as the workload across the integration space cannot be predicted a priori. We propose m-Cubes, a novel implementation of the well-known Vegas algorithm for execution on GPUs. Vegas transforms integration variables followed by calculation of a Monte Carlo integral estimate using adaptive partitioning of the resulting space. mCubes improves performance on GPUs by maintaining relatively uniform workload across the processors. As a result, our optimized Cuda implementation for Nvidia GPUs outperforms parallelization approaches proposed in past literature. We further demonstrate the efficiency of m-Cubes by evaluating a six-dimensional integral from a cosmology application, achieving significant speedup and greater precision than the Cuba library’s CPU implementation of Vegas. We also evaluate mCubes on a standard integrand test suite. m-Cubes outperforms the serial implementations of the Cuba and Gsl libraries by orders of magnitude speedup while maintaining comparable accuracy. Our approach yields a speedup of at least 10 when compared against publicly available Monte Carlo based GPU implementations. In summary, m-Cubes can solve integrals that are prohibitively expensive using standard libraries and custom implementations. A modern C++ interface header-only implementation makes m-Cubes portable, allowing its utilization in complicated pipelines with easy to define stateful integrals. Compatibility with non-Nvidia GPUs is achieved with our initial implementation of m-Cubes using the Kokkos framework.
Original Publication Citation
Sakiotis, I., Arumugam, K., Paterno, M., Ranjan, D., Terzić, B., & Zubair, M. (2022). m-CUBES: An efficient and portable implementation of multi-dimensional integration for GPUs. International Conference on High Performance Computing (ISC High Performance 2022), May 29- June 2, 2022 Hamburg, Germany. https://lss.fnal.gov/archive/2022/conf/fermilab-conf-22-043-ldrd-scd.pdf
Repository Citation
Sakiotis, I., Arumugam, K., Paterno, M., Ranjan, D., Terzić, B., & Zubair, M. (2022). m-CUBES: An efficient and portable implementation of multi-dimensional integration for GPUs. International Conference on High Performance Computing (ISC High Performance 2022), May 29- June 2, 2022 Hamburg, Germany. https://lss.fnal.gov/archive/2022/conf/fermilab-conf-22-043-ldrd-scd.pdf
ORCID
0000-0002-1988-0314 (Sakiotis), 0000-0002-9646-8155 (Terzić), 0000-0002-5449-1779 (Zubair)
Comments
"Fermilab reports are all published under a Creative Commons Attribution 4.0 International (CC BY 4.0) license. They may be freely distributed and copied and re-used for any purpose, but it is requested that in any subsequent use Fermilab be given appropriate acknowledgement."