File tree Expand file tree Collapse file tree 3 files changed +14
-15
lines changed Expand file tree Collapse file tree 3 files changed +14
-15
lines changed Original file line number Diff line number Diff line change 11from vpython import *
2- from time import clock
32
43# Hard-sphere gas.
54
4847p = []
4948apos = []
5049pavg = sqrt (2 * mass * 1.5 * k * T ) # average kinetic energy p**2/(2mass) = (3/2)kT
51-
50+
5251for i in range (Natoms ):
5352 x = L * random ()- L / 2
5453 y = L * random ()- L / 2
@@ -152,24 +151,24 @@ def checkCollisions():
152151
153152 # Update all positions
154153 for i in range (Natoms ): Atoms [i ].pos = apos [i ] = apos [i ] + (p [i ]/ mass )* dt
155-
154+
156155 # Check for collisions
157156 hitlist = checkCollisions ()
158157
159-
158+
160159 for i in range (Natoms ):
161160 loc = apos [i ]
162161 if abs (loc .x ) > L / 2 :
163162 if loc .x < 0 : p [i ].x = abs (p [i ].x )
164163 else : p [i ].x = - abs (p [i ].x )
165-
164+
166165 if abs (loc .y ) > L / 2 :
167166 if loc .y < 0 : p [i ].y = abs (p [i ].y )
168167 else : p [i ].y = - abs (p [i ].y )
169-
168+
170169 if abs (loc .z ) > L / 2 :
171170 if loc .z < 0 : p [i ].z = abs (p [i ].z )
172171 else : p [i ].z = - abs (p [i ].z )
173-
172+
174173timer = clock ()- timer
175174print (timer )
Original file line number Diff line number Diff line change 11import time
2- try :
3- _clock = time .perf_counter # time.clock is deprecated in Python 3.3, gone in 3.8
4- except :
5- _clock = time . clock
2+
3+ _clock = time .perf_counter
4+
5+
66_tick = 1 / 60
77
88#import platform
Original file line number Diff line number Diff line change 66from math import sqrt , tan , pi
77
88import time
9- try :
10- clock = time . perf_counter # time.clock is deprecated in Python 3.3, gone in 3.8
11- except :
12- clock = time . clock
9+
10+ # vpython provides clock in its namespace
11+ clock = time . perf_counter
12+
1313import sys
1414from . import __version__ , __gs_version__
1515from ._notebook_helpers import _isnotebook
You can’t perform that action at this time.
0 commit comments