[논문 리뷰] Certifying solutions to overdetermined and singular polynomial systems over Q

Tulay Ayyildiz Akoglu, Jonathan D. Hauenstein외 1명|arXiv (Cornell University)|2014. 08. 12.

Polynomial and algebraic computation참고 문헌 41인용 수 6

한 줄 요약

이 논문은 유리수 ℚ 위의 과잉결정 및 특이 다항식 계수계의 근사해를 정밀하게 검증하기 위해 기호적-수치적 혼합 방법을 제안한다. 수치적 근에서 정확한 유리수 일변수 표현(RUR)을 계산함으로써, 뉴턴 반복과 정확한 산술을 사용해 조건이 나쁜 계수계를 잘 조절된 RUR로 변환하고, α-이론을 통해 검증을 가능하게 한다. 이는 다항식 시간 복잡도를 가지며, RUR가 발견되면 조기 종료가 가능하다.

ABSTRACT

This paper is concerned with certifying that a given point is near an exact root of an overdetermined or singular polynomial system with rational coefficients. The difficulty lies in the fact that consistency of overdetermined systems is not a continuous property. Our certification is based on hybrid symbolic-numeric methods to compute the exact "rational univariate representation" (RUR) of a component of the input system from approximate roots. For overdetermined polynomial systems with simple roots, we compute an initial RUR from approximate roots. The accuracy of the RUR is increased via Newton iterations until the exact RUR is found, which we certify using exact arithmetic. Since the RUR is well-constrained, we can use it to certify the given approximate roots using alpha-theory. To certify isolated singular roots, we use a determinantal form of the "isosingular deflation", which adds new polynomials to the original system without introducing new variables. The resulting polynomial system is overdetermined, but the roots are now simple, thereby reducing the problem to the overdetermined case. We prove that our algorithms have complexity that are polynomial in the input plus the output size upon successful convergence, and we use worst case upper bounds for termination when our iteration does not converge to an exact RUR. Examples are included to demonstrate the approach.

연구 동기 및 목표

• 일관성이 연속적이지 않으며, 보편적인 하한값이 너무 낙관적이지 않은 과잉결정 및 특이 다항식 계수계에서 근사근을 검증하는 데 도전하는 것.
• 직접 근사근에서 정확한 유리수 일변수 표현(RUR)을 계산함으로써, 지나치게 보수적인 하한값을 피하는 실용적인 검증 방법을 개발하는 것.
• 정규화 및 반복적 정밀화를 통해 특이 근을 단순한 근으로 변환함으로써 α-이론 검증을 과잉결정 및 특이 계수계로 확장하는 것.
• 수치적 근 찾기와 정확한 산술을 조합하여 정확성과 효율성을 확보하고, 수렴 시 입력 및 출력 크기의 다항식 시간 복잡도를 달성하는 것.

제안 방법

• 과잉결정 또는 특이 계수계의 경우, 호모토피 연속법을 사용해 무작위로 잘 조절된 부분계수계의 근사근을 계산한다.
• 히우리스틱 필터링을 통해 고립점 또는 유리수 성분에 해당하는 후보 근사근을 선택한다.
• 선택된 근사근에서 라그랑주 보간법과 유리수 재구성 기법을 사용해 초도 RUR를 구성한다.
• 지역 또는 전역 뉴턴 반복을 사용해 RUR를 정밀화하여 정확도를 높이고, 정확한 산술을 통해 정확한 RUR에 도달할 때까지 반복한다.
• 행렬식 형태를 통한 등급일치 정규화를 적용하여 새로운 변수를 도입하지 않고도 다항식을 추가함으로써 특이 근을 단순한 근으로 변환한다.
• 최종 RUR를 α-이론 또는 헤르미트 행렬을 사용해 검증하여 원래의 근사근이 정확한 해 근처에 있음을 확인한다.

실험 결과

연구 질문

• RQ1 과잉결정 다항식 계수계의 근사근은 보편적인 하한값보다 정확한 RUR를 사용해 높은 신뢰도로 검증할 수 있는가?
• RQ2 다항식 계수계의 특이 근은 어떻게 단순한 근으로 변환되어 표준 검증 기법을 적용할 수 있는가?
• RQ3 근사근에서 정확한 RUR를 계산하는 데 드는 계산 복잡도는 무엇이며, 실세계 계수계에 대해 실용적으로 만들 수 있는가?
• RQ4 하이브리드 기호-수치적 방법은 잘 조절된 RUR에 대해 입력 및 출력 크기의 다항식 시간 수렴을 달성할 수 있는가?

주요 결과

• 제안된 방법은 올바른 근 후보가 선택된 경우, 뉴턴 반복을 통해 정확한 RUR를 성공적으로 계산하며 수렴이 보장된다.
• Caprasse 계수계의 경우, 각도가 2인 4개의 유리수 성분을 식별하였고, 분해체는 ℚ[√3]이며 8개의 중복근의 구조를 확인하였다.
• cyclic-4 계수계에서, 계수의 유리수를 가진 12차 RUR를 계산하였고, 8개의 임bedded 등급일치 점과 두 개의 기약 곡선을 정확히 식별하였다.
• 162개의 해를 가진 cyclic-9 계수계에서는 54개의 정상점과 54개의 이중점에 대해 최대 12차까지의 고차 RUR를 계산하였고, 아이디얼의 구조를 확인하였다.
• 시험 크기가 작을 경우 조기 종료가 이루어져 보편적인 하한값보다 실용성을 입증하였다.
• 성공적인 수렴 시 복잡도는 입력 및 출력 크기의 다항식 시간에 해당하며, 수렴하지 않는 경우의 최악의 경우 복잡도도 제시되었다.