diff --git a/tests/linalg_test.py b/tests/linalg_test.py
index 58c876ee..5af84b52 100644
--- a/tests/linalg_test.py
+++ b/tests/linalg_test.py
@@ -99,6 +99,8 @@ def test_jack_multi_dot():
         tt = pe.linalg.jack_matmul(my_array, my_array, my_array) - pe.linalg.matmul(my_array, my_array, my_array)
 
         for t, e in np.ndenumerate(tt):
+            e.gamma_method()
+            assert e.is_zero_within_error(0.01)
             assert e.is_zero(atol=1e-1), t
             assert np.isclose(e.value, 0.0)