diff --git a/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Fixed_border_iterative_parameterizer_3.h b/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Fixed_border_iterative_parameterizer_3.h
index 204ce03f211..2ade4d55fb3 100755
--- a/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Fixed_border_iterative_parameterizer_3.h
+++ b/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Fixed_border_iterative_parameterizer_3.h
@@ -9,7 +9,6 @@
 #define CGAL_SURFACE_MESH_PARAMETERIZATION_FIXED_BORDER_ITERATIVE_PARAMETERIZER_3_H
 
 #define DEBUG_L0 1
-#define DEBUG_L1 0
 
 #include <CGAL/license/Surface_mesh_parameterization.h>
 
@@ -239,11 +238,6 @@ public:
     //
     initialize_system_from_mesh_border(A, Bu, Bv, mesh, bhd, uvmap1, vimap1);
 
-    if(DEBUG_L1) {
-      printMatrix(vertices, vimap1, A, "A");
-      printVector(vertices, vimap1, Bu, "Bu");
-      printVector(vertices, vimap1, Bv, "Bv");
-    }
     // Fill the matrix for the inner vertices v_i: compute A's coefficient
     // w_ij for each neighbor j; then w_ii = - sum of w_ijs
     boost::unordered_set<vertex_descriptor> main_border;
@@ -298,12 +292,6 @@ public:
         }
       }
 
-      if(DEBUG_L1) {
-        print(A, Xu, Bu);
-        std::cout << std::endl;
-        print(A, Xv, Bv);
-      }
-
       // solve linear equations
       // Solve "A*Xu = Bu". On success, solution is (1/Du) * Xu.
       // Solve "A*Xv = Bv". On success, solution is (1/Dv) * Xv.
diff --git a/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Iterative_authalic_parameterizer_3.h b/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Iterative_authalic_parameterizer_3.h
index 87aab77bec7..3a313c43356 100755
--- a/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Iterative_authalic_parameterizer_3.h
+++ b/Surface_mesh_parameterization/include/CGAL/Surface_mesh_parameterization/Iterative_authalic_parameterizer_3.h
@@ -406,21 +406,12 @@ private:
   }
 
   double getvL2(const TriangleMesh& mesh, vertex_descriptor &vertex)  const {
-    bool debug = false;
-    if(vertex == 711111)
-      debug = true;
-    if(debug)
-      std::cout << "\t" << vertex << "\t" << std::flush;
     halfedge_descriptor hf = halfedge(vertex, mesh);
-    if(debug)
-      std::cout << hf << "\t" << std::flush;
     face_around_target_circulator f_j(hf, mesh), end_f_j = f_j;
     double varphi = 0.0;
     double localArea = 0.0;
     int i=0;
     CGAL_For_all(f_j, end_f_j)  {
-      if(debug)
-        std::cout << *f_j << "\t" << std::flush;
       if(*f_j > mesh.number_of_faces())
         continue;
       varphi += get(fL2Map,*f_j)*get(areaMap,*f_j);
@@ -428,11 +419,6 @@ private:
       i++;
     }
 
-    if(debug) {
-      std::cout << varphi << "\t" << std::flush;
-      std::cout << localArea << "\t" << std::flush;
-    }
-
     if(mesh.is_border(vertex) && mesh.degree(vertex) != i+1)
       std::cerr << std::endl;
     else if(!mesh.is_border(vertex) && mesh.degree(vertex) != i)