nrl-tracker 1.7.1__tar.gz → 1.7.3__tar.gz

{nrl_tracker-1.7.1/nrl_tracker.egg-info → nrl_tracker-1.7.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 1.7.1
+Version: 1.7.3
 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
@@ -63,7 +63,7 @@ Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualizati
 
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v1.7.1-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v1.7.2-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)
@@ -209,7 +209,46 @@ pytcl/
 
 - [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
 - [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
-- [Examples](examples/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
+
+## Examples & Tutorials
+
+The library includes 39 runnable code examples demonstrating all major features:
+
+### Examples (29 files in `/examples/`)
+
+Comprehensive demonstrations of library functionality:
+- **Tracking & Estimation**: Kalman filters, particle filters, smoothers
+- **Assignment**: Hungarian algorithm, k-best assignments, 3D assignment
+- **Coordinates**: Frame conversions, transformations, geodetic calculations
+- **Dynamics**: State models, motion models, dynamic systems
+- **Filtering**: Uncertainty visualization, multi-target tracking
+- **Astronomy**: Ephemerides, orbital mechanics, relativistic corrections
+- **Navigation**: INS/GNSS integration, geophysical modeling
+- **Signal Processing**: Detection, filtering, transforms
+- **Terrain & Atmosphere**: Elevation models, atmospheric properties
+
+**Status**: ✅ All 29 examples validated and passing (100% execution success)
+
+### Tutorials (10 modules in `/docs/tutorials/`)
+
+Interactive learning modules with visualizations:
+- Assignment algorithms and 3D assignment problems
+- Atmospheric and geophysical models
+- Dynamical systems and reference frames
+- Filtering and smoothing techniques
+- Sensor fusion and advanced filtering
+- Special functions and mathematical tools
+
+**Status**: ✅ All 10 tutorials validated and passing (100% execution success)
+
+## Documentation
+
+- [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
+- [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
 
 ## Comparison with Original MATLAB Library
 

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/README.md

@@ -1,6 +1,6 @@
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v1.7.1-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v1.7.2-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)
@@ -146,7 +146,46 @@ pytcl/
 
 - [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
 - [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
-- [Examples](examples/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
+
+## Examples & Tutorials
+
+The library includes 39 runnable code examples demonstrating all major features:
+
+### Examples (29 files in `/examples/`)
+
+Comprehensive demonstrations of library functionality:
+- **Tracking & Estimation**: Kalman filters, particle filters, smoothers
+- **Assignment**: Hungarian algorithm, k-best assignments, 3D assignment
+- **Coordinates**: Frame conversions, transformations, geodetic calculations
+- **Dynamics**: State models, motion models, dynamic systems
+- **Filtering**: Uncertainty visualization, multi-target tracking
+- **Astronomy**: Ephemerides, orbital mechanics, relativistic corrections
+- **Navigation**: INS/GNSS integration, geophysical modeling
+- **Signal Processing**: Detection, filtering, transforms
+- **Terrain & Atmosphere**: Elevation models, atmospheric properties
+
+**Status**: ✅ All 29 examples validated and passing (100% execution success)
+
+### Tutorials (10 modules in `/docs/tutorials/`)
+
+Interactive learning modules with visualizations:
+- Assignment algorithms and 3D assignment problems
+- Atmospheric and geophysical models
+- Dynamical systems and reference frames
+- Filtering and smoothing techniques
+- Sensor fusion and advanced filtering
+- Special functions and mathematical tools
+
+**Status**: ✅ All 10 tutorials validated and passing (100% execution success)
+
+## Documentation
+
+- [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
+- [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
 
 ## Comparison with Original MATLAB Library
 

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3/nrl_tracker.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 1.7.1
+Version: 1.7.3
 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
@@ -63,7 +63,7 @@ Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualizati
 
 # Tracker Component Library (Python)
 
-[![PyPI version](https://img.shields.io/badge/pypi-v1.7.1-blue.svg)](https://pypi.org/project/nrl-tracker/)
+[![PyPI version](https://img.shields.io/badge/pypi-v1.7.2-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)
@@ -209,7 +209,46 @@ pytcl/
 
 - [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
 - [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
-- [Examples](examples/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
+
+## Examples & Tutorials
+
+The library includes 39 runnable code examples demonstrating all major features:
+
+### Examples (29 files in `/examples/`)
+
+Comprehensive demonstrations of library functionality:
+- **Tracking & Estimation**: Kalman filters, particle filters, smoothers
+- **Assignment**: Hungarian algorithm, k-best assignments, 3D assignment
+- **Coordinates**: Frame conversions, transformations, geodetic calculations
+- **Dynamics**: State models, motion models, dynamic systems
+- **Filtering**: Uncertainty visualization, multi-target tracking
+- **Astronomy**: Ephemerides, orbital mechanics, relativistic corrections
+- **Navigation**: INS/GNSS integration, geophysical modeling
+- **Signal Processing**: Detection, filtering, transforms
+- **Terrain & Atmosphere**: Elevation models, atmospheric properties
+
+**Status**: ✅ All 29 examples validated and passing (100% execution success)
+
+### Tutorials (10 modules in `/docs/tutorials/`)
+
+Interactive learning modules with visualizations:
+- Assignment algorithms and 3D assignment problems
+- Atmospheric and geophysical models
+- Dynamical systems and reference frames
+- Filtering and smoothing techniques
+- Sensor fusion and advanced filtering
+- Special functions and mathematical tools
+
+**Status**: ✅ All 10 tutorials validated and passing (100% execution success)
+
+## Documentation
+
+- [API Reference](https://pytcl.readthedocs.io/en/latest/api/)
+- [User Guides](https://pytcl.readthedocs.io/en/latest/user_guide/)
+- [Examples](examples/) - 29 validated example scripts
+- [Tutorials](docs/tutorials/) - 10 interactive tutorial modules
 
 ## Comparison with Original MATLAB Library
 

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "nrl-tracker"
-version = "1.7.1"
+version = "1.7.3"
 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.7.1 → nrl_tracker-1.7.3}/pytcl/__init__.py

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

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/pytcl/astronomical/__init__.py

@@ -77,10 +77,10 @@ from pytcl.astronomical.orbital_mechanics import (
     vis_viva,
 )
 from pytcl.astronomical.reference_frames import (
-    earth_rotation_angle,  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
+    earth_rotation_angle,  # Time utilities; Precession; Nutation
 )
+from pytcl.astronomical.reference_frames import ecef_to_eci  # Time utilities
 from pytcl.astronomical.reference_frames import (
-    ecef_to_eci,
     eci_to_ecef,
     ecliptic_to_equatorial,
     equation_of_equinoxes,
@@ -116,7 +116,7 @@ from pytcl.astronomical.reference_frames import (
     true_obliquity,
 )
 from pytcl.astronomical.relativity import (
-    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation; Shapiro delay; Precession; PN effects; Range corrections
+    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation
 )
 from pytcl.astronomical.relativity import (
     G_GRAV,

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/pytcl/coordinate_systems/projections/__init__.py

@@ -27,7 +27,7 @@ Examples
 """
 
 from pytcl.coordinate_systems.projections.projections import (
-    WGS84_A,  # Constants; Result types; Azimuthal Equidistant; UTM; Lambert Conformal Conic; Mercator; Stereographic; Transverse Mercator
+    WGS84_A,  # Constants; Result types
 )
 from pytcl.coordinate_systems.projections.projections import (
     WGS84_B,

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/conftest.py

@@ -82,7 +82,9 @@ class NumpyTestCase:
 
     def assert_shape(self, arr, expected_shape):
         """Assert array has expected shape."""
-        assert arr.shape == expected_shape, f"Expected shape {expected_shape}, got {arr.shape}"
+        assert (
+            arr.shape == expected_shape
+        ), f"Expected shape {expected_shape}, got {arr.shape}"
 
     def assert_symmetric(self, arr, rtol=1e-10):
         """Assert matrix is symmetric."""

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_assignment_algorithms.py

@@ -3,7 +3,7 @@
 import numpy as np
 from numpy.testing import assert_allclose
 
-from pytcl.assignment_algorithms import (  # 2D Assignment; K-Best 2D Assignment; 3D Assignment; Gating; Data Association
+from pytcl.assignment_algorithms import (  # 2D Assignment
     Assignment3DResult,
     AssignmentResult,
     AssociationResult,

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_astronomical.py

@@ -150,7 +150,9 @@ class TestOrbitalElementConversions:
 
     def test_roundtrip_elements_state(self):
         """Test elements -> state -> elements roundtrip."""
-        elements_orig = OrbitalElements(a=8000, e=0.2, i=0.5, raan=0.3, omega=0.8, nu=1.0)
+        elements_orig = OrbitalElements(
+            a=8000, e=0.2, i=0.5, raan=0.3, omega=0.8, nu=1.0
+        )
         state = orbital_elements_to_state(elements_orig)
         elements_back = state_to_orbital_elements(state)
 

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_coordinate_systems.py

@@ -403,7 +403,9 @@ class TestSlerp:
 
         # Should be 45 deg rotation about z
         expected_angle = angle / 2
-        expected = np.array([np.cos(expected_angle / 2), 0, 0, np.sin(expected_angle / 2)])
+        expected = np.array(
+            [np.cos(expected_angle / 2), 0, 0, np.sin(expected_angle / 2)]
+        )
         np.testing.assert_allclose(result, expected, atol=1e-10)
 
 
@@ -1315,4 +1317,6 @@ class TestSEZConversions:
         # For northward target, the magnitude should be large and east component small
         total_distance = np.linalg.norm(sez)
         assert total_distance > 100000, "Distance should be significant"
-        assert abs(sez[1]) < abs(sez[0]), "East component should be smaller than meridional"
+        assert abs(sez[1]) < abs(
+            sez[0]
+        ), "East component should be smaller than meridional"

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_coverage_boost.py

@@ -438,7 +438,9 @@ class TestCoordinatedTurnProcessNoise:
         T = 0.1
         sigma_a = 1.0
         sigma_omega = 0.01
-        Q = q_coord_turn_2d(T, sigma_a, sigma_omega, state_type="position_velocity_omega")
+        Q = q_coord_turn_2d(
+            T, sigma_a, sigma_omega, state_type="position_velocity_omega"
+        )
         assert Q.shape == (5, 5)
         np.testing.assert_allclose(Q, Q.T)
 

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_coverage_boost_2.py

@@ -121,7 +121,9 @@ class TestParticleFilters:
         def likelihood_func(z, x):
             return gaussian_likelihood(z, x, R)
 
-        new_weights, log_lik = bootstrap_pf_update(particles, weights, z, likelihood_func)
+        new_weights, log_lik = bootstrap_pf_update(
+            particles, weights, z, likelihood_func
+        )
 
         assert new_weights.shape == weights.shape
         assert np.isclose(np.sum(new_weights), 1.0)
@@ -325,7 +327,9 @@ class TestParticleFilters:
 
         mean = particle_mean(particles, weights)
 
-        expected = 0.5 * np.array([0, 0]) + 0.3 * np.array([1, 1]) + 0.2 * np.array([2, 2])
+        expected = (
+            0.5 * np.array([0, 0]) + 0.3 * np.array([1, 1]) + 0.2 * np.array([2, 2])
+        )
         assert np.allclose(mean, expected)
 
     def test_particle_covariance(self):
@@ -991,7 +995,9 @@ class TestNavigationCoverage:
         tdoa12 = 0.0
         tdoa13 = 0.0
 
-        result1, result2 = great_circle_tdoa_loc(lat1, lon1, lat2, lon2, lat3, lon3, tdoa12, tdoa13)
+        result1, result2 = great_circle_tdoa_loc(
+            lat1, lon1, lat2, lon2, lat3, lon3, tdoa12, tdoa13
+        )
 
         # Should return at least one solution
         assert result1 is not None or result2 is not None

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_egm.py

@@ -98,7 +98,9 @@ class TestClenshawSumOrderDerivative:
         cos_theta = 0.5
         sin_theta = np.sqrt(1 - cos_theta**2)
 
-        result = clenshaw_sum_order_derivative(0, cos_theta, sin_theta, coef.C, coef.S, 10)
+        result = clenshaw_sum_order_derivative(
+            0, cos_theta, sin_theta, coef.C, coef.S, 10
+        )
 
         assert isinstance(result, tuple)
         assert len(result) == 4
@@ -140,8 +142,12 @@ class TestClenshawPotential:
         lat = np.radians(30)
         lon = 0
 
-        V_close = clenshaw_potential(lat, lon, coef.R, coef.C, coef.S, coef.R, coef.GM, 10)
-        V_far = clenshaw_potential(lat, lon, coef.R * 1.1, coef.C, coef.S, coef.R, coef.GM, 10)
+        V_close = clenshaw_potential(
+            lat, lon, coef.R, coef.C, coef.S, coef.R, coef.GM, 10
+        )
+        V_far = clenshaw_potential(
+            lat, lon, coef.R * 1.1, coef.C, coef.S, coef.R, coef.GM, 10
+        )
 
         # Potential magnitude should decrease with distance
         assert abs(V_far) < abs(V_close)
@@ -152,8 +158,12 @@ class TestClenshawPotential:
         lon = 0
         r = coef.R
 
-        V_north = clenshaw_potential(np.radians(45), lon, r, coef.C, coef.S, coef.R, coef.GM, 10)
-        V_south = clenshaw_potential(np.radians(-45), lon, r, coef.C, coef.S, coef.R, coef.GM, 10)
+        V_north = clenshaw_potential(
+            np.radians(45), lon, r, coef.C, coef.S, coef.R, coef.GM, 10
+        )
+        V_south = clenshaw_potential(
+            np.radians(-45), lon, r, coef.C, coef.S, coef.R, coef.GM, 10
+        )
 
         # For Earth's gravity field (with small asymmetries), should be close
         # The potential is dominated by the central term, so relative difference is small
@@ -185,7 +195,9 @@ class TestClenshawGravity:
         lon = np.radians(0)
         r = coef.R
 
-        g_r, g_lat, g_lon = clenshaw_gravity(lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10)
+        g_r, g_lat, g_lon = clenshaw_gravity(
+            lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10
+        )
 
         # Full gravity includes the central term GM/r^2 which is about 9.8 m/s^2
         # The sign depends on convention - we compute -dV/dr
@@ -245,7 +257,9 @@ class TestAssociatedLegendreScaled:
             for m in range(n + 1):
                 # exp is indexed by degree only (1D array)
                 P_reconstructed = P_scaled[n, m] * (10.0 ** exp[n])
-                assert_allclose(P_reconstructed, P_unscaled[n, m], rtol=1e-8, atol=1e-12)
+                assert_allclose(
+                    P_reconstructed, P_unscaled[n, m], rtol=1e-8, atol=1e-12
+                )
 
     def test_no_overflow_high_degree(self):
         """High degree computation doesn't overflow."""
@@ -306,7 +320,9 @@ class TestGeoidHeight:
         # Test several locations
         for lat in [0, 45, -45, 90]:
             for lon in [0, 90, -90, 180]:
-                N = geoid_height(np.radians(lat), np.radians(lon), coefficients=coef, n_max=36)
+                N = geoid_height(
+                    np.radians(lat), np.radians(lon), coefficients=coef, n_max=36
+                )
                 assert not np.isnan(N)
                 assert not np.isinf(N)
 
@@ -332,7 +348,9 @@ class TestGeoidHeights:
         lon = np.radians(-75)
 
         N_single = geoid_height(lat, lon, coefficients=coef, n_max=10)
-        N_batch = geoid_heights(np.array([lat]), np.array([lon]), coefficients=coef, n_max=10)
+        N_batch = geoid_heights(
+            np.array([lat]), np.array([lon]), coefficients=coef, n_max=10
+        )
 
         assert_allclose(N_batch[0], N_single)
 
@@ -362,7 +380,9 @@ class TestGravityDisturbance:
             np.radians(30), np.radians(60), h=0.0, coefficients=coef, n_max=10
         )
 
-        expected_mag = np.sqrt(result.delta_g_r**2 + result.delta_g_lat**2 + result.delta_g_lon**2)
+        expected_mag = np.sqrt(
+            result.delta_g_r**2 + result.delta_g_lat**2 + result.delta_g_lon**2
+        )
         assert_allclose(result.magnitude, expected_mag)
 
     def test_disturbance_decreases_with_altitude(self):
@@ -372,7 +392,9 @@ class TestGravityDisturbance:
         lon = np.radians(0)
 
         result_0 = gravity_disturbance(lat, lon, h=0.0, coefficients=coef, n_max=20)
-        result_10km = gravity_disturbance(lat, lon, h=10000.0, coefficients=coef, n_max=20)
+        result_10km = gravity_disturbance(
+            lat, lon, h=10000.0, coefficients=coef, n_max=20
+        )
 
         # Disturbance magnitude should decrease with altitude
         assert result_10km.magnitude < result_0.magnitude
@@ -385,7 +407,9 @@ class TestGravityAnomaly:
         """gravity_anomaly returns a float."""
         coef = create_test_coefficients(n_max=10)
 
-        anom = gravity_anomaly(np.radians(45), np.radians(-75), h=0.0, coefficients=coef, n_max=10)
+        anom = gravity_anomaly(
+            np.radians(45), np.radians(-75), h=0.0, coefficients=coef, n_max=10
+        )
 
         assert isinstance(anom, float)
         assert not np.isnan(anom)
@@ -394,7 +418,9 @@ class TestGravityAnomaly:
         """Gravity anomaly computation produces finite values."""
         coef = create_test_coefficients(n_max=36)
 
-        anom = gravity_anomaly(np.radians(45), np.radians(0), h=0.0, coefficients=coef, n_max=36)
+        anom = gravity_anomaly(
+            np.radians(45), np.radians(0), h=0.0, coefficients=coef, n_max=36
+        )
 
         # Note: With test coefficients (including full J2), values are large
         # because J2 represents the equatorial bulge which is part of the
@@ -423,7 +449,9 @@ class TestDeflectionOfVertical:
         """Deflection values are finite."""
         coef = create_test_coefficients(n_max=36)
 
-        xi, eta = deflection_of_vertical(np.radians(45), np.radians(0), coefficients=coef, n_max=36)
+        xi, eta = deflection_of_vertical(
+            np.radians(45), np.radians(0), coefficients=coef, n_max=36
+        )
 
         # Note: Test coefficients include full J2, so deflections are larger
         # than real-world values which are typically < 60 arcseconds.
@@ -493,12 +521,16 @@ class TestClenshawVsNaive:
         r = coef.R
 
         # Compute potential using Clenshaw
-        V_clenshaw = clenshaw_potential(lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10)
+        V_clenshaw = clenshaw_potential(
+            lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10
+        )
 
         # Compute using naive spherical harmonic sum
         # spherical_harmonic_sum takes (lat, lon, r, C, S, R, GM, n_max)
         # and returns (V, dV_r, dV_lat)
-        V_naive, _, _ = spherical_harmonic_sum(lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10)
+        V_naive, _, _ = spherical_harmonic_sum(
+            lat, lon, r, coef.C, coef.S, coef.R, coef.GM, 10
+        )
 
         # Should match to reasonable precision
         # (Clenshaw and naive may have different numerical properties)
@@ -546,8 +578,12 @@ class TestEdgeCases:
         """Negative longitude works correctly."""
         coef = create_test_coefficients(n_max=20)
 
-        N_pos = geoid_height(np.radians(45), np.radians(90), coefficients=coef, n_max=20)
-        N_neg = geoid_height(np.radians(45), np.radians(-270), coefficients=coef, n_max=20)
+        N_pos = geoid_height(
+            np.radians(45), np.radians(90), coefficients=coef, n_max=20
+        )
+        N_neg = geoid_height(
+            np.radians(45), np.radians(-270), coefficients=coef, n_max=20
+        )
 
         # 90° and -270° are the same longitude
         assert_allclose(N_pos, N_neg, rtol=1e-10)

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_emm.py

@@ -130,31 +130,41 @@ class TestEMMFunction:
 
     def test_total_intensity_reasonable(self, test_coefficients):
         """Total field intensity should be in expected range."""
-        result = emm(np.radians(45), np.radians(0), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(45), np.radians(0), 0, 2020.0, coefficients=test_coefficients
+        )
         # Field intensity typically 25,000-65,000 nT at mid-latitudes
         # but can exceed this at high latitudes, allow up to 100,000 nT
         assert 20000 < result.F < 100000
 
     def test_horizontal_intensity_formula(self, test_coefficients):
         """H = sqrt(X^2 + Y^2)."""
-        result = emm(np.radians(40), np.radians(-75), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(40), np.radians(-75), 0, 2020.0, coefficients=test_coefficients
+        )
         H_calc = np.sqrt(result.X**2 + result.Y**2)
         assert_allclose(result.H, H_calc, rtol=1e-10)
 
     def test_total_intensity_formula(self, test_coefficients):
         """F = sqrt(H^2 + Z^2)."""
-        result = emm(np.radians(40), np.radians(-75), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(40), np.radians(-75), 0, 2020.0, coefficients=test_coefficients
+        )
         F_calc = np.sqrt(result.H**2 + result.Z**2)
         assert_allclose(result.F, F_calc, rtol=1e-10)
 
     def test_declination_range(self, test_coefficients):
         """Declination should be within -180 to 180 degrees."""
-        result = emm(np.radians(45), np.radians(-75), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(45), np.radians(-75), 0, 2020.0, coefficients=test_coefficients
+        )
         assert -np.pi <= result.D <= np.pi
 
     def test_inclination_range(self, test_coefficients):
         """Inclination should be within -90 to 90 degrees."""
-        result = emm(np.radians(45), np.radians(-75), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(45), np.radians(-75), 0, 2020.0, coefficients=test_coefficients
+        )
         assert -np.pi / 2 <= result.I <= np.pi / 2
 
 
@@ -219,7 +229,9 @@ class TestWMMHR:
 
     def test_wmmhr_returns_result(self, test_coefficients):
         """WMMHR returns MagneticResult."""
-        result = wmmhr(np.radians(45), np.radians(-75), 0, 2025.0, coefficients=test_coefficients)
+        result = wmmhr(
+            np.radians(45), np.radians(-75), 0, 2025.0, coefficients=test_coefficients
+        )
         assert hasattr(result, "F")
         assert hasattr(result, "D")
         assert hasattr(result, "I")
@@ -259,7 +271,9 @@ class TestConvenienceFunctions:
 
     def test_emm_declination(self, test_coefficients):
         """emm_declination returns correct value."""
-        result = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
+        )
         D = emm_declination(
             np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
         )
@@ -267,7 +281,9 @@ class TestConvenienceFunctions:
 
     def test_emm_inclination(self, test_coefficients):
         """emm_inclination returns correct value."""
-        result = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
+        )
         incl = emm_inclination(
             np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
         )
@@ -275,7 +291,9 @@ class TestConvenienceFunctions:
 
     def test_emm_intensity(self, test_coefficients):
         """emm_intensity returns correct value."""
-        result = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients)
+        result = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
+        )
         F = emm_intensity(
             np.radians(40), np.radians(-105), 0, 2020.0, coefficients=test_coefficients
         )
@@ -313,8 +331,12 @@ class TestSecularVariation:
 
     def test_field_changes_with_time(self, test_coefficients):
         """Field should change slightly between years."""
-        result_2020 = emm(np.radians(45), np.radians(0), 0, 2020.0, coefficients=test_coefficients)
-        result_2022 = emm(np.radians(45), np.radians(0), 0, 2022.0, coefficients=test_coefficients)
+        result_2020 = emm(
+            np.radians(45), np.radians(0), 0, 2020.0, coefficients=test_coefficients
+        )
+        result_2022 = emm(
+            np.radians(45), np.radians(0), 0, 2022.0, coefficients=test_coefficients
+        )
         # Small but non-zero change expected
         # Secular variation is ~50-100 nT/year at mid-latitudes
         assert result_2020.F != result_2022.F
@@ -389,7 +411,9 @@ class TestHighDegreeEvaluation:
     def test_higher_degree_coefficients(self):
         """Test with higher degree coefficients (n_max=50)."""
         coef = create_emm_test_coefficients(n_max=50)
-        result = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=50)
+        result = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=50
+        )
         assert not np.isnan(result.F)
         assert result.F > 0
 
@@ -398,8 +422,12 @@ class TestHighDegreeEvaluation:
         coef = create_emm_test_coefficients(n_max=50)
 
         # Evaluate at different n_max values
-        result_12 = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=12)
-        result_50 = emm(np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=50)
+        result_12 = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=12
+        )
+        result_50 = emm(
+            np.radians(40), np.radians(-105), 0, 2020.0, coefficients=coef, n_max=50
+        )
 
         # Both should be valid
         assert result_12.F > 0

{nrl_tracker-1.7.1 → nrl_tracker-1.7.3}/tests/test_ephemerides.py

@@ -155,7 +155,9 @@ class TestMoonPosition:
     def test_moon_position_earth_centered(self):
         """Test Moon position relative to Earth."""
         jd = 2451545.0
-        r_earth_centered, v_earth_centered = self.eph.moon_position(jd, frame="earth_centered")
+        r_earth_centered, v_earth_centered = self.eph.moon_position(
+            jd, frame="earth_centered"
+        )
 
         # Moon is about 385,000 km = 0.00257 AU from Earth
         distance = np.linalg.norm(r_earth_centered)
@@ -163,7 +165,9 @@ class TestMoonPosition:
         distance_km = distance * au_to_km
 
         # Should be roughly 380,000-390,000 km
-        assert 370000 < distance_km < 400000, f"Moon distance {distance_km:.0f} km is unexpected"
+        assert (
+            370000 < distance_km < 400000
+        ), f"Moon distance {distance_km:.0f} km is unexpected"
 
     def test_moon_position_frames_consistency(self):
         """Test that Moon positions are consistent across frames."""