satkit 0.3.0__tar.gz → 0.3.1__tar.gz

{satkit-0.3.0 → satkit-0.3.1}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "satkit"
-version = "0.3.0"
+version = "0.3.1"
 edition = "2021"
 description = "Satellite Toolkit"
 readme = "README.md"
@@ -21,22 +21,22 @@ name = "satkit"
 
 [dependencies]
 nalgebra = "0.33.0"
-ndarray = "0.15.6"
-libc = "0.2"
+ndarray = "0.15.4"
+libc = "0.2.158"
 chrono = "0.4.19"
 cty = "0.2.2"
 num-traits = "0.2"
-thiserror = "1.0"
+thiserror = "1.0.63"
 simba = "0.9"
 once_cell = "1.6.0"
 numpy = { version = "0.21.0", optional = true }
 pyo3 = { version = "0.21", features = ["extension-module"], optional = true }
-ureq = "2.9.1"
+ureq = "2.10.1"
 json = "0.12.4"
 process_path = "0.1.4"
 
 [build-dependencies]
-cc = { version = "1.0.74", features = ["parallel"] }
+cc = { version = "1.1.15", features = ["parallel"] }
 chrono = "0.4.38"
 pyo3-build-config = { version = "0.21", optional = true }
 

{satkit-0.3.0/python/satkit.egg-info → satkit-0.3.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: satkit
-Version: 0.3.0
+Version: 0.3.1
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>

{satkit-0.3.0 → satkit-0.3.1}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "satkit"
-version = "0.3.0"
+version = "0.3.1"
 dependencies = ["numpy>=1.20.0"]
 requires-python = ">= 3.8"
 authors = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]

{satkit-0.3.0 → satkit-0.3.1}/python/satkit/satkit.pyi

@@ -1490,7 +1490,7 @@ class satstate:
         Args:
             time (satkit.time|satkit.duration): Time or duration from current time to which to propagate the state
 
-        Keyword Arguments:
+        Keyword Args:
             propsettings: satkit.satprop.propsettings object describing settings to use in the propagation. If omitted, default is used 
 
         Returns:

{satkit-0.3.0 → satkit-0.3.1}/python/satkit/satprop.pyi

