diff --git a/apps/sample/polygonmesh_quadsimpl/quadsimpl.cpp b/apps/sample/polygonmesh_quadsimpl/quadsimpl.cpp
index db32f966..786a102f 100644
--- a/apps/sample/polygonmesh_quadsimpl/quadsimpl.cpp
+++ b/apps/sample/polygonmesh_quadsimpl/quadsimpl.cpp
@@ -464,6 +464,9 @@ int main(int argc, char *argv[]) {
     // After loading check mesh consistency
     assert(vcg::tri::UpdateHalfEdges<MyPolyMesh>::CheckConsistency(pm));
 
+    HalfedgeQuadClean<MyPolyMesh>::remove_singlets(pm);
+    HalfedgeQuadClean<MyPolyMesh>::remove_doublets(pm);
+
     vcg::LocalOptimization<MyPolyMesh> loc(pm);
     init_heap(pm, loc, adaptive);