fix bug again

trixi-framework · Dec 5, 2023 · 4ce423c · 4ce423c
1 parent c28ad27
commit 4ce423c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/src/schemes/boundary/monaghan_kajtar/monaghan_kajtar.jl b/src/schemes/boundary/monaghan_kajtar/monaghan_kajtar.jl
@@ -77,9 +77,9 @@ function Base.show(io::IO, model::BoundaryModelMonaghanKajtar)
 end
 
 @inline function pressure_acceleration(pressure_correction, m_b, p_a, p_b,
-                                       boundary_model::BoundaryModelMonaghanKajtar,
                                        rho_a, rho_b, pos_diff::SVector{NDIMS},
                                        smoothing_length, grad_kernel,
+                                       boundary_model::BoundaryModelMonaghanKajtar,
                                        density_calculator) where {NDIMS}
     (; K, beta, boundary_particle_spacing) = boundary_model
 

diff --git a/src/schemes/solid/total_lagrangian_sph/rhs.jl b/src/schemes/solid/total_lagrangian_sph/rhs.jl
@@ -104,10 +104,9 @@ function interact!(dv, v_particle_system, u_particle_system,
         # Boundary forces
         # Note: neighbor and particle pressure are switched in this call
         #       and `pressure_correction` is set to `1.0` (no correction)
-        dv_boundary = pressure_acceleration(1.0, m_b, p_b, p_a, boundary_model,
-                                            rho_b, rho_a, pos_diff,
+        dv_boundary = pressure_acceleration(1.0, m_b, p_b, p_a, rho_b, rho_a, pos_diff,
                                             neighbor_system.smoothing_length,
-                                            grad_kernel, density_calculator)
+                                            grad_kernel, boundary_model, density_calculator)
         dv_particle = dv_boundary + dv_viscosity
 
         for i in 1:ndims(particle_system)