correct bug

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;
+          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);