  We present the results of a numerical benchmark study for the
MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two-
and three-dimensional compressible hydrodynamics problems relevant to stellar
interiors. MUSIC is an implicit large eddy simulation code that uses implicit
time integration, implemented as a Jacobian-free Newton Krylov method. A
physics based preconditioning technique which can be adjusted to target varying
physics is used to improve the performance of the solver. The problems used for
this benchmark study include the Rayleigh-Taylor and Kelvin-Helmholtz
instabilities, and the decay of the Taylor-Green vortex. Additionally we show a
test of hydrostatic equilibrium, in a stellar environment which is dominated by
radiative effects. In this setting the flexibility of the preconditioning
technique is demonstrated. This work aims to bridge the gap between the
hydrodynamic test problems typically used during development of numerical
methods and the complex flows of stellar interiors. A series of
multi-dimensional tests are performed and analysed. Each of these test cases is
analysed with a simple, scalar diagnostic, with the aim of enabling direct code
comparisons. As the tests performed do not have analytic solutions we verify
MUSIC by comparing to established codes including ATHENA and the PENCIL code.
MUSIC is able to both reproduce behaviour from established and widely-used
codes as well as results expected from theoretical predictions. This
benchmarking study concludes a series of papers describing the development of
the MUSIC code and provides confidence in the future applications.
