Since the appearance of computers, numerical methods for discontinuous solutions of quasi-linear hyperbolic systems of partial differential equations have been among the most important research subjects in numerical analysis.