code: def fixed_point_uniform(points, cells, *args, **kwargs):
      
  
      def get_new_points(mesh):
          return get_new_points_volume_averaged(mesh, mesh.cell_barycenters)
  
      mesh = MeshTri(points, cells)
      runner(get_new_points, mesh, *args, **kwargs)
      return mesh.node_coords, mesh.cells["nodes"]
docs: Idea:
      Move interior mesh points into the weighted averages of the centroids
      (barycenters) of their adjacent cells.