@pimms.value<EOL><INDENT>def face_centers(face_coordinates):<DEDENT>
return pimms.imm_array(np.sum(face_coordinates, axis=<NUM_LIT:0>) / <NUM_LIT>)<EOL>
mesh.face_centers is the (d x n) array of the centers of each triangle in the given mesh.