OMPEGAS: Optimized Relativistic Code for Multicore Architecture

Elena N. Akimova, Vladimir E. Misilov, Igor M. Kulikov, Igor G. Chernykh

Resultado de pesquisa: Articlerevisão de pares


The paper presents a new hydrodynamical code, OMPEGAS, for the 3D simulation of astrophysical flows on shared memory architectures. It provides a numerical method for solving the three-dimensional equations of the gravitational hydrodynamics based on Godunov's method for solving the Riemann problem and the piecewise parabolic approximation with a local stencil. It obtains a high order of accuracy and low dissipation of the solution. The code is implemented for multicore processors with vector instructions using the OpenMP technology, Intel SDLT library, and compiler auto-vectorization tools. The model problem of simulating a star explosion was used to study the developed code. The experiments show that the presented code reproduces the behavior of the explosion correctly. Experiments for the model problem with a grid size of 128 x 128 x 128 were performed on an 16-core Intel Core i9-12900K CPU to study the efficiency and performance of the developed code. By using the autovectorization, we achieved a 3.3-fold increase in speed in comparison with the non-vectorized program on the processor with AVX2 support. By using multithreading with OpenMP, we achieved an increase in speed of 2.6 times on a 16-core processor in comparison with the vectorized single-threaded program. The total increase in speed was up to ninefold.
Idioma originalEnglish
Número do artigo2546
Número de emissão14
Estado da publicaçãoPublished - 2022

ASJC Scopus subject areas

  • Computer Science (miscellaneous)
  • Engineering (miscellaneous)
  • Mathematics(all)

WoS ResearchAreas Categories

  • Mathematics

Impressão digital

Mergulhe nos tópicos de investigação de “OMPEGAS: Optimized Relativistic Code for Multicore Architecture“. Em conjunto formam uma impressão digital única.

Citar isto