Comparison of various numerical choices implemented in diffhydro in terms of fractional error in density through the center of a Sedov blast wave. We choose for our fiducial model RK2 integrator, HLLC Riemann solver, MUSCL reconstruction with a Vanleer flux limiter. We also show a reference ATHENA++ run with the same initial conditions for comparison. Top: choice of integrator used. Upper middle: choice of flux-reconstruction method used. Lower middle: choice of Reimann solver used. Bottom: flux limiter used.