Modeling power consumption of 3D MPDATA and the CG method on ARM and Intel multicore architectures

We propose an approach to estimate the power consumption of algorithms, as a function of the frequency and number of cores, using only a very reduced set of real power measures. In addition, we also provide the formulation of a method to select the voltage–frequency scaling–concurrency throttling configurations that should be tested in order to obtain accurate estimations of the power dissipation. The power models and selection methodology are verified using two real scientific application: the stencil-based 3D MPDATA algorithm and the conjugate gradient (CG) method for sparse linear systems. MPDATA is a crucial component of the EULAG model, which is widely used in weather forecast simulations. The CG algorithm is the keystone for iterative solution of sparse symmetric positive definite linear systems via Krylov subspace methods. The reliability of the method is confirmed for a variety of ARM and Intel architectures, where the estimated results correspond to the real measured values with the average error being slightly below 5% in all cases

Performance and Energy Optimization of the Iterative Solution of Sparse Linear Systems on Multicore Processors

En esta tesis doctoral se aborda la solución de sistemas dispersos de ecuaciones lineales utilizando métodos iterativos precondicionados basados en subespacios de Krylov. En concreto, se centra en ILUPACK, una biblioteca que implementa precondicionadores de tipo ILU multinivel para la solución eficiente de sistemas lineales dispersos. El incremento en el número de ecuaciones, y la aparición de nuevas arquitecturas, motiva el desarrollo de una versión paralela de ILUPACK que optimice tanto el tiempo de ejecución como el consumo energético en arquitecturas multinúcleo actuales y en clusters de nodos construidos con esta tecnología. El objetivo principal de la tesis es el diseño, implementación y valuación de resolutores paralelos energéticamente eficientes para sistemas lineales dispersos orientados a procesadores multinúcleo así como aceleradores hardware como el Intel Xeon Phi. Para lograr este objetivo, se aprovecha el paralelismo de tareas mediante OmpSs y MPI, y se desarrolla un entorno automático para detectar ineficiencias energéticas.In this dissertation we target the solution of large sparse systems of linear equations using preconditioned iterative methods based on Krylov subspaces. Specifically, we focus on ILUPACK, a library that offers multi-level ILU preconditioners for the effective solution of sparse linear systems. The increase of the number of equations and the introduction of new HPC architectures motivates us to develop a parallel version of ILUPACK which optimizes both execution time and energy consumption on current multicore architectures and clusters of nodes built from this type of technology. Thus, the main goal of this thesis is the design, implementation and evaluation of parallel and energy-efficient iterative sparse linear system solvers for multicore processors as well as recent manycore accelerators such as the Intel Xeon Phi. To fulfill the general objective, we optimize ILUPACK exploiting task parallelism via OmpSs and MPI, and also develope an automatic framework to detect energy inefficiencies

Salinity Transport in a Finite-Volume Sigma-Layer Three-Dimensional Model

The objective of this study was to develop a 3-D model for The Pontchartrain Estuary that was capable of long-term mass conservative simulation of salinities. This was accomplished in a multi-stage approach involving: a physical model of salinity exchange through a pass; a 3-D FVCOM model of the physical experiment; the development and testing of an FVCOM model for an idealized Pontchartrain Basin; and for the entire estuary. The data from the physical model tests were used to validate the performance of the FVCOM model with density-driven flows. These results showed that hydrostatic FVCOM captured the primary internal wave movement. The idealized basin simulations were used to evaluate several issues related to salinity transport, namely the relative importance of baroclinic forcing, tidal forcing and hydrology. The idealized domain also permitted the testing of sigma-gradients, spatial distribution of friction coefficients, wind stress and various boundary treatments. The results showed that the density-driven exchange of saltwater at the open boundary required a baroclinic boundary condition for salinity as well as a lateral filter at the boundary on each sigma layer. A new radiative baroclinic open boundary condition was developed for FVCOM. When tides and hydrology were included, the FVCOM model was shown to reproduce the seasonal salinity that has been observed for long-term periods. It was also found that the simulation of tides and salinity in FVCOM is very sensitive to the spatial distribution of the friction coefficient; relatively low friction was required in the open water regions and high friction was needed in the passes and waterways to reproduce the tides and salinity distribution. A variable friction coefficient option was coded on FVCOM. The findings from the idealized model were utilized to setup two models for the actual estuary. Both models extend from Lake Maurepas, one to the Chandeleurs Islands and the other to Mobile Bay. The baroclinic open boundary and variable friction were implemented in these models. They were calibrated for tides and salinity. The 2008 Bonnet Carré Spillway Opening was applied to the first model. A tidal pumping effect in Lake Pontchartrain was observed and captured by the model

Salinity Transport in a Finite-Volume Sigma-Layer Three-Dimensional Model

The objective of this study was to develop a 3-D model for The Pontchartrain Estuary that was capable of long-term mass conservative simulation of salinities. This was accomplished in a multi-stage approach involving: a physical model of salinity exchange through a pass; a 3-D FVCOM model of the physical experiment; the development and testing of an FVCOM model for an idealized Pontchartrain Basin; and for the entire estuary. The data from the physical model tests were used to validate the performance of the FVCOM model with density-driven flows. These results showed that hydrostatic FVCOM captured the primary internal wave movement. The idealized basin simulations were used to evaluate several issues related to salinity transport, namely the relative importance of baroclinic forcing, tidal forcing and hydrology. The idealized domain also permitted the testing of sigma-gradients, spatial distribution of friction coefficients, wind stress and various boundary treatments. The results showed that the density-driven exchange of saltwater at the open boundary required a baroclinic boundary condition for salinity as well as a lateral filter at the boundary on each sigma layer. A new radiative baroclinic open boundary condition was developed for FVCOM. When tides and hydrology were included, the FVCOM model was shown to reproduce the seasonal salinity that has been observed for long-term periods. It was also found that the simulation of tides and salinity in FVCOM is very sensitive to the spatial distribution of the friction coefficient; relatively low friction was required in the open water regions and high friction was needed in the passes and waterways to reproduce the tides and salinity distribution. A variable friction coefficient option was coded on FVCOM. The findings from the idealized model were utilized to setup two models for the actual estuary. Both models extend from Lake Maurepas, one to the Chandeleurs Islands and the other to Mobile Bay. The baroclinic open boundary and variable friction were implemented in these models. They were calibrated for tides and salinity. The 2008 Bonnet Carré Spillway Opening was applied to the first model. A tidal pumping effect in Lake Pontchartrain was observed and captured by the model