This function extracts a triangulation data structure
from an triangulation object created by tri.mesh.
The vertices in the returned matrix (let's denote it with
retval) are ordered
counterclockwise with the first vertex taken
to be the one with smallest index. Thus,
retval[i,"node2"] and retval[i,"node3"] are larger
than
retval[i,"node3"] and index adjacent neighbors of
node retval[i,"node1"]. The columns trx and
arcx, x=1,2,3 index the triangle and arc,
respectively, which are opposite (not shared
by) node nodex, with trix= 0 if
arcx indexes a boundary arc. Vertex
indexes range from 1 to N, triangle indexes
from 0 to NT, and, if included, arc indexes
from 1 to NA = NT+N-1. The triangles are
ordered on first (smallest) vertex indexes,
except that the sets of constraint triangles
(triangles contained in the closure of a constraint
region) follow the non-constraint
triangles.