diff --git a/meanas/test/test_fdfd.py b/meanas/test/test_fdfd.py
index e9f3a21..533ef91 100644
--- a/meanas/test/test_fdfd.py
+++ b/meanas/test/test_fdfd.py
@@ -96,7 +96,8 @@ def sim(request, shape, epsilon, dxes, j_distribution, omega, pec, pmc):
 
     j_vec = vec(j_distribution)
     eps_vec = vec(epsilon)
-    e_vec = fdfd.solvers.generic(J=j_vec, omega=omega, dxes=dxes, epsilon=eps_vec, matrix_solver_opts={'atol': 1e-12})
+    e_vec = fdfd.solvers.generic(J=j_vec, omega=omega, dxes=dxes, epsilon=eps_vec,
+                                 matrix_solver_opts={'atol': 1e-15, 'tol': 1e-10})
     e = unvec(e_vec, shape[1:])
 
     sim = SimResult(