satkit 0.3.4__tar.gz → 0.3.6__tar.gz

{satkit-0.3.4 → satkit-0.3.6}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "satkit"
-version = "0.3.3"
+version = "0.3.6"
 edition = "2021"
 description = "Satellite Toolkit"
 readme = "README.md"

{satkit-0.3.4/python/satkit.egg-info → satkit-0.3.6}/PKG-INFO

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

{satkit-0.3.4 → satkit-0.3.6}/pyproject.toml

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

{satkit-0.3.4 → satkit-0.3.6/python/satkit.egg-info}/PKG-INFO

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

{satkit-0.3.4 → satkit-0.3.6}/src/earthgravity.rs

@@ -6,7 +6,7 @@ use std::path::PathBuf;
 use nalgebra as na;
 type CoeffTable = na::DMatrix<f64>;
 
-type DivisorTable = na::SMatrix<f64, 20, 20>;
+type DivisorTable = na::SMatrix<f64, 44, 44>;
 
 use once_cell::sync::OnceCell;
 
@@ -171,7 +171,9 @@ impl Gravity {
     pub fn accel(&self, pos: &Vec3, order: usize) -> Vec3 {
         // This is tedious, but using generics allows for vectors to be
         // allocated on the stack, which is faster
-        if order == 2 {
+        if order == 1 {
+            self.accel_t::<1, 5>(pos)
+        } else if order == 2 {
             self.accel_t::<2, 6>(pos)
         } else if order == 3 {
             self.accel_t::<3, 7>(pos)
@@ -255,7 +257,9 @@ impl Gravity {
     pub fn accel_and_partials(&self, pos: &Vec3, order: usize) -> (Vec3, na::Matrix3<f64>) {
         // This is tedious, but using generics allows for vectors to be
         // allocated on the stack, which is faster
-        if order == 2 {
+        if order == 1 {
+            self.accel_and_partials_t::<1, 5>(pos)
+        } else if order == 2 {
             self.accel_and_partials_t::<2, 6>(pos)
         } else if order == 3 {
             self.accel_and_partials_t::<3, 7>(pos)
@@ -613,7 +617,7 @@ impl Gravity {
             coeffs: cs,
             divisor_table: {
                 let mut dt: DivisorTable = DivisorTable::zeros();
-                for m in 0..19 {
+                for m in 0..43 {
                     if m > 0 {
                         dt[(m, m)] = 2.0 * m as f64 - 1.0;
                     }