[ca] Els Models del Sistema Terrestre (ESMs) són crucials per a comprendre el canvi climàtic, ja que brinden informació valuosa sobre fenòmens i sistemes del món real, amb nombroses aplicacions en ciència i enginyeria. No obstant això, la comunitat científica requereix una major potència computacional per a reduir les incerteses en aquests models. S'estan aconseguint avanços significatius mitjançant la transició a sistemes informàtics heterogenis que combinen Unitats Centrals de Processament (CPUs) i Unitats de Processament Gràfic (GPUs). Aquests models es poden accelerar notablement mitjançant l'optimització dels components que consumeixen més temps, com la química, així reduint les incerteses i millorant la nostra comprensió de la dinàmica climàtica. L'objectiu principal d'aquesta tesi és millorar el rendiment dels solucionadors de química en models atmosfèrics mitjançant el desenvolupament de solucions eficients adaptades a arquitectures heterogènies basades en GPU. Específicament, la tesis es centra en maximitzar el rendiment de la GPU, millorar la precisió, garantir la portabilitat i minimitzar l'esforç de desenvolupament mitjançant la realització de modificacions menors al solucionador de química existent. La contribució fonamental d'aquesta tesi és el desenvolupament d'una versió CPU-GPU del solucionador de química CAMP per al model atmosfèric MONARCH, aconseguint una acceleració significativa. A més, el solver químic desenvolupat es pot adaptar fàcilment per a resoldre sistemes genèrics d'Equacions Diferencials Ordinàries (ODE) a la GPU, gràcies a la seva similitud amb la llibreria CVODE utilitzada per al solver CPU de CAMP. El codi font obert s'ha publicat com una actualització de CAMP per a facilitar la portabilitat. Les principals contribucions inclouen un mètode per integrar fàcilment codi GPU als solucionadors de química (Multi-cells), distribuir la càrrega computacional entre els processos de la GPU (Block-cells) i un algoritme d'equilibri de càrrega automàtic per a l'execució concurrent de la CPU i la GPU. Particularment, l'algoritme d'equilibri de càrrega automàtic està dissenyat per ésser simple, garantint la portabilitat i oferint una execució accelerada. A més, s'ha desenvolupat un solucionador lineal de GPU seguint l'enfocament Block-cells, la qual cosa permet una major exploració de Block-cells en solucionadors iteratius i facilita l'acoblament amb altres sistemes, una consideració important donat l'ús generalitzat de solucionadors lineals en les comunitats científiques i matemàtiques. Els resultats demostren una acceleració de 250x en comparació amb la versió de CPU d'un sol procés utilitzant quatre GPU en un sol node. En una comparació de node a node, utilitzant quatre GPU enfront de 80 nuclis de CPU, l'acceleració obtinguda és de 8,14x. Normalitzant els nuclis utilitzats enfront d'altres solucions d'última generació, l'acceleració és de 16x. L'acceleració substancial aconseguida amb canvis algorítmics mínims subratlla l'importància de les noves estratègies desenvolupades en aquesta tesi.

[en] Earth System Models (ESMs) are crucial for understanding climate change, providing valuable insights into real-world phenomena and systems, with numerous applications in science and engineering. However, the scientific community requires greater computational power to reduce uncertainties in these models. Significant progress is being made by transitioning to heterogeneous computing systems that combine Central Processing Units (CPUs) and Graphics Processing Units (GPUs). These models can be accelerated notably by optimizing the most time-consuming components, such as chemistry, thereby reducing uncertainties and enhancing our understanding of climate dynamics. The primary objective of this thesis is to improve the performance of chemistry solvers in atmospheric models by developing efficient solutions tailored to GPU-based heterogeneous architectures. Specifically, the focus is maximizing GPU performance, enhancing accuracy, ensuring portability, and minimizing development effort by minor modifications to the existing chemistry solver. This thesis's critical contribution is developing a CPU-GPU version of the CAMP chemistry solver for the atmospheric model MONARCH, achieving significant acceleration. Furthermore, the developed chemical solver can be easily adapted to solve generic systems of Ordinary Differential Equations (ODEs) on the GPU due to its similarity to the CVODE library used in the CAMP CPU solver. The open-source code has been released as an upgrade to CAMP to facilitate portability. The main contributions include a method for seamlessly integrating GPU code into chemistry solvers (Multi-cells), distributing computational loads across GPU threads (Block-cells), and an automatic load-balancing algorithm for concurrent CPU-GPU execution. Notably, the automatic load-balancing algorithm is designed to be simple, ensuring portability while delivering accelerated execution. Furthermore, a GPU linear solver has been developed following the Block-cells approach, enabling further exploration of Block-cells in iterative solvers and facilitating coupling with other systems—an important consideration given the widespread use of linear solvers in scientific and mathematical communities. The results demonstrate a 250x speedup compared to the single-threaded CPU version when using four GPUs in a single node. In a node-to-node comparison, utilizing four GPUs versus 80 CPU cores, the speedup achieved is 8.14x. When normalizing the cores used against other state-of-the-art solutions, the speedup is 16x. The substantial speedup achieved with minimal algorithmic changes underscores the significance of the new strategies developed in this thesis.