diff --git a/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/repair_polygon_soup.h b/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/repair_polygon_soup.h
index b6df34136cb..005eb344e3d 100644
--- a/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/repair_polygon_soup.h
+++ b/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/repair_polygon_soup.h
@@ -126,43 +126,21 @@ bool simplify_polygon(PointRange& points,
 {
   const std::size_t ini_polygon_size = polygon.size();
 
-  // Start at the last since if two points are identical, the second one gets removed.
-  // By starting at 'last', we ensure that 'to_remove' is ordered from closest to .begin()
-  // to closest to .end()
-  std::size_t last = ini_polygon_size - 1, i = last;
-  bool stop = false;
-  std::vector<std::size_t> to_remove;
-
-  do
+  for(std::size_t i=0; i<polygon.size(); ++i)
   {
-    std::size_t next_i = (i == last) ? 0 : i+1;
-    stop = (next_i == last);
+    const std::size_t s = polygon.size();
+    if(s == 1)
+      break;
 
-    while(polygon[i] == polygon[next_i] || // combinatorial equality
-          traits.equal_3_object()(points[polygon[i]], points[polygon[next_i]])) // geometric equality
+    const std::size_t ni = (i + 1) % s;
+    if(polygon[i] == polygon[ni] ||
+       traits.equal_3_object()(points[polygon[i]], points[polygon[ni]]))
     {
-      to_remove.push_back(next_i);
 #ifdef CGAL_PMP_REPAIR_POLYGON_SOUP_VERBOSE_PP
-      std::cout << "Duplicate point: polygon[" << next_i << "] = " << polygon[next_i] << std::endl;
+      std::cout << "Duplicate point: polygon[" << ni << "] = " << polygon[ni] << std::endl;
 #endif
-      next_i = (next_i == last) ? 0 : next_i+1;
-
-      // Every duplicate in front of 'last' (circularly-speaking) has already been cleared
-      if(next_i == last)
-      {
-        stop = true;
-        break;
-      }
+      polygon.erase(polygon.begin() + i--);
     }
-
-    i = next_i;
-  }
-  while(!stop);
-
-  while(!to_remove.empty())
-  {
-    polygon.erase(polygon.begin() + to_remove.back());
-    to_remove.pop_back();
   }
 
   const std::size_t removed_points_n = ini_polygon_size - polygon.size();