Function getCombinedPoint

src/utilities/geometry/Geometry.cpp:423–432

}

Point3d getCombinedPoint(const Point3d& point3d, std::vector<Point3d>& allPoints, double tol) {
  for (const Point3d& otherPoint : allPoints) {
    if (std::sqrt(std::pow(point3d.x() - otherPoint.x(), 2) + std::pow(point3d.y() - otherPoint.y(), 2) + std::pow(point3d.z() - otherPoint.z(), 2))
        < tol) {
      return otherPoint;
    }
  }
  allPoints.push_back(point3d);
  return point3d;
}

std::vector<std::vector<Point3d>> computeTriangulation(const Point3dVector& vertices, const std::vector<std::vector<Point3d>>& holes, double tol) {
  std::vector<std::vector<Point3d>> result;