... | ... | @@ -17,3 +17,23 @@ How to define accuracy? |
|
|
Run for 10 orbital times, find distance between initial and final position, and see how this changes with 'eta'
|
|
|
|
|
|
![eta_test_plot](/uploads/d5d23c553f45cd21bfcbee79cbd9f8e9/eta_test_plot.png)
|
|
|
|
|
|
Also calculated analytic position at the times of the snapshots, and computed orbit with scipy ODE solver (using 10^6 timesteps per orbit)
|
|
|
|
|
|
![orbit_plot_99](/uploads/cf391e9143a6bb3476c74933efe9c596/orbit_plot_99.png)
|
|
|
|
|
|
From this we can calculate the error.
|
|
|
|
|
|
![orbit_plot_error_99](/uploads/b2970e35495b58988566e4ca976d5586/orbit_plot_error_99.png)
|
|
|
|
|
|
Strange feature: If I start the orbit from periapsis rather that apoapsis, then the errors are much worse (seems like the orbit loses energy somehow)
|
|
|
|
|
|
![orbit_plot_101](/uploads/61e87602ccdc96f795e5a112a73537af/orbit_plot_101.png)
|
|
|
|
|
|
![orbit_plot_error_101](/uploads/d549b4383a7ebe744c844251b58237fe/orbit_plot_error_101.png)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|