--- /dev/null
+75c75
+< constexpr unsigned int n_global_refinements = 3;
+---
+> constexpr unsigned int n_global_refinements = 2;
+82c82
+< constexpr double final_time = testcase == 0 ? 10 : 2.0;
+---
+> constexpr double final_time = testcase == 0 ? 5 : 2.0;
+2196,2198d2195
+< const unsigned int n_vect_number = VectorizedArray<Number>::size();
+< const unsigned int n_vect_bits = 8 * sizeof(Number) * n_vect_number;
+<
+2202,2206d2198
+< pcout << "Vectorization over " << n_vect_number << ' '
+< << (std::is_same_v<Number, double> ? "doubles" : "floats") << " = "
+< << n_vect_bits << " bits ("
+< << Utilities::System::get_current_vectorization_level() << ')'
+< << std::endl;
+2284d2275
+< timer.print_wall_time_statistics(MPI_COMM_WORLD);
--- /dev/null
+Running with 1 MPI processes
+Number of degrees of freedom: 36864 ( = 4 [vars] x 256 [cells] x 36 [dofs/cell/var] )
+Time step size: 0.0137935, minimal h: 0.625, initial transport scaling: 0.205622
+
+Time: 0, dt: 0.014, error rho: 1.867e-05, rho * u: 7.97e-05, energy: 0.0001589
+Time: 1.01, dt: 0.014, error rho: 5.715e-05, rho * u: 0.0001503, energy: 0.0002395
+Time: 2, dt: 0.014, error rho: 7.119e-05, rho * u: 0.0001799, energy: 0.0002681
+Time: 3.01, dt: 0.014, error rho: 7.981e-05, rho * u: 0.0001989, energy: 0.0002824
+Time: 4, dt: 0.014, error rho: 9.041e-05, rho * u: 0.0002087, energy: 0.0002955
+Time: 5.01, dt: 0.014, error rho: 7.695e-05, rho * u: 0.000158, energy: 0.0002245
+
--- /dev/null
+67c67
+< constexpr unsigned int fe_degree = 5;
+---
+> constexpr unsigned int fe_degree = 3;
+89c89
+< constexpr double final_time = testcase == 0 ? 10 : 2.0;
+---
+> constexpr double final_time = testcase == 0 ? 10 : 0.5;
+1519,1521d1518
+< const unsigned int n_vect_number = VectorizedArrayType::size();
+< const unsigned int n_vect_bits = 8 * sizeof(Number) * n_vect_number;
+<
+1525,1529d1521
+< pcout << "Vectorization over " << n_vect_number << ' '
+< << (std::is_same_v<Number, double> ? "doubles" : "floats") << " = "
+< << n_vect_bits << " bits ("
+< << Utilities::System::get_current_vectorization_level() << ')'
+< << std::endl;
+1593d1584
+< timer.print_wall_time_statistics(MPI_COMM_WORLD);
--- /dev/null
+Running with 1 MPI processes
+Number of degrees of freedom: 102400 ( = 4 [vars] x 1600 [cells] x 16 [dofs/cell/var] )
+Time step size: 0.000636798, minimal h: 0.0075, initial transport scaling: 0.00441186
+
+Time: 0, dt: 0.00064, norm rho: 3.914e-16, rho * u: 1.431e-16, energy: 8.334e-16
+Time: 0.0503, dt: 0.00054, norm rho: 0.0207, rho * u: 0.03801, energy: 0.08745
+Time: 0.1, dt: 0.00051, norm rho: 0.02202, rho * u: 0.04509, energy: 0.08908
+Time: 0.15, dt: 0.00049, norm rho: 0.02256, rho * u: 0.04585, energy: 0.08944
+Time: 0.2, dt: 0.0005, norm rho: 0.02052, rho * u: 0.04354, energy: 0.08191
+Time: 0.25, dt: 0.00051, norm rho: 0.01688, rho * u: 0.04195, energy: 0.06839
+Time: 0.3, dt: 0.00051, norm rho: 0.01649, rho * u: 0.04002, energy: 0.0658
+Time: 0.35, dt: 0.00049, norm rho: 0.01768, rho * u: 0.04255, energy: 0.07025
+Time: 0.4, dt: 0.0005, norm rho: 0.01599, rho * u: 0.04491, energy: 0.06383
+Time: 0.45, dt: 0.0005, norm rho: 0.01523, rho * u: 0.04347, energy: 0.06164
+Time: 0.5, dt: 0.0005, norm rho: 0.01567, rho * u: 0.04187, energy: 0.06332
+