void add_neighbourhood_to_hullPoints(pointVector& hullPoints, const Vertex_const_handle& vert, const double& tapeSize) { Point_3 pnt = vert->point(); hullPoints.push_back(pnt); // Add vertex point Nef_polyhedron::SVertex_const_iterator svcIt = vert->svertices_begin(), svcItEND = vert->svertices_end(); CGAL_For_all(svcIt,svcItEND) { Vector_3 vecR(pnt,svcIt->target()->point()); Vector_3 vecRnew = vecR * tapeSize / std::sqrt(CGAL::to_double(vecR.squared_length())); if ((vecR.squared_length()-OVERLAP_DIST_THRESHOLD) > vecRnew.squared_length()) hullPoints.push_back(pnt+vecRnew); // Add svertex neighbourhood point (tapesize away from vertex) else hullPoints.push_back(svcIt->target()->point()); }