Периодическая 3d триангуляция - индексирование ячеек

Question

Я использую CGAL для выполнения триангуляции 3d delaunay для большого набора точек.

Моя цель состоит в том, чтобы сгенерировать выходной файл с форматом:

ICell Vertex1 Vertex2 Vertex3 Vertex4 NeigborCell1 NeighborCell2 .... 

Где все значения индексов, соответствующие конкретной ячейки и вершины.

Для этого я включил библиотеку «Карта» и использовал ее для обработки вершинных дескрипторов и ручек Cell для привязки их к индексу.
Этот метод отлично работает для вершин, но в случае ячеек он выводит все ячейки, которые были сгенерированы в триангуляции (включая все периодические изображения той же ячейки).
Кроме того, поскольку изображения ячеек имеют разные ячейки, им присваивается другой индекс (так много соседей тоже являются изображениями)

Есть ли способ перебора между ячейками один раз? Как я могу индексировать их таким образом, чтобы каждый из них получал уникальный индекс?

Заранее за вашу помощь.

Answer 1

Я боюсь, что вы не можете проходить так, как хотите, по крайней мере, в текущей версии пакета. Вы можете сделать это только на тетраэдрах (геометрическое вложение ячеек) с помощью Periodic_tetrahedron_iterator с соответствующим типом Iterator_type. Возможно, вы можете взломать код Periodic_tetrahedron_iterator, чтобы сделать то, что вам нужно ...