[논문 리뷰] Deflation and Certified Isolation of Singular Zeros of Polynomial Systems

Angelos Mantzaflaris, Bernard Mourrain|arXiv (Cornell University)|2011. 01. 17.

Polynomial and algebraic computation참고 문헌 12인용 수 4

한 줄 요약

이 논문은 국소 환의 구조와 인증된 수치 계산을 활용하여 다항식 계수계의 특이 고립근을 인증하는 기호-수치 알고리즘을 제시한다. 중간 선형 시스템 크기를 줄이는 개선된 역계산 시스템 계산과 미분을 이용한 한 번의 단계에서의 정규화 기법을 도입하여 뉴턴 해법에서의 이차 수렴을 가능하게 하고, 간격 산술을 통한 존재성과 유일성의 인증을 보장한다.

ABSTRACT

We develop a new symbolic-numeric algorithm for the certification of singular isolated points, using their associated local ring structure and certified numerical computations. An improvement of an existing method to compute inverse systems is presented, which avoids redundant computation and reduces the size of the intermediate linear systems to solve. We derive a one-step deflation technique, from the description of the multiplicity structure in terms of differentials. The deflated system can be used in Newton-based iterative schemes with quadratic convergence. Starting from a polynomial system and a small-enough neighborhood, we obtain a criterion for the existence and uniqueness of a singular root of a given multiplicity structure, applying a well-chosen symbolic perturbation. Standard verification methods, based eg. on interval arithmetic and a fixed point theorem, are employed to certify that there exists a unique perturbed system with a singular root in the domain. Applications to topological degree computation and to the analysis of real branches of an implicit curve illustrate the method.

연구 동기 및 목표

• 계수값이 근사적이거나 조건이 나쁜 경우에도 특이 고립근을 인증하는 데 도전하는 문제를 해결하기 위해.
• 다중근에 대해 뉴턴 유형 해법에서의 이차 수렴을 복원하기 위해 다중근의 구조를 이용해 계수계를 정규화하는 방법을 개발하기 위해.
• 기호적 편향과 간격 산술을 통해 존재성과 유일성을 보장하는 특이근의 인증된 검증을 제공하기 위해.
• 중복 계산을 피하고 중간 선형 시스템 크기를 최소화하여 역계산 시스템 계산의 계산 비용을 줄이기 위해.
• 특이점에서의 은닉 곡선과 곡면의 위상적 분석을 가능하게 하기 위해 특이점에서의 실수 분지 수를 계산하기 위해.

제안 방법

• 이 방법은 다중근에서의 원시-쌍대 구조를 개선된 역계산 시스템 알고리즘을 사용하여 계산하며, 중복 작업을 피한다.
• 다중근의 구조에 기반하여 쌍대 공간의 미분 함수를 원래 방정식에 적용하여 한 번의 단계에서 정규화 시스템을 유도한다.
• 모르랭의 통합 기법을 통해 쌍대 기저의 인증된 수치 계산을 통해 정규화된 시스템을 구성하여 정확성을 보장한다.
• 기호적 편향을 적용하여 특이근을 편향된 시스템의 정확한 다중근으로 변환함으로써 인증을 가능하게 한다.
• 존재성과 유일성의 검증은 간격 산술과 고정점 정리에 의해 수행된다.
• 쌍대 공간에서의 곱셈 연산자 트레이스를 사용하여 특이점에서의 위상 차수와 실수 분지 수를 계산한다.

실험 결과

연구 질문

• RQ1 어떻게 효율적으로 다항식 계수계의 특이 고립근의 다중근 구조를 계산할 수 있는가?
• RQ2 다중근에 대해 뉴턴 방법에서 이차 수렴을 복원하는 데 가장 효과적인 정규화 시스템을 구성하는 방법은 무엇인가?
• RQ3 오직 수치적 근사와 기호적 편향만을 사용하여 특이근의 존재성과 유일성을 어떻게 인증할 수 있는가?
• RQ4 역계산 시스템을 계산하기 위해 필요한 최소한의 중간 선형 시스템 크기는 무엇이며, 이를 어떻게 줄일 수 있는가?
• RQ5 대수적 불변량을 사용하여 특이점에서의 은닉 곡선의 실수 분지 수를 어떻게 계산할 수 있는가?

주요 결과

• 제안된 역계산 시스템 계산은 매크로일의 방법에 비해 중간 선형 시스템 크기를 줄이며, 예제 6.2에서는 9×9 행렬을 사용한 반면 매크로일의 방법은 30×15 행렬을 사용한다.
• 예제 7.1에서 이 방법은 9×8 행렬을 사용하여 21개 방정식의 정규화 시스템을 계산하지만, 매크로일의 방법은 30×15 행렬이 필요하다.
• 다중근 131을 가진 DZ1 시스템에 대해 이 방법은 1450×524 행렬을 사용하며, 매크로일의 방법의 4004×1365 행렬보다 훨씬 작다.
• 이 방법은 특이점에서의 실수 분지 수를 성공적으로 계산한다: 예제 6.2에서 Br(f, (0,0)) = 6이며, 이는 tdeg = 3에서 유도된다.
• 예제 7.2에서 이 방법은 206×45 행렬을 사용하여 다중근 18을 계산하지만, 매크로일의 방법은 360×165 행렬을 필요로 하여 상당한 감소가 이루어졌다.
• 이 방법은 낮은 깊이의 구조에서 정규화를 중단함으로써 부분 정규화를 가능하게 하여, 예제 7.2에서 세 개의 방정식만을 가진 정규화 시스템을 얻었으며, 64인 비영인 야코비안 행렬식으로 검증되었다.