在对三角网格模型进行切分并且对切分模型三角化之后，可能会遇到分块数 > 2的情况，对于这种情况我们就需要得到各个分块的数量，并且将他们重组成独立的三角网格模型。
思路如下：
建立一个存放三角网格面的堆栈，遍历模型的所有的三角网格面，对每一个三角形都搜索与之相邻的三个三角形（每个边对应一个三角形），将这些三角形放入堆栈，直到堆栈为空表明已经找到属于同一个块的三角网格。再遍历以上过程，直到所有的三角网格都有属于自己的块。
在CGAL的Polyhedron_3.h头文件的Facet类里加入一个int型标志位mask，并且在构造器里将其初始化为0。接下来是具体的搜索过程，trilist是一个Facet_handle类型的堆栈。
int noc = 0; Facet_handle t , s; //所有的facet的mask都为0 for(Facet_iterator it = P1.facets_begin() , end = P1.facets_end() ; it != end;it++) { if(!it->mask) { noc ++ ; trilist.push(it); it->mask = noc ; while(trilist.size()) { t = trilist.top(); trilist.pop(); if((s=t->halfedge()->opposite()->facet()) != NULL && (!s->mask)) {trilist.push(s);s->mask=noc;} if((s=t->halfedge()->next()->opposite()->facet()) != NULL && (!s->mask)){trilist.push(s);s->mask=noc;} if((s=t->halfedge()->prev()->opposite()->facet()) != NULL && (!s->mask)){trilist.push(s);s->mask=noc;} } } //std::cout << "Mask of facet is "<< it->mask <<std::endl; } std::cout << "Number of components is "<< noc <<std::endl;
原文：http://www.cnblogs.com/jast/p/4169977.html