If matter sources are included, like fluids, then one might need
to consider flux conservative schema for the whole system of
equations. But probably it would be much better to use
relativistic dissipative fluids -whose global existence for small
data has been proven recently [30]-to dispense of shocks altogether when considering weak
data.