void MeshDenoisingBase::getFaceNeighbor(TriMesh &mesh, TriMesh::FaceHandle fh, FaceNeighborType face_neighbor_type, std::vector<TriMesh::FaceHandle> &face_neighbor) { face_neighbor.clear(); if(face_neighbor_type == kEdgeBased) { for(TriMesh::FaceFaceIter ff_it = mesh.ff_iter(fh); ff_it.is_valid(); ff_it++) face_neighbor.push_back(*ff_it); } else if(face_neighbor_type == kVertexBased) { std::set<int> neighbor_face_index; neighbor_face_index.clear(); for(TriMesh::FaceVertexIter fv_it = mesh.fv_begin(fh); fv_it.is_valid(); fv_it++) { for(TriMesh::VertexFaceIter vf_it = mesh.vf_iter(*fv_it); vf_it.is_valid(); vf_it++) { if((*vf_it) != fh) neighbor_face_index.insert(vf_it->idx()); } } for(std::set<int>::iterator iter = neighbor_face_index.begin(); iter != neighbor_face_index.end(); ++ iter) { face_neighbor.push_back(TriMesh::FaceHandle(*iter)); } } }