“…In recent years there have been many works on hyperbolic conservation laws with a spatially discontinuous flux function, providing a great number of results relating to existence, uniqueness, stability, and numerical approximations of entropy solutions [1,2,3,4,5,3,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,34,36,37,38,39,41,42,43,44,45,46,47]. Herein we are interested in numerical methods for the initial value problem u t + F(x, u) x = 0 for (x, t) ∈ Π T := R × (0, T ), u(x, 0) = u 0 (x) for x ∈ R, (1.1) F(x, u) := H(−x)g(u) + H(x)f (u), where H(x) is the Heaviside function.…”