diff --git a/src/equations/compressible_euler_multicomponent_1d.jl b/src/equations/compressible_euler_multicomponent_1d.jl index da434579f7..bcfec89e5c 100644 --- a/src/equations/compressible_euler_multicomponent_1d.jl +++ b/src/equations/compressible_euler_multicomponent_1d.jl @@ -202,16 +202,15 @@ function initial_condition_weak_blast_wave(x, t, r = abs(x_norm) cos_phi = x_norm > 0 ? 1 : -1 - prim_rho = SVector{ncomponents(equations), real(equations)}(r > 0.5f0 ? - 2^(i - 1) * (1 - 2) / - (1 - - 2^ncomponents(equations)) * - one(RealT) : - 2^(i - 1) * (1 - 2) / - (1 - - 2^ncomponents(equations)) * - convert(RealT, 1.1691) - for i in eachcomponent(equations)) + prim_rho = SVector{ncomponents(equations), RealT}(r > 0.5f0 ? + 2^(i - 1) * (1 - 2) / + (1 - + 2^ncomponents(equations)) : + 2^(i - 1) * (1 - 2) / + (1 - + 2^ncomponents(equations)) * + 1.1691 + for i in eachcomponent(equations)) v1 = r > 0.5f0 ? zero(RealT) : convert(RealT, 0.1882) * cos_phi p = r > 0.5f0 ? one(RealT) : convert(RealT, 1.245) diff --git a/src/equations/compressible_euler_multicomponent_2d.jl b/src/equations/compressible_euler_multicomponent_2d.jl index 2424ad0301..050f9dc899 100644 --- a/src/equations/compressible_euler_multicomponent_2d.jl +++ b/src/equations/compressible_euler_multicomponent_2d.jl @@ -220,16 +220,15 @@ function initial_condition_weak_blast_wave(x, t, phi = atan(y_norm, x_norm) sin_phi, cos_phi = sincos(phi) - prim_rho = SVector{ncomponents(equations), real(equations)}(r > 0.5f0 ? - 2^(i - 1) * (1 - 2) / - (1 - - 2^ncomponents(equations)) * - one(RealT) : - 2^(i - 1) * (1 - 2) / - (1 - - 2^ncomponents(equations)) * - convert(RealT, 1.1691) - for i in eachcomponent(equations)) + prim_rho = SVector{ncomponents(equations), RealT}(r > 0.5f0 ? + 2^(i - 1) * (1 - 2) / + (1 - + 2^ncomponents(equations)) : + 2^(i - 1) * (1 - 2) / + (1 - + 2^ncomponents(equations)) * + 1.1691 + for i in eachcomponent(equations)) v1 = r > 0.5f0 ? zero(RealT) : convert(RealT, 0.1882) * cos_phi v2 = r > 0.5f0 ? zero(RealT) : convert(RealT, 0.1882) * sin_phi