Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Fixed typos in relativity test scripts
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Nov 18, 2022
1 parent 848cecc commit 60cf506
Show file tree
Hide file tree
Showing 2 changed files with 891 additions and 13 deletions.
10 changes: 5 additions & 5 deletions examples/helio_gr_test/swiftest_relativity.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@
"source": [
"%%capture\n",
"tstep_out = 10.0\n",
"sim_gr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"whm\",general_relativity=True)"
"sim_gr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"helio\",general_relativity=True)"
]
},
{
Expand All @@ -51,7 +51,7 @@
"outputs": [],
"source": [
"%%capture\n",
"sim_nogr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"whm\",general_relativity=False)"
"sim_nogr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"helio\",general_relativity=False)"
]
},
{
Expand Down Expand Up @@ -124,12 +124,12 @@
"fig, ax = plt.subplots()\n",
"\n",
"ax.plot(t, varpi_obs, label=\"JPL Horizons\",linewidth=2.5)\n",
"ax.plot(tsim, varpisim_gr, label=\"Swiftest WHM GR\",linewidth=1.5)\n",
"ax.plot(tsim, varpisim_nogr, label=\"Swiftest WHM No GR\",linewidth=1.5)\n",
"ax.plot(tsim, varpisim_gr, label=\"Swiftest helio GR\",linewidth=1.5)\n",
"ax.plot(tsim, varpisim_nogr, label=\"Swiftest helio No GR\",linewidth=1.5)\n",
"ax.set_xlabel('Time (y)')\n",
"ax.set_ylabel('Mercury $\\\\varpi$ (deg)')\n",
"ax.legend()\n",
"plt.savefig(\"whm_gr_mercury_precession.png\",dpi=300)\n",
"plt.savefig(\"helio_gr_mercury_precession.png\",dpi=300)\n",
"print('Mean precession rate for Mercury long. peri. (arcsec/100 y)')\n",
"print(f'JPL Horizons : {np.mean(dvarpi_obs)}')\n",
"print(f'Swiftest No GR : {np.mean(dvarpi_nogr)}')\n",
Expand Down
Loading

0 comments on commit 60cf506

Please sign in to comment.