CGAL seam_mesh: как включить шов в halfedge_descriptor границы?

Я параметризую простую сетку (см. рисунок ниже) с помощью 8cpp-example.html" rel="nofollow noreferrer">seam_Polyhedron_3. cpp пример CGAL 4.14.

Я определил вершины для «разрезания», которые я записал в файл шва (они лежат вдоль одной полилинии, соединяющей нижнюю и верхнюю границу):

73 75 80 120

Тем не менее, я смущен тем, что линия

halfedge_descriptor bhd = CGAL::Polygon_mesh_processing::longest_border(mesh).first;

возвращает: кажется, что он перебирает только нижнюю границу. Это я мог проверить, посчитав вершины на границе,

typedef typename boost::graph_traits<Mesh>::vertex_descriptor Mesh_vertex_descriptor;
Mesh_vertex_descriptor beg = source(bhd, mesh);
int i=0;
do
{
    i++;
    bhd = next(bhd, mesh);
}
while(source(bhd, mesh) != beg);
std::cout << "There were " << i << " vertices on the boundary." << std::endl;

который вернул 34 (что, в свою очередь, является количеством вершин на нижней границе).

Вопрос. Как изменить конструкцию bhd, чтобы она охватывала всю (виртуальную) границу? То есть итерация, которую я только что показал, должна проходить через нижнюю границу, шов вверх, верхнюю границу и шов вниз и возвращать 32 + 3 + 32 + 3 = 70?

До сих пор я пытался найти альтернативу longest_boundary и посмотрите, есть ли возможность изменить поведение через именованные параметры . Я также пробовал добавлять вершины на нижней и верхней границах сетки, но это ничего не меняет, как указано в документации Seam_mesh. Изменение порядка вершин шва также не дало видимых изменений.

Заранее спасибо!