Fix rattle!(...) uniform random in radial component

JuliaMolSim · Oct 23, 2024 · 76995b0 · 76995b0
1 parent c9f37f9
commit 76995b0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/utils.jl b/src/utils.jl
@@ -117,7 +117,7 @@ function rattle!(at::FlexibleSystem, r::Quantity)
       𝐫ᵢ = p.position 
       T = typeof(ustrip(𝐫ᵢ[1]))
       ui = randn(Vec3{T})
-      p_new = _set_position(p, 𝐫ᵢ + r * ui / norm(ui))
+      p_new = _set_position(p, 𝐫ᵢ + rand() * r * ui / norm(ui))
       at.particles[i] = p_new
    end
    return at