use PDL::VectorValued::Utils;
# $pdl is ncoords=3,nvertices=3,ntriangles
my $uniqv = $pdl->uniqvec; # ncoords=3,nunique
my $face_indices = vsearchvec($pdl, $uniqv); # nvertices=3,ntriangles