fixing tests with logpolar

Author: dkazanc

zenodo-tests/test_recon/test_algorithm.py

@@ -1,3 +1,4 @@
+
 import cupy as cp
 import numpy as np
 import pytest
@@ -70,7 +71,7 @@ def test_reconstruct_FBP3d_tomobar_i12_dataset1(i12_dataset1):
     assert recon_data.shape == (2560, 50, 2560)
 
 
-def test_reconstruct_LPRec_i12_dataset1(i12_dataset1):
+def test_reconstruct_LPRec3d_tomobar_i12_dataset1(i12_dataset1):
     force_clean_gpu_memory()
     projdata = i12_dataset1[0]
     angles = i12_dataset1[1]
@@ -92,12 +93,14 @@ def test_reconstruct_LPRec_i12_dataset1(i12_dataset1):
     )
     assert recon_data.flags.c_contiguous
     recon_data = recon_data.get()
-    assert_allclose(np.sum(recon_data), 8973.761, rtol=1e-07, atol=1e-6)
+    assert isclose(np.sum(recon_data), 8973.755, abs_tol=10**-3)
+    assert isclose(np.max(recon_data), 0.006396, abs_tol=10**-6)
+    assert isclose(np.min(recon_data), -0.006159, abs_tol=10**-6)
     assert recon_data.dtype == np.float32
     assert recon_data.shape == (2560, 3, 2560)
 
 
-def test_reconstruct_LP_REC_i13_dataset1(i13_dataset1):
+def test_reconstruct_LPRec_tomobar_i13_dataset1(i13_dataset1):
     force_clean_gpu_memory()
     projdata = i13_dataset1[0]
     angles = i13_dataset1[1]
@@ -134,7 +137,9 @@ def test_reconstruct_LP_REC_i13_dataset1(i13_dataset1):
 
     assert recon_data.flags.c_contiguous
     recon_data = recon_data.get()
-    assert isclose(np.sum(recon_data), 1241.859, abs_tol=10**-3)
+    assert isclose(np.sum(recon_data), 1241.8612, abs_tol=10**-4)
+    assert isclose(np.max(recon_data), 0.008238, abs_tol=10**-6)
+    assert isclose(np.min(recon_data), -0.006565, abs_tol=10**-6)
     assert recon_data.dtype == np.float32
     assert recon_data.shape == (4646, 1, 4646)
 
@@ -200,7 +205,9 @@ def test_reconstruct_LPRec3d_tomobar_i13_dataset2(i13_dataset2):
     assert recon_data.flags.c_contiguous
     recon_data = recon_data.get()
 
-    assert isclose(np.sum(recon_data), 4095.577, abs_tol=10**-3)
+    assert isclose(np.sum(recon_data), 4095.6272, abs_tol=10**-4)
+    assert isclose(np.max(recon_data), 0.010567, abs_tol=10**-6)
+    assert isclose(np.min(recon_data), -0.008396, abs_tol=10**-6)    
     assert recon_data.dtype == np.float32
     assert recon_data.shape == (2560, 10, 2560)