diff --git a/sectionproperties/analysis/section.py b/sectionproperties/analysis/section.py
index 195b6ff9..d3d36668 100644
--- a/sectionproperties/analysis/section.py
+++ b/sectionproperties/analysis/section.py
@@ -193,7 +193,6 @@ def __init__(
             Polygon(self.geometry.mesh["vertices"][tri][0:3])
             for tri in self.geometry.mesh["triangles"]
         ]
-        self.poly_mesh_idx = dict((id(poly), i) for i, poly in enumerate(p_mesh))
         self.mesh_search_tree = STRtree(p_mesh)
 
         # initialise class storing section properties
@@ -2259,11 +2258,7 @@ def get_stress_at_points(
 
         for pt in pts:
             query_geom = Point(pt)
-            tri_ids = [
-                self.poly_mesh_idx[id(poly)]
-                for poly in self.mesh_search_tree.query(query_geom)  # TODO
-                if poly.intersects(query_geom)
-            ]
+            tri_ids = self.mesh_search_tree.query(query_geom, predicate="intersects")
             if len(tri_ids) == 0:
                 sig = None
             elif len(tri_ids) == 1: