satkit 0.9.0__tar.gz → 0.9.2__tar.gz

{satkit-0.9.0 → satkit-0.9.2}/.github/workflows/build.yml

@@ -64,7 +64,7 @@ jobs:
         run: |
           python -m pip install requests
           python python/test/download_data.py ${{matrix.data_path}}
-        
+
 
       # Try to recover satkit test vectors from cache
       - name: Cache Satkit Test Vectors
@@ -79,7 +79,7 @@ jobs:
       - name: Download Satkit Test Vectors
         if: steps.cache-satkit-testvecs.outputs.cache-hit != 'true'
         run: |
-          python -m pip install requests 
+          python -m pip install requests
           python python/test/download_testvecs.py ${{matrix.testvec_path}}
 
       - name: Cargo Cache
@@ -91,16 +91,16 @@ jobs:
             ~/.cargo/registry/index/
             ~/.cargo/registry/cache/
             ~/.cargo/git/db/
-            target/            
+            target/
           key: ${{ runner.os }}-cargo-${{ hashFiles('**/Cargo.lock') }}
           restore-keys: ${{ runner.os }}-cargo-
 
       - name: Build
         run: cargo build --release
-      
+
       - name: Test
-        run: cargo test
+        run: cargo test --features chrono
 
       - name: Doc
-        run: cargo doc --no-deps 
+        run: cargo doc --no-deps
 

satkit-0.9.2/.github/workflows/cargo-publish.yml

@@ -0,0 +1,35 @@
+name: Publish Rust crate
+
+on:
+  push:
+    tags:
+      - "rust_*"
+
+permissions:
+  contents: read
+
+jobs:
+  publish:
+    runs-on: ubuntu-latest
+    environment:
+      name: crates-io
+    steps:
+      - name: Checkout
+        uses: actions/checkout@v4
+
+      - name: Set up Rust
+        uses: dtolnay/rust-toolchain@stable
+
+      - name: Verify tag matches Cargo.toml
+        run: |
+          TAG_VERSION="${GITHUB_REF_NAME#rust_}"
+          CARGO_VERSION=$(cargo metadata --format-version 1 --no-deps | jq -r '.packages[0].version')
+          if [ "$TAG_VERSION" != "$CARGO_VERSION" ]; then
+            echo "Tag version ($TAG_VERSION) does not match Cargo.toml version ($CARGO_VERSION)."
+            exit 1
+          fi
+
+      - name: Publish to crates.io
+        env:
+          CARGO_REGISTRY_TOKEN: ${{ secrets.CARGO_REGISTRY_TOKEN }}
+        run: cargo publish

{satkit-0.9.0 → satkit-0.9.2}/.github/workflows/wheels.yml

@@ -32,7 +32,7 @@ jobs:
             platform: linux
             testvec_path: satkit-testvecs
             data_path: astro-data
-            activatestr: source .venv/bin/activate          
+            activatestr: source .venv/bin/activate
           - os: windows-latest
             testvec_path: satkit-testvecs
             data_path: astro-data
@@ -63,7 +63,7 @@ jobs:
           key: satkit-data-cache
 
       # Download data if it is not in the cache
-      - name: Download Satkit Data 
+      - name: Download Satkit Data
         if: steps.cache-satkit-data.outputs.cache-hit != 'true'
         run: |
           python -m venv .venv
@@ -79,11 +79,11 @@ jobs:
         with:
           path: ${{matrix.testvec_path}}
           enableCrossOsArchive: True
-          key: satkit-testvecs-cache        
+          key: satkit-testvecs-cache
 
       # Download test vectors if they are not in the cache
-      - name: Download Satkit Test Vectors 
-        if: steps.cache-satkit-testvecs.outputs.cache-hit != 'true' 
+      - name: Download Satkit Test Vectors
+        if: steps.cache-satkit-testvecs.outputs.cache-hit != 'true'
         run: |
           python -m venv .venv
           ${{matrix.activatestr}}
@@ -106,7 +106,7 @@ jobs:
           source .venv/bin/activate
           python -m pip install build
           python -m build --sdist . -o wheelhouse
-          
+
 
       - name: Upload Artifacts
         uses: actions/upload-artifact@v4
