Increase atol

gsrp_smart_util.py

@@ -127,7 +127,7 @@ def add(mem1, mem2, cc, t_max, id1, id2, n_ind, mem_limit=np.infty):
     return out_val, out_tij
 
 
-def mask_val(mem, val, atol=5e-9):
+def mask_val(mem, val, atol=1e-8):
     mask = mem[0] >= val - atol
     return mem[0][mask], np.compress(mask, mem[1], axis=1)