Currently, mesh adaptation is performed within ODESolverBase::steady_state(). To make the functions of ODESolver and MeshAdaptation separate, it would be better to move mesh adaptation to flow solver in future and call MeshAdaptation::adapt_mesh() after ODESolverBase::steady_state(). Some tests (such as the S-shock) do not use flow solver yet and these need to be modified to use the flow solver framework.
Currently, mesh adaptation is performed within ODESolverBase::steady_state(). To make the functions of ODESolver and MeshAdaptation separate, it would be better to move mesh adaptation to flow solver in future and call MeshAdaptation::adapt_mesh() after ODESolverBase::steady_state(). Some tests (such as the S-shock) do not use flow solver yet and these need to be modified to use the flow solver framework.