[논문 리뷰] Exact Symbolic-Numeric Computation of Planar Algebraic Curves

Eric Berberich, Pavel Emeliyanenko외 2명|arXiv (Cornell University)|2012. 01. 07.

Polynomial and algebraic computation참고 문헌 5인용 수 5

한 줄 요약

이 논문은 기호적 및 수치적 기법을 융합하여 평면 대수 곡선 배열을 인증 가능하고 완전하며 효율적으로 계산하는 알고리즘을 제안한다. 이는 두 가지 핵심 구성요소—이원수 시스템의 실근을 고립하는 데 쓰이는 Bisolve와 단일 곡선을 분석하는 데 쓰이는 GeoTop—을 도입하며, 기호 연산으로는 결과식(resultants)과 최대공약수(GCDs)만을 사용하고, GPU 계산을 통해 이를 가속화하며 좌표 변환을 피함으로써, 기존 방법에 비해 도전적인 사례에서 뛰어난 성능을 발휘한다.

ABSTRACT

We present a novel certified and complete algorithm to compute arrangements of real planar algebraic curves. It 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. From a high-level perspective, the overall method splits into two main subroutines, namely an algorithm denoted Bisolve to isolate the real solutions of a zero-dimensional bivariate system, and an algorithm denoted GeoTop to analyze a single algebraic curve. Compared to existing approaches based on elimination techniques, we considerably improve the corresponding lifting steps in both subroutines. As a result, generic position of the input system is never assumed, and thus our algorithm never demands for any change of coordinates. In addition, we significantly limit the types of involved exact operations, that is, we only use resultant and gcd computations as purely symbolic operations. The latter results are achieved by combining techniques from different fields such as (modular) symbolic computation, numerical analysis and algebraic geometry. We have implemented our algorithms as prototypical contributions to the C++-project CGAL. They exploit graphics hardware to expedite the symbolic computations. We have also compared our implementation with the current reference implementations, that is, LGP and Maple's Isolate for polynomial system solving, and CGAL's bivariate algebraic kernel for analyses and arrangement computations of algebraic curves. For various series of challenging instances, our exhaustive experiments show that the new implementations outperform the existing ones.

연구 동기 및 목표

• 좌표 변환을 요구하지 않고 실평면 대수 곡선 배열을 완전하고 인증 가능한 알고리즘으로 계산하는 것.
• 기호 연산을 결과식과 최대공약수 연산으로만 제한하여 고비용 대수적 변환에 의존도를 줄이는 것.
• GPU 가속 산술을 활용하여 기하학적으로 복잡한 곡선(특이점이 있거나 희박한 곡선 포함)의 효율성과 강건성을 향상시키는 것.
• 비표준 입력(예: 공유 좌표가 있는 경우)이더라도 원래 좌표계에서 정확하고 완전한 기하-위상적 분석을 가능하게 하는 것.
• CGAL에 통합된 실용적이고 고성능 구현을 제공하여 현재 기준 구현보다 뛰어난 성능을 발휘하는 것.

제안 방법

• 알고리즘은 배열 계산을 두 가지 주요 서브루틴으로 분해한다: 이원수 시스템의 실근을 고립하는 데 쓰이는 Bisolve와 단일 대수 곡선의 위상을 분석하는 데 쓰이는 GeoTop.
• Bisolve는 결과식 계산을 통해 x 및 y 방향으로 해를 투영한 후, 결과식 데이터와 호모토피적 고려를 융합한 새로운 포함 조건(predicate)을 적용하여 좌표 변경 없이 고립 상자(box)를 인증한다.
• 기호 연산으로는 이원수 다항식에 대한 결과식 계산과 일변수 다항식에 대한 최대공약수 계산 두 가지뿐이 사용되며, 이는 기호 연산 오버헤드를 최소화한다.
• 결과식과 최대공약수의 GPU 가속 계산을 통해 CPU 기반 계산 대비 100배 이상의 속도 향상을 달성한다.
• 모든 좌표 변환을 피함으로써 비표준 입력(예: 공유 좌표가 있는 경우)에서도 원래 좌표계에서 직접 계산이 가능하다.
• 필터링 기법과 인증된 근사 산술을 사용하여 고비용 기호 연산의 수를 줄이고, 핵심 단계에 집중한 계산을 가능하게 한다.

실험 결과

연구 질문

• RQ1 기호-수치 혼합 접근법이 기호 연산의 수를 최소화하면서도 완전성과 인증 가능성을 확보하고 평면 대수 곡선 배열을 계산할 수 있는가?
• RQ2 결과식과 최대공약수의 GPU 가속 계산이 대수 곡선을 위한 기호-수치 알고리즘의 성능을 어떻게 향상시킬 수 있는가?
• RQ3 제거 방법이 좌표 변환 없이도 강건성이나 효율성을 손상시키지 않고 얼마나 효과적으로 작용할 수 있는가?
• RQ4 제안된 방법이 특이점이 있거나 희박하거나 교차하는 등의 기하학적으로 복잡한 곡선을 이전 방법보다 더 효율적으로 처리할 수 있는가?
• RQ5 새로운 알고리즘이 Lgp, Maple의 Isolate, CGAL의 이원수 커널과 같은 기준 구현 대비 도전적인 사례에서 얼마나 높은 성능 향상을 보이는가?

주요 결과

• 제안된 알고리즘인 GeoTopAK_2는 특수 곡선의 조합이나 투영된 교차 곡선 등의 어려운 사례에서 기존 구현보다 뛰어난 성능을 보이며, 일부 벤치마크에서 최대 10배의 속도 향상을 기록했다.
• 특수 곡선 조합 10개에 대해 GeoTopAK_2는 1201.34초 만에 완료되었고, Ak_2는 4000초 후 타임아웃에 걸려 성능 면에서 뛰어난 확장성을 보였다.
• 8개 이상의 동일한 유형 곡선에 대해서는 Ak_2가 타임아웃에 걸렸지만, GeoTopAK_2는 1200초 이내에 완료되어 복잡성에 대한 뛰어난 내성성을 입증했다.
• 결과식과 최대공약수의 GPU 가속 계산으로 기호 연산의 성능 저하 문제를 해결하여 더 이상 알고리즘의 성능 저하 요인이 되지 않았다.
• 비표준 입력(예: 공유 좌표가 있는 경우)에서도 좌표 변환 없이 완전한 인증과 완전성을 확보했으며, 이는 강건성을 입증한다.
• 새로운 알고리즘을 웹 데모에 통합한 결과, 광범위한 사용이 이루어졌으며, 실용성과 사용성에 대한 실질적 의미를 보였다.