[논문 리뷰] An implementation of CAD in Maple utilising McCallum projection

Matthew England|arXiv (Cornell University)|2013. 02. 26.

Advanced Numerical Analysis Techniques참고 문헌 23인용 수 6

한 줄 요약

이 논문은 매클럴의 CADW 알고리즘을 매플에서 완전히 구현한 것으로, 분할 다항식을 사용한 그의 개선된 투영 연산자를 활용하여 순서 불변성의 원통형 대수적 분해를 생성한다. 이는 CADW의 유일한 완전한 구현이며, Qepcad와 같은 기존 도구들이 영향을 받는 nullification 문제로 인해 실패하는 문제들에 대해서도 해결할 수 있게 해주며, 특정 논리 문제를 위한 최소 분해를 지원한다.

ABSTRACT

Cylindrical algebraic decomposition (CAD) is an important tool for the investigation of semi-algebraic sets. Originally introduced by Collins in the 1970s for use in quantifier elimination it has since found numerous applications within algebraic geometry and beyond. Following from his original work in 1988, McCallum presented an improved algorithm, CADW, which offered a huge increase in the practical utility of CAD. In 2009 a team based at the University of Western Ontario presented a new and quite separate algorithm for CAD, which was implemented and included in the computer algebra system Maple. As part of a wider project at Bath investigating CAD and its applications, Collins and McCallum's CAD algorithms have been implemented in Maple. This report details these implementations and compares them to Qepcad and the Ontario algorithm. The implementations were originally undertaken to facilitate research into the connections between the algorithms. However, the ability of the code to guarantee order-invariant output has led to its use in new research on CADs which are minimal for certain problems. In addition, the implementation described here is of interest as the only full implementation of CADW, (since Qepcad does not currently make use of McCallum's delineating polynomials), and hence can solve problems not admissible to other CAD implementations.

연구 동기 및 목표

• McCallum의 CADW 알고리즘을 매플에 구현하여, 분할 다항식을 사용한 그의 개선된 투영 연산자의 완전하고 작동 가능한 버전을 제공한다.
• 다양한 CAD 계산 접근 방식 간의 구조적 및 알고리즘적 연관성, 특히 투영-리프팅 방법과 정규 체인 방법 간의 연관성에 대한 연구를 가능하게 한다.
• 특정 논리 문제를 위한 최소 분해가 필요한 새로운 응용 분야를 지원하기 위해 순서 불변성의 CAD를 제공한다.
• Qepcad와 같은 기존 도구의 한계를 해결하며, 이는 McCallum의 분할 다항식을 구현하지 않아 일부 입력에 대해 실패할 수 있기 때문이다.
• 연구 및 실용적 용도를 위해 사용자 정의 변수 순서와 출력 제어 기능을 지원하는 신뢰성 있고 확장 가능한 CAD 프레임워크를 매플에 제공한다.

제안 방법

• McCallum의 1988년 CADW 알고리즘을 매플에 맞게 변형하여, 분할 가능성 이론과 최소 분할 다항식을 기반으로 한 동일한 투영 연산자를 사용한다.
• 브라운의 최소 분할 다항식을 통합하여 최종 CAD 출력에서 순서 불변성을 보장한다. 이는 고급 응용 분야에서 핵심적인 기능이다.
• 변수 순서, 셀 생성, 각 셀 위의 스택 형성 기능을 지원하는 매플 내에서 완전한 투영-리프팅 파이프라인을 구현한다.
• 출력 유형에 대한 사용자 제어를 가능하게 하며, 입력 다항식이 반드시 필요로 하지 않더라도 순서 불변성을 강제 적용할 수 있도록 한다.
• 성능 및 정확성을 평가하기 위해 벤치마크 문제 세트를 대상으로 이 구현체(PCAD)를 Qepcad(QCAD) 및 RegularChains 기반 CAD(TCAD)와 비교한다.
• 리프팅 과정에서의 nullification 상황을 처리하기 위해 계산을 계속 진행할 수 있도록 하여, 안전한 경우에 한해 조기 실패를 방지한다.

실험 결과

연구 질문

• RQ1 McCallum의 CADW 알고리즘이 분할 다항식을 사용하여 순서 불변성을 보장하는 방식으로 매플에 완전하고 정확하게 구현될 수 있는가?
• RQ2 Qepcad가 nullification 문제로 실패하는 문제들에 대해서도 새로운 매플 기반 CADW 구현이 성공적으로 해결할 수 있는가?
• RQ3 이 구현체의 순서 불변 출력이 최소 CAD와 같은 새로운 응용 분야, 예를 들어 특정 논리 문제를 위한 최소 분해에 얼마나 기여하는가?
• RQ4 매플 기반 CADW 구현체의 성능은 Qepcad 및 RegularChains 기반 CAD와 비교해 셀 수와 계산 시간 측면에서 어떻게 나타나는가?
• RQ5 이 구현체는 브라운-맥컬럼 연산자와 같은 향후 개선 사항을 통합하여 셀 수를 줄이고 효율성을 향상시킬 수 있는가?

주요 결과

• 이 구현체는 Qepcad가 그의 분할 다항식을 구현하지 않기 때문에, McCallum의 CADW 알고리즘에 대한 유일한 완전하고 작동 가능한 버전이다.
• 위트니의 우산과 x축 타원 문제에서, Qepcad가 nullification로 인해 실패하는 동안 이 구현체는 성공적으로 CAD를 계산하여 강건성을 입증했다.
• x축 타원 문제에서 매플 구현체는 20,225개의 셀을 생성했으며, Qepcad의 출력과 정확히 일치하지만, Qepcad는 nullification로 인해 실패했다고 보고했다.
• 순서 불변성을 강제로 적용했을 경우, $f = zy - x^2$ 예제의 경우 23개의 셀을 정확히 계산했으며, sign-invariant 모드에서는 21개의 셀을 생성했다.
• 평균적으로 Qepcad는 속도와 셀 수 측면에서 매플 구현체를 앞서지만, Qepcad가 실패하는 문제들에 대해서는 매플 버전이 필수적이다.
• 이 구현체는 분할 다항식을 올바르게 처리하고 순서 불변성이 보장되므로, 최근 연구에서 요구하는 특정 문제를 위한 최소 CAD 생성을 지원한다.