Bump pkg versions [skip ci]

Author: islent

Project.toml

@@ -16,7 +16,6 @@ DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
 Optim = "429524aa-4258-5aef-a3af-852621145aeb"
 PhysicalConstants = "5ad8b20f-a522-5ce9-bfc9-ddf1d5bda6ab"
 PhysicalParticles = "78724250-4eb2-11e9-0f34-13a7aab864c5"
-Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 PrecompileTools = "aea7be01-6a6a-4083-8856-8a6e6704d82a"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
@@ -37,11 +36,10 @@ DocStringExtensions = "0.8, 0.9"
 Optim = "1"
 PhysicalConstants = "0.2"
 PhysicalParticles = "1"
-Plots = "1"
 PrecompileTools = "1"
 QuadGK = "2"
 Reexport = "1"
-StructArrays = "0.6"
+StructArrays = "0.6, 0.7"
 Unitful = "1"
 UnitfulAstro = "1"
 julia = "1.6"

src/Tools.jl

@@ -87,23 +87,23 @@ $(TYPEDSIGNATURES)
 Rotate along the positive z-axis.
 - `ratio`: rotational motion v.s. random motion. If equals `1`, only rotational component
 """
-function rotational_velocity_acc(x, y, z, a, ratio = 1.0)
+function rotational_velocity_acc(x, y, z, a, rot_ratio = 1.0)
     u = unit(x)
     r = sqrt(x^2 + y^2 + z^2)
     v = sqrt(a * r)
-    v_vec = -normalize(PVector(ustrip(u,x), ustrip(u,y), 0.0) × PVector(0.0, 0.0, 1.0)) * v
-    v_vec = ratio * v_vec + (1-ratio) * randn(PVector{Float64}) * v # this will change the average velocity
-    v_vec = v_vec + (1-ratio) * randn(PVector{Float64}) * v
+    v_vec_rot = -normalize(PVector(ustrip(u,x), ustrip(u,y), 0.0) × PVector(0.0, 0.0, 1.0)) * v
+    v_vec = rot_ratio * v_vec_rot + (1-rot_ratio) * normalize(randn(PVector{Float64})) * v 
+    # v_vec = v_vec_rot + (1-rot_ratio) * randn(PVector{Float64}) * v # this will change the average velocity
 end
 
 """
 $(TYPEDSIGNATURES)
 Rotate along the positive z-axis.
-- `ratio`: rotational motion v.s. random motion. If equals `1`, only rotational component
+- `rot_ratio`: rotational motion v.s. random motion. If equals `1`, only rotational component
 """
-function rotational_velocity(x, y, v, ratio = 1.0)
+function rotational_velocity(x, y, v, rot_ratio = 1.0)
     u = unit(x)
-    v_vec = -normalize(PVector(ustrip(u,x), ustrip(u,y), 0.0) × PVector(0.0, 0.0, 1.0)) * v
-    # v_vec = ratio * v_vec + (1-ratio) * randn(PVector{Float64}) * v # this will change the average velocity
-    v_vec = v_vec + (1-ratio) * randn(PVector{Float64}) * v
+    v_vec_rot = -normalize(PVector(ustrip(u,x), ustrip(u,y), 0.0) × PVector(0.0, 0.0, 1.0)) * v
+    v_vec = rot_ratio * v_vec_rot + (1-rot_ratio) * normalize(randn(PVector{Float64})) * v 
+    # v_vec = v_vec_rot + (1-rot_ratio) * randn(PVector{Float64}) * v # this will change the average velocity
 end