A Simflowny-based finite-difference code for high-performance computing in numerical relativity

The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processors and advanced computational strategies are adopted. To achieve this goal we have developed Simflowny, an open platform for scientific dynamical models which automatically generates parallel code for different simulation frameworks, allowing the use of HPC infrastructures to nonspecialist scientists. One of these frameworks is SAMRAI, a mature patch-based structured adaptive mesh refinement infrastructure, capable of reaching exascale in some specific problems. Here we present the numerical techniques that we have implemented on this framework by using Simflowny in order to perform fast, efficient, accurate and highly-scalable simulations. These techniques involve high-order schemes for smooth and non-smooth solutions, Adaptive Mesh Refinement with arbitrary resolution ratios and an optimal strategy for the sub-cycling in time. We validate the automatically generated codes for the SAMRAI infrastructure with some simple test examples (i.e. wave equation and Newtonian MHD) and finally with the Einstein equations.