효율적인 가시다각형 계산을 위한 구현과 실험

초록

본 논문은 CGAL에 포함될 세 가지 가시다각형 알고리즘(조·심프슨 O(n), 아사노 O(n log n), 새로운 삼각형 확장 O(n²))의 정확한 구현과 광범위한 실험을 제시한다. 특히 삼각형 확장 알고리즘은 최악의 경우 O(n²) 복잡도를 갖지만 실제 데이터에서는 두 자릿수 정도 빠른 성능을 보이며, 모든 구현은 정확한 수학 연산과 퇴화 케이스 처리를 지원한다.

상세 분석

논문은 먼저 기존의 가시다각형 알고리즘을 정리하고, CGAL에 부재하던 정확한 구현을 메우는 목표를 제시한다. 조·심프슨 알고리즘은 단순 다각형에 대해 경계 스택을 순차적으로 스캔하면서 O(n) 시간·공간 복잡도를 달성한다. 구현 시 다각형이 360° 이상 감싸는 경우를 처리하기 위해 winding counter를 도입했으며, 스택 조작을 상태 머신 형태로 구조화해 코드 가독성과 디버깅을 용이하게 했다. 아사노 알고리즘은 회전 스위프 방식을 사용한다. 모든 정점을 q에 대한 극각 순으로 정렬하고, 현재 레이와 교차하는 선분들을 균형 이진 탐색 트리(T)로 관리한다. T의 최솟값이 변할 때마다 가시다각형의 새로운 꼭짓점을 출력한다. 핵심은 레이와 선분 사이의 순서를 교차점을 실제로 계산하지 않고, 최대 다섯 개의 방향 판정(orientation)만으로 결정하는 기법이다. 이는 정확한 부동소수점 연산 없이도 정수 연산만으로 비교를 가능하게 하여 CGAL의 정확한 커널과 잘 맞는다. 새로운 삼각형 확장 알고리즘은 사전 삼각분할(Constrained Delaunay)을 이용한다. 질의점 q가 포함된 삼각형을 찾은 뒤, 그 삼각형의 경계와 다각형 외곽을 따라 재귀적으로 시야를 확장한다. 각 재귀 단계에서는 현재 시야를 제한하는 두 ‘반사점’(reflex vertices)을 유지하고, 새로운 삼각형을 탐색하면서 필요한 경우 시야를 다시 분할한다. 최악의 경우 모든 삼각형을 두 번씩 방문해 O(n²) 시간이지만, 실제 입력에서는 시야에 포함되는 삼각형만 탐색하므로 사실상 출력 민감적이다. 특히 구멍이 없는 단순 다각형에서는 O(n)에 수렴한다. 구현은 스택 기반 DFS와 삼각형 인접 정보를 활용해 O(n) 메모리를 유지한다. 실험에서는 노르웨이(≈20k 정점)와 대성당(≈1.2k 정점) 같은 실제 데이터에서 삼각형 확장 알고리즘이 평균 0.18 ms·쿼리로 가장 빠른 성능을 보였으며, 인위적인 최악 사례에서는 O(n log n) 아사노 알고리즘이 결국 앞선다. 전체적으로 구현은 정확한 연산, 퇴화 처리, 안테나(1차원 선분) 포함 여부 옵션을 제공해 실무 적용성을 높였다.