diff --git a/Polygon_mesh_processing/examples/Polygon_mesh_processing/CMakeLists.txt b/Polygon_mesh_processing/examples/Polygon_mesh_processing/CMakeLists.txt
index 9d70695db33..39c1bde74dd 100644
--- a/Polygon_mesh_processing/examples/Polygon_mesh_processing/CMakeLists.txt
+++ b/Polygon_mesh_processing/examples/Polygon_mesh_processing/CMakeLists.txt
@@ -74,6 +74,7 @@ include( CGAL_CreateSingleSourceCGALProgram )
 
 create_single_source_cgal_program( "self_intersections_example.cpp" )
 create_single_source_cgal_program( "hole_filling_example.cpp" )
+create_single_source_cgal_program( "hole_filling_example_SM.cpp" )
 create_single_source_cgal_program( "stitch_borders_example.cpp" )
 create_single_source_cgal_program( "compute_normals_example.cpp" )
 create_single_source_cgal_program( "point_inside_example.cpp")
diff --git a/Polygon_mesh_processing/examples/Polygon_mesh_processing/hole_filling_example_SM.cpp b/Polygon_mesh_processing/examples/Polygon_mesh_processing/hole_filling_example_SM.cpp
new file mode 100644
index 00000000000..818942d9d9f
--- /dev/null
+++ b/Polygon_mesh_processing/examples/Polygon_mesh_processing/hole_filling_example_SM.cpp
@@ -0,0 +1,62 @@
+#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
+#include <CGAL/Surface_mesh.h>
+#include <CGAL/boost/graph/graph_traits_Surface_mesh.h>
+
+#include <CGAL/Polygon_mesh_processing/triangulate_hole.h>
+
+#include <iostream>
+#include <fstream>
+#include <vector>
+#include <boost/foreach.hpp>
+
+typedef CGAL::Exact_predicates_inexact_constructions_kernel Kernel;
+typedef Kernel::Point_3 Point;
+typedef CGAL::Surface_mesh<Point> Mesh;
+
+typedef boost::graph_traits<Mesh>::halfedge_descriptor    halfedge_descriptor;
+typedef boost::graph_traits<Mesh>::face_descriptor       face_descriptor;
+typedef boost::graph_traits<Mesh>::vertex_descriptor      vertex_descriptor;
+
+int main(int argc, char* argv[])
+{
+  const char* filename = (argc > 1) ? argv[1] : "data/mech-holes-shark.off";
+  std::ifstream input(filename);
+
+  Mesh mesh;
+  if ( !input || !(input >> mesh) ) {
+    std::cerr << "Not a valid off file." << std::endl;
+    return 1;
+  }
+
+  // Incrementally fill the holes
+  unsigned int nb_holes = 0;
+  BOOST_FOREACH(halfedge_descriptor h, halfedges(mesh))
+  {
+    if(is_border(h,mesh))
+    {
+      std::vector<face_descriptor>  patch_facets;
+      std::vector<vertex_descriptor> patch_vertices;
+      bool success = CGAL::cpp11::get<0>(
+        CGAL::Polygon_mesh_processing::triangulate_refine_and_fair_hole(
+                  mesh,
+                  h,
+                  std::back_inserter(patch_facets),
+                  std::back_inserter(patch_vertices),
+     CGAL::Polygon_mesh_processing::parameters::vertex_point_map(get(CGAL::vertex_point, mesh)).
+                  geom_traits(Kernel())) );
+
+      std::cout << "* Number of facets in constructed patch: " << patch_facets.size() << std::endl;
+      std::cout << "  Number of vertices in constructed patch: " << patch_vertices.size() << std::endl;
+      std::cout << "  Is fairing successful: " << success << std::endl;
+      nb_holes++;
+    }
+  }
+
+  std::cout << std::endl;
+  std::cout << nb_holes << " holes have been filled" << std::endl;
+  
+  std::ofstream out("filled_SM.off");
+  out.precision(17);
+  out << mesh << std::endl;
+  return 0;
+}