diff --git a/src/mesh/MeshSmoother.cpp b/src/mesh/MeshSmoother.cpp
index 5b77dcdf7cf44245bb8112feb2342773dbd52718..8d81332b797a5825c420b530c85eacc402bb3129 100644
--- a/src/mesh/MeshSmoother.cpp
+++ b/src/mesh/MeshSmoother.cpp
@@ -173,10 +173,10 @@ class MeshSmootherHandler::MeshSmoother
       m_given_mesh.numberOfNodes(), PUGS_LAMBDA(const NodeId node_id) {
         const auto& node_cell_list = node_to_cell_matrix[node_id];
 
-        Rd mean_position(zero);
-        size_t number_of_neighbours = 0;
         for (size_t i_cell = 0; i_cell < node_cell_list.size(); ++i_cell) {
           const size_t i_cell_node = node_number_in_their_cells(node_id, i_cell);
+          Rd mean_position(zero);
+          size_t number_of_neighbours = 0;
 
           const CellId cell_id       = node_cell_list[i_cell];
           const auto& cell_node_list = cell_to_node_matrix[cell_id];
@@ -187,13 +187,13 @@ class MeshSmootherHandler::MeshSmoother
               number_of_neighbours++;
             }
           }
+          mean_position    = 1. / number_of_neighbours * mean_position;
+          shift_r[node_id] = mean_position - given_xr[node_id];
         }
-        mean_position    = 1 / number_of_neighbours * mean_position;
-        shift_r[node_id] = mean_position - given_xr[node_id];
-        double nshift    = sqrt(dot(shift_r[node_id], shift_r[node_id]));
-        if (nshift > max_delta_xr[node_id]) {
-          shift_r[node_id] = max_delta_xr[node_id] / nshift * shift_r[node_id];
-        }
+        // double nshift    = sqrt(dot(shift_r[node_id], shift_r[node_id]));
+        // if (nshift > max_delta_xr[node_id]) {
+        //   shift_r[node_id] = max_delta_xr[node_id] / nshift * shift_r[node_id];
+        // }
       });
 
     this->_applyBC(shift_r);