[논문 리뷰] Arrangement Computation for Planar Algebraic Curves

Eric Berberich, Pavel Emeliyanenko외 2명|arXiv (Cornell University)|2011. 03. 24.

Advanced Numerical Analysis Techniques참고 문헌 26인용 수 4

한 줄 요약

이 논문은 결과식(resultant)과 GCD만을 기호 연산으로 사용하고, 인증된 복소수 근 해소기로 하이브리드 기호-수치 리프팅 단계를 수행하는, 인증되고 완전한 평면 대수 곡선 배열을 계산하는 알고리즘을 제시한다. 이 방법은 좌표 변환 없이 GPU 가속 산술을 활용하고, 좌표 변형을 피함으로써, 특히 고차수 또는 특이 곡선인 경우 기존 최고 수준의 성능(CGAL의 Ak_2) 대비 최대 120배 빠른 성능 향상을 달성한다.

ABSTRACT

We present a new certified and complete algorithm to compute arrangements of real planar algebraic curves. Our algorithm provides a geometric-topological analysis of the decomposition of the plane induced by a finite number of algebraic curves in terms of a cylindrical algebraic decomposition of the plane. Compared to previous approaches, we improve in two main aspects: Firstly, we significantly reduce the amount of exact operations, that is, our algorithms only uses resultant and gcd as purely symbolic operations. Secondly, we introduce a new hybrid method in the lifting step of our algorithm which combines the usage of a certified numerical complex root solver and information derived from the resultant computation. Additionally, we never consider any coordinate transformation and the output is also given with respect to the initial coordinate system. We implemented our algorithm as a prototypical package of the C++-library CGAL. Our implementation exploits graphics hardware to expedite the resultant and gcd computation. We also compared our implementation with the current reference implementation, that is, CGAL's curve analysis and arrangement for algebraic curves. For various series of challenging instances, our experiments show that the new implementation outperforms the existing one.

연구 동기 및 목표

• 좌표 변환 없이 실수 평면 대수 곡선 배열을 계산하기 위한 인증되고 완전한 알고리즘을 개발하기.
• 곡선 위상 계산에서 비용이 많이 드는 기호 연산인 부분결과식(subresultants)과 그로버 기저(Gröbner bases)에 대한 의존도를 줄이기.
• 이전 방법이 계산 블로킹 문제를 겪는 고차수 및 특이 곡선의 효율성을 향상시키기.
• 결과식과 GCD 계산의 GPU 가속을 통해 실용적인 성능 향상을 이룰 수 있도록 하기.
• 기호 계산을 최소화하고 변형을 피하면서 기하학적 정확성과 완전성을 유지하기.

제안 방법

• 알고리즘은 평면의 원통형 대수적 분할을 사용하여 대수 곡선의 배열을 계산하며, 단일 및 쌍별 곡선을 기하적-위상적 분할을 통해 분석한다.
• 기호 연산으로 결과식과 GCD만을 사용하고, 이를 그래픽스 하드웨어로 이관하여 가속화한다.
• 인증된 수치적 복소수 근 해소기와 테이시에르(Teissier)의 복소수 근 개수에 대한 이론적 경계를 조합하여 실근을 효율적으로 고립하는 하이브리드 리프팅 단계를 사용한다.
• 기본 근 고립 방법으로 기호-수치 필터인 FastLift를 사용하며, 희귀한 경우를 제외하고는 실패하지 않는다.
• FastLift가 실패하는 경우, [4]에 기반한 후행 방법을 사용하여 완전성을 보장한다.
• 전체 계산이 원래 좌표계에서 수행되어 입력의 희소성과 변형으로 인한 계수의 과도한 증가를 방지한다.

실험 결과

연구 질문

• RQ1 대수 곡선 배열 계산에서 기호 연산의 계산 비용을 정확성과 완전성을 유지하면서 크게 줄일 수 있는가?
• RQ2 기호-수치 접근법이 곡선 위상 계산의 리프팅 단계에서 근 고립에 얼마나 효과적인가?
• RQ3 GPU 가속이 대수 곡선의 결과식과 GCD와 같은 기호 계산 작업에서 성능 향상에 어느 정도 기여하는가?
• RQ4 특이 곡선 또는 비일반 곡선의 배열에서 좌표 변형의 부재가 성능과 수치 안정성에 어떤 영향을 미치는가?
• RQ5 인증되고 완전한 알고리즘이 고차수 또는 투영된 교차 곡선과 같은 어려운 벤치마크 인스턴스에서 현재 기준 구현(CGAL의 Ak_2)을 초월할 수 있는가?

주요 결과

• 쉬운 인스턴스인 일반 위치의 비특이 곡선에 대해, 새로운 구현인 FastKernel은 GPU 가속된 Ak_2보다 평균적으로 두 배 빠르다.
• 표면의 투영된 교차 곡선이나 특수한 특이성을 가진 곡선에 대해서는 성능 향상이 5배에서 120배에 이르며, 곡선의 차수 증가에 따라 성능 향상이 커진다.
• 가장 도전적인 테스트 케이스(유형 X의 11개 곡선)에서 FastKernel은 2696.15초 만에 완료되었고, Ak_2는 4000초 후 타임아웃되었다.
• 대부분의 경우에서 기호-수치 필터인 FastLift는 모든 피브어(fibers)에서 성공했으며, 희귀한 경우에만 후행 방법인 Lift이 호출되었다.
• 성능 향상은 고차수 곡선에서 가장 두드러지며, 이 경우 Ak_2의 부분결과식 계산이 주요 블로킹 요소가 된다.
• 변형을 피함으로써 계수의 희소성을 유지하고, 이전 방법에서 성능 저하를 초래하는 계수 비트 크기의 과도한 증가를 방지한다.