set minimum eps based on tol

evalf · Jan 5, 2021 · d09c087 · d09c087
1 parent 240c81f
commit d09c087
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/nutils/topology.py b/nutils/topology.py
@@ -536,7 +536,7 @@ def locate(self, geom, coords, *, tol, eps=0, maxiter=100, arguments=None, weigh
               xi += numpy.linalg.solve(J_xi, coord - coord_xi)
             except numpy.linalg.LinAlgError:
               break
-          if converged and ref.inside(xi, eps=eps):
+          if converged and ref.inside(xi, eps=max(eps, *tol/abs(numpy.linalg.eigvals(J_xi)))):
             ielems[ipoint] = ielem
             xis[ipoint] = xi
             break