Test SRMHD wave speeds and quartic solver, add shocktube.c