added Atol and Btol in scipy lsqr for all tests (#486)

* atol botl fixed to 1e-8 in basic op test

* atol botl fixed to 1e-8 in test_convolve

* atol botl fixed to 1e-8 in test_diagonal

* atol botl fixed to 1e-8 in test_dwts

* atol botl fixed to 1e-8 in test_ffts

* atol botl fixed to 1e-8 in test_fredholm

* atol botl fixed to 1e-8 in test_functionoperator

* atol botl fixed to 1e-8 in test_kronecker

* atol botl fixed to 1e-8 in test_nonstatconvolve

* atol botl fixed to 1e-8 in test_shift

* atol botl fixed to 1e-8 in test_smoothing

---------

Co-authored-by: aniket singh rawat <[email protected]>

Author: dikwickley

pytests/test_basicoperators.py

@@ -63,7 +63,7 @@ def test_LinearRegression(par):
     assert dottest(LRop, par["ny"], 2)
 
     x = np.array([1.0, 2.0], dtype=np.float32)
-    xlsqr = lsqr(LRop, LRop * x, damp=1e-10, iter_lim=300, show=0)[0]
+    xlsqr = lsqr(LRop, LRop * x, damp=1e-10, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=3)
 
     y = LRop * x
@@ -82,7 +82,7 @@ def test_MatrixMult(par):
     assert dottest(Gop, par["ny"], par["nx"], complexflag=0 if par["imag"] == 0 else 3)
 
     x = np.ones(par["nx"]) + par["imag"] * np.ones(par["nx"])
-    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, show=0)[0]
+    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=4)
 
 
@@ -100,7 +100,7 @@ def test_MatrixMult_sparse(par):
     assert dottest(Gop, par["ny"], par["nx"], complexflag=0 if par["imag"] == 1 else 3)
 
     x = np.ones(par["nx"]) + par["imag"] * np.ones(par["nx"])
-    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, show=0)[0]
+    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=4)
 
 
@@ -133,7 +133,7 @@ def test_MatrixMult_repeated(par):
     )
 
     x = (np.ones((par["nx"], 5)) + par["imag"] * np.ones((par["nx"], 5))).ravel()
-    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, show=0)[0]
+    xlsqr = lsqr(Gop, Gop * x, damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=4)
 
 

pytests/test_convolve.py

@@ -133,7 +133,7 @@ def test_Convolve1D(par):
 
         x = np.zeros((par["nx"]))
         x[par["nx"] // 2] = 1.0
-        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, show=0)[0]
+        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, atol=1e-8, btol=1e-8, show=0)[0]
         assert_array_almost_equal(x, xlsqr, decimal=1)
 
     # 1D on 2D
@@ -153,7 +153,7 @@ def test_Convolve1D(par):
             int(par["nx"] / 2 - 3) : int(par["nx"] / 2 + 3),
         ] = 1.0
         x = x.ravel()
-        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, show=0)[0]
+        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, atol=1e-8, btol=1e-8, show=0)[0]
         assert_array_almost_equal(x, xlsqr, decimal=1)
 
     # 1D on 3D
@@ -175,7 +175,7 @@ def test_Convolve1D(par):
         int(par["nx"] / 2 - 3) : int(par["nx"] / 2 + 3),
     ] = 1.0
     x = x.ravel()
-    xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, show=0)[0]
+    xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=1)
 
 
@@ -200,7 +200,7 @@ def test_Convolve2D(par):
             int(par["nx"] / 2 - 3) : int(par["nx"] / 2 + 3),
         ] = 1.0
         x = x.ravel()
-        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, show=0)[0]
+        xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, atol=1e-8, btol=1e-8, show=0)[0]
         assert_array_almost_equal(x, xlsqr, decimal=1)
 
     # 2D on 3D
@@ -224,7 +224,7 @@ def test_Convolve2D(par):
         int(par["nx"] / 2 - 3) : int(par["nx"] / 2 + 3),
     ] = 1.0
     x = x.ravel()
