--- a/src/Grids/mapped_grid.jl	Wed Aug 21 19:23:42 2024 +0200
+++ b/src/Grids/mapped_grid.jl	Thu Aug 22 08:14:04 2024 +0200
@@ -92,6 +92,28 @@
     return ms
 end
 
+function min_spacing(g::MappedGrid{T,2} where T)
+    n, m = size(g)
+
+    ms = Inf
+    for i ∈ 1:n-1, j ∈ 1:m-1 # loop over each cell of the grid
+
+        ms = min(
+            ms,
+            norm(g[i+1,j]-g[i,j]),
+            norm(g[i,j+1]-g[i,j]),
+
+            norm(g[i+1,j]-g[i+1,j+1]),
+            norm(g[i,j+1]-g[i+1,j+1]),
+
+            norm(g[i+1,j+1]-g[i,j]),
+            norm(g[i+1,j]-g[i,j+1]),
+        )
+        # NOTE: This could be optimized to avoid checking all interior edges twice.
+    end
+
+    return ms
+end
 
 """
     normal(g::MappedGrid, boundary)