diff --git a/pyproject.toml b/pyproject.toml
index ea6c234b84..ffa81e3062 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -17,7 +17,7 @@ authors = [
 keywords = ["geophysics", "inverse problem"]
 
 dependencies = [
-    "discretize>=0.10.0",
+    "discretize>=0.11.0",
     "empymod>=2.0.0",
     "geoana>=0.5.0",
     "geoh5py>=0.11.0a1, <0.12.dev",
diff --git a/simpeg/maps.py b/simpeg/maps.py
index 3361fe725a..e2ebdd6699 100644
--- a/simpeg/maps.py
+++ b/simpeg/maps.py
@@ -6238,7 +6238,7 @@ def projection(self):
         Set the projection matrix with partial volumes
         """
         if getattr(self, "_projection", None) is None:
-            in_local = self.local_mesh._get_containing_cell_indexes(
+            in_local = self.local_mesh.get_containing_cells(
                 self.global_mesh.cell_centers
             )
 
@@ -6257,7 +6257,7 @@ def projection(self):
 
             if self.enforce_active:
                 self.local_active[
-                    self.local_mesh._get_containing_cell_indexes(
+                    self.local_mesh.get_containing_cells(
                         self.global_mesh.cell_centers[~self.global_active, :]
                     )
                 ] = False