nrl-tracker 1.1.3__tar.gz → 1.3.0__tar.gz

{nrl_tracker-1.1.3/nrl_tracker.egg-info → nrl_tracker-1.3.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 1.1.3
+Version: 1.3.0
 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

{nrl_tracker-1.1.3 → nrl_tracker-1.3.0/nrl_tracker.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 1.1.3
+Version: 1.3.0
 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

{nrl_tracker-1.1.3 → nrl_tracker-1.3.0}/nrl_tracker.egg-info/SOURCES.txt

@@ -27,6 +27,7 @@ pytcl/astronomical/reference_frames.py
 pytcl/astronomical/relativity.py
 pytcl/astronomical/time_systems.py
 pytcl/atmosphere/__init__.py
+pytcl/atmosphere/ionosphere.py
 pytcl/atmosphere/models.py
 pytcl/clustering/__init__.py
 pytcl/clustering/dbscan.py
@@ -34,6 +35,7 @@ pytcl/clustering/gaussian_mixture.py
 pytcl/clustering/hierarchical.py
 pytcl/clustering/kmeans.py
 pytcl/containers/__init__.py
+pytcl/containers/base.py
 pytcl/containers/cluster_set.py
 pytcl/containers/covertree.py
 pytcl/containers/kd_tree.py
@@ -64,6 +66,8 @@ pytcl/dynamic_estimation/kalman/__init__.py
 pytcl/dynamic_estimation/kalman/extended.py
 pytcl/dynamic_estimation/kalman/linear.py
 pytcl/dynamic_estimation/kalman/square_root.py
+pytcl/dynamic_estimation/kalman/sr_ukf.py
+pytcl/dynamic_estimation/kalman/ud_filter.py
 pytcl/dynamic_estimation/kalman/unscented.py
 pytcl/dynamic_estimation/measurement_update/__init__.py
 pytcl/dynamic_estimation/particle_filters/__init__.py
@@ -184,6 +188,7 @@ tests/test_relativity.py
 tests/test_rhumb.py
 tests/test_signal_processing.py
 tests/test_smoothers.py
+tests/test_spatial_containers_parametrized.py
 tests/test_spatial_structures.py
 tests/test_special_functions_phase12.py
 tests/test_static_estimation.py
@@ -195,4 +200,5 @@ tests/test_tracking_containers.py
 tests/test_transforms.py
 tests/test_v030_comprehensive.py
 tests/test_v030_features.py
+tests/test_validation.py
 tests/unit/test_core.py

{nrl_tracker-1.1.3 → nrl_tracker-1.3.0}/pyproject.toml

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

{nrl_tracker-1.1.3 → nrl_tracker-1.3.0}/pytcl/__init__.py

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

{nrl_tracker-1.1.3 → nrl_tracker-1.3.0}/pytcl/astronomical/reference_frames.py

@@ -21,6 +21,8 @@ References
        A&A, 2003.
 """
 
+import logging
+from functools import lru_cache
 from typing import Tuple
 
 import numpy as np
@@ -28,6 +30,22 @@ from numpy.typing import NDArray
 
 from pytcl.astronomical.time_systems import JD_J2000
 
+# Module logger
+_logger = logging.getLogger("pytcl.astronomical.reference_frames")
+
+# Cache configuration
+_CACHE_JD_DECIMALS = 6  # ~86ms precision for JD quantization
+_CACHE_MAXSIZE = 128  # Max cached epochs
+
+
+def _quantize_jd(jd: float) -> float:
+    """Quantize Julian date for cache key compatibility.
+
+    Rounds to _CACHE_JD_DECIMALS decimal places (~86ms precision).
+    This enables cache hits for nearly identical epochs.
+    """
+    return round(jd, _CACHE_JD_DECIMALS)
+
 
 def julian_centuries_j2000(jd: float) -> float:
     """
@@ -78,6 +96,37 @@ def precession_angles_iau76(T: float) -> Tuple[float, float, float]:
     )
 
 
+@lru_cache(maxsize=_CACHE_MAXSIZE)
+def _precession_matrix_cached(jd_quantized: float) -> tuple:
+    """Cached precession matrix computation (internal).
+
+    Returns tuple of tuples for hashability.
+    """
+    T = julian_centuries_j2000(jd_quantized)
+    zeta, theta, z = precession_angles_iau76(T)
+
+    cos_zeta = np.cos(zeta)
+    sin_zeta = np.sin(zeta)
+    cos_theta = np.cos(theta)
+    sin_theta = np.sin(theta)
+    cos_z = np.cos(z)
+    sin_z = np.sin(z)
+
+    return (
+        (
+            cos_zeta * cos_theta * cos_z - sin_zeta * sin_z,
+            -sin_zeta * cos_theta * cos_z - cos_zeta * sin_z,
+            -sin_theta * cos_z,
+        ),
+        (
+            cos_zeta * cos_theta * sin_z + sin_zeta * cos_z,
+            -sin_zeta * cos_theta * sin_z + cos_zeta * cos_z,
+            -sin_theta * sin_z,
+        ),
+        (cos_zeta * sin_theta, -sin_zeta * sin_theta, cos_theta),
+    )
+
+
 def precession_matrix_iau76(jd: float) -> NDArray[np.floating]:
     """
     Compute IAU 1976 precession matrix from J2000 to date.
@@ -92,34 +141,15 @@ def precession_matrix_iau76(jd: float) -> NDArray[np.floating]:
     P : ndarray
         Precession rotation matrix (3x3).
         Transforms from J2000 (GCRF) to mean of date.
-    """
-    T = julian_centuries_j2000(jd)
-    zeta, theta, z = precession_angles_iau76(T)
 
-    cos_zeta = np.cos(zeta)
-    sin_zeta = np.sin(zeta)
-    cos_theta = np.cos(theta)
-    sin_theta = np.sin(theta)
-    cos_z = np.cos(z)
-    sin_z = np.sin(z)
-
-    P = np.array(
-        [
-            [
-                cos_zeta * cos_theta * cos_z - sin_zeta * sin_z,
-                -sin_zeta * cos_theta * cos_z - cos_zeta * sin_z,
-                -sin_theta * cos_z,
-            ],
-            [
-                cos_zeta * cos_theta * sin_z + sin_zeta * cos_z,
-                -sin_zeta * cos_theta * sin_z + cos_zeta * cos_z,
-                -sin_theta * sin_z,
-            ],
-            [cos_zeta * sin_theta, -sin_zeta * sin_theta, cos_theta],
-        ]
-    )
-
-    return P
+    Notes
+    -----
+    Results are cached for repeated queries at the same epoch.
+    Cache key is quantized to ~86ms precision.
+    """
+    jd_q = _quantize_jd(jd)
+    cached = _precession_matrix_cached(jd_q)
+    return np.array(cached)
 
 
 def nutation_angles_iau80(jd: float) -> Tuple[float, float]:
@@ -203,6 +233,38 @@ def mean_obliquity_iau80(jd: float) -> float:
     return eps0_arcsec * np.pi / (180 * 3600)
 
 
+@lru_cache(maxsize=_CACHE_MAXSIZE)
+def _nutation_matrix_cached(jd_quantized: float) -> tuple:
+    """Cached nutation matrix computation (internal).
+
+    Returns tuple of tuples for hashability.
+    """
+    dpsi, deps = nutation_angles_iau80(jd_quantized)
+    eps0 = mean_obliquity_iau80(jd_quantized)
+    eps = eps0 + deps
+
+    cos_eps0 = np.cos(eps0)
+    sin_eps0 = np.sin(eps0)
+    cos_eps = np.cos(eps)
+    sin_eps = np.sin(eps)
+    cos_dpsi = np.cos(dpsi)
+    sin_dpsi = np.sin(dpsi)
+
+    return (
+        (cos_dpsi, -sin_dpsi * cos_eps0, -sin_dpsi * sin_eps0),
+        (
+            sin_dpsi * cos_eps,
+            cos_dpsi * cos_eps0 * cos_eps + sin_eps0 * sin_eps,
+            cos_dpsi * sin_eps0 * cos_eps - cos_eps0 * sin_eps,
+        ),
+        (
+            sin_dpsi * sin_eps,
+            cos_dpsi * cos_eps0 * sin_eps - sin_eps0 * cos_eps,
+            cos_dpsi * sin_eps0 * sin_eps + cos_eps0 * cos_eps,
+        ),
+    )
+
+
 def nutation_matrix(jd: float) -> NDArray[np.floating]:
     """
     Compute nutation matrix.
@@ -217,35 +279,15 @@ def nutation_matrix(jd: float) -> NDArray[np.floating]:
     N : ndarray
         Nutation rotation matrix (3x3).
         Transforms from mean of date to true of date.
-    """
-    dpsi, deps = nutation_angles_iau80(jd)
-    eps0 = mean_obliquity_iau80(jd)
-    eps = eps0 + deps
 
-    cos_eps0 = np.cos(eps0)
-    sin_eps0 = np.sin(eps0)
-    cos_eps = np.cos(eps)
-    sin_eps = np.sin(eps)
-    cos_dpsi = np.cos(dpsi)
-    sin_dpsi = np.sin(dpsi)
-
-    N = np.array(
-        [
-            [cos_dpsi, -sin_dpsi * cos_eps0, -sin_dpsi * sin_eps0],
-            [
-                sin_dpsi * cos_eps,
-                cos_dpsi * cos_eps0 * cos_eps + sin_eps0 * sin_eps,
-                cos_dpsi * sin_eps0 * cos_eps - cos_eps0 * sin_eps,
-            ],
-            [
-                sin_dpsi * sin_eps,
-                cos_dpsi * cos_eps0 * sin_eps - sin_eps0 * cos_eps,
-                cos_dpsi * sin_eps0 * sin_eps + cos_eps0 * cos_eps,
-            ],
-        ]
-    )
-
-    return N
+    Notes
+    -----
+    Results are cached for repeated queries at the same epoch.
+    Cache key is quantized to ~86ms precision.
+    """
+    jd_q = _quantize_jd(jd)
+    cached = _nutation_matrix_cached(jd_q)
+    return np.array(cached)
 
 
 def earth_rotation_angle(jd_ut1: float) -> float:
@@ -647,6 +689,33 @@ def equatorial_to_ecliptic(
     return R @ r_eq
 
 
+def clear_transformation_cache() -> None:
+    """Clear cached transformation matrices.
+
+    Call this function to clear all cached precession and nutation
+    matrices. Useful when memory is constrained or after processing
+    a batch of observations at different epochs.
+    """
+    _precession_matrix_cached.cache_clear()
+    _nutation_matrix_cached.cache_clear()
+    _logger.debug("Transformation matrix cache cleared")
+
+
+def get_cache_info() -> dict:
+    """Get cache statistics for transformation matrices.
+
+    Returns
+    -------
+    dict
+        Dictionary with 'precession' and 'nutation' keys, each containing
+        CacheInfo namedtuple with hits, misses, maxsize, currsize.
+    """
+    return {
+        "precession": _precession_matrix_cached.cache_info(),
+        "nutation": _nutation_matrix_cached.cache_info(),
+    }
+
+
 __all__ = [
     # Time utilities
     "julian_centuries_j2000",
@@ -674,4 +743,7 @@ __all__ = [
     # Ecliptic/equatorial
     "ecliptic_to_equatorial",
     "equatorial_to_ecliptic",
+    # Cache management
+    "clear_transformation_cache",
+    "get_cache_info",
 ]

nrl_tracker-1.3.0/pytcl/atmosphere/__init__.py

@@ -0,0 +1,68 @@
+"""
+Atmospheric models for tracking applications.
+
+This module provides standard atmosphere models used for computing
+temperature, pressure, density, and other properties at various altitudes.
+
+Submodules
+----------
+models : Standard atmosphere models (US76, ISA)
+ionosphere : Ionospheric models for GPS/GNSS corrections
+"""
+
+from pytcl.atmosphere.ionosphere import (
+    DEFAULT_KLOBUCHAR,
+    F_L1,
+    F_L2,
+    IonosphereState,
+    KlobucharCoefficients,
+    dual_frequency_tec,
+    ionospheric_delay_from_tec,
+    klobuchar_delay,
+    magnetic_latitude,
+    scintillation_index,
+    simple_iri,
+)
+from pytcl.atmosphere.models import G0  # Constants
+from pytcl.atmosphere.models import (
+    GAMMA,
+    P0,
+    RHO0,
+    T0,
+    AtmosphereState,
+    R,
+    altitude_from_pressure,
+    isa_atmosphere,
+    mach_number,
+    true_airspeed_from_mach,
+    us_standard_atmosphere_1976,
+)
+
+__all__ = [
+    # Atmosphere state and models
+    "AtmosphereState",
+    "us_standard_atmosphere_1976",
+    "isa_atmosphere",
+    "altitude_from_pressure",
+    "mach_number",
+    "true_airspeed_from_mach",
+    # Atmosphere constants
+    "T0",
+    "P0",
+    "RHO0",
+    "G0",
+    "R",
+    "GAMMA",
+    # Ionosphere
+    "IonosphereState",
+    "KlobucharCoefficients",
+    "DEFAULT_KLOBUCHAR",
+    "klobuchar_delay",
+    "dual_frequency_tec",
+    "ionospheric_delay_from_tec",
+    "simple_iri",
+    "magnetic_latitude",
+    "scintillation_index",
+    "F_L1",
+    "F_L2",
+]