diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/cb.swiftest.in b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/cb.swiftest.in index e343dee10..2ce982a28 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/cb.swiftest.in +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/cb.swiftest.in @@ -1,4 +1,4 @@ -39.47692640889762629 +0.00029591220819207776388 0.0046504672609621575315 4.7535806948127356533e-12 -2.2473967953572827815e-18 diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/config.swiftest.in b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/config.swiftest.in index 7a75d8951..ce9a380aa 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/config.swiftest.in +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/config.swiftest.in @@ -1,7 +1,7 @@ ! Swiftest input file generated using init_cond.py T0 0 -TSTOP 1.0 -DT 0.0027378507871321013 +TSTOP 365.25 +DT 1.0 CB_IN cb.swiftest.in PL_IN pl.swiftest.in TP_IN tp.in @@ -10,7 +10,7 @@ ISTEP_OUT 1 ISTEP_DUMP 1 BIN_OUT bin.swiftest.dat OUT_TYPE REAL8 -OUT_FORM XV +OUT_FORM EL OUT_STAT REPLACE CHK_CLOSE yes CHK_RMIN 0.004650467260962157 @@ -26,4 +26,4 @@ ROTATION no GR no MU2KG 1.988409870698051e+30 DU2M 149597870700.0 -TU2S 31557600.0 +TU2S 86400.0 diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/init_cond.py b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/init_cond.py index fc20896c6..d38f942e6 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/init_cond.py +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/init_cond.py @@ -20,13 +20,13 @@ MU2KG = np.longdouble(GMSunSI / GC) #Conversion from mass unit to kg DU2M = np.longdouble(AU2M) #Conversion from radius unit to centimeters -TU2S = np.longdouble(year) #Conversion from time unit to seconds +TU2S = np.longdouble(JD) #Conversion from time unit to seconds GU = np.longdouble(GC / (DU2M**3 / (MU2KG * TU2S**2))) GMSun = np.longdouble(GMSunSI / (DU2M**3 / TU2S**2)) # Solar oblatenes values: From Mecheri et al. (2004), using Corbard (b) 2002 values (Table II) -J2 = 0.0 #np.longdouble(2.198e-7) * (Rsun / DU2M)**2 -J4 = 0.0 #np.longdouble(-4.805e-9) * (Rsun / DU2M)**4 +J2 = np.longdouble(2.198e-7) * (Rsun / DU2M)**2 +J4 = np.longdouble(-4.805e-9) * (Rsun / DU2M)**4 npl = 9 ntp = 2 * npl @@ -190,7 +190,7 @@ print(f'ISTEP_DUMP {iout:d}') print(f'BIN_OUT {swifter_bin}') print(f'OUT_TYPE REAL8') -print(f'OUT_FORM XV') +print(f'OUT_FORM EL') print(f'OUT_STAT NEW') print(f'J2 {J2}') print(f'J4 {J4}') @@ -221,7 +221,7 @@ print(f'ISTEP_DUMP {iout:d}') print(f'BIN_OUT {swiftest_bin}') print(f'OUT_TYPE REAL8') -print(f'OUT_FORM XV') +print(f'OUT_FORM EL') print(f'OUT_STAT REPLACE') print(f'CHK_CLOSE yes') print(f'CHK_RMIN {rmin}') diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/param.swifter.in b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/param.swifter.in index b6827600e..1f400dfb2 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/param.swifter.in +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/param.swifter.in @@ -1,7 +1,7 @@ ! Swifter input file generated using init_cond.py T0 0 -TSTOP 1.0 -DT 0.0027378507871321013 +TSTOP 365.25 +DT 1.0 PL_IN pl.swifter.in TP_IN tp.in IN_TYPE ASCII @@ -9,7 +9,7 @@ ISTEP_OUT 1 ISTEP_DUMP 1 BIN_OUT bin.swifter.dat OUT_TYPE REAL8 -OUT_FORM XV +OUT_FORM EL OUT_STAT NEW J2 4.7535806948127355e-12 J4 -2.2473967953572827e-18 diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swifter.in b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swifter.in index 55a416e66..c1ddfc2bf 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swifter.in +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swifter.in @@ -1,40 +1,40 @@ 10 -1 39.47692640889762629 +1 0.00029591220819207776388 0.0 0.0 0.0 0.0 0.0 0.0 -2 6.553709809565314e-06 0.0014751258227142052 ! mercury +2 4.9125474498983625e-11 0.0014751258227142052 ! mercury 1.6306381826061646e-05 0.008059842448018334 -0.4616051037329109 -0.03846017738329229 0.02248719132054853 0.001934639213990692 -0.001904656977422976 -3 9.663313399581539e-05 0.006759134232034941 ! venus +3 7.243452483873647e-10 0.006759134232034941 ! venus 4.0453784346544176e-05 -0.5115875215389065 0.5030818749037324 0.03642547299277956 -0.01425515725454357 -0.01452868630179309 0.0006232072038298823 -4 0.00012002693582795246 0.010044625087011915 ! earthmoon +4 8.997011382166019e-10 0.010044625087011915 ! earthmoon 4.25875607065041e-05 -0.1090020607540907 -1.009893805009766 4.823302918632528e-05 0.01682491922568941 -0.001910549762056979 3.992660742687128e-08 -5 1.2739802010675942e-05 0.0072467897902424765 ! mars +5 9.549535102761465e-11 0.0072467897902424765 ! mars 2.2657408050928896e-05 -1.342897929331636 0.9778655112682739 0.05343398538723887 -0.007712315645393206 -0.01011917844182223 -2.287744801261131e-05 -6 0.037692251088985676 0.3552720805286442 ! jupiter +6 2.825345908631355e-07 0.3552720805286442 ! jupiter 0.0004673261703049093 3.923184193414315 -3.168419770483168 -0.0746147877972047 0.004655552638985802 0.006232623300954468 -0.0001300429201057457 -7 0.011285899820091273 0.4376460836930155 ! saturn +7 8.459715183006416e-08 0.4376460836930155 ! saturn 0.00038925687730393614 6.185794462795267 -7.804174837804826 -0.110498432926239 0.004066833203985018 0.003458637040736611 -0.0002219310939327014 -8 0.001723658947826773 0.46946272948265794 ! uranus +8 1.2920249163736674e-08 0.46946272948265794 ! uranus 0.00016953449859497232 14.9290976575471 12.92949673572929 -0.1454099139559955 -0.002599557960646664 0.002795888198858545 4.391864857782088e-05 -9 0.0020336100526728304 0.78119478483336 ! neptune +9 1.5243589003230834e-08 0.78119478483336 ! neptune 0.00016458790412449367 29.54416169025338 -4.716921603714237 -0.5838030174427992 0.0004792636209523189 0.00312573757291745 -7.53264045199501e-05 -10 2.924216771029454e-07 0.05379680851617536 ! plutocharon +10 2.1919422829042796e-12 0.05379680851617536 ! plutocharon 7.943294877391593e-06 14.54448346259197 -31.05223519593471 -0.8828000265625595 0.002923077617691739 0.0006625916902153526 -0.0009142553677224461 diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swiftest.in b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swiftest.in index c0d273ca7..02a8f9bef 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swiftest.in +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/pl.swiftest.in @@ -1,37 +1,37 @@ 9 -2 6.553709809565314e-06 ! mercury +2 4.9125474498983625e-11 ! mercury 1.6306381826061646e-05 0.008059842448018334 -0.4616051037329109 -0.03846017738329229 0.02248719132054853 0.001934639213990692 -0.001904656977422976 -3 9.663313399581539e-05 ! venus +3 7.243452483873647e-10 ! venus 4.0453784346544176e-05 -0.5115875215389065 0.5030818749037324 0.03642547299277956 -0.01425515725454357 -0.01452868630179309 0.0006232072038298823 -4 0.00012002693582795246 ! earthmoon +4 8.997011382166019e-10 ! earthmoon 4.25875607065041e-05 -0.1090020607540907 -1.009893805009766 4.823302918632528e-05 0.01682491922568941 -0.001910549762056979 3.992660742687128e-08 -5 1.2739802010675942e-05 ! mars +5 9.549535102761465e-11 ! mars 2.2657408050928896e-05 -1.342897929331636 0.9778655112682739 0.05343398538723887 -0.007712315645393206 -0.01011917844182223 -2.287744801261131e-05 -6 0.037692251088985676 ! jupiter +6 2.825345908631355e-07 ! jupiter 0.0004673261703049093 3.923184193414315 -3.168419770483168 -0.0746147877972047 0.004655552638985802 0.006232623300954468 -0.0001300429201057457 -7 0.011285899820091273 ! saturn +7 8.459715183006416e-08 ! saturn 0.00038925687730393614 6.185794462795267 -7.804174837804826 -0.110498432926239 0.004066833203985018 0.003458637040736611 -0.0002219310939327014 -8 0.001723658947826773 ! uranus +8 1.2920249163736674e-08 ! uranus 0.00016953449859497232 14.9290976575471 12.92949673572929 -0.1454099139559955 -0.002599557960646664 0.002795888198858545 4.391864857782088e-05 -9 0.0020336100526728304 ! neptune +9 1.5243589003230834e-08 ! neptune 0.00016458790412449367 29.54416169025338 -4.716921603714237 -0.5838030174427992 0.0004792636209523189 0.00312573757291745 -7.53264045199501e-05 -10 2.924216771029454e-07 ! plutocharon +10 2.1919422829042796e-12 ! plutocharon 7.943294877391593e-06 14.54448346259197 -31.05223519593471 -0.8828000265625595 0.002923077617691739 0.0006625916902153526 -0.0009142553677224461 diff --git a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb index b8e418539..ed3e55e8a 100644 --- a/examples/rmvs_swifter_comparison/9pl_18tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb +++ b/examples/rmvs_swifter_comparison/9pl_18tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb @@ -67,6 +67,508 @@ "swiftdiff = swiftdiff.rename({'time' : 'time (d)'})" ] }, + { + "cell_type": "code", + "execution_count": 25, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
<xarray.DataArray 'e' (time: 366)>\n", + "array([0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999247, 0.99999247, 0.99999247,\n", + " 0.99999247, 0.99999247, 0.99999247, 0.99999247, 0.99999248,\n", + " 0.99999249, 0.99999264, 0.9999928 , 0.99999281, 0.99999282,\n", + " 0.99999282, 0.99999283, 0.99999283, 0.99999283, 0.99999283,\n", + " 0.99999283, 0.99999283, 0.99999283, 0.99999283, 0.99999283,\n", + " 0.99999283, 0.99999283, 0.99999283, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + "...\n", + " 0.99999318, 0.99999318, 0.99999318, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.9999932 , 0.99999319, 0.99999319, 0.9999932 ,\n", + " 0.9999932 , 0.9999932 , 0.99999321, 0.99999321, 0.99999321,\n", + " 0.99999321, 0.99999322, 0.99999322, 0.99999322, 0.99999322,\n", + " 0.99999322, 0.99999322, 0.99999322, 0.99999322, 0.99999322,\n", + " 0.99999322])\n", + "Coordinates:\n", + " id int64 3\n", + " * time (time) float64 0.0 0.002738 0.005476 ... 0.9938 0.9966 0.9993
array([0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999245,\n", + " 0.99999245, 0.99999245, 0.99999245, 0.99999245, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999246, 0.99999246, 0.99999246,\n", + " 0.99999246, 0.99999246, 0.99999247, 0.99999247, 0.99999247,\n", + " 0.99999247, 0.99999247, 0.99999247, 0.99999247, 0.99999248,\n", + " 0.99999249, 0.99999264, 0.9999928 , 0.99999281, 0.99999282,\n", + " 0.99999282, 0.99999283, 0.99999283, 0.99999283, 0.99999283,\n", + " 0.99999283, 0.99999283, 0.99999283, 0.99999283, 0.99999283,\n", + " 0.99999283, 0.99999283, 0.99999283, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + " 0.99999284, 0.99999284, 0.99999284, 0.99999284, 0.99999284,\n", + "...\n", + " 0.99999318, 0.99999318, 0.99999318, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.99999319, 0.99999319, 0.99999319, 0.99999319,\n", + " 0.99999319, 0.9999932 , 0.99999319, 0.99999319, 0.9999932 ,\n", + " 0.9999932 , 0.9999932 , 0.99999321, 0.99999321, 0.99999321,\n", + " 0.99999321, 0.99999322, 0.99999322, 0.99999322, 0.99999322,\n", + " 0.99999322, 0.99999322, 0.99999322, 0.99999322, 0.99999322,\n", + " 0.99999322])
array(3)
array([0. , 0.002738, 0.005476, ..., 0.99384 , 0.996578, 0.999316])