If you want to get the faces of the triangulation of your non-convex
polygon, you can just start from a face that you know is inside
(which can be found by using 'locate()' on a point you know is
inside the polygon), and then using a spreading algorithm to "color"
faces as long as they are not separated by a constrained edge.