Member model::Polyhedron< dim, order >::finalize (): This function uses ConvexHull and therefore will yield a wrong face boundary for non-convex faces. In order to support non-convex faces we need to have a structure with face edges.