Merge pull request #62 from Gromwud/main

Sparsity bugfix

Author: Gromwud

epde/operators/common/sparsity.py

@@ -73,8 +73,8 @@ def fit(self, X, y, sample_weights):
                 rho = np.dot(X[:, j], residual) + old_coef * norm_sq
 
                 # Use CV-based Thresholding
-                # threshold = cv[j] * y_sq_sum
-                threshold = cv[j] * self.n_samples
+                threshold = cv[j] * y_sq_sum
+                # threshold = cv[j] * self.n_samples
                 # threshold = cv[j] * norm_sq * abs(old_coef)
                 new_coef = self._soft_threshold(rho, threshold) / norm_sq