diff --git a/src/cellcomplex/property_topomesh/example_topomesh.py b/src/cellcomplex/property_topomesh/example_topomesh.py
index ab475b9b7fa675121e07e7000230491da0cf25a0..e2d84c7ced6a2a45b45efb12fece4d71b915b38f 100644
--- a/src/cellcomplex/property_topomesh/example_topomesh.py
+++ b/src/cellcomplex/property_topomesh/example_topomesh.py
@@ -742,9 +742,9 @@ def circle_voronoi_topomesh(size = 1,cell_size = 1.,circle_size = 100.,z_coef =
         iterations = n_cells//2
 
     for iteration in range(iterations):
-        cell_to_cell_vectors = cell_points[np.newaxis] - cell_points[:, np.newaxis]
+        cell_to_cell_vectors = cell_points[:, np.newaxis] - cell_points[np.newaxis]
         cell_to_cell_distances = np.linalg.norm(cell_to_cell_vectors,axis=2)/radius
-        cell_to_circle_vectors = np.array([[p-q for q in circle_points] for p in cell_points])
+        cell_to_circle_vectors = cell_points[:, np.newaxis] - circle_points[np.newaxis]
         cell_to_circle_distances = np.linalg.norm(cell_to_circle_vectors,axis=2)/radius
 
         deformation_force = np.zeros_like(cell_points)