diff --git a/examples/tree_1d_dgsem/elixir_advection_PERK2.jl b/examples/tree_1d_dgsem/elixir_advection_PERK2.jl
index f8c20adbdb9..86936209360 100644
--- a/examples/tree_1d_dgsem/elixir_advection_PERK2.jl
+++ b/examples/tree_1d_dgsem/elixir_advection_PERK2.jl
@@ -60,5 +60,3 @@ sol = Trixi.solve(ode, ode_algorithm,
 # Print the timer summary
 summary_callback()
 
-# Show the error of the elixir
-analysis_callback(sol)
diff --git a/src/time_integration/paired_explicit_runge_kutta/methods_PERK2.jl b/src/time_integration/paired_explicit_runge_kutta/methods_PERK2.jl
index 6dd9957a9b2..d4a66af7fdd 100644
--- a/src/time_integration/paired_explicit_runge_kutta/methods_PERK2.jl
+++ b/src/time_integration/paired_explicit_runge_kutta/methods_PERK2.jl
@@ -280,7 +280,7 @@ function solve!(integrator::PERK2Integrator)
             end
 
             # Higher stages
-            for stage in 3:(alg.num_stages)
+            for stage in 3:alg.num_stages
                 # Construct current state
                 @threaded for i in eachindex(integrator.du)
                     integrator.u_tmp[i] = integrator.u[i] +
diff --git a/test/test_tree_1d_advection.jl b/test/test_tree_1d_advection.jl
index d00bc82deeb..39ead7707e3 100644
--- a/test/test_tree_1d_advection.jl
+++ b/test/test_tree_1d_advection.jl
@@ -84,8 +84,8 @@ end
 
 @trixi_testset "elixir_advection_PERK2.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_PERK2.jl"),
-                        l2=[0.011662300515980219],
-                        linf=[0.01647256923710194])
+                        l2=[0.0006911538643023453],
+                        linf=[0.0009859411767757509])
     # Ensure that we do not have excessive memory allocations
     # (e.g., from type instabilities)
     let