Porting OpenMP to GPGPU to accelerate Krylov space computations in coupled geodynamics simulator

Understanding the causality between the events leading upto and post fault slipand the earthquake recording is important for seismic design and monitoring ofunderground structures, bridges and reinforced concrete buildings as well as climatemitigation projects like carbon sequestration and energy technologies like enhancedgeothermal systems or oilfield wastewater disposal. While the events leading uptofault slip are typically governed by poroelastostatics, the events post fault slip caneasily transition into poroelastodynamics territory due to runaway fault slip velocities.An understanding of expected fault slip velocities is critical apriori, as an algorithmwhich can seamlessly transition from time marching in poroelastostatics realm toporoelastodynamics realm and vice-versa is extremely difficult to achieve. That beingsaid, every effort in the direction of accelerating the computations on the flow sideare a necessary step forward in rendering a fast coupled geodynamics simulator.In this document, we present a framework in which we study the porting of theOpenMP parallelism of the flow simulator to a GPGPU