minor: fixed some error messages

Author: mrava87

pylops/avo/poststack.py

@@ -94,14 +94,12 @@ def _PoststackLinearModelling(
             D = ncp.diag(0.5 * ncp.ones(nt0 - 1, dtype=dtype), k=1) - ncp.diag(
                 0.5 * ncp.ones(nt0 - 1, dtype=dtype), -1
             )
-            # D[0] = D[-1] = 0.
             D = inplace_set(ncp.array(0.0), D, 0)
             D = inplace_set(ncp.array(0.0), D, -1)
         else:
             D = ncp.diag(ncp.ones(nt0 - 1, dtype=dtype), k=1) - ncp.diag(
                 ncp.ones(nt0, dtype=dtype), k=0
             )
-            # D[-1] = 0.
             D = inplace_set(ncp.array(0.0), D, -1)
 
         # Create wavelet operator

pylops/jaxoperator.py

@@ -50,19 +50,19 @@ def _rmatvecad(self, x: JaxType, y: JaxType) -> JaxType:
         Parameters
         ----------
         x : :obj:`jaxlib.xla_extension.ArrayImpl`
-            Input array
+            Input array for forward
         y : :obj:`jaxlib.xla_extension.ArrayImpl`
-            Output array (where to store the
-            Vector-Jacobian product)
+            Input array for adjoint
 
         Returns
-        ----------
-        y : :obj:`jaxlib.xla_extension.ArrayImpl`
+        -------
+        xadj : :obj:`jaxlib.xla_extension.ArrayImpl`
             Output array
 
         """
         _, f_vjp = jax.vjp(self._matvec, x)
-        return jax.jit(f_vjp)(y)[0]
+        xadj = jax.jit(f_vjp)(y)[0]
+        return xadj
 
     def rmatvecad(self, x: JaxType, y: JaxType) -> JaxType:
         """Adjoint matrix-vector multiplication with AD
@@ -84,7 +84,9 @@ def rmatvecad(self, x: JaxType, y: JaxType) -> JaxType:
         M, N = self.shape
 
         if x.shape != (M,) and x.shape != (M, 1):
-            raise ValueError("dimension mismatch")
+            raise ValueError(
+                f"Dimension mismatch. Got {x.shape}, but expected {(M, 1)} or {(M,)}."
+            )
 
         y = self._rmatvecad(x, y)
 
@@ -93,5 +95,8 @@ def rmatvecad(self, x: JaxType, y: JaxType) -> JaxType:
         elif x.ndim == 2:
             y = y.reshape(N, 1)
         else:
-            raise ValueError("invalid shape returned by user-defined rmatvecad()")
+            raise ValueError(
+                f"Invalid shape returned by user-defined rmatvecad(). "
+                f"Expected 2-d ndarray or matrix, not {x.ndim}-d ndarray"
+            )
         return y

pylops/linearoperator.py

@@ -508,7 +508,9 @@ def matvec(self, x: NDArray) -> NDArray:
         M, N = self.shape
 
         if x.shape != (N,) and x.shape != (N, 1):
-            raise ValueError("dimension mismatch")
+            raise ValueError(
+                f"Dimension mismatch. Got {x.shape}, but expected {(M, 1)} or {(M,)}."
+            )
 
         y = self._matvec(x)
 
@@ -517,7 +519,7 @@ def matvec(self, x: NDArray) -> NDArray:
         elif x.ndim == 2:
             y = y.reshape(M, 1)
         else:
-            raise ValueError("invalid shape returned by user-defined matvec()")
+            raise ValueError("Invalid shape returned by user-defined matvec()")
         return y
 
     @count(forward=False)
@@ -542,7 +544,9 @@ def rmatvec(self, x: NDArray) -> NDArray:
         M, N = self.shape
 
         if x.shape != (M,) and x.shape != (M, 1):
-            raise ValueError("dimension mismatch")
+            raise ValueError(
+                f"Dimension mismatch. Got {x.shape}, but expected {(M, 1)} or {(M,)}."
+            )
 
         y = self._rmatvec(x)
 
@@ -551,7 +555,7 @@ def rmatvec(self, x: NDArray) -> NDArray:
         elif x.ndim == 2:
             y = y.reshape(N, 1)
         else:
-            raise ValueError("invalid shape returned by user-defined rmatvec()")
+            raise ValueError("Invalid shape returned by user-defined rmatvec()")
         return y
 
     @count(forward=True, matmat=True)
@@ -574,9 +578,9 @@ def matmat(self, X: NDArray) -> NDArray:
 
         """
         if X.ndim != 2:
-            raise ValueError("expected 2-d ndarray or matrix, " "not %d-d" % X.ndim)
+            raise ValueError(f"Expected 2-d ndarray or matrix, not {X.ndim}-d ndarray")
         if X.shape[0] != self.shape[1]:
-            raise ValueError("dimension mismatch: %r, %r" % (self.shape, X.shape))
+            raise ValueError(f"Dimension mismatch: {self.shape}, {X.shape}")
         Y = self._matmat(X)
         return Y
 
@@ -600,9 +604,9 @@ def rmatmat(self, X: NDArray) -> NDArray:
 
         """
         if X.ndim != 2:
-            raise ValueError("expected 2-d ndarray or matrix, " "not %d-d" % X.ndim)
+            raise ValueError(f"Expected 2-d ndarray or matrix, not {X.ndim}-d ndarray")
         if X.shape[0] != self.shape[0]:
-            raise ValueError("dimension mismatch: %r, %r" % (self.shape, X.shape))
+            raise ValueError(f"Dimension mismatch: {self.shape}, {X.shape}")
         Y = self._rmatmat(X)
         return Y