diff --git a/src/HRSC/dragger.jl b/src/HRSC/dragger.jl new file mode 100644 index 0000000000000000000000000000000000000000..9deca293327223562d95b6ded02f34db78b10e36 --- /dev/null +++ b/src/HRSC/dragger.jl @@ -0,0 +1,19 @@ +####################################################################### +# Type # +####################################################################### + + +struct DraggedArtificialViscosity{T_AV<:AbstractArtificialViscosity} <: AbstractArtificialViscosity + av::T_AV + drag::Float64 +end + + +####################################################################### +# Methods # +####################################################################### + + +function drag(mu, mu_m1, av::DraggedArtificialViscosity) + @. mu = (1-drag)*mu + drag*mu_m1 +end diff --git a/src/HRSC/weno.jl b/src/HRSC/weno.jl index 000712f9648d30f7220abcb12f2c2a1362431e90..9dbf38a08d21d3c4c9c8990c0ba8eeba2df9fb31 100644 --- a/src/HRSC/weno.jl +++ b/src/HRSC/weno.jl @@ -220,7 +220,7 @@ function reconstruct!(u, flags, weno::WENOReconstruction; isperiodic::Bool=false # normalize weights, (3.71) weights ./= sum(weights, dims=1) - + # reconstruct solution, (3.72) fill!(tmp_uweno, 0.0) for idx = 1:w+1