[논문 리뷰] Efficient Computation of Newton Polytopes of Specialized Resultants

Ioannis Z. Emiris, Vissarion Fisikopoulos외 2명|arXiv (Cornell University)|2011. 08. 30.

Advanced Numerical Analysis Techniques인용 수 1

한 줄 요약

이 논문은 정점 및 면 쿼리에 대한 오라클을 사용하여, 비용이 많이 드는 다면체 이동을 피하면서도 전문화된 결과 다항식의 뉴턴 다면체를 효율적이고 출력에 민감하게 계산하는 알고리즘을 제시한다. 높은 성능을 달성하여 5-, 6-, 7차원 다면체를 각각 최대 35,000개의 정점으로 2시간 이내에 처리하며, 해싱을 통해 최대 100배 빠르게 계산을 가속화하여 저차원에서 타원기하 도구를 능가하고, 빠르고 정확한 내부/외부 부피 근사치를 가능하게 한다.

ABSTRACT

We design an algorithm to compute the Newton polytope of the resultant, known as resultant polytope, or its orthogonal projection along a given direction. The resultant is fundamental in algebraic elimination, optimization, and geometric modeling. Our algorithm exactly computes vertex- and halfspace-representations of the polytope using an oracle producing resultant vertices in a given direction, thus avoiding walking on the polytope whose dimension is alpha-n-1, where the input consists of alpha points in Z^n. Our approach is output-sensitive as it makes one oracle call per vertex and facet. It extends to any polytope whose oracle-based definition is advantageous, such as the secondary and discriminant polytopes. Our publicly available implementation uses the experimental CGAL package triangulation. Our method computes 5-, 6- and 7-dimensional polytopes with 35K, 23K and 500 vertices, respectively, within 2hrs, and the Newton polytopes of many important surface equations encountered in geometric modeling in <1sec, whereas the corresponding secondary polytopes are intractable. It is faster than tropical geometry software up to dimension 5 or 6. Hashing determinantal predicates accelerates execution up to 100 times. One variant computes inner and outer approximations with, respectively, 90% and 105% of the true volume, up to 25 times faster.

연구 동기 및 목표

• 대수기하학 및 기하 모델링에서 결과 다항식의 고차원 뉴턴 다면체의 계산 불가능성 문제를 해결하기 위해.
• 차원 수와 점의 수가 증가할수록 성능이 급격히 떨어지는 다면체 이동 알고리즘의 한계를 극복하기 위해.
• 모든 정점과 면을 완전히 열거하지 않고도 오라클 호출 수를 최소화하는 출력에 민감한 방법을 개발하기 위해.
• 이러한 접근을 정차 다면체 및 판별 다면체와 같은 다른 다면체로 일반화하여 오라클 기반 정의를 갖는 경우에 적용하기 위해.
• 기하 모델링 및 최적화 분야의 실세계 응용을 위해 실용적이고 고성능의 구현을 제공하기 위해.

제안 방법

• 주어진 방향에서 결과 다면체의 정점과 면을 반환하는 오라클을 사용하여, 전체 다면체의 구조를 거쳐가는 것을 피하면서도 정확한 계산을 가능하게 한다.
• 정점 및 반공간 표현 전략을 활용하여 각 정점과 면마다 정확히 한 번의 오라클 호출만을 수행함으로써 출력에 민감한 성능을 확보한다.
• 실험적 기하 계산을 위해 고차원 다면체 구축을 지원하는 CGAL 삼각분할 패키지를 활용한다.
• 행렬식 조건자의 해싱을 도입하여 오라클 쿼리의 성능을 가속화하고, 실행 시간을 최대 100배 빠르게 한다.
• 내부 및 외부 부피 근사치를 계산하는 변형을 구현하여 각각 진짜 부피의 90%와 105%를 제공함으로써 성능 향상을 달성한다.

실험 결과

연구 질문

• RQ1 결과 다면체의 뉴턴 다면체는 전체 다면체의 구조를 거쳐가는 것 없이 효율적으로 계산될 수 있는가?
• RQ2 고차원 결과 다항식에 대해 기존의 다면체 이동 기반 접근 방식과 비교했을 때 오라클 기반 접근 방식은 확장성과 성능 면에서 어떻게 다른가?
• RQ3 기하 조건자의 해싱은 결과 다면체의 계산을 얼마나 가속화할 수 있는가?
• RQ4 이 방법은 정차 다면체 및 판별 다면체와 같이 오라클 기반 정의를 갖는 다른 다면체로 일반화될 수 있는가?
• RQ5 정확한 계산과 비교했을 때 내부 및 외부 부피 근사치의 정확도와 효율성은 어떠한가?

주요 결과

• 알고리즘은 각각 35,000개, 23,000개, 500개의 정점을 가진 5-, 6-, 7차원 결과 다면체를 2시간 이내에 계산한다.
• 기하 모델링에서 중요한 표면 방정식의 뉴턴 다면체는 1초 이내에 계산되지만, 그 정차 다면체는 여전히 계산이 불가능하다.
• 차원 수가 5 또는 6 이하일 경우 기존의 타원기하 소프트웨어를 능가하는 성능을 보인다.
• 행렬식 조건자의 해싱은 계산 속도를 최대 100배 빠르게 한다.
• 내부 및 외부 부피 근사치 변형은 정확한 계산 대비 최대 25배 빠르게 동작하며, 각각 진짜 부피의 90%와 105%를 제공한다.
• 이러한 접근은 정차 다면체 및 판별 다면체와 같이 오라클 기반 정의를 갖는 다른 다면체로 일반화 가능하다.