i've lost track of all the changes.  they include a number of modifications to EulerEquationOfState to try to discern where the numerical inaccuracy comes from when I analytically provide the flux values rather than recomputing them from the eigenbasis, and from when I analytically provide the eigenvector inverses rather than numerically computing them from the eigenvectors.