@pimms.value<EOL><INDENT>def edge_centers(edge_coordinates):<DEDENT>
return pimms.imm_array(<NUM_LIT:0.5> * np.sum(edge_coordinates, axis=<NUM_LIT:0>))<EOL>
mesh.edge_centers is the (d x n) array of the centers of each edge in the given mesh.