diff --git a/src/geo/GModelParametrize.cpp b/src/geo/GModelParametrize.cpp
index 3824ecbb8a05e9068be8b00399448f9d24d1cfa4..d83443ea7d78922a9a846a75d4919fd9d2c8d89e 100644
--- a/src/geo/GModelParametrize.cpp
+++ b/src/geo/GModelParametrize.cpp
@@ -732,9 +732,7 @@ static int isTriangulationParametrizable(const std::vector<MTriangle *> &t,
 
   std::set<MVertex *> v;
   std::map<MEdge, int, MEdgeLessThan> e;
-  double surf = 0;
   for(std::size_t i = 0; i < t.size(); ++i) {
-    surf += t[i]->getVolume();
     for(int j = 0; j < 3; j++) {
       v.insert(t[i]->getVertex(j));
       auto it = e.find(t[i]->getEdge(j));
diff --git a/src/mesh/delaunay3d.cpp b/src/mesh/delaunay3d.cpp
index 8467afcd9248cac7a66671876e3670ad644c7e96..b45906328388776f8ccbe789d9e510c355ad7e7b 100644
--- a/src/mesh/delaunay3d.cpp
+++ b/src/mesh/delaunay3d.cpp
@@ -993,7 +993,6 @@ void delaunayTrgl(const std::size_t numThreads, const std::size_t NPTS_AT_ONCE,
 #endif
 
     double totSearch = 0;
-    double totCavity = 0;
     std::vector<std::size_t> _negatives;
     std::vector<cavityContainer> cavity(NPTS_AT_ONCE);
     std::vector<connContainer> bnd(NPTS_AT_ONCE);
@@ -1076,7 +1075,6 @@ void delaunayTrgl(const std::size_t numThreads, const std::size_t NPTS_AT_ONCE,
           faceToTet.clear();
           const std::size_t cSize = cavityK.size();
           const std::size_t bSize = bndK.size();
-          totCavity += cSize;
           Choice[K] = cavityK[0];
           for(std::size_t i = 0; i < bSize; i++) {
             // reuse memory slots of invalid elements
diff --git a/src/mesh/meshGFaceDelaunayInsertion.cpp b/src/mesh/meshGFaceDelaunayInsertion.cpp
index afbd7b0741f12162c0baf2d3acad1f048d42b385..81ee34328bbca7a2ff30bbb64941b88067b1492e 100644
--- a/src/mesh/meshGFaceDelaunayInsertion.cpp
+++ b/src/mesh/meshGFaceDelaunayInsertion.cpp
@@ -1826,9 +1826,6 @@ void delaunayMeshIn2D(std::vector<MVertex *> &v,
   MVertex *box[4];
   initialSquare(v, box, t);
 
-  double AVG_ITER = 0;
-  double AVG_CAVSIZE = 0;
-
   if(hilbertSort) SortHilbert(v);
 
   for(size_t i = 0; i < v.size(); i++) {
@@ -1836,7 +1833,6 @@ void delaunayMeshIn2D(std::vector<MVertex *> &v,
 
     int NITER = 0;
     MTri3 *found = getTriToBreak(pv, t, NITER);
-    AVG_ITER += (double)NITER;
     if(!found) {
       Msg::Error("Cannot insert a point in 2D Delaunay");
       continue;
@@ -1845,7 +1841,6 @@ void delaunayMeshIn2D(std::vector<MVertex *> &v,
     cavity.clear();
 
     recurFindCavity(shell, cavity, pv, found);
-    AVG_CAVSIZE += (double)cavity.size();
 
     std::vector<MTri3 *> extended_cavity;
     for(std::size_t count = 0; count < shell.size(); count++) {