Tibeuleu/KozaiLidov
1
0
Fork 0
Code Activity

correct step incrementation

Browse Source
This commit is contained in:
Thibault Barnouin
2021-11-10 13:52:24 +01:00
parent ca9ca60b92
commit e461e96d21
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
lib/objects.py
View File

@@ -132,13 +132,13 @@ class System:
E[j] = self.Eval() E[j] = self.Eval()
L[j] = self.Lval() L[j] = self.Lval()
if display and (j+1)%100==0: if display and j%100==0:
# display progression # display progression
q_array = self.get_positions() q_array = self.get_positions()
if len(self.bodylist) == 1: if len(self.bodylist) == 1:
d.on_running(q_array[0], q_array[1], q_array[2], step=j, label="step {0:d}/{1:d}".format(j+1,N)) d.on_running(q_array[0], q_array[1], q_array[2], step=j, label="step {0:d}/{1:d}".format(j,N))
else: else:
d.on_running(q_array[:,0], q_array[:,1], q_array[:,2], step=j, label="step {0:d}/{1:d}".format(j+1,N)) d.on_running(q_array[:,0], q_array[:,1], q_array[:,2], step=j, label="step {0:d}/{1:d}".format(j,N))
if display: if display:
system("convert -delay 5 -loop 0 tmp/??????.png tmp/temp.gif && rm tmp/??????.png") system("convert -delay 5 -loop 0 tmp/??????.png tmp/temp.gif && rm tmp/??????.png")
system("convert tmp/temp.gif -fuzz 30% -layers Optimize plots/dynsyst.gif && rm tmp/temp.gif") system("convert tmp/temp.gif -fuzz 30% -layers Optimize plots/dynsyst.gif && rm tmp/temp.gif")
@@ -228,13 +228,13 @@ class System:
E[j] = self.Eval() E[j] = self.Eval()
L[j] = self.Lval() L[j] = self.Lval()
if display and (j+1)%100==0: if display and j%100==0:
# display progression # display progression
q_array = self.get_positions() q_array = self.get_positions()
if len(self.bodylist) == 1: if len(self.bodylist) == 1:
d.on_running(q_array[0], q_array[1], q_array[2], step=j, label="step {0:d}/{1:d}".format(j+1,N)) d.on_running(q_array[0], q_array[1], q_array[2], step=j, label="step {0:d}/{1:d}".format(j,N))
else: else:
d.on_running(q_array[:,0], q_array[:,1], q_array[:,2], step=j, label="step {0:d}/{1:d}".format(j+1,N)) d.on_running(q_array[:,0], q_array[:,1], q_array[:,2], step=j, label="step {0:d}/{1:d}".format(j,N))
if display: if display:
system("convert -delay 5 -loop 0 tmp/??????.png tmp/temp.gif && rm tmp/??????.png") system("convert -delay 5 -loop 0 tmp/??????.png tmp/temp.gif && rm tmp/??????.png")
system("convert tmp/temp.gif -fuzz 30% -layers Optimize plots/dynsyst.gif && rm tmp/temp.gif") system("convert tmp/temp.gif -fuzz 30% -layers Optimize plots/dynsyst.gif && rm tmp/temp.gif")

BIN
plots/dynsyst.gif Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 379 KiB