-    xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, show=0)[0]
+    xlsqr = lsqr(Cop, Cop * x, damp=1e-20, iter_lim=200, atol=1e-8, btol=1e-8, show=0)[0]
     # due to ringing in solution we cannot use assert_array_almost_equal
     assert np.linalg.norm(xlsqr - x) / np.linalg.norm(xlsqr) < 2e-1
 
@@ -251,7 +251,7 @@ def test_Convolve3D(par):
     ] = 1.0
     x = x.ravel()
     y = Cop * x
-    xlsqr = lsqr(Cop, y, damp=1e-20, iter_lim=400, show=0)[0]
+    xlsqr = lsqr(Cop, y, damp=1e-20, iter_lim=400, atol=1e-8, btol=1e-8, show=0)[0]
     # due to ringing in solution we cannot use assert_array_almost_equal
     assert np.linalg.norm(xlsqr - x) / np.linalg.norm(xlsqr) < 2e-1
 

pytests/test_diagonal.py

@@ -23,7 +23,7 @@ def test_Diagonal_1dsignal(par):
         assert dottest(Dop, ddim, ddim, complexflag=0 if par["imag"] == 0 else 3)
 
         x = np.ones(ddim) + par["imag"] * np.ones(ddim)
-        xlsqr = sp_lsqr(Dop, Dop * x, damp=1e-20, iter_lim=300, show=0)[0]
+        xlsqr = sp_lsqr(Dop, Dop * x, damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x, xlsqr, decimal=4)
 
@@ -45,7 +45,7 @@ def test_Diagonal_2dsignal(par):
         x = np.ones((par["nx"], par["nt"])) + par["imag"] * np.ones(
             (par["nx"], par["nt"])
         )
-        xlsqr = sp_lsqr(Dop, Dop * x.ravel(), damp=1e-20, iter_lim=300, show=0)[0]
+        xlsqr = sp_lsqr(Dop, Dop * x.ravel(), damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x.ravel(), xlsqr.ravel(), decimal=4)
 
@@ -69,7 +69,7 @@ def test_Diagonal_3dsignal(par):
         x = np.ones((par["ny"], par["nx"], par["nt"])) + par["imag"] * np.ones(
             (par["ny"], par["nx"], par["nt"])
         )
-        xlsqr = sp_lsqr(Dop, Dop * x.ravel(), damp=1e-20, iter_lim=300, show=0)[0]
+        xlsqr = sp_lsqr(Dop, Dop * x.ravel(), damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x.ravel(), xlsqr.ravel(), decimal=4)
 
@@ -91,7 +91,7 @@ def test_Diagonal_2dsignal_unflattened(par):
         x = np.ones((par["nx"], par["nt"])) + par["imag"] * np.ones(
             (par["nx"], par["nt"])
         )
-        xlsqr = lsqr(Dop, Dop * x, damp=1e-20, niter=300, show=0)[0]
+        xlsqr = lsqr(Dop, Dop * x, damp=1e-20, niter=300, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x, xlsqr, decimal=4)
 
@@ -115,6 +115,6 @@ def test_Diagonal_3dsignal_unflattened(par):
         x = np.ones((par["ny"], par["nx"], par["nt"])) + par["imag"] * np.ones(
             (par["ny"], par["nx"], par["nt"])
         )
-        xlsqr = lsqr(Dop, Dop * x, damp=1e-20, niter=300, show=0)[0]
+        xlsqr = lsqr(Dop, Dop * x, damp=1e-20, niter=300, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x, xlsqr, decimal=4)

pytests/test_dwts.py

@@ -33,7 +33,7 @@ def test_DWT_1dsignal(par):
 
     y = DWTop * x
     xadj = DWTop.H * y  # adjoint is same as inverse for dwt
-    xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+    xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     assert_array_almost_equal(x, xadj, decimal=8)
     assert_array_almost_equal(x, xinv, decimal=8)
@@ -57,7 +57,7 @@ def test_DWT_2dsignal(par):
 
         y = DWTop * x.ravel()
         xadj = DWTop.H * y  # adjoint is same as inverse for dwt
-        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x.ravel(), xadj, decimal=8)
         assert_array_almost_equal(x.ravel(), xinv, decimal=8)