@@ -228,17 +228,15 @@ def propagate(
         vel (npt.ArrayLike[float]): 3-element numpy array representing satellite GCRF velocity in m/s
         tm (satkit.time): satkit.time object representing instant at which satellite is at "pos" & "vel"
 
-    Keyword Arguments:
+    Keyword Args: 
         stop_time (satkit.time): satkit.time object representing instant at which new position and velocity will be computed
         duration_secs (float): duration in seconds from "tm" for at which new position and velocity will be computed.
         duration_days (float): duration in days from "tm" at which new position and velocity will be computed.
         duration (satkit.duration): duration from "tm" at which new position & velocity will be computed.
         output_phi (bool): Output 6x6 state transition matrix between "starttime" and "stoptime" (and at intervals, if specified)
         propsettings (propsettings): "propsettings" object with input settings for the propagation. if left out, default will be used.
-        satproperties (satproperties_static): "SatPropertiesStatic" object with drag and radiation pressure succeptibility of satellite. 
-                                            If left out, drag and radiation pressure are neglected
-        output_dense: boolean indicting whether or not dense output should be recorded.  Default is false.  If true, this will allow for calling the "interp" function
-                      to query states at arbitrary times between the start time and the stop time
+        satproperties (satproperties_static): "SatPropertiesStatic" object with drag and radiation pressure succeptibility of satellite. If left out, drag and radiation pressure are neglected
+        output_dense: boolean indicting whether or not dense output should be recorded.  Default is false.  If true, this will allow for calling the "interp" function to query states at arbitrary times between the start time and the stop time
 
     Returns:
         (propresult): Propagation result object holding state outputs, statistics, and dense output if requested

{satkit-0.3.0 → satkit-0.3.1}/python/satkit/utils.pyi

@@ -12,7 +12,7 @@ import satkit
 def update_datafiles(**kwargs):
     """Download & store data files needed for "satkit" computations
 
-    Keyword Arguments:
+    Keyword Args:
 
       overwrite (bool):  Download and overwrite files if they already exist
       dir(string): Target directory for files.  Uses existing data directory if not specified. (see "datadir" function)

{satkit-0.3.0 → satkit-0.3.1/python/satkit.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: satkit
-Version: 0.3.0
+Version: 0.3.1
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>

{satkit-0.3.0 → satkit-0.3.1}/python/test/orbitprop_gps_fit.py

@@ -25,26 +25,20 @@ import os
 basedir = os.getenv(
     "SATKIT_TESTVEC_ROOT", default="../.." + os.path.sep + "satkit-testvecs"
 )
+fname = (
+    basedir + os.path.sep + "orbitprop" + os.path.sep
+    + "ESA0OPSFIN_20233640000_01D_05M_ORB.SP3"
+)
+
+[pitrf, timearr] = read_sp3file(fname)
+pgcrf = np.stack(
+    np.fromiter(
+        (q * p for q, p in zip(sk.frametransform.qitrf2gcrf(timearr), pitrf)),
+        list,
+    ),
+    axis=0,
+)
 
-sp3names = ['mid22630.sp3', 'mid22631.sp3', 'mid22632.sp3',
-          'mid22633.sp3', 'mid22634.sp3', 'mid22635.sp3',
-          'mid22636.sp3']
-
-sp3names = ['mid22630.sp3', 'mid22631.sp3', 'mid22632.sp3']
-pitrf = np.zeros((0,3), np.float64)
-timearr = np.array([])
-
-for sp3name in sp3names:
-    fname = (
-        basedir
-        + os.path.sep
-        + "orbitprop"
-        + os.path.sep
-        + sp3name
-    )
-    [pitrf1, timearray1] = read_sp3file(fname)
-    pitrf = np.append(pitrf, pitrf1, axis=0)
-    timearr = np.append(timearr, timearray1)
 
 print(pitrf.shape)
 print(timearr.shape)
@@ -62,14 +56,15 @@ vgcrf = (pgcrf[1, :] - pgcrf[0, :]) / (timearr[1] - timearr[0]).seconds()
 
 # Initial state for non-linear least squares is initial velocity
 # and susceptibility to radiation pressuer : Cr A / m
-v0 = np.array([vgcrf[0], vgcrf[1], vgcrf[2], 0.01])
-
+#v0 = np.array([vgcrf[0], vgcrf[1], vgcrf[2], 0.02])
+v0 = np.array([2.47130555e03, 2.94682777e03, -5.34171918e02, 2.13018578e-02])
 
 def minfunc(v):
     tstart = timearr[0]
     tend = timearr[-1]
     settings = sk.satprop.propsettings()
     settings.use_jplephem = False
+    settings.gravity_order = 10
     satprops = sk.satprop.satproperties_static()
     satprops.craoverm = v[3]
 
@@ -78,12 +73,15 @@ def minfunc(v):
         v[0:3],
         tstart,
         stoptime=tend,
-        dt=sk.duration.from_minutes(5),
         propsettings=settings,
-        satproperties=satprops,
+        satproperties=satprops,            
+        output_dense=True,
     )
+    pest = np.zeros((len(timearr), 3))
+    for i in range(len(timearr)):
+        pest[i, :] = res.interp(timearr[i])[0:3]
 
-    return np.sum(np.sum((res["pos"] - pgcrf) ** 2, axis=1), axis=0)
+    return np.sum(np.sum((pest - pgcrf) ** 2, axis=1), axis=0)
 
 
 # Minimize difference between true satellite state and propagated state
@@ -91,6 +89,10 @@ def minfunc(v):
 r = minimize(minfunc, v0, method="Nelder-Mead")
 print(r)
 
+# %%
+print(v0)
+print(r.x)
+
 # %%
 import plotly.graph_objects as go
 
@@ -100,11 +102,15 @@ satprops = sk.satprop.satproperties_static()
 satprops.craoverm = r.x[3]
 res = sk.satprop.propagate(
     pgcrf[0,:], r.x[0:3], timearr[0],
-    stoptime=timearr[-1], dt=sk.duration(minutes=5),
-    propsettings=settings, satproperties=satprops
+    stoptime=timearr[-1],
+    propsettings=settings, satproperties=satprops,
+    output_dense=True
 )
-perr = res['pos'] - pgcrf
+perr = np.zeros((len(timearr), 3))
+for i in range(len(timearr)):
+    perr[i,:] = res.interp(timearr[i])[0:3] - pgcrf[i,:]
 fig = go.Figure()
+
 fig.add_trace(go.Scatter(x=[t.datetime() for t in timearr], y=perr[:,0], mode='lines', name='X'))
 fig.add_trace(go.Scatter(x=[t.datetime() for t in timearr], y=perr[:,1], mode='lines', name='Y'))
 fig.add_trace(go.Scatter(x=[t.datetime() for t in timearr], y=perr[:,2], mode='lines', name='Z'))

{satkit-0.3.0 → satkit-0.3.1}/python/test/test.py

@@ -598,13 +598,19 @@ class TestHighPrecisionPropagation:
         settings = sk.satprop.propsettings()
         settings.use_jplephem = False
 
+        # Determined by orbitprop_gps_fit.py
+        fitparam = np.array(
+            [ 2.47130562e+03,  2.94682753e+03, -5.34172176e+02,  2.32565692e-02]
+        )
+
+
         # Values for craoverm and velocity come from orbitprop_gps_fit.py
         satprops = sk.satprop.satproperties_static()
-        satprops.craoverm = 2.13018578e-2
-        print(satprops)
+        satprops.craoverm = fitparam[3]
+
         res = sk.satprop.propagate(
             pgcrf[0, :],
-            np.array([2.47130555e03, 2.94682777e03, -5.34171918e02]),
+            fitparam[0:3],
             timearr[0],
             stoptime=timearr[-1],
             propsettings=settings,
@@ -613,12 +619,12 @@ class TestHighPrecisionPropagation:
         )
         print(res)
 
-        # See if propagator is accurate to < 10 meters over 1 day on
+        # See if propagator is accurate to < 8 meters over 1 day on
         # each Cartesian axis
         for iv in range(pgcrf.shape[0]-5):            
             state = res.interp(timearr[iv])
             for ix in range(0, 3):
-                assert m.fabs(state[ix] - pgcrf[iv, ix]) < 10
+                assert m.fabs(state[ix] - pgcrf[iv, ix]) < 8
 
 class TestSatState:
     def test_lvlh(self):

{satkit-0.3.0 → satkit-0.3.1}/src/earthgravity.rs

@@ -159,7 +159,7 @@ type Mat3 = na::Matrix3<f64>;
 /// * `pos` - Position as ITRF coordinate (satkit.itrfcoord) or numpy
 ///                3-vector representing ITRF position in meters
 ///
-/// * `order` - Order of the gravity model
+/// * `order` - Order of the gravity model, up to 40
 ///
 /// # References
 ///
@@ -199,8 +199,56 @@ impl Gravity {
             self.accel_t::<14, 18>(pos)
         } else if order == 15 {
             self.accel_t::<15, 19>(pos)
-        } else {
+        } else if order == 16 {
             self.accel_t::<16, 20>(pos)
+        } else if order == 17 {
+            self.accel_t::<17, 21>(pos)
+        } else if order == 18 {
+            self.accel_t::<18, 22>(pos)
+        } else if order == 19 {
+            self.accel_t::<19, 23>(pos)
+        } else if order == 20 {
+            self.accel_t::<20, 24>(pos)
+        } else if order == 21 {
+            self.accel_t::<21, 25>(pos)
+        } else if order == 22 {
+            self.accel_t::<22, 26>(pos)
+        } else if order == 23 {
+            self.accel_t::<23, 27>(pos)
+        } else if order == 24 {
+            self.accel_t::<24, 28>(pos)
+        } else if order == 25 {
+            self.accel_t::<25, 29>(pos)
+        } else if order == 26 {
+            self.accel_t::<26, 30>(pos)
+        } else if order == 27 {
+            self.accel_t::<27, 31>(pos)
+        } else if order == 28 {
+            self.accel_t::<28, 32>(pos)
+        } else if order == 29 {
+            self.accel_t::<29, 33>(pos)
+        } else if order == 30 {
+            self.accel_t::<30, 34>(pos)
+        } else if order == 31 {
+            self.accel_t::<31, 35>(pos)
+        } else if order == 32 {
+            self.accel_t::<32, 36>(pos)
+        } else if order == 33 {
+            self.accel_t::<33, 37>(pos)
+        } else if order == 34 {
+            self.accel_t::<34, 38>(pos)
+        } else if order == 35 {
+            self.accel_t::<35, 39>(pos)
+        } else if order == 36 {
+            self.accel_t::<36, 40>(pos)
+        } else if order == 37 {
+            self.accel_t::<37, 41>(pos)
+        } else if order == 37 {
+            self.accel_t::<38, 42>(pos)
+        } else if order == 39 {
+            self.accel_t::<39, 43>(pos)
+        } else {
+            self.accel_t::<40, 44>(pos)
         }
     }
 
@@ -235,8 +283,56 @@ impl Gravity {
             self.accel_and_partials_t::<14, 18>(pos)
         } else if order == 15 {
             self.accel_and_partials_t::<15, 19>(pos)
-        } else {
+        } else if order == 16 {
             self.accel_and_partials_t::<16, 20>(pos)
+        } else if order == 17 {
+            self.accel_and_partials_t::<17, 21>(pos)
+        } else if order == 18 {
+            self.accel_and_partials_t::<18, 22>(pos)
+        } else if order == 19 {
+            self.accel_and_partials_t::<19, 23>(pos)
+        } else if order == 20 {
+            self.accel_and_partials_t::<20, 24>(pos)
+        } else if order == 21 {
+            self.accel_and_partials_t::<21, 25>(pos)
+        } else if order == 22 {
+            self.accel_and_partials_t::<22, 26>(pos)
+        } else if order == 23 {
+            self.accel_and_partials_t::<23, 27>(pos)
+        } else if order == 24 {
+            self.accel_and_partials_t::<24, 28>(pos)
+        } else if order == 25 {
+            self.accel_and_partials_t::<25, 29>(pos)
+        } else if order == 26 {
+            self.accel_and_partials_t::<26, 30>(pos)
+        } else if order == 27 {
+            self.accel_and_partials_t::<27, 31>(pos)
+        } else if order == 28 {
+            self.accel_and_partials_t::<28, 32>(pos)
+        } else if order == 29 {
+            self.accel_and_partials_t::<29, 33>(pos)
+        } else if order == 30 {
+            self.accel_and_partials_t::<30, 34>(pos)
+        } else if order == 31 {
+            self.accel_and_partials_t::<31, 35>(pos)
+        } else if order == 32 {
+            self.accel_and_partials_t::<32, 36>(pos)
+        } else if order == 33 {
+            self.accel_and_partials_t::<33, 37>(pos)
+        } else if order == 34 {
+            self.accel_and_partials_t::<34, 38>(pos)
+        } else if order == 35 {
+            self.accel_and_partials_t::<35, 39>(pos)
+        } else if order == 36 {
+            self.accel_and_partials_t::<36, 40>(pos)
+        } else if order == 37 {
+            self.accel_and_partials_t::<37, 15>(pos)
+        } else if order == 38 {
+            self.accel_and_partials_t::<38, 16>(pos)
+        } else if order == 39 {
+            self.accel_and_partials_t::<39, 17>(pos)
+        } else {
+            self.accel_and_partials_t::<40, 20>(pos)
         }
     }
 

{satkit-0.3.0 → satkit-0.3.1}/src/ode/rk_adaptive.rs

@@ -18,21 +18,33 @@ pub trait RKAdaptive<const N: usize, const NI: usize> {
     /// (first compute of next iteration is same as last compute of last iteration)
     const FSAL: bool;
 
+    /// Interpolate densely calculated solution onto
+    /// values that are evenly spaced in "x"
+    ///
     fn interpolate<S: ODEState>(xinterp: f64, sol: &ODESolution<S>) -> ODEResult<S> {
         if sol.dense.is_none() {
             return Err(Box::new(ODEError::NoDenseOutputInSolution));
         }
+        if sol.x < xinterp {
+            return Err(Box::new(ODEError::InterpExceedsSolutionBounds));
+        }
         let dense = sol.dense.as_ref().unwrap();
-        if xinterp < dense.x[0] || xinterp > sol.x {
+        if sol.x < dense.x[0] {
             return Err(Box::new(ODEError::InterpExceedsSolutionBounds));
         }
+
+        // We know indices are monotonically increasing, so only search from
+        // last found position in the array forward
         let mut idx = match dense.x.iter().position(|x| *x >= xinterp) {
             Some(v) => v,
-            None => dense.x.len() - 1,
+            None => dense.x.len(),
         };
         if idx > 0 {
             idx -= 1;
         }
+
+        // t is fractional distance beween x at idx and idx+1
+        // and is in range [0,1]
         let t = (xinterp - dense.x[idx]) / dense.h[idx];
 
         // Compute interpolant coefficient as funciton of t
@@ -47,6 +59,7 @@ pub trait RKAdaptive<const N: usize, const NI: usize> {
         let bi: Vec<f64> = Self::BI
             .iter()
             .map(|biarr| {
+                // Coefficients multiply increasing powers of t
                 let mut tj = 1.0;
                 biarr.iter().fold(0.0, |acc, bij| {
                     tj = tj * t;
@@ -60,14 +73,15 @@ pub trait RKAdaptive<const N: usize, const NI: usize> {
         //
         // This is equation(5) of:
         // https://link.springer.com/article/10.1023/A:1021190918665
-        let yarr = dense.yprime[idx]
+        //
+        let mut y = dense.yprime[idx]
             .iter()
             .enumerate()
             .fold(dense.y[idx].clone() / dense.h[idx], |acc, (ix, k)| {
                 acc + k.clone() * bi[ix]
             });
-
-        Ok(yarr * dense.h[idx])
+        y = y * dense.h[idx];
+        Ok(y)
     }
 
     ///

{satkit-0.3.0 → satkit-0.3.1}/src/ode/rkv87.rs

@@ -77,8 +77,8 @@ mod tests {
         (0..100).for_each(|idx| {
             let x = idx as f64 * PI / testcount as f64;
             let interp = RKV87::interpolate(x, &sol).unwrap();
-            assert!((interp[0] - x.cos()).abs() < 1e-10);
-            assert!((interp[1] + x.sin()).abs() < 1e-10);
+            assert!((interp[0] - x.cos()).abs() < 1e-12);
+            assert!((interp[1] + x.sin()).abs() < 1e-12);
         });
 
         Ok(())

{satkit-0.3.0 → satkit-0.3.1}/src/ode/rkv98.rs

@@ -78,8 +78,8 @@ mod tests {
         (0..100).for_each(|idx| {
             let x = idx as f64 * PI / testcount as f64;
             let interp = RKV98::interpolate(x, &sol).unwrap();
-            assert!((interp[0] - x.cos()).abs() < 1e-11);
-            assert!((interp[1] + x.sin()).abs() < 1e-11);
+            assert!((interp[0] - x.cos()).abs() < 1e-12);
+            assert!((interp[1] + x.sin()).abs() < 1e-12);
         });
 
         Ok(())

{satkit-0.3.0 → satkit-0.3.1}/src/orbitprop/propagator.rs

@@ -752,10 +752,15 @@ mod tests {
             })
             .collect();
 
-        let v0 = na::vector![2.47130555e03, 2.94682777e03, -5.34171918e02, 2.13018578e-02];
+        let v0 = na::vector![
+            2.47130562e+03,
+            2.94682753e+03,
+            -5.34172176e+02,
+            2.32565692e-02
+        ];
 
         let state0 = na::vector![pgcrf[0][0], pgcrf[0][1], pgcrf[0][2], v0[0], v0[1], v0[2]];
-        let satprops: SatPropertiesStatic = SatPropertiesStatic::new(0.0, 2.13018578e-2);
+        let satprops: SatPropertiesStatic = SatPropertiesStatic::new(0.0, v0[3]);
         let settings = PropSettings::default();
 
         let res = propagate(
@@ -769,10 +774,10 @@ mod tests {
 
         // We've propagated over a day; assert that the difference in position on all three coordinate axes
         // is less than 10 meters for all 5-minute intervals
-        for iv in 0..(pgcrf.len() - 5) {
+        for iv in 0..(pgcrf.len()) {
             let interp_state = res.interp(&times[iv])?;
             for ix in 0..3 {
-                assert!((pgcrf[iv][ix] - interp_state[ix]).abs() < 10.0);
+                assert!((pgcrf[iv][ix] - interp_state[ix]).abs() < 8.0);
             }
         }
 

{satkit-0.3.0 → satkit-0.3.1}/src/pybindings/pypropagate.rs

@@ -77,7 +77,6 @@ use pyo3::types::{PyDict, PyString};
 ///         and possibly dense output for interpolation
 ///
 ///
-///
 #[pyfunction(signature=(pos, vel, start, **kwargs))]
 pub fn propagate(
     pos: &Bound<'_, np::PyArray1<f64>>,

{satkit-0.3.0 → satkit-0.3.1}/src/pybindings/pysgp4.rs

@@ -4,8 +4,8 @@ use pyo3::types::{PyDict, PyList};
 use super::pyastrotime::ToTimeVec;
 use super::pytle::PyTLE;
 use crate::sgp4 as psgp4;
-use numpy::{PyArray, PyArray1};
 use numpy::PyArrayMethods;
+use numpy::{PyArray, PyArray1};
 
 // Thin Python wrapper around SGP4 Error
 #[allow(non_camel_case_types)]
@@ -71,8 +71,6 @@ impl From<psgp4::SGP4Error> for PySGP4Error {
     }
 }
 
-
-
 /// """SGP-4 propagator for TLE
 ///    
 /// Note:
@@ -95,7 +93,7 @@ impl From<psgp4::SGP4Error> for PySGP4Error {
 /// Returns:
 ///     tuple[npt.ArrayLike[np.float64], npt.ArrayLike[np.float64]]: position and velocity in meters and meters/second, respectively, in the TEME frame at each of the "Ntime" input times and each of the "Ntle" tles
 ///
-/// 
+///
 /// Example:
 ///
 ///
@@ -117,13 +115,17 @@ impl From<psgp4::SGP4Error> for PySGP4Error {
 /// >>>
 /// >>> # convert to ITRF coordinate object
 /// >>> coord = satkit.itrfcoord.from_vector(pitrf)
-/// >>> 
+/// >>>
 /// >>> # Print ITRF coordinate object location
 /// >>> print(coord)
 /// ITRFCoord(lat:  -0.0363 deg, lon:  -2.2438 deg, hae: 35799.51 km)
 #[pyfunction]
 #[pyo3(signature=(tle, time, **kwds))]
-pub fn sgp4(tle: &Bound<'_, PyAny>, time: &Bound<'_, PyAny>, kwds: Option<&Bound<'_, PyDict>>) -> PyResult<PyObject> {
+pub fn sgp4(
+    tle: &Bound<'_, PyAny>,
+    time: &Bound<'_, PyAny>,
+    kwds: Option<&Bound<'_, PyDict>>,
+) -> PyResult<PyObject> {
     let mut output_err = false;
     let mut opsmode: OpsMode = OpsMode::afspc;
     let mut gravconst: GravConst = GravConst::wgs72;
@@ -161,11 +163,13 @@ pub fn sgp4(tle: &Bound<'_, PyAny>, time: &Bound<'_, PyAny>, kwds: Option<&Bound
             // hence the switch
             if output_err == false {
                 Ok((
-                    PyArray1::from_slice_bound(py, r.data.as_slice()).as_gil_ref()
+                    PyArray1::from_slice_bound(py, r.data.as_slice())
+                        .as_gil_ref()
                         .reshape(dims.clone())
                         .unwrap()
                         .to_object(py),
-                    PyArray1::from_slice_bound(py, v.data.as_slice()).as_gil_ref()
+                    PyArray1::from_slice_bound(py, v.data.as_slice())
+                        .as_gil_ref()
                         .reshape(dims)
                         .unwrap()
                         .to_object(py),
@@ -173,10 +177,12 @@ pub fn sgp4(tle: &Bound<'_, PyAny>, time: &Bound<'_, PyAny>, kwds: Option<&Bound
                     .to_object(py))
             } else {
                 Ok((
-                    PyArray1::from_slice_bound(py, r.data.as_slice()).as_gil_ref()
+                    PyArray1::from_slice_bound(py, r.data.as_slice())
+                        .as_gil_ref()
                         .reshape(dims.clone())
                         .unwrap(),
-                    PyArray1::from_slice_bound(py, v.data.as_slice()).as_gil_ref()
+                    PyArray1::from_slice_bound(py, v.data.as_slice())
+                        .as_gil_ref()
                         .reshape(dims)
                         .unwrap(),
                     PyArray::from_owned_object_array_bound(
@@ -200,9 +206,9 @@ pub fn sgp4(tle: &Bound<'_, PyAny>, time: &Bound<'_, PyAny>, kwds: Option<&Bound
 
         pyo3::Python::with_gil(|py| -> PyResult<PyObject> {
             let n = tles.len() * tmarray.len() * 3;
-            
+
             let parr = PyArray1::zeros_bound(py, [n], false);
-            let varr= PyArray1::zeros_bound(py, [n], false);
+            let varr = PyArray1::zeros_bound(py, [n], false);
             let ntimes = tmarray.len();
 
             // I'd prefer to create this uninitialized, which would probably be a bit faster,