  We present a novel methodology based on Physics-Informed Neural Networks
(PINNs) for solving systems of partial differential equations admitting
discontinuous solutions. Our method, called Gradient-Annihilated PINNs
(GA-PINNs), introduces a modified loss function that requires the model to
partially ignore high-gradients in the physical variables, achieved by
introducing a suitable weighting function. The method relies on a set of
hyperparameters that control how gradients are treated in the physical loss and
how the activation functions of the neural model are dynamically accounted for.
The performance of our GA-PINN model is demonstrated by solving Riemann
problems in special relativistic hydrodynamics, extending earlier studies with
PINNs in the context of the classical Euler equations. The solutions obtained
with our GA-PINN model correctly describe the propagation speeds of
discontinuities and sharply capture the associated jumps. We use the relative
$l^{2}$ error to compare our results with the exact solution of special
relativistic Riemann problems, used as the reference ``ground truth'', and with
the error obtained with a second-order, central, shock-capturing scheme. In all
problems investigated, the accuracy reached by our GA-PINN model is comparable
to that obtained with a shock-capturing scheme and significantly higher than
that achieved by a baseline PINN algorithm. An additional benefit worth
stressing is that our PINN-based approach sidesteps the costly recovery of the
primitive variables from the state vector of conserved ones, a well-known
drawback of grid-based solutions of the relativistic hydrodynamics equations.
Due to its inherent generality and its ability to handle steep gradients, the
GA-PINN method discussed could be a valuable tool to model relativistic flows
in astrophysics and particle physics, characterized by the prevalence of
discontinuous solutions.
