Fix float type of radial rand()

JuliaMolSim · Oct 23, 2024 · 9d6c1a4 · 9d6c1a4
1 parent 76995b0
commit 9d6c1a4
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, 𝐫ᵢ + rand() * r * ui / norm(ui))
+      p_new = _set_position(p, 𝐫ᵢ + rand(T) * r * ui / norm(ui))
       at.particles[i] = p_new
    end
    return at