This paper introduces a code designed for the three-dimensional relativistic hydrodynamic simulation of astrophysical flows on massive parallel architectures. The code utilizes numerical techniques based on Godunov’s method and the piecewise parabolic approximation with a local stencil to solve the equations of gravitational hydrodynamics. The implementation of the code leverages the hybrid use of MPI, OpenMP, and vectorization technologies. The collision of relativistic jets serves as a case study to evaluate the code. Experiments were conducted to assess its efficiency, performance, and scalability. © 2023, The Author(s), under exclusive license to Springer Nature Switzerland AG.