@@ -135,3 +135,31 @@ jobs:
         uses: pypa/gh-action-pypi-publish@release/v1
         with:
           packages-dir: dist/
+
+  readthedocs:
+      runs-on: ubuntu-latest
+      needs: publish_wheels
+      environment: readthedocs
+      steps:
+        - name: Trigger Read the Docs builds (stable + latest)
+          env:
+            RTD_TOKEN: ${{ secrets.RTD_TOKEN }}   # store this in the 'readthedocs' Environment secrets
+            RTD_HOST: app.readthedocs.org
+            RTD_PROJECT: satellite-toolkit
+          run: |
+            set -euo pipefail
+
+            # Give PyPI a moment to propagate the new release
+            sleep 20
+
+            # Optional: sync versions/tags
+            curl -fsS -X POST "https://${RTD_HOST}/api/v3/projects/${RTD_PROJECT}/sync-versions/" \
+              -H "Authorization: Token ${RTD_TOKEN}"
+
+            # Build stable
+            curl -fsS -X POST "https://${RTD_HOST}/api/v3/projects/${RTD_PROJECT}/versions/stable/builds/" \
+              -H "Authorization: Token ${RTD_TOKEN}"
+
+            # Build latest
+            curl -fsS -X POST "https://${RTD_HOST}/api/v3/projects/${RTD_PROJECT}/versions/latest/builds/" \
+              -H "Authorization: Token ${RTD_TOKEN}"

{satkit-0.9.0 → satkit-0.9.2}/CHANGES.md

@@ -155,3 +155,15 @@
   - add "as_datetime" function in python for satkit.time (for consistent nomenclature)
   - Rename time-interval boundary nomenclature from `start/stop` to `begin/end` across Rust + Python APIs (including propagation functions and related settings/results).
 
+
+## Rust 0.9.1
+  - Functions that accept time as input now accept structs that implement the new `TimeLike` trait.
+  - Add a `chrono` feature that enables interoperability with the chrono crate.  implement `TimeLike` for `chrono::DateTime`
+
+## Rust 0.9.2, Python 0.9.1
+  - bugfix: handle high-precision propagation without error when duration is zero
+  - Fix Rust documentation for SGP4 to accurately represent sgp4 source
+  - Add a "zero" static function for duration to represent zero duration
+
+## Rust 0.9.3, Python 0.9.2
+  - bugfix: Fix (and document in python) the `to_enu` and `to_ned` functions in ITRFCoord

{satkit-0.9.0 → satkit-0.9.2}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "satkit"
-version = "0.9.0"
+version = "0.9.3"
 edition = "2021"
 description = "Satellite Toolkit"
 readme = "README.md"
@@ -39,11 +39,11 @@ serde-pickle = "1.2.0"
 itertools = "0.14.0"
 anyhow = "1"
 rmpfit = "0.3.0"
-
+chrono = { version = "0.4", optional = true }
 
 [build-dependencies]
 cc = { version = "1.2", features = ["parallel"] }
-chrono = "0.4.38"
+chrono = "0.4"
 pyo3-build-config = { version = "0.27.1", optional = true }
 
 [dev-dependencies]
@@ -52,7 +52,8 @@ approx = "0.5"
 rand_distr = "0.5.1"
 
 [features]
-pybindings = ["pyo3", "numpy", "pyo3-build-config"]
+pybindings = ["dep:pyo3", "dep:numpy", "dep:pyo3-build-config"]
+chrono = ["dep:chrono"]
 
 [profile.test]
 opt-level = 3

