diff --git a/src/solvers/dgsem_tree/indicators.jl b/src/solvers/dgsem_tree/indicators.jl index 04b3bc98aec..61ef682e8e9 100644 --- a/src/solvers/dgsem_tree/indicators.jl +++ b/src/solvers/dgsem_tree/indicators.jl @@ -116,7 +116,7 @@ function (indicator_hg::IndicatorHennemannGassner)(u, mesh, equations, dg::DGSEM # magic parameters # TODO: Are there better values for Float32? RealT = real(dg) - threshold = 0.5f0 * 10^(convert(RealT, -1.8) * nnodes(dg)^convert(RealT, 0.25)) + threshold = 0.5f0 * 10^(convert(RealT, -1.8) * nnodes(dg)^0.25f0) o_0001 = convert(RealT, 0.0001) parameter_s = log((1 - o_0001) / o_0001)