Skip to content

Commit

Permalink
remove kinetic energy output
Browse files Browse the repository at this point in the history
  • Loading branch information
Chenggong committed Nov 8, 2023
1 parent cf719a0 commit 620ac6b
Showing 1 changed file with 1 addition and 6 deletions.
7 changes: 1 addition & 6 deletions Sfilter/util/Langevin.py
Original file line number Diff line number Diff line change
Expand Up @@ -66,13 +66,8 @@ def baoab(potential, max_steps, dt, gamma, kBT, initial_position, initial_veloci
v = velocity_update(v, force, dt)

if step_number % save_frequency == 0 and step_number > 0:
kinetic = .5 * v * v
e_total = kinetic + potential_E

positions.append(x)
# velocities.append(v)
# potential_energies.append(potential_E)
# save_times.append(t)


t = t + dt
step_number = step_number + 1
Expand Down

0 comments on commit 620ac6b

Please sign in to comment.