Toward a GPU-enabled billionaire SVD in pyLOM

We develop and implement an accelerated high-performance and open-source computing environment for model order reduction in fluid dynamics called pyLOM. It contains singular value decompositionbased algorithms implemented for massively parallel GPU architectures. The library is profiled in detail un...

Full description

Bibliographic Details
Authors: Miró Jané, Arnau|||0000-0002-2772-6050; Eiximeno Franch, Benet; Silva, Lucas Gasparino Ferreira da; Kutz, Nathan; Rodríguez Pérez, Ivette María|||0000-0002-3749-277X; Lehmkuhl Barba, Oriol|||0000-0002-2670-1871
Format: Article
Publication date:2026
Country:España
Institution:Universitat Politècnica de Catalunya (UPC)
Repository:UPCommons. Portal del coneixement obert de la UPC
Language:English
OAI Identifier:oai:upcommons.upc.edu:2117/452499
Online Access:https://hdl.handle.net/2117/452499
https://dx.doi.org/10.1007/s00707-025-04621-1
Access Level:Open access
Keyword:Data reduction
Energy efficiency
Environmental impact
Graphics processing unit
Matrix factorization
Parallel architectures
Principal component analysis
Program processors
Supercomputers
Àrees temàtiques de la UPC::Informàtica::Arquitectura de computadors::Arquitectures paral·leles
Àrees temàtiques de la UPC::Física::Física de fluids
Description
Summary:We develop and implement an accelerated high-performance and open-source computing environment for model order reduction in fluid dynamics called pyLOM. It contains singular value decompositionbased algorithms implemented for massively parallel GPU architectures. The library is profiled in detail under the MareNostrum V supercomputer. The largest case has been computed under 20s with 100 GPUs and consisted of a billion nodes by a thousand snapshots matrix. A hybrid CPU-GPU parallel randomized QR factorization has been found to be able to leverage such large matrices. The largest speedup factor of 83 has been found on the QR factorization, while the matrix–matrix multiplication has shown a speedup factor of about 2. Additionally, two examples of application are provided in the flow around a cylinder and the Windsor body, whose PODiscomputed under 3s with 100 GPUs. This showcases the efficiency of GPUs, resulting in a 97% reduction in energy to solution and a reduction of 0.11 kg of CO2 emissions. The scalability and efficiency achieved suggest that this framework can play a key role in reducing the energy demands and environmental impact of large-scale data analysis and model order reduction across a wide range of applications.