Remove self.w

superdsm/dsm.py

@@ -337,7 +337,7 @@ def grad(self, params):
         grad = np.sum([term1 * q for q in self.q], axis=1)
         term1_sparse = coo_matrix(term1).transpose(copy=False)
         if self.smooth_mat.shape[1] > 0:
-            grad2  = (self.w.reshape(-1)[None, :] @ self.smooth_mat.multiply(term1_sparse)).reshape(-1)
+            grad2  = self.smooth_mat.multiply(term1_sparse).sum(axis=0).reshape(-1)
             grad2 += self.alpha * (params.ξ / self.term2)
             grad   = np.concatenate([grad, grad2])
         return grad