SNC was absorbing the Cr changes!

Author: ChristopherRabotin

src/dynamics/spacecraft.rs

@@ -419,11 +419,8 @@ impl Dynamics for SpacecraftDynamics {
             // Add this force model's estimation if applicable.
             if let Some(idx) = model.estimation_index() {
                 for j in 0..3 {
-                    grad[(idx, j)] += model_grad[(3, j)] / total_mass;
+                    grad[(j + 3, idx)] += model_grad[(3, j)] / total_mass;
                 }
-                // grad[(idx, j)] += model_grad[(3, j)] / total_mass;
-                // Cr = 1.4999976638233514 +/-0.20000001212152405
-                // RMAG error = 0.815728 m VMAG error = 1.346671 mm/s
             }
         }
 

tests/orbit_determination/spacecraft.rs

@@ -403,9 +403,9 @@ fn od_val_sc_srp_estimation(
     let initial_estimate = sc.to_estimate().unwrap();
 
     // let ckf = KF::no_snc(initial_estimate);
-    let ckf = KF::new(
+    let ckf = KF::no_snc(
         initial_estimate,
-        SNC3::from_diagonal(2 * Unit::Minute, &[1e-12, 1e-12, 1e-12]),
+        // SNC3::from_diagonal(2 * Unit::Minute, &[1e-14, 1e-14, 1e-14]),
     );
 
     let mut odp = ODProcess::ekf(