satkit 0.16.0__tar.gz → 0.16.2__tar.gz

{satkit-0.16.0 → satkit-0.16.2}/Cargo.toml

@@ -3,7 +3,7 @@ members = [".", "python"]
 
 [package]
 name = "satkit"
-version = "0.16.0"
+version = "0.16.2"
 edition = "2021"
 description = "Satellite Toolkit"
 readme = "README.md"
@@ -26,15 +26,11 @@ numeris = { version = "0.5.7", features = ["serde", "estimate", "ode"] }
 num-traits = "0.2.19"
 thiserror = "2.0"
 ureq = "3.1.2"
-json = "0.12.4"
 process_path = "0.1.4"
 serde = { version = "1.0", features = ["derive"] }
 serde_json = "1.0"
 quick-xml = { version = "0.38", features = ["serialize"], optional = true }
-serde-pickle = "1.2.0"
-itertools = "0.14.0"
 anyhow = "1"
-rmpfit = "0.3.0"
 chrono = { version = "0.4", optional = true }
 
 [build-dependencies]

{satkit-0.16.0/python/satkit.egg-info → satkit-0.16.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: satkit
-Version: 0.16.0
+Version: 0.16.2
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>
@@ -43,7 +43,7 @@ Satkit is a high-performance orbital mechanics library written in Rust with comp
 **[Documentation and tutorials](https://satkit.dev/)** (Python examples, but the concepts and API apply equally to Rust) | **[Rust API reference](https://docs.rs/satkit/)**
 
 > [!NOTE]
-> **Version 0.14.1** replaces `nalgebra` with [`numeris`](https://crates.io/crates/numeris) for all linear algebra (matrices, quaternions, ODE solvers). `numeris` is a standalone library that unifies these in a single crate, making them available to projects beyond satkit. Performance is equivalent, and all tests pass with matching residuals. If your code uses `nalgebra` types, enable the `nalgebra` feature on `numeris` for zero-cost conversions. Please open an issue if you hit any problems.
+> **Version 0.16.0** introduces several breaking changes: `Frame::RIC` is renamed to the canonical `Frame::RTN` (with `RIC` / `RSW` remaining as aliases, so existing code still compiles); a new `Frame::NTW` (velocity-aligned) is added; `LVLH` is now accepted as a maneuver/thrust frame; the uncertainty API is unified into `set_pos_uncertainty(sigma, frame)` / `set_vel_uncertainty(sigma, frame)` (the four old per-frame methods are removed, not deprecated); `PropSettings::default()` now uses `GravityModel::EGM96` instead of `JGM3`; the **Gauss-Jackson 8** fixed-step multistep integrator is available for long-duration propagation; and `PropSettings::max_steps` is now configurable. See `CHANGELOG.md` for the full list.
 
 ## Installation
 
@@ -133,7 +133,7 @@ let (pos, vel) = jplephem::geocentric_state(SolarSystem::Moon, &time)?;
 
 ### Coordinate Frames
 
-Full IAU-2006/2000 reduction with Earth orientation parameters:
+Full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation) with Earth orientation parameters:
 
 | Frame | Description |
 |-------|-------------|
@@ -154,7 +154,7 @@ Plus ENU, NED, and geodesic distance (Vincenty) utilities.
 
 ### Orbit Maneuvers
 
-- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (normal/tangent/cross-track — natural for prograde burns on eccentric orbits), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
+- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (velocity-aligned — natural for prograde burns on eccentric orbits, where a pure +T delta-v adds exactly Δv to |v|), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
 - **Continuous thrust** -- Constant-acceleration thrust arcs over time windows in any of the frames above, integrated directly into the force model
 - **Automatic segmentation** -- Propagation through maneuver sequences is handled transparently, including backward propagation
 
@@ -180,7 +180,7 @@ Seamless conversion between UTC, TAI, TT, TDB, UT1, and GPS time scales with ful
 SatKit uses [numeris](https://crates.io/crates/numeris) for all linear algebra (vectors, matrices, quaternions, ODE integration). If you also use nalgebra in your project, enable the `nalgebra` feature on numeris for zero-cost `From`/`Into` conversions between types:
 
 ```toml
-numeris = { version = "0.5.5", features = ["nalgebra"] }
+numeris = { version = "0.5.7", features = ["nalgebra"] }
 ```
 
 ### Cargo Features
@@ -207,7 +207,7 @@ The library is validated against:
 - **ICGEM** reference values for gravity field calculations
 - **GPS SP3** precise ephemerides for multi-day numerical propagation
 
-142 tests (106 unit + 36 doc-tests) run on every commit across Linux, macOS, and Windows.
+157 Rust tests and 81 Python tests run on every commit across Linux, macOS, and Windows.
 
 ### Running Tests Locally
 

{satkit-0.16.0 → satkit-0.16.2}/README.md

@@ -19,7 +19,7 @@ Satkit is a high-performance orbital mechanics library written in Rust with comp
 **[Documentation and tutorials](https://satkit.dev/)** (Python examples, but the concepts and API apply equally to Rust) | **[Rust API reference](https://docs.rs/satkit/)**
 
 > [!NOTE]
-> **Version 0.14.1** replaces `nalgebra` with [`numeris`](https://crates.io/crates/numeris) for all linear algebra (matrices, quaternions, ODE solvers). `numeris` is a standalone library that unifies these in a single crate, making them available to projects beyond satkit. Performance is equivalent, and all tests pass with matching residuals. If your code uses `nalgebra` types, enable the `nalgebra` feature on `numeris` for zero-cost conversions. Please open an issue if you hit any problems.
+> **Version 0.16.0** introduces several breaking changes: `Frame::RIC` is renamed to the canonical `Frame::RTN` (with `RIC` / `RSW` remaining as aliases, so existing code still compiles); a new `Frame::NTW` (velocity-aligned) is added; `LVLH` is now accepted as a maneuver/thrust frame; the uncertainty API is unified into `set_pos_uncertainty(sigma, frame)` / `set_vel_uncertainty(sigma, frame)` (the four old per-frame methods are removed, not deprecated); `PropSettings::default()` now uses `GravityModel::EGM96` instead of `JGM3`; the **Gauss-Jackson 8** fixed-step multistep integrator is available for long-duration propagation; and `PropSettings::max_steps` is now configurable. See `CHANGELOG.md` for the full list.
 
 ## Installation
 
@@ -109,7 +109,7 @@ let (pos, vel) = jplephem::geocentric_state(SolarSystem::Moon, &time)?;
 
 ### Coordinate Frames
 
-Full IAU-2006/2000 reduction with Earth orientation parameters:
+Full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation) with Earth orientation parameters:
 
 | Frame | Description |
 |-------|-------------|
@@ -130,7 +130,7 @@ Plus ENU, NED, and geodesic distance (Vincenty) utilities.
 
 ### Orbit Maneuvers
 
-- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (normal/tangent/cross-track — natural for prograde burns on eccentric orbits), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
+- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (velocity-aligned — natural for prograde burns on eccentric orbits, where a pure +T delta-v adds exactly Δv to |v|), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
 - **Continuous thrust** -- Constant-acceleration thrust arcs over time windows in any of the frames above, integrated directly into the force model
 - **Automatic segmentation** -- Propagation through maneuver sequences is handled transparently, including backward propagation
 
@@ -156,7 +156,7 @@ Seamless conversion between UTC, TAI, TT, TDB, UT1, and GPS time scales with ful
 SatKit uses [numeris](https://crates.io/crates/numeris) for all linear algebra (vectors, matrices, quaternions, ODE integration). If you also use nalgebra in your project, enable the `nalgebra` feature on numeris for zero-cost `From`/`Into` conversions between types:
 
 ```toml
-numeris = { version = "0.5.5", features = ["nalgebra"] }
+numeris = { version = "0.5.7", features = ["nalgebra"] }
 ```
 
 ### Cargo Features
@@ -183,7 +183,7 @@ The library is validated against:
 - **ICGEM** reference values for gravity field calculations
 - **GPS SP3** precise ephemerides for multi-day numerical propagation
 
-142 tests (106 unit + 36 doc-tests) run on every commit across Linux, macOS, and Windows.
+157 Rust tests and 81 Python tests run on every commit across Linux, macOS, and Windows.
 
 ### Running Tests Locally
 

{satkit-0.16.0 → satkit-0.16.2}/pyproject.toml

@@ -9,7 +9,7 @@ requires-python = ">= 3.10"
 authors = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 maintainers = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 readme = "README.md"
-version = "0.16.0"
+version = "0.16.2"
 license = "MIT"
 description = "Satellite Orbital Dynamics Toolkit"
 keywords = [

{satkit-0.16.0 → satkit-0.16.2}/python/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "satkit-python"
-version = "0.16.0"
+version = "0.16.2"
 edition = "2021"
 publish = false
 
@@ -16,7 +16,6 @@ numeris = { version = "0.5.6", features = ["serde", "estimate"] }
 anyhow = "1"
 serde = { version = "1.0", features = ["derive"] }
 serde-pickle = "1.2.0"
-rmpfit = "0.3.0"
 process_path = "0.1.4"
 
 [build-dependencies]

{satkit-0.16.0 → satkit-0.16.2}/python/satkit/__init__.pyi

@@ -48,5 +48,6 @@ __all__ = [
     "propresult",
     "propstats",
     "omm_from_url",
+    "tlefitstatus",
     "__version__",
 ]

{satkit-0.16.0 → satkit-0.16.2}/python/satkit/frametransform.pyi

@@ -491,7 +491,7 @@ def qgcrf2itrf(
     """Quaternion representing rotation from the Geocentric Celestial Reference Frame (GCRF) to the International Terrestrial Reference Frame (ITRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -529,7 +529,7 @@ def qgcrf2itrf(
     """Quaternion representing rotation from the Geocentric Celestial Reference Frame (GCRF) to the International Terrestrial Reference Frame (ITRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -552,7 +552,7 @@ def qgcrf2itrf(*args, **kwargs):
     """Quaternion representing rotation from the Geocentric Celestial Reference Frame (GCRF) to the International Terrestrial Reference Frame (ITRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -590,7 +590,7 @@ def qitrf2gcrf(
     """Quaternion representing rotation from the International Terrestrial Reference Frame (ITRF) to the Geocentric Celestial Reference Frame (GCRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -621,7 +621,7 @@ def qitrf2gcrf(
     """Quaternion representing rotation from the International Terrestrial Reference Frame (ITRF) to the Geocentric Celestial Reference Frame (GCRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -644,7 +644,7 @@ def qitrf2gcrf(*args, **kwargs):
     """Quaternion representing rotation from the International Terrestrial Reference Frame (ITRF) to the Geocentric Celestial Reference Frame (GCRF)
 
     Notes:
-        - Uses full IAU2010 Reduction
+        - Uses full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation)
         - See IERS Technical Note 36, Chapter 5
         - Does not include solid tides, ocean tides
         - Very computationally expensive
@@ -873,16 +873,21 @@ def to_gcrf(
 def itrf_to_gcrf_state(
     pos_itrf: npt.ArrayLike,
     vel_itrf: npt.ArrayLike,
-    time: time,
+    time: time | list[time] | npt.NDArray,
 ) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
     """Transform a satellite state (position + velocity) from ITRF to GCRF.
 
+    Accepts either a single state (``pos``/``vel`` are length-3 vectors and
+    ``time`` is a single ``satkit.time``) or a batch of ``N`` states
+    (``pos``/``vel`` are shape ``(N, 3)`` arrays and ``time`` is a length-``N``
+    array/list of times). The output shape matches the input.
+
     Unlike the raw :func:`qitrf2gcrf` quaternion, this function correctly
     handles the Earth-rotation contribution to velocity. A point at rest
     on Earth's surface has zero velocity in ITRF but ~465 m/s in GCRF at
     the equator, and this function accounts for that term.
 
-    The IAU 2010 ITRF → GCRF reduction decomposes into three stages:
+    The IERS 2010 ITRF → GCRF reduction decomposes into three stages:
     polar motion (ITRF → TIRS), Earth rotation about the CIO polar axis
     (TIRS → CIRS), and precession-nutation (CIRS → GCRF). The
     Earth-rotation sweep term ``omega_earth x r`` is computed in
@@ -897,9 +902,9 @@ def itrf_to_gcrf_state(
     1. Rotate ``pos_itrf`` and ``vel_itrf`` into TIRS via polar motion.
     2. Add ``omega_earth x r_tirs`` to the velocity in TIRS, where
        ``omega_earth = (0, 0, OMEGA_EARTH)`` exactly.
-    3. Rotate TIRS → CIRS → GCRF via the full IAU 2010 chain.
+    3. Rotate TIRS → CIRS → GCRF via the full IERS 2010 chain.
 
-    Uses the full IAU 2010 reduction (polar motion + Earth rotation +
+    Uses the full IERS 2010 reduction (polar motion + Earth rotation +
     precession-nutation with dX/dY corrections from Earth orientation
     parameters).
 
@@ -932,7 +937,7 @@ def itrf_to_gcrf_state(
 def gcrf_to_itrf_state(
     pos_gcrf: npt.ArrayLike,
     vel_gcrf: npt.ArrayLike,
-    time: time,
+    time: time | list[time] | npt.NDArray,
 ) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
     """Transform a satellite state (position + velocity) from GCRF to ITRF.
 
@@ -941,7 +946,12 @@ def gcrf_to_itrf_state(
     term **in TIRS** (where Earth's rotation axis is exactly along +z),
     then applies inverse polar motion to reach ITRF. A geostationary
     satellite (whose GCRF velocity is pure orbital motion) produces
-    zero velocity in ITRF. Uses the full IAU 2010 reduction.
+    zero velocity in ITRF. Uses the full IERS 2010 reduction.
+
+    Accepts either a single state or a batch of ``N`` states: when
+    ``pos``/``vel`` are shape ``(N, 3)`` arrays, ``time`` must be a
+    length-``N`` array/list of times, and the returned arrays have
+    shape ``(N, 3)``.
 
     Args:
         pos_gcrf: 3-element position vector in GCRF [m]
@@ -954,6 +964,32 @@ def gcrf_to_itrf_state(
     """
     ...
 
+def itrf_to_gcrf_state_approx(
+    pos_itrf: npt.ArrayLike,
+    vel_itrf: npt.ArrayLike,
+    time: time | list[time] | npt.NDArray,
+) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
+    """Approximate ITRF → GCRF state transform using the IAU-76/FK5
+    reduction (accurate to ~1 arcsec on position).
+
+    Faster alternative to :func:`itrf_to_gcrf_state` when the full IERS
+    2010 precision is not required. Neglects polar motion, so the
+    Earth-rotation sweep ``omega_earth x r`` is evaluated in ITRF directly.
+    Accepts scalar or batched inputs like :func:`itrf_to_gcrf_state`.
+    """
+    ...
+
+def gcrf_to_itrf_state_approx(
+    pos_gcrf: npt.ArrayLike,
+    vel_gcrf: npt.ArrayLike,
+    time: time | list[time] | npt.NDArray,
+) -> tuple[npt.NDArray[np.float64], npt.NDArray[np.float64]]:
+    """Approximate GCRF → ITRF state transform using the IAU-76/FK5
+    reduction. Inverse of :func:`itrf_to_gcrf_state_approx`; accurate to
+    ~1 arcsec on position. Accepts scalar or batched inputs.
+    """
+    ...
+
 def from_gcrf(
     frame: frame,
     pos: npt.ArrayLike,

{satkit-0.16.0 → satkit-0.16.2}/python/satkit/satkit.pyi

@@ -295,14 +295,15 @@ class TLE:
             Input GCRF states are rotated into TEME frame used by SGP4.
             First and second derivatives of mean motion are ignored, as they are not used by SGP4.
 
-            Non-linear Levenberg-Marquardt optimization (via the Rust ``rmpfit`` crate)
-            is performed to fit inclination, eccentricity, RAAN, argument of perigee,
-            mean anomaly, mean motion, and drag (bstar) to the provided states.
+            Non-linear Levenberg-Marquardt optimization is performed to fit
+            inclination, eccentricity, RAAN, argument of perigee, mean anomaly,
+            mean motion, and drag (bstar) to the provided states. The solver
+            is built on top of the ``numeris`` linear algebra crate.
 
             The results dictionary includes the following keys:
-            ``success``, ``best_norm``, ``orig_norm``, ``n_iter``, ``n_fev``,
-            ``n_par``, ``n_free``, ``n_pegged``, ``n_func``, ``resid``,
-            ``xerror``, ``covar``.
+            ``status`` (a :class:`tlefitstatus`), ``converged`` (bool),
+            ``orig_norm``, ``best_norm``, ``grad_norm``, ``n_iter``,
+            ``n_res_evals``.
 
         Example:
             ```python
@@ -314,7 +315,7 @@ class TLE:
             epoch = satkit.time(2024, 1, 1)
 
             tle, results = satkit.TLE.fit_from_states(states, times, epoch)
-            if results["success"] == satkit.mpsuccess.MP_OK_CHI:
+            if results["converged"]:
                 print("Fit successful")
             ```
         """
@@ -605,51 +606,38 @@ class weekday:
     Saturday: ClassVar[weekday]
     """Saturday"""
 
-class mpsuccess:
+class tlefitstatus:
     """
-    State of Levenberg-Marquardt optimization from the `rmpfit` rust library
-
-    For details see: <https://docs.rs/rmpfit/latest/rmpfit/>
+    Termination status of the TLE non-linear least-squares fit performed by
+    :meth:`TLE.fit_from_states`.
 
     Values:
 
-    - `NotDone`: Not finished iterations
-    - `Chi`: Convergence in chi-square value
-    - `Par`: Convergence in parameter value
-    - `Both`: Convergence in both chi-square and parameter values
-    - `Dir`: Convergence in orthogonality
-    - `MaxIter`: Maximum iterations reached
-    - `Ftol`: ftol is too small; no further improvement
-    - `Xtol`: xtol is too small; no further improvement
-    - `Gtol`: gtol is too small; no further improvement
+    - ``GradientConverged``: converged on gradient norm tolerance
+    - ``StepConverged``: converged on relative step size tolerance
+    - ``CostConverged``: converged on relative cost change tolerance
+    - ``MaxIterations``: maximum number of iterations reached
+    - ``DampingSaturated``: Levenberg-Marquardt damping parameter saturated
     """
 
-    NotDone: ClassVar[mpsuccess]
-    """Not finished iterations"""
-
-    Chi: ClassVar[mpsuccess]
-    """Convergence in chi-square value"""
-
-    Par: ClassVar[mpsuccess]
-    """Convergence in parameter value"""
+    GradientConverged: ClassVar[tlefitstatus]
+    """Converged on gradient norm tolerance"""
 
-    Both: ClassVar[mpsuccess]
-    """Convergence in both chi-square and parameter values"""
+    StepConverged: ClassVar[tlefitstatus]
+    """Converged on relative step size tolerance"""
 
-    Dir: ClassVar[mpsuccess]
-    """Convergence in orthogonality"""
+    CostConverged: ClassVar[tlefitstatus]
+    """Converged on relative cost change tolerance"""
 
-    MaxIter: ClassVar[mpsuccess]
-    """Maximum iterations reached"""
+    MaxIterations: ClassVar[tlefitstatus]
+    """Maximum number of iterations reached"""
 
-    Ftol: ClassVar[mpsuccess]
-    """ftol is too small; no further improvement"""
+    DampingSaturated: ClassVar[tlefitstatus]
+    """Levenberg-Marquardt damping parameter saturated"""
 
-    Xtol: ClassVar[mpsuccess]
-    """xtol is too small; no further improvement"""
-
-    Gtol: ClassVar[mpsuccess]
-    """gtol is too small; no further improvement"""
+    def converged(self) -> bool:
+        """True if the fit converged successfully."""
+        ...
 
 class timescale:
     """

{satkit-0.16.0 → satkit-0.16.2/python/satkit.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: satkit
-Version: 0.16.0
+Version: 0.16.2
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>
@@ -43,7 +43,7 @@ Satkit is a high-performance orbital mechanics library written in Rust with comp
 **[Documentation and tutorials](https://satkit.dev/)** (Python examples, but the concepts and API apply equally to Rust) | **[Rust API reference](https://docs.rs/satkit/)**
 
 > [!NOTE]
-> **Version 0.14.1** replaces `nalgebra` with [`numeris`](https://crates.io/crates/numeris) for all linear algebra (matrices, quaternions, ODE solvers). `numeris` is a standalone library that unifies these in a single crate, making them available to projects beyond satkit. Performance is equivalent, and all tests pass with matching residuals. If your code uses `nalgebra` types, enable the `nalgebra` feature on `numeris` for zero-cost conversions. Please open an issue if you hit any problems.
+> **Version 0.16.0** introduces several breaking changes: `Frame::RIC` is renamed to the canonical `Frame::RTN` (with `RIC` / `RSW` remaining as aliases, so existing code still compiles); a new `Frame::NTW` (velocity-aligned) is added; `LVLH` is now accepted as a maneuver/thrust frame; the uncertainty API is unified into `set_pos_uncertainty(sigma, frame)` / `set_vel_uncertainty(sigma, frame)` (the four old per-frame methods are removed, not deprecated); `PropSettings::default()` now uses `GravityModel::EGM96` instead of `JGM3`; the **Gauss-Jackson 8** fixed-step multistep integrator is available for long-duration propagation; and `PropSettings::max_steps` is now configurable. See `CHANGELOG.md` for the full list.
 
 ## Installation
 
@@ -133,7 +133,7 @@ let (pos, vel) = jplephem::geocentric_state(SolarSystem::Moon, &time)?;
 
 ### Coordinate Frames
 
-Full IAU-2006/2000 reduction with Earth orientation parameters:
+Full IERS 2010 Conventions reduction (IAU 2006/2000A precession-nutation) with Earth orientation parameters:
 
 | Frame | Description |
 |-------|-------------|
@@ -154,7 +154,7 @@ Plus ENU, NED, and geodesic distance (Vincenty) utilities.
 
 ### Orbit Maneuvers
 
-- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (normal/tangent/cross-track — natural for prograde burns on eccentric orbits), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
+- **Impulsive maneuvers** -- Instantaneous delta-v applied at a scheduled time during propagation. Supported frames: GCRF (inertial), RTN (radial/tangential/normal — the CCSDS OEM convention, also exposed as `RSW` and `RIC` aliases), NTW (velocity-aligned — natural for prograde burns on eccentric orbits, where a pure +T delta-v adds exactly Δv to |v|), and LVLH (Local Vertical / Local Horizontal). Ergonomic helpers `add_prograde` / `add_retrograde` / `add_radial` / `add_normal` for common scalar-magnitude burns.
 - **Continuous thrust** -- Constant-acceleration thrust arcs over time windows in any of the frames above, integrated directly into the force model
 - **Automatic segmentation** -- Propagation through maneuver sequences is handled transparently, including backward propagation
 
@@ -180,7 +180,7 @@ Seamless conversion between UTC, TAI, TT, TDB, UT1, and GPS time scales with ful
 SatKit uses [numeris](https://crates.io/crates/numeris) for all linear algebra (vectors, matrices, quaternions, ODE integration). If you also use nalgebra in your project, enable the `nalgebra` feature on numeris for zero-cost `From`/`Into` conversions between types:
 
 ```toml
-numeris = { version = "0.5.5", features = ["nalgebra"] }
+numeris = { version = "0.5.7", features = ["nalgebra"] }
 ```
 
 ### Cargo Features
@@ -207,7 +207,7 @@ The library is validated against:
 - **ICGEM** reference values for gravity field calculations
 - **GPS SP3** precise ephemerides for multi-day numerical propagation
 
-142 tests (106 unit + 36 doc-tests) run on every commit across Linux, macOS, and Windows.
+157 Rust tests and 81 Python tests run on every commit across Linux, macOS, and Windows.
 
 ### Running Tests Locally
 

{satkit-0.16.0 → satkit-0.16.2}/python/satkit.egg-info/SOURCES.txt

@@ -37,7 +37,6 @@ python/src/pylambert.rs
 python/src/pylpephem_moon.rs
 python/src/pylpephem_planets.rs
 python/src/pylpephem_sun.rs
-python/src/pympsuccess.rs
 python/src/pynrlmsise.rs
 python/src/pyomm.rs
 python/src/pypropagate.rs
@@ -50,6 +49,7 @@ python/src/pysgp4.rs
 python/src/pysolarsystem.rs
 python/src/pythrust.rs
 python/src/pytle.rs
+python/src/pytlefitstatus.rs
 python/src/pyukf.rs
 python/src/pyutils.rs
 src/consts.rs

{satkit-0.16.0 → satkit-0.16.2}/python/src/lib.rs

@@ -15,7 +15,7 @@ mod pykepler;
 mod pylpephem_moon;
 mod pylpephem_planets;
 mod pylpephem_sun;
-mod pympsuccess;
+mod pytlefitstatus;
 mod pynrlmsise;
 mod pypropresult;
 mod pyquaternion;
@@ -134,6 +134,10 @@ fn frametransform(_py: Python, m: &Bound<'_, PyModule>) -> PyResult<()> {
         .unwrap();
     m.add_function(wrap_pyfunction!(pyft::gcrf_to_itrf_state, m)?)
         .unwrap();
+    m.add_function(wrap_pyfunction!(pyft::itrf_to_gcrf_state_approx, m)?)
+        .unwrap();
+    m.add_function(wrap_pyfunction!(pyft::gcrf_to_itrf_state_approx, m)?)
+        .unwrap();
     m.add_function(wrap_pyfunction!(pyft::qmod2gcrf, m)?)
         .unwrap();
     m.add_function(wrap_pyfunction!(pyft::qtod2mod_approx, m)?)
@@ -168,7 +172,7 @@ pub fn satkit(_py: Python, m: &Bound<'_, PyModule>) -> PyResult<()> {
     m.add_class::<pyconsts::Consts>()?;
     m.add_class::<SolarSystem>()?;
     m.add_class::<pytle::PyTLE>()?;
-    m.add_class::<pympsuccess::PyMPSuccess>()?;
+    m.add_class::<pytlefitstatus::PyTleFitStatus>()?;
 
     m.add_class::<PyGeodet>()?;
     m.add_class::<PyITRFCoord>()?;