Multithreaded simulations freeze sometimes

I've had multithreaded simulations freeze multiple times lately. Interrupting the simulation works immediately, though.
The interruption error always looks like this:
```
ERROR: LoadError: InterruptException:
Stacktrace:
[1] try_yieldto(undo::typeof(Base.ensure_rescheduled))
@ Base .\task.jl:710
[2] wait
@ .\task.jl:769 [inlined]
[3] yield()
@ Base .\task.jl:662
[4] wait
@ C:\Users\Erik\.julia\packages\ThreadingUtilities\pkz6e\src\threadtasks.jl:62 [inlined]
[5] wait
@ C:\Users\Erik\.julia\packages\ThreadingUtilities\pkz6e\src\threadtasks.jl:57 [inlined]
[6] macro expansion
@ C:\Users\Erik\.julia\packages\Polyester\0DPCU\src\batch.jl:89 [inlined]
[7] _batch_no_reserve
@ C:\Users\Erik\.julia\packages\Polyester\0DPCU\src\batch.jl:53 [inlined]
[8] batch
@ C:\Users\Erik\.julia\packages\Polyester\0DPCU\src\batch.jl:195 [inlined]
[9] macro expansion
@ C:\Users\Erik\.julia\packages\Polyester\0DPCU\src\closure.jl:164 [inlined]
[10] macro expansion
@ c:\Users\Erik\git\Trixi.jl\src\auxiliary\auxiliary.jl:181 [inlined]
[11] calc_volume_integral!(du::StrideArraysCore.PtrArray{Tuple{Static.StaticInt{5}, Static.StaticInt{4}, Static.StaticInt{4}, Static.StaticInt{4}, Int64}, (true, true, true, true, true), Float64, 5, 1, 0, (1, 2, 3, 4, 5), Tuple{Static.StaticInt{8}, Static.StaticInt{40}, Static.StaticInt{160}, Static.StaticInt{640}, Static.StaticInt{2560}}, NTuple{5, Static.StaticInt{1}}}, u::StrideArraysCore.PtrArray{Tuple{Static.StaticInt{5}, Static.StaticInt{4}, Static.StaticInt{4}, Static.StaticInt{4}, Int64}, (true, true, true, true, true), Float64, 5, 1, 0, (1, 2, 3, 4, 5), Tuple{Static.StaticInt{8}, Static.StaticInt{40}, Static.StaticInt{160}, Static.StaticInt{640}, Static.StaticInt{2560}}, NTuple{5, Static.StaticInt{1}}}, mesh::P4estMesh{3, Float64, Ptr{P4est.LibP4est.p8est}, 5, 4}, nonconservative_terms::Val{false}, equations::CompressibleEulerEquations3D{Float64}, volume_integral::VolumeIntegralShockCapturingHG{typeof(flux_chandrashekar), FluxLaxFriedrichs{typeof(max_abs_speed_naive)}, IndicatorHennemannGassner{Float64, typeof(density_pressure), NamedTuple{(:alpha, :alpha_tmp, :indicator_threaded, :modal_threaded, :modal_tmp1_threaded, :modal_tmp2_threaded), Tuple{Vector{Float64}, Vector{Float64}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}}}}}, dg::DGSEM{LobattoLegendreBasis{Float64, 4, SVector{4, Float64},
Matrix{Float64}, Matrix{Float64}, Matrix{Float64}}, Trixi.LobattoLegendreMortarL2{Float64, 4, Matrix{Float64}, Matrix{Float64}}, SurfaceIntegralWeakForm{FluxLaxFriedrichs{typeof(max_abs_speed_naive)}}, VolumeIntegralShockCapturingHG{typeof(flux_chandrashekar), FluxLaxFriedrichs{typeof(max_abs_speed_naive)}, IndicatorHennemannGassner{Float64, typeof(density_pressure), NamedTuple{(:alpha, :alpha_tmp, :indicator_threaded, :modal_threaded, :modal_tmp1_threaded, :modal_tmp2_threaded), Tuple{Vector{Float64}, Vector{Float64}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}}}}}}, cache::NamedTuple{(:elements, :interfaces, :boundaries, :mortars, :element_ids_dg, :element_ids_dgfv, :fstar1_threaded, :fstar2_threaded, :fstar3_threaded, :fstar_threaded, :fstar_tmp_threaded, :u_threaded), Tuple{Trixi.P4estElementContainer{3, Float64, Float64, 4, 5, 6}, Trixi.P4estInterfaceContainer{3, Float64, 5}, Trixi.P4estBoundaryContainer{3, Float64, 4}, Trixi.P4estMortarContainer{3, Float64, 4, 6}, Vector{Int64}, Vector{Int64}, Vector{Array{Float64, 4}}, Vector{Array{Float64, 4}}, Vector{Array{Float64, 4}}, Vector{Array{Float64, 4}}, Vector{Array{Float64, 3}}, Vector{Array{Float64, 3}}}})
@ Trixi c:\Users\Erik\git\Trixi.jl\src\solvers\dgsem_tree\dg_3d.jl:434
```
It's always the `@threaded` block of the used volume integral.

Has this happened to anyone else before? Is there a way to debug stuff like this?