@@ -83,7 +83,7 @@ def test_DWT_3dsignal(par):
 
         y = DWTop * x.ravel()
         xadj = DWTop.H * y  # adjoint is same as inverse for dwt
-        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x.ravel(), xadj, decimal=8)
         assert_array_almost_equal(x.ravel(), xinv, decimal=8)
@@ -103,7 +103,7 @@ def test_DWT2D_2dsignal(par):
 
     y = DWTop * x.ravel()
     xadj = DWTop.H * y  # adjoint is same as inverse for dwt
-    xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+    xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     assert_array_almost_equal(x.ravel(), xadj, decimal=8)
     assert_array_almost_equal(x.ravel(), xinv, decimal=8)
@@ -129,7 +129,7 @@ def test_DWT2D_3dsignal(par):
 
         y = DWTop * x.ravel()
         xadj = DWTop.H * y  # adjoint is same as inverse for dwt
-        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+        xinv = lsqr(DWTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         assert_array_almost_equal(x.ravel(), xadj, decimal=8)
         assert_array_almost_equal(x.ravel(), xinv, decimal=8)

pytests/test_ffts.py

@@ -269,7 +269,7 @@ def test_FFT_random_real(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -403,7 +403,7 @@ def test_FFT_random_complex(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -464,7 +464,7 @@ def test_FFT2D_random_real(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -544,7 +544,7 @@ def test_FFT2D_random_complex(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -606,7 +606,7 @@ def test_FFTND_random_real(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -686,7 +686,7 @@ def test_FFTND_random_complex(par):
 
     # Ensure inverse and adjoint recover x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=0, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xadj, decimal=decimal)
     assert_array_almost_equal(x, xinv, decimal=decimal)
 
@@ -850,7 +850,7 @@ def test_FFT_1dsignal(par):
 
     y = FFTop * x
     xadj = FFTop.H * y  # adjoint is same as inverse for fft
-    xinv = lsqr(FFTop, y, damp=1e-10, iter_lim=10, show=0)[0]
+    xinv = lsqr(FFTop, y, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     # check all signal if nt>nfft and only up to nfft if nfft<nt
     imax = par["nt"] if par["nfft"] is None else min([par["nt"], par["nfft"]])
@@ -874,7 +874,7 @@ def test_FFT_1dsignal(par):
         assert_array_almost_equal(y_fftshift, np.fft.fftshift(y))
 
         xadj = FFTop_fftshift.H * y_fftshift  # adjoint is same as inverse for fft
-        xinv = lsqr(FFTop_fftshift, y_fftshift, damp=1e-10, iter_lim=10, show=0)[0]
+        xinv = lsqr(FFTop_fftshift, y_fftshift, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
         assert_array_almost_equal(x[:imax], xadj[:imax], decimal=decimal)
         assert_array_almost_equal(x[:imax], xinv[:imax], decimal=decimal)
 
@@ -930,7 +930,7 @@ def test_FFT_2dsignal(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is same as inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj.reshape(nt, nx))
     dinv = np.real(dinv.reshape(nt, nx))
@@ -958,7 +958,7 @@ def test_FFT_2dsignal(par):
         assert_array_almost_equal(D_fftshift, D2)
 
         dadj = FFTop_fftshift.H * D_fftshift  # adjoint is same as inverse for fft
-        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, show=0)[0]
+        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         dadj = np.real(dadj.reshape(nt, nx))
         dinv = np.real(dinv.reshape(nt, nx))
@@ -988,7 +988,7 @@ def test_FFT_2dsignal(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj.reshape(nt, nx))
     dinv = np.real(dinv.reshape(nt, nx))
@@ -1016,7 +1016,7 @@ def test_FFT_2dsignal(par):
         assert_array_almost_equal(D_fftshift, D2)
 
         dadj = FFTop_fftshift.H * D_fftshift  # adjoint is same as inverse for fft
-        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, show=0)[0]
+        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         dadj = np.real(dadj.reshape(nt, nx))
         dinv = np.real(dinv.reshape(nt, nx))
@@ -1085,7 +1085,7 @@ def test_FFT_3dsignal(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is same as inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj.reshape(nt, nx, ny))
     dinv = np.real(dinv.reshape(nt, nx, ny))
@@ -1125,7 +1125,7 @@ def test_FFT_3dsignal(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj.reshape(nt, nx, ny))
     dinv = np.real(dinv.reshape(nt, nx, ny))
@@ -1165,7 +1165,7 @@ def test_FFT_3dsignal(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj.reshape(nt, nx, ny))
     dinv = np.real(dinv.reshape(nt, nx, ny))
@@ -1193,7 +1193,7 @@ def test_FFT_3dsignal(par):
         assert_array_almost_equal(D_fftshift, D2)
 
         dadj = FFTop_fftshift.H * D_fftshift  # adjoint is same as inverse for fft
-        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, show=0)[0]
+        dinv = lsqr(FFTop_fftshift, D_fftshift, damp=1e-10, iter_lim=10, atol=1e-8, btol=1e-8, show=0)[0]
 
         dadj = np.real(dadj.reshape(nt, nx, ny))
         dinv = np.real(dinv.reshape(nt, nx, ny))
@@ -1250,7 +1250,7 @@ def test_FFT2D(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj).reshape(par["nt"], par["nx"])
     dinv = np.real(dinv).reshape(par["nt"], par["nx"])
@@ -1296,7 +1296,7 @@ def test_FFT2D(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj).reshape(par["nt"], par["nx"])
     dinv = np.real(dinv).reshape(par["nt"], par["nx"])
@@ -1356,7 +1356,7 @@ def test_FFT3D(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj).reshape(par["nt"], par["nx"], par["ny"])
     dinv = np.real(dinv).reshape(par["nt"], par["nx"], par["ny"])
@@ -1407,7 +1407,7 @@ def test_FFT3D(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj).reshape(par["nt"], par["nx"], par["ny"])
     dinv = np.real(dinv).reshape(par["nt"], par["nx"], par["ny"])
@@ -1454,7 +1454,7 @@ def test_FFT3D(par):
 
     D = FFTop * d.ravel()
     dadj = FFTop.H * D  # adjoint is inverse for fft
-    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, show=0)[0]
+    dinv = lsqr(FFTop, D, damp=1e-10, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
 
     dadj = np.real(dadj).reshape(par["nt"], par["nx"], par["ny"])
     dinv = np.real(dinv).reshape(par["nt"], par["nx"], par["ny"])

pytests/test_fredholm.py

@@ -96,6 +96,6 @@ def test_Fredholm1(par):
         par["nsl"] * par["ny"] * par["nz"],
         complexflag=0 if par["imag"] == 0 else 3,
     )
-    xlsqr = lsqr(Fop, Fop * x.ravel(), damp=1e-20, iter_lim=30, show=0)[0]
+    xlsqr = lsqr(Fop, Fop * x.ravel(), damp=1e-20, iter_lim=30, atol=1e-8, btol=1e-8, show=0)[0]
     xlsqr = xlsqr.reshape(par["nsl"], par["ny"], par["nz"])
     assert_array_almost_equal(x, xlsqr, decimal=3)

pytests/test_functionoperator.py

@@ -77,7 +77,7 @@ def adjoint_f(y):
 
     # Only test inversion for square or overdetermined systems
     if par["nc"] <= par["nr"]:
-        xlsqr = lsqr(Fop, F_x, damp=0, iter_lim=100, show=0)[0]
+        xlsqr = lsqr(Fop, F_x, damp=0, iter_lim=100, atol=1e-8, btol=1e-8, show=0)[0]
         assert_array_almost_equal(x, xlsqr, decimal=4)
 
 

pytests/test_kronecker.py

@@ -29,7 +29,7 @@ def test_Kroneker(par):
         Kop, par["ny"] ** 2, par["nx"] ** 2, complexflag=0 if par["imag"] == 0 else 3
     )
 
-    xlsqr = lsqr(Kop, Kop * x, damp=1e-20, iter_lim=300, show=0)[0]
+    xlsqr = lsqr(Kop, Kop * x, damp=1e-20, iter_lim=300, atol=1e-8, btol=1e-8, show=0)[0]
     assert_array_almost_equal(x, xlsqr, decimal=2)
 
     # Comparison with numpy