Added Roes Riemann solver

We only ever seem to use `nvr - nvl`; compute this directly in lieu of `nvl` and `nvr` in a single `${pyfr.dot(...})` call.