Tibeuleu/KozaiLidov
1
0
Fork 0
Code Activity

debugg gif savegit status

Browse Source
This commit is contained in:
Thibault Barnouin
2022-01-07 13:26:07 +01:00
parent d64d947f44
commit 5728ff39f1
4 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
lib/LeapFrog.py
View File

@@ -67,7 +67,7 @@ def leapfrog(dyn_syst, bin_syst, duration, dt, recover_param=False, display=Fals
d.on_running(dyn_syst, step=step, label="{0:.2f} years".format(j*dt/yr)) d.on_running(dyn_syst, step=step, label="{0:.2f} years".format(j*dt/yr))
if display: if display:
d.close() d.close()
if not savename is None: if gif:
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 10% -layers Optimize plots/{0:s}_dynsyst.gif".format(savename)) system("convert tmp/temp.gif -fuzz 10% -layers Optimize plots/{0:s}_dynsyst.gif".format(savename))

2
lib/hermite.py
View File

@@ -114,7 +114,7 @@ def hermite(dyn_syst, bin_syst, duration, dt, recover_param=False, display=False
d.on_running(dyn_syst, step=step, label="{0:.2f} years".format(j*dt/yr)) d.on_running(dyn_syst, step=step, label="{0:.2f} years".format(j*dt/yr))
if display: if display:
d.close() d.close()
if not savename is None: if gif:
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 10% -layers Optimize plots/{0:s}_dynsyst.gif".format(savename)) system("convert tmp/temp.gif -fuzz 10% -layers Optimize plots/{0:s}_dynsyst.gif".format(savename))

4
main.py
View File

@@ -36,7 +36,7 @@ def main():
savename = "{0:d}bodies_{1:s}".format(n_bodies, integrator) savename = "{0:d}bodies_{1:s}".format(n_bodies, integrator)
#simulation start #simulation start
E, L = [], [] E, L, sma, ecc = [], [], [], []
for i,step0 in enumerate(step): for i,step0 in enumerate(step):
bodylist = [] bodylist = []
for j in range(n_bodies): for j in range(n_bodies):
@@ -52,6 +52,8 @@ def main():
E0, L0, sma, ecc = hermite(dyn_syst, bin_syst, duration, step0, recover_param=True, display=display, savename=savename, gif=gif) E0, L0, sma, ecc = hermite(dyn_syst, bin_syst, duration, step0, recover_param=True, display=display, savename=savename, gif=gif)
E.append(E0) E.append(E0)
L.append(L0) L.append(L0)
sma.append(sma0)
ecc.append(ecc0)
parameters = [duration, step, dyn_syst, integrator] parameters = [duration, step, dyn_syst, integrator]
display_parameters(E, L, sma, ecc, parameters=parameters, savename=savename) display_parameters(E, L, sma, ecc, parameters=parameters, savename=savename)

2
main_concurrent.py
View File

@@ -29,7 +29,7 @@ def main():
v = np.array([v1, v2, v3],dtype=np.longdouble) v = np.array([v1, v2, v3],dtype=np.longdouble)
#integration parameters #integration parameters
duration, step = 5000.*yr, np.array([10./(365.25)],dtype=np.longdouble)*yr #integration time and step in seconds duration, step = 100.*yr, np.array([1./(365.25)],dtype=np.longdouble)*yr #integration time and step in seconds
step = np.sort(step)[::-1] step = np.sort(step)[::-1]
integrator = "leapfrog" integrator = "leapfrog"
n_bodies = 3 n_bodies = 3