From a5ddabe6a99e96f5df73e109d0600034b8564f82 Mon Sep 17 00:00:00 2001 From: Daniel_Doehring Date: Wed, 11 Oct 2023 08:50:47 +0200 Subject: [PATCH] Add comments --- examples/tree_2d_dgsem/elixir_navierstokes_shearlayer_amr.jl | 1 + src/solvers/dgsem_tree/dg_2d_parabolic.jl | 2 ++ src/solvers/dgsem_tree/dg_3d_parabolic.jl | 2 ++ 3 files changed, 5 insertions(+) diff --git a/examples/tree_2d_dgsem/elixir_navierstokes_shearlayer_amr.jl b/examples/tree_2d_dgsem/elixir_navierstokes_shearlayer_amr.jl index 82c85aa16ba..2df05fbf210 100644 --- a/examples/tree_2d_dgsem/elixir_navierstokes_shearlayer_amr.jl +++ b/examples/tree_2d_dgsem/elixir_navierstokes_shearlayer_amr.jl @@ -65,6 +65,7 @@ analysis_callback = AnalysisCallback(semi, interval=analysis_interval) alive_callback = AliveCallback(analysis_interval=analysis_interval,) +# This uses velocity-based AMR amr_indicator = IndicatorLöhner(semi, variable=v1) amr_controller = ControllerThreeLevel(semi, amr_indicator, base_level = 3, diff --git a/src/solvers/dgsem_tree/dg_2d_parabolic.jl b/src/solvers/dgsem_tree/dg_2d_parabolic.jl index 220ad75fbe6..e1c2daf8f08 100644 --- a/src/solvers/dgsem_tree/dg_2d_parabolic.jl +++ b/src/solvers/dgsem_tree/dg_2d_parabolic.jl @@ -515,6 +515,8 @@ function calc_boundary_flux_by_direction_divergence!(surface_flux_values::Abstra return nothing end +# `cache` is the hyperbolic cache, i.e., in particular not `cache_parabolic`. +# This is because mortar handling is done in the (hyperbolic) `cache`. function prolong2mortars!(cache, flux_viscous::Vector{Array{uEltype, 4}}, mesh::TreeMesh{2}, equations_parabolic::AbstractEquationsParabolic, diff --git a/src/solvers/dgsem_tree/dg_3d_parabolic.jl b/src/solvers/dgsem_tree/dg_3d_parabolic.jl index 4bcbdb0b8c5..34c75070d82 100644 --- a/src/solvers/dgsem_tree/dg_3d_parabolic.jl +++ b/src/solvers/dgsem_tree/dg_3d_parabolic.jl @@ -598,6 +598,8 @@ function calc_boundary_flux_by_direction_divergence!(surface_flux_values::Abstra return nothing end +# `cache` is the hyperbolic cache, i.e., in particular not `cache_parabolic`. +# This is because mortar handling is done in the (hyperbolic) `cache`. function prolong2mortars!(cache, flux_viscous::Vector{Array{uEltype, 5}}, mesh::TreeMesh{3},