Following https://iopscience.iop.org/article/10.1088/0004-637X/806/2/203

At the same time, we should also follow the Riemann solvers / diffusion coefficients and make the stopping times template specialized classes