  We present an extension of the HLLC approximate Riemann solver by Toro,
Spruce and Speares to the relativistic equations of fluid dynamics. The solver
retains the simplicity of the original two-wave formulation proposed by Harten,
Lax and van Leer (HLL) but it restores the missing contact wave in the solution
of the Riemann problem. The resulting numerical scheme is computationally
efficient, robust and positively conservative. The performance of the new
solver is evaluated through numerical testing in one and two dimensions.
