From 9ef606687986fc039a33bdf8dd25dadfd67b4b45 Mon Sep 17 00:00:00 2001 From: Charles Kawczynski Date: Sat, 26 Oct 2024 01:09:34 -0400 Subject: [PATCH] Add support for make_q_bounds with new datalayouts --- src/Limiters/quasimonotone.jl | 14 ++++++++++++++ 1 file changed, 14 insertions(+) diff --git a/src/Limiters/quasimonotone.jl b/src/Limiters/quasimonotone.jl index 9f41645948..b038cd91de 100644 --- a/src/Limiters/quasimonotone.jl +++ b/src/Limiters/quasimonotone.jl @@ -74,6 +74,13 @@ function make_q_bounds( _, _, _, _, Nh = size(ρq) return DataLayouts.IFH{S, 2}(similar(parent(ρq), (2, Nf, Nh))) end +function make_q_bounds( + ρq::Union{DataLayouts.IHF{S}, DataLayouts.IJHF{S}}, +) where {S} + Nf = DataLayouts.ncomponents(ρq) + _, _, _, _, Nh = size(ρq) + return DataLayouts.IHF{S, 2}(similar(parent(ρq), (2, Nh, Nf))) +end function make_q_bounds( ρq::Union{DataLayouts.VIFH{S}, DataLayouts.VIJFH{S}}, ) where {S} @@ -81,6 +88,13 @@ function make_q_bounds( _, _, _, Nv, Nh = size(ρq) return DataLayouts.VIFH{S, Nv, 2}(similar(parent(ρq), (Nv, 2, Nf, Nh))) end +function make_q_bounds( + ρq::Union{DataLayouts.VIHF{S}, DataLayouts.VIJHF{S}}, +) where {S} + Nf = DataLayouts.ncomponents(ρq) + _, _, _, Nv, Nh = size(ρq) + return DataLayouts.VIHF{S, Nv, 2}(similar(parent(ρq), (Nv, 2, Nh, Nf))) +end """