diff --git a/src/solvers/dgsem_tree/dg_1d.jl b/src/solvers/dgsem_tree/dg_1d.jl
index a06498abe1c..58664b2ea98 100644
--- a/src/solvers/dgsem_tree/dg_1d.jl
+++ b/src/solvers/dgsem_tree/dg_1d.jl
@@ -398,7 +398,7 @@ end
 
 function prolong2interfaces!(cache, u_or_flux_viscous,
                              mesh::TreeMesh{1}, equations,
-                             dg::DG{<:GaussLegendreBasis})
+                             dg::DGSEM{<:GaussLegendreBasis})
     @unpack interfaces = cache
     @unpack neighbor_ids = interfaces
     @unpack boundary_interpolation = dg.basis
@@ -535,7 +535,7 @@ end
 
 function prolong2boundaries!(cache, u_or_flux_viscous,
                              mesh::TreeMesh{1}, equations,
-                             dg::DG{<:GaussLegendreBasis})
+                             dg::DGSEM{<:GaussLegendreBasis})
     @unpack boundaries = cache
     @unpack neighbor_sides = boundaries
     @unpack boundary_interpolation = dg.basis
@@ -698,7 +698,7 @@ end
 
 function calc_surface_integral!(du, u, mesh::Union{TreeMesh{1}, StructuredMesh{1}},
                                 equations, surface_integral::SurfaceIntegralWeakForm,
-                                dg::DG{<:GaussLegendreBasis}, cache)
+                                dg::DGSEM{<:GaussLegendreBasis}, cache)
     @unpack boundary_interpolation_inverse_weights = dg.basis
     @unpack surface_flux_values = cache.elements
 
diff --git a/src/solvers/dgsem_tree/dg_2d.jl b/src/solvers/dgsem_tree/dg_2d.jl
index e8b464649c2..a944fdca4b6 100644
--- a/src/solvers/dgsem_tree/dg_2d.jl
+++ b/src/solvers/dgsem_tree/dg_2d.jl
@@ -524,7 +524,7 @@ function prolong2interfaces!(cache, u, mesh::TreeMesh{2}, equations, dg::DG)
 end
 
 function prolong2interfaces!(cache, u, mesh::TreeMesh{2}, equations,
-                             dg::DG{<:GaussLegendreBasis})
+                             dg::DGSEM{<:GaussLegendreBasis})
     @unpack interfaces = cache
     @unpack orientations, neighbor_ids = interfaces
     @unpack boundary_interpolation = dg.basis
@@ -701,7 +701,7 @@ end
 
 function prolong2boundaries!(cache, u,
                              mesh::TreeMesh{2}, equations,
-                             dg::DG{<:GaussLegendreBasis})
+                             dg::DGSEM{<:GaussLegendreBasis})
     @unpack boundaries = cache
     @unpack orientations, neighbor_sides = boundaries
     @unpack boundary_interpolation = dg.basis
@@ -1252,7 +1252,7 @@ function calc_surface_integral!(du, u,
                                 mesh::Union{TreeMesh{2}, StructuredMesh{2},
                                             StructuredMeshView{2}},
                                 equations, surface_integral::SurfaceIntegralWeakForm,
-                                dg::DG{<:GaussLegendreBasis}, cache)
+                                dg::DGSEM{<:GaussLegendreBasis}, cache)
     @unpack boundary_interpolation_inverse_weights = dg.basis
     @unpack surface_flux_values = cache.elements
 
diff --git a/test/test_parabolic_1d.jl b/test/test_parabolic_1d.jl
index 75af0c07c60..a5ee4b6ab25 100644
--- a/test/test_parabolic_1d.jl
+++ b/test/test_parabolic_1d.jl
@@ -289,6 +289,28 @@ end
     @test_allocations(Trixi.rhs_parabolic!, semi, sol, 1000)
 end
 
+@trixi_testset "TreeMesh1D: elixir_navierstokes_convergence_walls.jl (Gauss-Legendre)" begin
+    @test_trixi_include(joinpath(EXAMPLES_DIR, "tree_1d_dgsem",
+                                 "elixir_navierstokes_convergence_walls.jl"),
+                        solver=DGSEM(polydeg = 3, surface_flux = flux_hll,
+                                     basis_type = GaussLegendreBasis),
+                        l2=[
+                            4.201407214726233e-5,
+                            9.758795848409044e-5,
+                            0.00042000319972803127
+                        ],
+                        linf=[
+                            0.00015354693449309842,
+                            0.00041974539800903807,
+                            0.0016951216452660844
+                        ],
+                        atol=1e-10)
+    # Ensure that we do not have excessive memory allocations
+    # (e.g., from type instabilities)
+    @test_allocations(Trixi.rhs!, semi, sol, 1000)
+    @test_allocations(Trixi.rhs_parabolic!, semi, sol, 1000)
+end
+
 @trixi_testset "TreeMesh1D: elixir_navierstokes_convergence_walls_amr.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "tree_1d_dgsem",
                                  "elixir_navierstokes_convergence_walls_amr.jl"),
@@ -355,7 +377,7 @@ end
     @test_allocations(Trixi.rhs_parabolic!, semi, sol, 1000)
 end
 
-@trixi_testset "TreeMesh1D: elixir_navierstokes_viscous_shock.jl" begin
+@trixi_testset "TreeMesh1D: elixir_navierstokes_viscous_shock.jl (Gauss-Legendre)" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "tree_1d_dgsem",
                                  "elixir_navierstokes_viscous_shock.jl"),
                         solver=DGSEM(polydeg = 3, surface_flux = flux_hlle,