Add HLLC riemann solver for nrmhd, move source to riemann.c