Commit 51eb665a authored by Olaf's avatar Olaf

Show calculation progress

parent 4f43b32c
......@@ -402,6 +402,8 @@ def dynamics(bodies, D, h, t_max, barnes_hut, theta, approximation, j, startsave
force1, U = force_barneshut(positions, mass, theta)
velo = velocity(velo, force0, force1, h)
force0 = force1
if i%round((T+1)/10) == 0:
print(int(100*round(i/(T+1), 1)), "%")
if i%j == 0 and i>=startsave:
index += 1
all_pos[index] = positions
......@@ -414,6 +416,8 @@ def dynamics(bodies, D, h, t_max, barnes_hut, theta, approximation, j, startsave
force1, U = force_approx(mass, positions, D, n)
velo = velocity(velo, force0, force1, h)
force0 = force1
if i%round((T+1)/10) == 0:
print(int(100*round(i/(T+1), 1)), "%")
if i%j == 0 and i>=startsave:
index += 1
all_pos[index] = positions
......@@ -426,6 +430,8 @@ def dynamics(bodies, D, h, t_max, barnes_hut, theta, approximation, j, startsave
force1, U = force_calc(mass, positions, D)
velo = velocity(velo, force0, force1, h)
force0 = force1
if i%round((T+1)/10) == 0:
print(int(100*round(i/(T+1), 1)), "%")
if i%j == 0 and i>=startsave:
index += 1
all_pos[index] = positions
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment