bug: fix bug when applying power of Op to CuPy arrays

Author: mrava87

pylops/linearoperator.py

@@ -1474,7 +1474,7 @@ def __init__(self, A: LinearOperator, p: int) -> None:
         self.args = (A, p)
 
     def _power(self, fun: Callable, x: NDArray) -> NDArray:
-        res = np.array(x, copy=True)
+        res = x.copy()
         for _ in range(self.args[1]):
             res = fun(res)
         return res