def nbody_debug(bodies, time_step): """Run simulation with visualization""" x_max = 500 y_max = 500 z_max = 500 galaxy = random_galaxy(x_max, y_max, z_max, bodies) gfx = gfx_init(x_max, y_max, z_max) dt = 100.0 # One hundred year timesteps ensures that we see movement start = time.time() for _ in range(time_step): move(galaxy, dt) show(gfx, galaxy) stop = time.time() print 'Simulated ' + str(bodies) + ' bodies for ' + str(time_step) \ + ' timesteps in ' + str(stop - start) + ' seconds'
def nbody_benchmark(bodies_list, time_step): """Run benchmark simulation without visualization""" x_max = 1e18 y_max = 1e18 z_max = 1e18 dt = 1e12 for i in range(3): for bodies in bodies_list: solarsystem, asteroids = random_system(x_max, y_max, z_max, 10, bodies) start = time.time() for _ in range(time_step): move(solarsystem, asteroids, dt) stop = time.time() print "Simulated " + str(bodies) + " bodies for " + str(time_step) + " timesteps in " + str( stop - start ) + " seconds"
def nbody_benchmark(bodies_list, time_step): """Run benchmark simulation without visualization""" x_max = 500 y_max = 500 z_max = 500 dt = 1.0 # One year timesteps for better accuracy for bodies in bodies_list: galaxy = random_galaxy(x_max, y_max, z_max, bodies) start = time.time() for _ in range(time_step): move(galaxy, dt) stop = time.time() print 'Simulated ' + str(bodies) + ' bodies for ' \ + str(time_step) + ' timesteps in ' + str(stop - start) \ + ' seconds'
def nbody_debug(n, bodies, time_step): """Run simulation with visualization""" x_max = 1e18 y_max = 1e18 z_max = 1e18 solarsystem, asteroids= random_system(x_max, y_max, z_max, n, bodies) P3 = gfx_init(x_max, y_max, z_max) dt = 1e12 # One hundred year timesteps ensures that we see movement start = time.time() for step in range(time_step): if step%1 == 0: show(P3, solarsystem, asteroids) move(solarsystem, asteroids, dt) print step stop = time.time() print 'Simulated ' + str(bodies) + ' bodies for ' + str(time_step) \ + ' timesteps in ' + str(stop - start) + ' seconds'