grss 4.2.0__tar.gz → 4.3.4__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- {grss-4.2.0/grss.egg-info → grss-4.3.4}/PKG-INFO +1 -1
- {grss-4.2.0 → grss-4.3.4}/grss/fit/fit_ades.py +1 -0
- {grss-4.2.0 → grss-4.3.4}/grss/fit/fit_optical.py +2 -1
- {grss-4.2.0 → grss-4.3.4}/grss/fit/fit_radar.py +0 -9
- {grss-4.2.0 → grss-4.3.4}/grss/fit/fit_simulation.py +290 -118
- {grss-4.2.0 → grss-4.3.4}/grss/kernels/get_kernels.py +10 -10
- {grss-4.2.0 → grss-4.3.4}/grss/prop/prop_utils.py +25 -7
- {grss-4.2.0 → grss-4.3.4}/grss/utils.py +13 -0
- grss-4.3.4/grss/version.txt +1 -0
- {grss-4.2.0 → grss-4.3.4/grss.egg-info}/PKG-INFO +1 -1
- {grss-4.2.0 → grss-4.3.4}/include/approach.h +6 -0
- {grss-4.2.0 → grss-4.3.4}/include/gr15.h +18 -6
- {grss-4.2.0 → grss-4.3.4}/include/grss.h +5 -2
- {grss-4.2.0 → grss-4.3.4}/include/simulation.h +69 -25
- {grss-4.2.0 → grss-4.3.4}/include/stm.h +6 -0
- {grss-4.2.0 → grss-4.3.4}/include/utilities.h +10 -0
- grss-4.3.4/src/approach.cpp +1097 -0
- {grss-4.2.0 → grss-4.3.4}/src/force.cpp +143 -0
- {grss-4.2.0 → grss-4.3.4}/src/gr15.cpp +102 -42
- {grss-4.2.0 → grss-4.3.4}/src/grss.cpp +131 -58
- {grss-4.2.0 → grss-4.3.4}/src/interpolate.cpp +20 -28
- {grss-4.2.0 → grss-4.3.4}/src/observe.cpp +101 -66
- {grss-4.2.0 → grss-4.3.4}/src/pck.cpp +15 -3
- {grss-4.2.0 → grss-4.3.4}/src/simulation.cpp +207 -47
- {grss-4.2.0 → grss-4.3.4}/src/stm.cpp +51 -0
- {grss-4.2.0 → grss-4.3.4}/tests/cpp/prop/pck_map.cpp +5 -5
- grss-4.3.4/tests/python/fit/chesley.ipynb +275 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/fit/eggl.ipynb +17 -17
- grss-4.3.4/tests/python/fit/farnocchia.ipynb +275 -0
- grss-4.3.4/tests/python/fit/shantanunaidu.ipynb +277 -0
- grss-4.2.0/grss/version.txt +0 -1
- grss-4.2.0/src/approach.cpp +0 -984
- grss-4.2.0/tests/python/fit/chesley.ipynb +0 -275
- grss-4.2.0/tests/python/fit/farnocchia.ipynb +0 -281
- grss-4.2.0/tests/python/fit/shantanunaidu.ipynb +0 -277
- {grss-4.2.0 → grss-4.3.4}/CMakeLists.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/LICENSE +0 -0
- {grss-4.2.0 → grss-4.3.4}/MANIFEST.in +0 -0
- {grss-4.2.0 → grss-4.3.4}/README.md +0 -0
- {grss-4.2.0 → grss-4.3.4}/build_cpp.sh +0 -0
- {grss-4.2.0 → grss-4.3.4}/docs/source/conf.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/extern/get_cspice.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/__init__.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/debias/get_debiasing_data.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/fit/__init__.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/fit/fit_utils.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/prop/__init__.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/prop/prop_parallel.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss/prop/prop_unscented.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss.egg-info/SOURCES.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss.egg-info/dependency_links.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss.egg-info/requires.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/grss.egg-info/top_level.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/elements.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/force.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/interpolate.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/observe.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/parallel.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/pck.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/spk.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/include/timeconvert.h +0 -0
- {grss-4.2.0 → grss-4.3.4}/pyproject.toml +0 -0
- {grss-4.2.0 → grss-4.3.4}/requirements.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/setup.cfg +0 -0
- {grss-4.2.0 → grss-4.3.4}/setup.py +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/CMakeLists.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/elements.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/parallel.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/spk.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/timeconvert.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/src/utilities.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/cpp/prop/CMakeLists.txt +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/cpp/prop/apophis.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/cpp/prop/didymos.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/cpp/prop/spk_map.cpp +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/prop/apophis.ipynb +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/prop/didymos.ipynb +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/prop/icarus.ipynb +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/prop/moshup.ipynb +0 -0
- {grss-4.2.0 → grss-4.3.4}/tests/python/run_tests.sh +0 -0
|
@@ -479,7 +479,8 @@ def apply_debiasing_scheme(obs_df, lowres, verbose):
|
|
|
479
479
|
unbias_counter += len(group)
|
|
480
480
|
continue
|
|
481
481
|
if cat_code not in biased_catalogs:
|
|
482
|
-
|
|
482
|
+
if verbose:
|
|
483
|
+
print(f"\tUnknown star catalog: {cat}")
|
|
483
484
|
obs_df.loc[group.index, 'biasRA'] = 0.0
|
|
484
485
|
obs_df.loc[group.index, 'biasDec'] = 0.0
|
|
485
486
|
continue
|
|
@@ -71,7 +71,6 @@ def add_radar_obs(obs_df, t_min_tdb=None, t_max_tdb=None, verbose=False):
|
|
|
71
71
|
print(f"Read in {num_obs} radar observations from JPL radar API.")
|
|
72
72
|
data = raw_data['data']
|
|
73
73
|
time_range_count = 0
|
|
74
|
-
pre1972_count = 0
|
|
75
74
|
for i in range(num_obs):
|
|
76
75
|
obs = data[num_obs-i-1]
|
|
77
76
|
date = Time(obs[1], format='iso', scale='utc')
|
|
@@ -105,14 +104,6 @@ def add_radar_obs(obs_df, t_min_tdb=None, t_max_tdb=None, verbose=False):
|
|
|
105
104
|
obs_df.loc[idx,'frq'] = freq
|
|
106
105
|
obs_df.loc[idx,'sigDelay'] = obs_sigma if delay else np.nan
|
|
107
106
|
obs_df.loc[idx,'sigDoppler'] = obs_sigma if doppler else np.nan
|
|
108
|
-
if date.utc.mjd < 41317:
|
|
109
|
-
pre1972_count += 1
|
|
110
|
-
obs_df.loc[idx,'selAst'] = 'd'
|
|
111
|
-
if pre1972_count > 0:
|
|
112
|
-
print(f"\tWARNING: {pre1972_count} radar observations were taken before 1972.\n"
|
|
113
|
-
"\t\tThese residuals cannot be reliably computed because high-accuracy Earth "
|
|
114
|
-
"orientation kernels are not available for this time period.\n"
|
|
115
|
-
"\t\tForce deleting them...")
|
|
116
107
|
if verbose:
|
|
117
108
|
print(f"\tFiltered to {num_obs-time_range_count} observations that satisfy the "
|
|
118
109
|
"time range constraints.")
|