It could be useful to have some testing workflow for the Riemann solvers. I was thinking about writing some Pytest-based regression tests for the Python and f2py-compiled Fortran solvers. The idea would be to compare, for a set of problems, the output of each solver (waves, speeds, normal and transverse fluctuations) with some reference values stored in dictionaries, instead of focusing on the Riemann solution itself. Would something like this make sense?