✨ Developing PerturbationMatrix class (#59)

* ✨ Developing PerturbationMatrix class

* ✅ fixing flake8

* ✨ added a 3D matrix

* 🎨 fairly good sketch

* Update src/psfmachine/perturbation.py

Co-authored-by: Jorge Martínez-Palomera <[email protected]>

* Update src/psfmachine/perturbation.py

* fixing identation and _clean_vectors when segments=True

* ⚡ faster!

* 🎨 adding pixel masking

* 📝 added docstrings

* ✨ Added PCA method

* 🎨 added fbpca explicitly as dependency

* 📝 docs

* 🐛 fixing pca

* 🐛 fix @jorgemarpa bug

* 📝 update docstrings

* Update src/psfmachine/perturbation.py

Co-authored-by: Jorge Martínez-Palomera <[email protected]>

* ✨ now pca has smoothed components

* 🐛 update pca

* 🎨 flake8

* 🎨 no comment block

* 📝 docstrings

* 🎨 updating sane defaults, adding smooth PCA components

* 🎨 sane defaults

* 📝 documentation

Co-authored-by: Jorge Martínez-Palomera <[email protected]>
Co-authored-by: Jorge <[email protected]>

Author: christinahedges

poetry.lock

@@ -887,7 +887,7 @@ def build_time_model(self, plot=False, downsample=False, split_time_model=False)
             dy,
             n_knots=self.n_time_knots,
             radius=self.time_radius,
-            spacing=self.cartesian_knot_spacing,
+            knot_spacing_type=self.cartesian_knot_spacing,
         )
         A2 = sparse.vstack([A_c] * time_binned.shape[0], format="csr")
         # Cartesian spline with time dependence
@@ -1049,7 +1049,7 @@ def plot_time_model(self, segment=0):
             dy,
             n_knots=self.n_time_knots,
             radius=self.time_radius,
-            spacing=self.cartesian_knot_spacing,
+            knot_spacing_type=self.cartesian_knot_spacing,
         )
         # if self.seg_splits.shape[0] == 2 and segment == 0:
         #     seg_mask = np.ones(time_binned.shape[0], dtype=bool)
@@ -1576,7 +1576,7 @@ def fit_model(self, fit_va=False):
                 dy,
                 n_knots=self.n_time_knots,
                 radius=self.time_radius,
-                spacing=self.cartesian_knot_spacing,
+                knot_spacing_type=self.cartesian_knot_spacing,
             )
             A_cp3 = sparse.hstack([A_cp, A_cp, A_cp, A_cp], format="csr")