fix: change BlendingContinuous to work with pylops solvers

BlendingContinuous matvec/rmatvec used * and .H * for some
internal calls to operators, which does not allow running this
operator with pylops solvers as they disable ndarray_multiplication.
Changing this to matvec/rmatvec calls followed by reshape enables
the use of this operator with pylops solvers.

Author: mrava87

pylops/waveeqprocessing/blending.py

@@ -160,7 +160,11 @@ def _matvec_largerecs(self, x: NDArray) -> NDArray:
             if self.ShiftOps[i] is None:
                 blended_data[:, shift_int : shift_int + self.nt] += x[i, :, :]
             else:
-                shifted_data = self.ShiftOps[i] * self.PadOp * x[i, :, :]
+                shifted_data = (
+                    self.ShiftOps[i]
+                    .matvec(self.PadOp.matvec(x[i, :, :].ravel()))
+                    .reshape(self.ShiftOps[i].dimsd)
+                )
                 blended_data[:, shift_int : shift_int + self.nt + 1] += shifted_data
         return blended_data
 
@@ -172,11 +176,11 @@ def _rmatvec_largerecs(self, x: NDArray) -> NDArray:
             if self.ShiftOps[i] is None:
                 deblended_data[i, :, :] = x[:, shift_int : shift_int + self.nt]
             else:
-                shifted_data = (
-                    self.PadOp.H
-                    * self.ShiftOps[i].H
-                    * x[:, shift_int : shift_int + self.nt + 1]
-                )
+                shifted_data = self.PadOp.rmatvec(
+                    self.ShiftOps[i].rmatvec(
+                        x[:, shift_int : shift_int + self.nt + 1].ravel()
+                    )
+                ).reshape(self.PadOp.dims)
                 deblended_data[i, :, :] = shifted_data
         return deblended_data