{satkit-0.9.0 → satkit-0.9.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: satkit
-Version: 0.9.0
+Version: 0.9.2
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>
@@ -77,13 +77,18 @@ Satkit provides robust, high-performance satellite orbital mechanics calculation
   - Python versions 3.8 through 3.14 supported
   - Full documentation at <https://satellite-toolkit.readthedocs.io/latest/>
 
+## Optional Features
+
+- **chrono**: Enables interoperability with `chrono::DateTime` by implementing the `TimeLike` trait.
+  Activate with Cargo feature `chrono`.
+
 ## Key Capabilities
 
 ### Coordinate Frame Transformations
 High-precision conversions between multiple reference frames with full support for time-varying Earth orientation:
   - **ITRF** - International Terrestrial Reference Frame (Earth-fixed)
   - **GCRF** - Geocentric Celestial Reference Frame using IAU-2000/2006 reduction (inertial)
-  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation
+  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation (almost inertial)
   - **CIRF** - Celestial Intermediate Reference Frame (IAU-2006 intermediate)
   - **TIRF** - Terrestrial Intermediate Reference Frame (Earth-rotation intermediate)
   - **Geodetic** - Latitude, longitude, altitude with WGS-84 ellipsoid
@@ -98,7 +103,7 @@ Multiple propagation methods optimized for different accuracy and performance re
     - Full AFSPC and improved mode support
     - TLE fitting from high-precision states with drag estimation
     - Batch processing for multiple satellites
-    - Support TLE and Orbital Mean-Element Messages
+    - Support TLE and Orbital Mean-Element Messages inputs
   - **Keplerian**: Fast analytical two-body propagation for preliminary analysis
 
 ### Numerical Integration Force Models
@@ -125,21 +130,21 @@ Access to high-precision solar system body positions:
 Comprehensive support for all standard astronomical time scales:
   - **UTC** - Coordinated Universal Time with leap second handling
   - **TAI** - International Atomic Time
-  - **TT** - Terrestrial Time
+  - **TT**  - Terrestrial Time
   - **TDB** - Barycentric Dynamical Time
   - **UT1** - Universal Time with Earth orientation corrections
   - **GPS** - GPS Time
   - Automatic conversion between all time scales with microsecond precision
 
 ### Geodetic Utilities
-- **Geodesic Calculations**: Accurate distance and azimuth between ground locations using Vincenty's formulae
-- **Coordinate Conversions**: ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
-- **Elevation/Azimuth**: Topocentric coordinate transformations for ground station analysis
+- **Geodesic Calculations** - Accurate distance and azimuth between ground locations using Vincenty's formulae
+- **Coordinate Conversions** - ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
+- **Elevation/Azimuth** - Topocentric coordinate transformations for ground station analysis
 
 ### Sun / Moon Calculations
-- **Sun rise / set**: Compute sun rise / set times as function of day & location
-- **Moon Phase**: Phase of moon and fraction illuminated
-- **Ephemeris**: Fast low-precision ephemeris for sun & moon
+- **Sun rise / set** - Compute sun rise / set times as function of day & location
+- **Moon Phase** - Phase of moon and fraction illuminated
+- **Ephemeris** - Fast low-precision ephemeris for sun & moon
 
 ## Technical Details
 

{satkit-0.9.0 → satkit-0.9.2}/README.md

@@ -30,13 +30,18 @@ Satkit provides robust, high-performance satellite orbital mechanics calculation
   - Python versions 3.8 through 3.14 supported
   - Full documentation at <https://satellite-toolkit.readthedocs.io/latest/>
 
+## Optional Features
+
+- **chrono**: Enables interoperability with `chrono::DateTime` by implementing the `TimeLike` trait.
+  Activate with Cargo feature `chrono`.
+
 ## Key Capabilities
 
 ### Coordinate Frame Transformations
 High-precision conversions between multiple reference frames with full support for time-varying Earth orientation:
   - **ITRF** - International Terrestrial Reference Frame (Earth-fixed)
   - **GCRF** - Geocentric Celestial Reference Frame using IAU-2000/2006 reduction (inertial)
-  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation
+  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation (almost inertial)
   - **CIRF** - Celestial Intermediate Reference Frame (IAU-2006 intermediate)
   - **TIRF** - Terrestrial Intermediate Reference Frame (Earth-rotation intermediate)
   - **Geodetic** - Latitude, longitude, altitude with WGS-84 ellipsoid
@@ -51,7 +56,7 @@ Multiple propagation methods optimized for different accuracy and performance re
     - Full AFSPC and improved mode support
     - TLE fitting from high-precision states with drag estimation
     - Batch processing for multiple satellites
-    - Support TLE and Orbital Mean-Element Messages
+    - Support TLE and Orbital Mean-Element Messages inputs
   - **Keplerian**: Fast analytical two-body propagation for preliminary analysis
 
 ### Numerical Integration Force Models
@@ -78,21 +83,21 @@ Access to high-precision solar system body positions:
 Comprehensive support for all standard astronomical time scales:
   - **UTC** - Coordinated Universal Time with leap second handling
   - **TAI** - International Atomic Time
-  - **TT** - Terrestrial Time
+  - **TT**  - Terrestrial Time
   - **TDB** - Barycentric Dynamical Time
   - **UT1** - Universal Time with Earth orientation corrections
   - **GPS** - GPS Time
   - Automatic conversion between all time scales with microsecond precision
 
 ### Geodetic Utilities
-- **Geodesic Calculations**: Accurate distance and azimuth between ground locations using Vincenty's formulae
-- **Coordinate Conversions**: ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
-- **Elevation/Azimuth**: Topocentric coordinate transformations for ground station analysis
+- **Geodesic Calculations** - Accurate distance and azimuth between ground locations using Vincenty's formulae
+- **Coordinate Conversions** - ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
+- **Elevation/Azimuth** - Topocentric coordinate transformations for ground station analysis
 
 ### Sun / Moon Calculations
-- **Sun rise / set**: Compute sun rise / set times as function of day & location
-- **Moon Phase**: Phase of moon and fraction illuminated
-- **Ephemeris**: Fast low-precision ephemeris for sun & moon
+- **Sun rise / set** - Compute sun rise / set times as function of day & location
+- **Moon Phase** - Phase of moon and fraction illuminated
+- **Ephemeris** - Fast low-precision ephemeris for sun & moon
 
 ## Technical Details
 

{satkit-0.9.0 → satkit-0.9.2}/pyproject.toml

@@ -9,7 +9,7 @@ requires-python = ">= 3.8"
 authors = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 maintainers = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 readme = "README.md"
-version = "0.9.0"
+version = "0.9.2"
 license = { file = "LICENSE" }
 description = "Satellite Orbital Dynamics Toolkit"
 keywords = [
@@ -40,6 +40,8 @@ classifiers = [
 [project.optional-dependencies]
 test = ["pytest", "xmltodict"]
 
+[tool.setuptools]
+license-files = ["LICENSE"]
 
 [tool.setuptools.packages.find]
 where = ["python"]

{satkit-0.9.0 → satkit-0.9.2}/python/requirements.txt

@@ -103,7 +103,7 @@ rfc3986-validator==0.1.1; python_version >= '2.7' and python_version not in '3.0
 rpds-py==0.22.3; python_version >= '3.9'
 rsa==4.9; python_version >= '3.6'
 satkit-data==0.8.0; python_version >= '3.8'
-satkit==0.9.0; python_version >= '3.8'
+satkit; python_version >= '3.8'
 scipy==1.16.2; python_version >= '3.10'
 send2trash==1.8.3; python_version >= '2.7' and python_version not in '3.0, 3.1, 3.2, 3.3, 3.4, 3.5'
 setuptools==75.8.0; python_version >= '3.9'

{satkit-0.9.0 → satkit-0.9.2}/python/satkit/_version.py

@@ -14,5 +14,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = "0.9.0"
-__version_tuple__ = version_tuple = (0, 9, 0)
+__version__ = version = "0.9.2"
+__version_tuple__ = version_tuple = (0, 9, 2)

{satkit-0.9.0 → satkit-0.9.2}/python/satkit/satkit.pyi

@@ -2066,6 +2066,34 @@ class itrfcoord:
             satkit.quaternion: Quaternion representiong rotation from East-North-Up (ENU) to ITRF at this location
         """
 
+    def to_enu(self, refcoord: itrfcoord) -> npt.NDArray[np.float64]:
+        """Return vector from reference coordinate to this coordinate in East-North-Up (ENU) frame of the reference coordinate
+
+        Args:
+            refcoord (itrfcoord): Reference ITRF coordinate representing origin of ENU frame
+
+        Returns:
+            npt.NDArray[np.float64]: 3-element numpy array representing vector from reference coordinate to this coordinate in East-North-Up (ENU) frame of reference at the reference coordinate
+
+        Note:
+            * This is equivalent to calling: refcoord.qenu2itrf.conj * (self - refcoord)
+
+        """
+
+    def to_ned(self, refcoord: itrfcoord) -> npt.NDArray[np.float64]:
+        """Return vector from reference coordinate to this coordinate in North-East-Down (NED) at the reference coordinate
+
+        Args:
+            refcoord (itrfcoord): Reference ITRF coordinate representing origin of NED frame
+
+        Returns:
+            npt.NDArray[np.float64]: 3-element numpy array representing vector from reference coordinate to this coordinate in North-East-Down (NED) frame of reference at the reference coordinate
+
+        Note:
+            * This is equivalent to calling: refcoord.qned2itrf.conj * (self - refcoord)
+
+        """
+
     def __sub__(self, other: itrfcoord) -> npt.NDArray[np.float64]:
         """Subtract another ITRF coordinate from this one
 

{satkit-0.9.0 → satkit-0.9.2}/python/satkit.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: satkit
-Version: 0.9.0
+Version: 0.9.2
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>
@@ -77,13 +77,18 @@ Satkit provides robust, high-performance satellite orbital mechanics calculation
   - Python versions 3.8 through 3.14 supported
   - Full documentation at <https://satellite-toolkit.readthedocs.io/latest/>
 
+## Optional Features
+
+- **chrono**: Enables interoperability with `chrono::DateTime` by implementing the `TimeLike` trait.
+  Activate with Cargo feature `chrono`.
+
 ## Key Capabilities
 
 ### Coordinate Frame Transformations
 High-precision conversions between multiple reference frames with full support for time-varying Earth orientation:
   - **ITRF** - International Terrestrial Reference Frame (Earth-fixed)
   - **GCRF** - Geocentric Celestial Reference Frame using IAU-2000/2006 reduction (inertial)
-  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation
+  - **TEME** - True Equinox Mean Equator frame used in SGP4 propagation (almost inertial)
   - **CIRF** - Celestial Intermediate Reference Frame (IAU-2006 intermediate)
   - **TIRF** - Terrestrial Intermediate Reference Frame (Earth-rotation intermediate)
   - **Geodetic** - Latitude, longitude, altitude with WGS-84 ellipsoid
@@ -98,7 +103,7 @@ Multiple propagation methods optimized for different accuracy and performance re
     - Full AFSPC and improved mode support
     - TLE fitting from high-precision states with drag estimation
     - Batch processing for multiple satellites
-    - Support TLE and Orbital Mean-Element Messages
+    - Support TLE and Orbital Mean-Element Messages inputs
   - **Keplerian**: Fast analytical two-body propagation for preliminary analysis
 
 ### Numerical Integration Force Models
@@ -125,21 +130,21 @@ Access to high-precision solar system body positions:
 Comprehensive support for all standard astronomical time scales:
   - **UTC** - Coordinated Universal Time with leap second handling
   - **TAI** - International Atomic Time
-  - **TT** - Terrestrial Time
+  - **TT**  - Terrestrial Time
   - **TDB** - Barycentric Dynamical Time
   - **UT1** - Universal Time with Earth orientation corrections
   - **GPS** - GPS Time
   - Automatic conversion between all time scales with microsecond precision
 
 ### Geodetic Utilities
-- **Geodesic Calculations**: Accurate distance and azimuth between ground locations using Vincenty's formulae
-- **Coordinate Conversions**: ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
-- **Elevation/Azimuth**: Topocentric coordinate transformations for ground station analysis
+- **Geodesic Calculations** - Accurate distance and azimuth between ground locations using Vincenty's formulae
+- **Coordinate Conversions** - ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
+- **Elevation/Azimuth** - Topocentric coordinate transformations for ground station analysis
 
 ### Sun / Moon Calculations
-- **Sun rise / set**: Compute sun rise / set times as function of day & location
-- **Moon Phase**: Phase of moon and fraction illuminated
-- **Ephemeris**: Fast low-precision ephemeris for sun & moon
+- **Sun rise / set** - Compute sun rise / set times as function of day & location
+- **Moon Phase** - Phase of moon and fraction illuminated
+- **Ephemeris** - Fast low-precision ephemeris for sun & moon
 
 ## Technical Details
 

{satkit-0.9.0 → satkit-0.9.2}/python/satkit.egg-info/SOURCES.txt

@@ -9,6 +9,7 @@ README.md
 build.rs
 pyproject.toml
 .github/workflows/build.yml
+.github/workflows/cargo-publish.yml
 .github/workflows/wheels.yml
 extern/nrlmsise/DOCUMENTATION
 extern/nrlmsise/nrlmsise-00.c
@@ -169,6 +170,7 @@ src/sgp4/satrec.rs
 src/sgp4/sgp4_impl.rs
 src/sgp4/sgp4_lowlevel.rs
 src/sgp4/sgp4init.rs
+src/time/chrono.rs
 src/time/duration.rs
 src/time/instant.rs
 src/time/instant_err.rs
@@ -176,6 +178,7 @@ src/time/instant_ops.rs
 src/time/instantparse.rs
 src/time/mod.rs
 src/time/tests.rs
+src/time/timelike.rs
 src/time/timescale.rs
 src/time/weekday.rs
 src/tle/fitting.rs

{satkit-0.9.0 → satkit-0.9.2}/python/test/test.py

@@ -492,6 +492,85 @@ class TestITRFCoord:
         assert itrf.longitude_deg == pytest.approx(46.4464)
 
 
+    def test_ned_enu(self):
+
+        """
+        Test NED and ENU conversions
+        """
+        refcoord = sk.itrfcoord(
+            latitude_deg=30.0, longitude_deg=-90.0, altitude=0.0
+        )
+        testcoord = sk.itrfcoord(
+            latitude_deg=30.0, longitude_deg=-90.0, altitude=100.0
+        )
+
+        ned = testcoord.to_ned(refcoord)
+        enu = testcoord.to_enu(refcoord)
+
+        # Check NED values
+        # North component
+        assert ned[0] == pytest.approx(0, abs=1e-8)
+        # East component
+        assert ned[1] == pytest.approx(0, abs=1e-8)
+        # Down component
+        assert ned[2] == pytest.approx(-100.0, rel=1e-8)
+
+        # Check ENU values
+        # East component
+        assert enu[0] == pytest.approx(0, abs=1e-8)
+        # North component
+        assert enu[1] == pytest.approx(0, abs=1e-8)
+        # Up component
+        assert enu[2] == pytest.approx(100.0, abs=1e-8)
+
+        coord1 = sk.itrfcoord(latitude_deg=42.466, longitude_deg=-71.1516, altitude=10.0)
+
+        # Go east 10 meters and check
+        coord2 = sk.itrfcoord(coord1.vector + coord1.qenu2itrf * np.array([10.0, 0.0, 0.0]))
+        enu = coord2.to_enu(coord1)
+        assert enu[0] == pytest.approx(10.0, abs=1e-8)
+        assert enu[1] == pytest.approx(0.0, abs=1e-8)
+        assert enu[2] == pytest.approx(0.0, abs=1e-8)
+
+        # Go north 10 meters and check
+        coord2 = sk.itrfcoord(coord1.vector + coord1.qenu2itrf * np.array([0.0, 10.0, 0.0]))
+        enu = coord2.to_enu(coord1)
+        assert enu[0] == pytest.approx(0.0, abs=1e-8)
+        assert enu[1] == pytest.approx(10.0, abs=1e-8)
+        assert enu[2] == pytest.approx(0.0, abs=1e-8)
+
+        # Go up 10 meters and check
+        coord2 = sk.itrfcoord(coord1.vector + coord1.qenu2itrf * np.array([0.0, 0.0, 10.0]))
+        enu = coord2.to_enu(coord1)
+        assert enu[0] == pytest.approx(0.0, abs=1e-8)
+        assert enu[1] == pytest.approx(0.0, abs=1e-8)
+        assert enu[2] == pytest.approx(10.0, abs=1e-8)
+
+        for ix in range(50):
+            # Create random coordinates
+            lat1 = np.random.uniform(-90.0, 90.0)
+            lon1 = np.random.uniform(-180.0, 180.0)
+            alt1 = np.random.uniform(100.0, 40000.0)
+            lat2 = np.random.uniform(-90.0, 90.0)
+            lon2 = np.random.uniform(-180.0, 180.0)
+            alt2 = np.random.uniform(100.0, 40000.0)
+            coord1 = sk.itrfcoord(latitude_deg=lat1, longitude_deg=lon1, altitude=alt1)
+            coord2 = sk.itrfcoord(latitude_deg=lat2, longitude_deg=lon2, altitude=alt2)
+
+            # Check to_ned, to_enu against manually computed values
+            ned = coord2.to_ned(coord1)
+            enu = coord2.to_enu(coord1)
+            ned2 = coord1.qned2itrf.conj * (coord2-coord1)
+            enu2 = coord1.qenu2itrf.conj * (coord2-coord1)
+            assert ned[0] == pytest.approx(ned2[0], rel=1e-8)
+            assert ned[1] == pytest.approx(ned2[1], rel=1e-8)
+            assert ned[2] == pytest.approx(ned2[2], rel=1e-8)
+            assert enu[0] == pytest.approx(enu2[0], rel=1e-8)
+            assert enu[1] == pytest.approx(enu2[1], rel=1e-8)
+            assert enu[2] == pytest.approx(enu2[2], rel=1e-8)
+
+
+
 class TestMoon:
     def test_moonpos(self):
         """

{satkit-0.9.0 → satkit-0.9.2}/src/earth_orientation_params.rs

@@ -178,7 +178,7 @@ pub fn update() -> Result<()> {
     let d = datadir()?;
     if d.metadata()?.permissions().readonly() {
         bail!(
-            r#"Data directory is read-only. 
+            r#"Data directory is read-only.
              Try setting the environment variable SATKIT_DATA
              to a writeable directory and re-starting or explicitly set
              data directory"#
@@ -236,7 +236,7 @@ pub fn eop_from_mjd_utc(mjd_utc: f64) -> Option<[f64; 6]> {
 
                     This warning will only be shown once for each library load.
 
-                    To disable this warning, run `satkit.frametransform.disable_eop_time_warning()`                                                        
+                    To disable this warning, run `satkit.frametransform.disable_eop_time_warning()`
                     ",
                     mjd_utc
                 );
@@ -251,7 +251,7 @@ pub fn eop_from_mjd_utc(mjd_utc: f64) -> Option<[f64; 6]> {
 
                     This warning will only be shown once for each library load.
 
-                    To disable this warning, run `satkit::earth_orientation_params::disable_eop_time_warning()`                                                        
+                    To disable this warning, run `satkit::earth_orientation_params::disable_eop_time_warning()`
                     ",
                     mjd_utc
                 );
@@ -269,10 +269,10 @@ pub fn eop_from_mjd_utc(mjd_utc: f64) -> Option<[f64; 6]> {
                     eprintln!(
                         r"
                     Warning: EOP data not available for MJD UTC = {} (too early).
-                   
+
                     This warning will only be shown once for each library load.
 
-                    To disable this warning, run `satkit.frametransform.disable_eop_time_warning()`                                                        
+                    To disable this warning, run `satkit.frametransform.disable_eop_time_warning()`
                     ",
                         mjd_utc
                     );
@@ -280,11 +280,11 @@ pub fn eop_from_mjd_utc(mjd_utc: f64) -> Option<[f64; 6]> {
                     #[cfg(not(feature = "pybindings"))]
                     eprintln!(
                 "
-                    Warning: EOP data not available for MJD UTC = {} (too early).          
+                    Warning: EOP data not available for MJD UTC = {} (too early).
 
                     This warning will only be shown once for each library load.
 
-                    To disable this warning, run `satkit::earth_orientation_params::disable_eop_time_warning()`                                                        
+                    To disable this warning, run `satkit::earth_orientation_params::disable_eop_time_warning()`
                     ",
                     mjd_utc
                 );
@@ -336,7 +336,7 @@ pub fn eop_from_mjd_utc(mjd_utc: f64) -> Option<[f64; 6]> {
 /// ```
 ///
 #[inline]
-pub fn get(tm: &crate::Instant) -> Option<[f64; 6]> {
+pub fn get<T: crate::TimeLike>(tm: &T) -> Option<[f64; 6]> {
     eop_from_mjd_utc(tm.as_mjd_with_scale(crate::TimeScale::UTC))
 }