nrl-tracker 0.21.5__tar.gz → 0.22.1__tar.gz

{nrl_tracker-0.21.5/nrl_tracker.egg-info → nrl_tracker-0.22.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 0.21.5
+Version: 0.22.1
 Summary: Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms
 Author: Original: David F. Crouse, Naval Research Laboratory
 Maintainer: Python Port Contributors
@@ -60,7 +60,7 @@ Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualizati
 
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v0.21.5-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v0.22.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/README.md

@@ -1,6 +1,6 @@
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v0.21.5-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v0.22.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1/nrl_tracker.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 0.21.5
+Version: 0.22.1
 Summary: Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms
 Author: Original: David F. Crouse, Naval Research Laboratory
 Maintainer: Python Port Contributors
@@ -60,7 +60,7 @@ Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualizati
 
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v0.21.5-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v0.22.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/nrl_tracker.egg-info/SOURCES.txt

@@ -19,9 +19,11 @@ pytcl/assignment_algorithms/two_dimensional/__init__.py
 pytcl/assignment_algorithms/two_dimensional/assignment.py
 pytcl/assignment_algorithms/two_dimensional/kbest.py
 pytcl/astronomical/__init__.py
+pytcl/astronomical/ephemerides.py
 pytcl/astronomical/lambert.py
 pytcl/astronomical/orbital_mechanics.py
 pytcl/astronomical/reference_frames.py
+pytcl/astronomical/relativity.py
 pytcl/astronomical/time_systems.py
 pytcl/atmosphere/__init__.py
 pytcl/atmosphere/models.py
@@ -163,6 +165,7 @@ tests/test_coverage_boost_2.py
 tests/test_dynamic_models.py
 tests/test_egm.py
 tests/test_emm.py
+tests/test_ephemerides.py
 tests/test_gaussian_mixtures.py
 tests/test_geophysical.py
 tests/test_great_circle.py
@@ -176,6 +179,7 @@ tests/test_performance_evaluation.py
 tests/test_phase6_specialized.py
 tests/test_plotting.py
 tests/test_projections.py
+tests/test_relativity.py
 tests/test_rhumb.py
 tests/test_signal_processing.py
 tests/test_smoothers.py

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "nrl-tracker"
-version = "0.21.5"
+version = "0.22.1"
 description = "Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms"
 readme = "README.md"
 authors = [

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/__init__.py

@@ -20,7 +20,7 @@ References
        no. 5, pp. 18-27, May 2017.
 """
 
-__version__ = "0.21.5"
+__version__ = "0.22.1"
 __author__ = "Python Port Contributors"
 __original_author__ = "David F. Crouse, Naval Research Laboratory"
 

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/assignment_algorithms/data_association.py

@@ -10,13 +10,8 @@ from typing import List, NamedTuple, Optional
 import numpy as np
 from numpy.typing import ArrayLike, NDArray
 
-from pytcl.assignment_algorithms.gating import (
-    mahalanobis_batch,
-    mahalanobis_distance,
-)
-from pytcl.assignment_algorithms.two_dimensional import (
-    assign2d,
-)
+from pytcl.assignment_algorithms.gating import mahalanobis_batch, mahalanobis_distance
+from pytcl.assignment_algorithms.two_dimensional import assign2d
 
 
 class AssociationResult(NamedTuple):

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/assignment_algorithms/gating.py

@@ -19,9 +19,9 @@ def _mahalanobis_distance_2d(
     S_inv: np.ndarray,
 ) -> float:
     """JIT-compiled Mahalanobis distance for 2D innovations."""
-    return innovation[0] * (
-        S_inv[0, 0] * innovation[0] + S_inv[0, 1] * innovation[1]
-    ) + innovation[1] * (S_inv[1, 0] * innovation[0] + S_inv[1, 1] * innovation[1])
+    return innovation[0] * (S_inv[0, 0] * innovation[0] + S_inv[0, 1] * innovation[1]) + innovation[
+        1
+    ] * (S_inv[1, 0] * innovation[0] + S_inv[1, 1] * innovation[1])
 
 
 @njit(cache=True, fastmath=True)

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/assignment_algorithms/jpda.py

@@ -94,9 +94,7 @@ def compute_measurement_likelihood(
         return 0.0
 
     mahal_sq = innovation @ np.linalg.solve(innovation_cov, innovation)
-    likelihood = (
-        detection_prob * np.exp(-0.5 * mahal_sq) / np.sqrt((2 * np.pi) ** m * det_S)
-    )
+    likelihood = detection_prob * np.exp(-0.5 * mahal_sq) / np.sqrt((2 * np.pi) ** m * det_S)
 
     return likelihood
 
@@ -210,9 +208,7 @@ def jpda_probabilities(
     if n_tracks <= 5 and n_meas <= 5:
         beta = _jpda_exact(likelihood_matrix, gated, detection_prob, clutter_density)
     else:
-        beta = _jpda_approximate(
-            likelihood_matrix, gated, detection_prob, clutter_density
-        )
+        beta = _jpda_approximate(likelihood_matrix, gated, detection_prob, clutter_density)
 
     return beta
 
@@ -258,9 +254,7 @@ def _jpda_exact(
             if gated[track_idx, meas_idx] and track_idx not in used_tracks:
                 current_assignment.append(track_idx)
                 used_tracks.add(track_idx)
-                yield from generate_hypotheses(
-                    meas_idx + 1, current_assignment, used_tracks
-                )
+                yield from generate_hypotheses(meas_idx + 1, current_assignment, used_tracks)
                 used_tracks.remove(track_idx)
                 current_assignment.pop()
 
@@ -298,9 +292,7 @@ def _jpda_exact(
         hypothesis_probs = [p / total_prob for p in hypothesis_probs]
 
     # Compute marginal association probabilities
-    for h_idx, (assignment, prob) in enumerate(
-        zip(hypothesis_assignments, hypothesis_probs)
-    ):
+    for h_idx, (assignment, prob) in enumerate(zip(hypothesis_assignments, hypothesis_probs)):
         detected_tracks = set()
         for j, track_idx in enumerate(assignment):
             if track_idx >= 0:

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/assignment_algorithms/two_dimensional/kbest.py

@@ -335,9 +335,7 @@ def _partition_solution(
 
     for i in range(start_idx, n_assigned):
         # Require assignments 0..i-1, forbid assignment i
-        new_required = required + [
-            (row_ind[j], col_ind[j]) for j in range(start_idx, i)
-        ]
+        new_required = required + [(row_ind[j], col_ind[j]) for j in range(start_idx, i)]
         new_forbidden = forbidden + [(row_ind[i], col_ind[i])]
 
         # Solve constrained problem

{nrl_tracker-0.21.5 → nrl_tracker-0.22.1}/pytcl/astronomical/__init__.py

@@ -2,7 +2,8 @@
 Astronomical calculations for target tracking.
 
 This module provides time system conversions, orbital mechanics,
-Lambert problem solvers, and reference frame transformations.
+Lambert problem solvers, reference frame transformations, and high-precision
+ephemerides for celestial bodies.
 
 Examples
 --------
@@ -18,8 +19,19 @@ Examples
 >>> r1 = np.array([5000, 10000, 2100])
 >>> r2 = np.array([-14600, 2500, 7000])
 >>> sol = lambert_universal(r1, r2, 3600)
+
+>>> # Query Sun position with high precision
+>>> from pytcl.astronomical import sun_position
+>>> r_sun, v_sun = sun_position(2451545.0)  # J2000.0
 """
 
+from pytcl.astronomical.ephemerides import (
+    DEEphemeris,
+    barycenter_position,
+    moon_position,
+    planet_position,
+    sun_position,
+)
 from pytcl.astronomical.lambert import (
     LambertSolution,
     bi_elliptic_transfer,
@@ -28,8 +40,10 @@ from pytcl.astronomical.lambert import (
     lambert_universal,
     minimum_energy_transfer,
 )
-from pytcl.astronomical.orbital_mechanics import (  # Constants; Types; Anomaly conversions; Element conversions; Propagation; Orbital quantities
-    GM_EARTH,
+from pytcl.astronomical.orbital_mechanics import (
+    GM_EARTH,  # Constants; Types; Anomaly conversions; Element conversions; Propagation; Orbital quantities
+)
+from pytcl.astronomical.orbital_mechanics import (
     GM_JUPITER,
     GM_MARS,
     GM_MOON,
@@ -62,8 +76,10 @@ from pytcl.astronomical.orbital_mechanics import (  # Constants; Types; Anomaly
     true_to_mean_anomaly,
     vis_viva,
 )
-from pytcl.astronomical.reference_frames import (  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
-    earth_rotation_angle,
+from pytcl.astronomical.reference_frames import (
+    earth_rotation_angle,  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
+)
+from pytcl.astronomical.reference_frames import (
     ecef_to_eci,
     eci_to_ecef,
     ecliptic_to_equatorial,
@@ -83,8 +99,25 @@ from pytcl.astronomical.reference_frames import (  # Time utilities; Precession;
     sidereal_rotation_matrix,
     true_obliquity,
 )
-from pytcl.astronomical.time_systems import (  # Julian dates; Time scales; Unix time; GPS week; Sidereal time; Leap seconds; Constants
-    JD_GPS_EPOCH,
+from pytcl.astronomical.relativity import (
+    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation; Shapiro delay; Precession; PN effects; Range corrections
+)
+from pytcl.astronomical.relativity import (
+    G_GRAV,
+    geodetic_precession,
+    gravitational_time_dilation,
+    lense_thirring_precession,
+    post_newtonian_acceleration,
+    proper_time_rate,
+    relativistic_range_correction,
+    schwarzschild_precession_per_orbit,
+    schwarzschild_radius,
+    shapiro_delay,
+)
+from pytcl.astronomical.time_systems import (
+    JD_GPS_EPOCH,  # Julian dates; Time scales; Unix time; GPS week; Sidereal time; Leap seconds; Constants
+)
+from pytcl.astronomical.time_systems import (
     JD_J2000,
     JD_UNIX_EPOCH,
     MJD_OFFSET,
@@ -218,4 +251,24 @@ __all__ = [
     # Reference frames - Ecliptic/equatorial
     "ecliptic_to_equatorial",
     "equatorial_to_ecliptic",
+    # Ephemerides - Classes
+    "DEEphemeris",
+    # Ephemerides - Functions
+    "sun_position",
+    "moon_position",
+    "planet_position",
+    "barycenter_position",
+    # Relativity - Constants
+    "C_LIGHT",
+    "G_GRAV",
+    # Relativity - Functions
+    "schwarzschild_radius",
+    "gravitational_time_dilation",
+    "proper_time_rate",
+    "shapiro_delay",
+    "schwarzschild_precession_per_orbit",
+    "post_newtonian_acceleration",
+    "geodetic_precession",
+    "lense_thirring_precession",
+    "relativistic_range_correction",
 ]