  Multidimensional shock-capturing numerical schemes for special relativistic
hydrodynamics (RHD) are computationally more expensive than their correspondent
Euler versions, due to the nonlinear relations between conservative and
primitive variables and to the consequent complexity of the Jacobian matrices
(needed for the spectral decomposition in most of the approximate Riemann
solvers of common use). Here an efficient and easy-to-implement
three-dimensional (3-D) shock-capturing scheme for ideal RHD is presented.
Based on the algorithms developed by P. Londrillo and L. Del Zanna ({\em
Astrophys. J.} 530, 508-524, 2000) for the non-relativistic magnetohydrodynamic
(MHD) case, and having in mind its relativistic MHD extension (to appear in a
forthcoming paper), the scheme uses high order (third) Convex Essentially
Non-Oscillatory (CENO) finite difference interpolation routines and
central-type averaged Riemann solvers, which do not make use of time-consuming
characteristic decomposition. The scheme is very efficient and robust, and it
gives results comparable to those obtained with more sophisticated algorithms,
even in ultrarelativistic multidimensional test problems.
