Add convenience function for per-face data

[ffreyer]

This allows

mesh(..., attribute_name = per_face(data, faces))
# or alternatively:
mesh(old_mesh, attribute_name = per_face(data, old_mesh))