A certified numerical algorithm for the topology of resultant and discriminant curves

详细信息    查看全文

• 作者：Ré ; mi Imbach^a ; ^b ; ^{Remi.Imbach@inria.fr} ; Guillaume Moroz^a ; ^b ; ^{Guillaume.Moroz@inria.fr} ; Marc Pouget^a ; ^b ; ^{Marc.Pouget@inria.fr}
• 关键词：Topology of algebraic curves ; Resultant ; Discriminant ; Subresultant ; Numerical algorithm ; Singularities ; Interval arithmetic ; Node and cusp singularities
• 刊名：Journal of Symbolic Computation
• 出版年：2017
• 出版时间：May-June 2017
• 年：2017
• 卷：80
• 期：part_P2
• 页码：285-306
• 全文大小：748 K
• 卷排序：80

文摘

Let CC be a real plane algebraic curve defined by the resultant of two polynomials (resp. by the discriminant of a polynomial). Geometrically such a curve is the projection of the intersection of the surfaces P(x,y,z)=Q(x,y,z)=0P(x,y,z)=Q(x,y,z)=0 (resp. P(x,y,z)=∂P∂z(x,y,z)=0), and generically its singularities are nodes (resp. nodes and ordinary cusps). State-of-the-art numerical algorithms compute the topology of smooth curves but usually fail to certify the topology of singular ones. The main challenge is to find practical numerical criteria that guarantee the existence and the uniqueness of a singularity inside a given box B, while ensuring that B   does not contain any closed loop of CC. We solve this problem by first providing a square deflation system, based on subresultants, that can be used to certify numerically whether B contains a unique singularity p   or not. Then we introduce a numeric adaptive separation criterion based on interval arithmetic to ensure that the topology of CC in B is homeomorphic to the local topology at p. Our algorithms are implemented and experiments show their efficiency compared to state-of-the-art symbolic or homotopic methods.