Propagate inbounds into viscous_velocity functions

Author: efaulhaber

src/schemes/boundary/wall_boundary/system.jl

@@ -145,7 +145,8 @@ end
     return zero(SVector{ndims(system), eltype(system)})
 end
 
-@inline function viscous_velocity(v, system::WallBoundarySystem, particle, v_particle)
+@propagate_inbounds function viscous_velocity(v, system::WallBoundarySystem,
+                                              particle, v_particle)
     return viscous_velocity(v, system.boundary_model.viscosity, system,
                             particle, v_particle)
 end
@@ -155,7 +156,7 @@ end
     return v_particle
 end
 
-@inline function viscous_velocity(v, viscosity, system, particle, v_particle)
+@propagate_inbounds function viscous_velocity(v, viscosity, system, particle, v_particle)
     # Wall velocity in the viscosity calculation contains the physical wall velocity
     # and an interpolated velocity when a wall viscosity (no-slip BC) is used.
     return extract_svector(system.boundary_model.cache.wall_velocity, system, particle)

src/schemes/structure/rigid_body/system.jl

@@ -506,7 +506,8 @@ end
     return neighbor_system.boundary_model.viscosity
 end
 
-@inline function viscous_velocity(v, system::RigidBodySystem, particle, v_particle)
+@propagate_inbounds function viscous_velocity(v, system::RigidBodySystem,
+                                              particle, v_particle)
     # This function is only used in fluid-structure interaction,
     # so it is never called when `boundary_model` is `nothing`
     return viscous_velocity(v, system.boundary_model.viscosity, system,

src/schemes/structure/total_lagrangian_sph/system.jl

@@ -333,7 +333,8 @@ end
     return current_density(v, system.boundary_model, system)
 end
 
-@inline function viscous_velocity(v, system::TotalLagrangianSPHSystem, particle, v_particle)
+@propagate_inbounds function viscous_velocity(v, system::TotalLagrangianSPHSystem,
+                                              particle, v_particle)
     # This function is only used in fluid-structure interaction,
     # so it is never called when `boundary_model` is `nothing`
     return viscous_velocity(v, system.boundary_model.viscosity, system,