  In this article, we present entropy stable discontinuous Galerkin numerical
schemes for equations of special relativistic hydrodynamics with the ideal
equation of state. The numerical schemes use the summation by parts (SBP)
property of Gauss-Lobatto quadrature rules. To achieve entropy stability for
the scheme, we use two-point entropy conservative numerical flux inside the
cells and a suitable entropy stable numerical flux at the cell interfaces. The
resulting semi-discrete scheme is then shown to entropy stable. Time
discretization is performed using SSP Runge-Kutta methods. Several numerical
test cases are presented to validate the accuracy and stability of the proposed
schemes.
