diff --git a/pylops/signalprocessing/_fourierradon2d_numba.py b/pylops/signalprocessing/_fourierradon2d_numba.py
index e2139c6b..4657cd92 100755
--- a/pylops/signalprocessing/_fourierradon2d_numba.py
+++ b/pylops/signalprocessing/_fourierradon2d_numba.py
@@ -16,7 +16,6 @@ def _radon_inner_2d(X, Y, f, px, h, flim0, flim1, npx, nh):
         for ifr in range(flim0, flim1):
             for ipx in range(npx):
                 Y[ih, ifr] += X[ipx, ifr] * exp(-1j * 2 * pi * f[ifr] * px[ipx] * h[ih])
-    return Y
 
 
 @jit(nopython=True, parallel=parallel, nogil=True, cache=True, fastmath=True)
@@ -25,4 +24,3 @@ def _aradon_inner_2d(X, Y, f, px, h, flim0, flim1, npx, nh):
         for ifr in range(flim0, flim1):
             for ih in range(nh):
                 X[ipx, ifr] += Y[ih, ifr] * exp(1j * 2 * pi * f[ifr] * px[ipx] * h[ih])
-    return X
diff --git a/pylops/signalprocessing/_fourierradon3d_numba.py b/pylops/signalprocessing/_fourierradon3d_numba.py
index 43571db5..1ba1ef3a 100755
--- a/pylops/signalprocessing/_fourierradon3d_numba.py
+++ b/pylops/signalprocessing/_fourierradon3d_numba.py
@@ -24,7 +24,6 @@ def _radon_inner_3d(X, Y, f, py, px, hy, hx, flim0, flim1, npy, npx, nhy, nhx):
                             * f[ifr]
                             * (py[ipy] * hy[ihy] + px[ipx] * hx[ihx])
                         )
-    return Y
 
 
 @jit(nopython=True, parallel=parallel, nogil=True, cache=True, fastmath=True)
@@ -41,4 +40,3 @@ def _aradon_inner_3d(X, Y, f, py, px, hy, hx, flim0, flim1, npy, npx, nhy, nhx):
                             * f[ifr]
                             * (py[ipy] * hy[ihy] + px[ipx] * hx[ihx])
                         )
-    return X
diff --git a/pylops/signalprocessing/fourierradon2d.py b/pylops/signalprocessing/fourierradon2d.py
index 1d7b524c..83bc1d94 100755
--- a/pylops/signalprocessing/fourierradon2d.py
+++ b/pylops/signalprocessing/fourierradon2d.py
@@ -262,7 +262,7 @@ def _matvec_numba(self, x: NDArray) -> NDArray:
         y = np.zeros((self.nh, self.nfft2), dtype=self.cdtype)
 
         x = sp.fft.rfft(x, n=self.nfft, axis=-1)
-        y = _radon_inner_2d(
+        _radon_inner_2d(
             x,
             y,
             self.f,
@@ -281,7 +281,7 @@ def _rmatvec_numba(self, y: NDArray) -> NDArray:
         x = np.zeros((self.npx, self.nfft2), dtype=self.cdtype)
 
         y = sp.fft.rfft(y, n=self.nfft, axis=-1)
-        x = _aradon_inner_2d(
+        _aradon_inner_2d(
             x,
             y,
             self.f,
diff --git a/pylops/signalprocessing/fourierradon3d.py b/pylops/signalprocessing/fourierradon3d.py
index 8b1e9fea..135dece2 100755
--- a/pylops/signalprocessing/fourierradon3d.py
+++ b/pylops/signalprocessing/fourierradon3d.py
@@ -309,7 +309,7 @@ def _matvec_numba(self, x: NDArray) -> NDArray:
         y = np.zeros((self.nhy, self.nhx, self.nfft2), dtype=self.cdtype)
 
         x = sp.fft.rfft(x, n=self.nfft, axis=-1)
-        y = _radon_inner_3d(
+        _radon_inner_3d(
             x,
             y,
             self.f,
@@ -332,7 +332,7 @@ def _rmatvec_numba(self, y: NDArray) -> NDArray:
         x = np.zeros((self.npy, self.npx, self.nfft2), dtype=self.cdtype)
 
         y = sp.fft.rfft(y, n=self.nfft, axis=-1)
-        x = _aradon_inner_3d(
+        _aradon_inner_3d(
             x,
             y,
             self.f,