diff --git a/Bhat_Moussouni_Numerical_Simulation_Report.pdf b/Bhat_Moussouni_Numerical_Simulation_Report.pdf new file mode 100644 index 0000000..d532d46 Binary files /dev/null and b/Bhat_Moussouni_Numerical_Simulation_Report.pdf differ diff --git a/Figs/dt_vs_cpu_time_loglog.pdf b/Figs/dt_vs_cpu_time_loglog.pdf index 4d4859c..2127acc 100644 Binary files a/Figs/dt_vs_cpu_time_loglog.pdf and b/Figs/dt_vs_cpu_time_loglog.pdf differ diff --git a/Figs/orbit_dt.pdf b/Figs/orbit_dt.pdf index 69239df..c4684bb 100644 Binary files a/Figs/orbit_dt.pdf and b/Figs/orbit_dt.pdf differ diff --git a/Figs/orbit_dt_0.0003.pdf b/Figs/orbit_dt_0.0003.pdf index 8dfe130..c4cad32 100644 Binary files a/Figs/orbit_dt_0.0003.pdf and b/Figs/orbit_dt_0.0003.pdf differ diff --git a/Figs/orbit_dt_0.0004.pdf b/Figs/orbit_dt_0.0004.pdf index 9dd6b21..d527842 100644 Binary files a/Figs/orbit_dt_0.0004.pdf and b/Figs/orbit_dt_0.0004.pdf differ diff --git a/Figs/orbit_dt_0.0006.pdf b/Figs/orbit_dt_0.0006.pdf index e50fbdb..5003247 100644 Binary files a/Figs/orbit_dt_0.0006.pdf and b/Figs/orbit_dt_0.0006.pdf differ diff --git a/Figs/orbit_dt_0.0008.pdf b/Figs/orbit_dt_0.0008.pdf index 72ea070..5f370d4 100644 Binary files a/Figs/orbit_dt_0.0008.pdf and b/Figs/orbit_dt_0.0008.pdf differ diff --git a/Figs/orbit_dt_0.0010.pdf b/Figs/orbit_dt_0.0010.pdf index 91ca1a8..0d83565 100644 Binary files a/Figs/orbit_dt_0.0010.pdf and b/Figs/orbit_dt_0.0010.pdf differ diff --git a/Figs/orbit_dt_0.0012.pdf b/Figs/orbit_dt_0.0012.pdf index 806401d..e3bc0f9 100644 Binary files a/Figs/orbit_dt_0.0012.pdf and b/Figs/orbit_dt_0.0012.pdf differ diff --git a/Figs/orbit_dt_0.0016.pdf b/Figs/orbit_dt_0.0016.pdf index 653cd01..9a1deed 100644 Binary files a/Figs/orbit_dt_0.0016.pdf and b/Figs/orbit_dt_0.0016.pdf differ diff --git a/Figs/orbit_dt_0.0022.pdf b/Figs/orbit_dt_0.0022.pdf index 3d90abf..9a7e73e 100644 Binary files a/Figs/orbit_dt_0.0022.pdf and b/Figs/orbit_dt_0.0022.pdf differ diff --git a/Figs/orbit_dt_0.0031.pdf b/Figs/orbit_dt_0.0031.pdf index c681d7c..fa003fa 100644 Binary files a/Figs/orbit_dt_0.0031.pdf and b/Figs/orbit_dt_0.0031.pdf differ diff --git a/Figs/orbit_dt_0.0043.pdf b/Figs/orbit_dt_0.0043.pdf index f287a2b..0c48ca3 100644 Binary files a/Figs/orbit_dt_0.0043.pdf and b/Figs/orbit_dt_0.0043.pdf differ diff --git a/Figs/orbit_dt_0.0060.pdf b/Figs/orbit_dt_0.0060.pdf index 2bd3f7c..ba0e313 100644 Binary files a/Figs/orbit_dt_0.0060.pdf and b/Figs/orbit_dt_0.0060.pdf differ diff --git a/Figs/orbit_dt_0.0083.pdf b/Figs/orbit_dt_0.0083.pdf index bf949d8..dc01994 100644 Binary files a/Figs/orbit_dt_0.0083.pdf and b/Figs/orbit_dt_0.0083.pdf differ diff --git a/Figs/orbit_dt_0.0114.pdf b/Figs/orbit_dt_0.0114.pdf index 3dd2d67..8322001 100644 Binary files a/Figs/orbit_dt_0.0114.pdf and b/Figs/orbit_dt_0.0114.pdf differ diff --git a/Figs/orbit_dt_0.0158.pdf b/Figs/orbit_dt_0.0158.pdf index 4a82bee..47e06e2 100644 Binary files a/Figs/orbit_dt_0.0158.pdf and b/Figs/orbit_dt_0.0158.pdf differ diff --git a/Figs/orbit_dt_0.0220.pdf b/Figs/orbit_dt_0.0220.pdf index e45e3c1..029458a 100644 Binary files a/Figs/orbit_dt_0.0220.pdf and b/Figs/orbit_dt_0.0220.pdf differ diff --git a/Figs/orbit_dt_0.0304.pdf b/Figs/orbit_dt_0.0304.pdf index 440da26..e904619 100644 Binary files a/Figs/orbit_dt_0.0304.pdf and b/Figs/orbit_dt_0.0304.pdf differ diff --git a/Figs/orbit_dt_0.0422.pdf b/Figs/orbit_dt_0.0422.pdf index 74095e8..9c1fc71 100644 Binary files a/Figs/orbit_dt_0.0422.pdf and b/Figs/orbit_dt_0.0422.pdf differ diff --git a/Figs/orbit_dt_0.0584.pdf b/Figs/orbit_dt_0.0584.pdf index 3055e5b..16a3dea 100644 Binary files a/Figs/orbit_dt_0.0584.pdf and b/Figs/orbit_dt_0.0584.pdf differ diff --git a/Figs/orbit_dt_0.0810.pdf b/Figs/orbit_dt_0.0810.pdf index 56cf490..48ca9df 100644 Binary files a/Figs/orbit_dt_0.0810.pdf and b/Figs/orbit_dt_0.0810.pdf differ diff --git a/Figs/orbit_dt_0.1122.pdf b/Figs/orbit_dt_0.1122.pdf index 58f5df9..6cd8c64 100644 Binary files a/Figs/orbit_dt_0.1122.pdf and b/Figs/orbit_dt_0.1122.pdf differ diff --git a/Figs/orbit_dt_0.1555.pdf b/Figs/orbit_dt_0.1555.pdf index a8c2201..f19ad76 100644 Binary files a/Figs/orbit_dt_0.1555.pdf and b/Figs/orbit_dt_0.1555.pdf differ diff --git a/Figs/orbit_dt_0.2154.pdf b/Figs/orbit_dt_0.2154.pdf index 066d4dd..7fc2482 100644 Binary files a/Figs/orbit_dt_0.2154.pdf and b/Figs/orbit_dt_0.2154.pdf differ diff --git a/Figs/orbit_dt_0.2985.pdf b/Figs/orbit_dt_0.2985.pdf index fa9601c..6ead517 100644 Binary files a/Figs/orbit_dt_0.2985.pdf and b/Figs/orbit_dt_0.2985.pdf differ diff --git a/Figs/orbit_dt_0.4137.pdf b/Figs/orbit_dt_0.4137.pdf index fe71974..667fac0 100644 Binary files a/Figs/orbit_dt_0.4137.pdf and b/Figs/orbit_dt_0.4137.pdf differ diff --git a/Figs/orbit_dt_0.5732.pdf b/Figs/orbit_dt_0.5732.pdf index d4220dc..4521af0 100644 Binary files a/Figs/orbit_dt_0.5732.pdf and b/Figs/orbit_dt_0.5732.pdf differ diff --git a/Figs/orbit_dt_0.7943.pdf b/Figs/orbit_dt_0.7943.pdf index 62630c9..d297693 100644 Binary files a/Figs/orbit_dt_0.7943.pdf and b/Figs/orbit_dt_0.7943.pdf differ diff --git a/Figs/pot_hh.png b/Figs/pot_hh.png index 59b7f73..d6727ea 100644 Binary files a/Figs/pot_hh.png and b/Figs/pot_hh.png differ diff --git a/Figs/timestep_vs_final_energy_error_loglog1.pdf b/Figs/timestep_vs_final_energy_error_loglog1.pdf index aec10b4..57263a5 100644 Binary files a/Figs/timestep_vs_final_energy_error_loglog1.pdf and b/Figs/timestep_vs_final_energy_error_loglog1.pdf differ diff --git a/Source/test_integrators.py b/Source/test_integrators.py index b0df030..8e6ae10 100644 --- a/Source/test_integrators.py +++ b/Source/test_integrators.py @@ -140,7 +140,7 @@ for h in h_range: color=colors['RK4'], label="RK4") - ax.set_title("$\\Var{{t}} = {:.4f}$".format(h)) + ax.set_title("$h = {:.4f}$".format(h)) ax.set_xlabel("$x$") ax.set_ylabel("$y$") ax.set_aspect("equal") @@ -177,7 +177,7 @@ for h in h_range: axs[i].set_aspect("equal") - fig.suptitle("$\\Var{{t}} = {:.4f}$".format(h)) + #fig.suptitle("$\\Var{{t}} = {:.4f}$".format(h)) axs[0].set_xlabel("$x$") axs[0].set_ylabel("$y$") axs[0].legend(loc="upper left") diff --git a/Source/test_potentials.py b/Source/test_potentials.py index 59d3f9a..2c537d0 100644 --- a/Source/test_potentials.py +++ b/Source/test_potentials.py @@ -64,7 +64,7 @@ def hh(W): Y = W[1] POT = pot.hh_potential(W, position_only=True) fig, ax = plt.subplots(1) - ax.set_title("Hénon–Heils potential") + ax.set_title("Hénon--Heils potential") pcm = ax.pcolormesh(X, Y, POT) fig.colorbar(pcm, label="potential") ax.set_aspect("equal") @@ -73,7 +73,6 @@ def hh(W): fig.savefig("Figs/pot_hh.png") return 0 - if __name__ == "__main__": kepler(W) hh(W)