element_centers Mass centers of given elements.
x0 = element_centers(msh, indices_to_mesh_elements)
x0 = element_centers(msh, element_definitions_by_node_indices)