Skip to content

Commit

Permalink
Cosmetic change
Browse files Browse the repository at this point in the history
  • Loading branch information
Kolaru committed Feb 14, 2024
1 parent 5c1d591 commit d56421b
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions src/NormalModes.jl
Original file line number Diff line number Diff line change
Expand Up @@ -162,12 +162,12 @@ momentum.
"""
function StatsBase.sample(rng::AbstractRNG, nm::NormalDecomposition, n_samples)
hbar = 1 # Atomic units
X = nm.M * nm.U
Δx_dist = MvNormal(Diagonal(1/2 * (hbar ./ nm.ωs)))
Δp_dist = MvNormal(Diagonal(1/2 * hbar * nm.ωs))
MU = nm.M * nm.U

Δx = MU * rand(rng, Δx_dist, n_samples)
Δp = inv(nm.M)^2 * MU * rand(rng, Δp_dist, n_samples)
Δx = X * rand(rng, Δx_dist, n_samples)
Δp = inv(nm.M)^2 * X * rand(rng, Δp_dist, n_samples)

return Δx * aunit(u"m"), Δp * aunit(u"kg*m/s")
end
Expand Down

0 comments on commit d56421b

Please sign in to comment.