Using dspharm_full instead of access by indices

corbass/inversion.py

@@ -242,7 +242,7 @@ class Inversion:
         self.cor_obs_full = np.copy(self.cov_NEZ_obs[:, self.idx_NEZ_full])
 
         self.mu_NEZ_full = np.dot(self.g_bar,
-                                  self.dspharm_obs[:, self.idx_NEZ_full])
+                                  dspharm_full)
 
         self.mu_NEZ_obs = np.dot(self.g_bar, self.dspharm_obs) \
             + np.linalg.multi_dot([self.cor_obs_full,
@@ -254,7 +254,7 @@ class Inversion:
                                                  self.cor_obs_full.T])
 
         self.R_obs = self.dspharm_obs \
-            - np.linalg.multi_dot((self.dspharm_obs[:, self.idx_NEZ_full],
+            - np.linalg.multi_dot((dspharm_full,
                                    self.prc_NEZ_full,
                                    self.cor_obs_full.T))