@pimms.param<EOL><INDENT>def coordinates(crds):<DEDENT>
crds = np.asarray(crds)<EOL>if crds.shape[<NUM_LIT:0>] != <NUM_LIT:2> and crds.shape[<NUM_LIT:0>] != <NUM_LIT:3>:<EOL><INDENT>crds = crds.T<EOL>if crds.shape[<NUM_LIT:0>] != <NUM_LIT:2> and crds.shape[<NUM_LIT:0>] != <NUM_LIT:3>:<EOL><INDENT>raise ValueError('<STR_LIT>')<EOL><DEDENT><DEDENT>return pimms.imm_array(crds)<EOL>
mesh.coordinates is a read-only numpy array of size (d x n) where d is the number of
  dimensions and n is the number of vertices in the mesh.