Visualizing Planar and Space Implicit Real Algebraic Curves with Singularities

This paper presents a new method for visualizing implicit real algebraic curves inside a bounding box in the 2-D or 3-D ambient space based on numerical continuation and critical point methods. The underlying techniques work also for tracing space curve in higher-dimensional space. Since the topology of a curve near a singular point of it is not numerically stable, the authors trace only the curve outside neighborhoods of singular points and replace each neighborhood simply by a point, which produces a polygonal approximation that is e-close to the curve. Such an approximation is more stable for defining the numerical connectedness of the complement of the projection of the curve in ℝ 2 , which is important for applications such as solving bi-parametric polynomial systems. The algorithm starts by computing three types of key points of the curve, namely the intersection of the curve with small spheres centered at singular points, regular critical points of every connected components of the curve, as well as intersection points of the curve with the given bounding box. It then traces the curve starting with and in the order of the above three types of points. This basic scheme is further enhanced by several optimizations, such as grouping singular points in natural clusters, tracing the curve by a try-and-resume strategy and handling “pseudo singular points”. The effectiveness of the algorithm is illustrated by numerous examples. This manuscript extends the proposed preliminary results that appeared in CASC 2018.

[1]  Guoliang Xu,et al.  Piecewise Rational Approximations of Real Algebraic Curves , 1997 .

[2]  Changbo Chen,et al.  A Continuation Method for Visualizing Planar Real Algebraic Curves with Singularities , 2018, CASC.

[3]  Alexandre Goldsztejn,et al.  Certified Parallelotope Continuation for One-Manifolds , 2013, SIAM J. Numer. Anal..

[4]  Michael Sagraloff,et al.  Visualizing Arcs of Implicit Algebraic Curves, Exactly and Fast , 2009, ISVC.

[5]  G. Stewart Perturbation theory for the singular value decomposition , 1990 .

[6]  Chee-Keng Yap,et al.  Planar Minimization Diagrams via Subdivision with Applications to Anisotropic Voronoi Diagrams , 2016, Comput. Graph. Forum.

[7]  Kai Jin,et al.  Isotopic epsilon-meshing of real algebraic space curves , 2014, SNC.

[8]  Abel J. P. Gomes A continuation algorithm for planar implicit curves with singularities , 2014, Comput. Graph..

[9]  Changbo Chen,et al.  Full Rank Representation of Real Algebraic Sets and Applications , 2017, CASC.

[10]  Richard E. Chandler,et al.  A tracking algorithm for implicitly defined curves , 1988, IEEE Computer Graphics and Applications.

[11]  Marc Moreno Maza,et al.  On Solving Parametric Polynomial Systems , 2012, Mathematics in Computer Science.

[12]  Changbo Chen,et al.  Penalty Function Based Critical Point Approach to Compute Real Witness Solution Points of Polynomial Systems , 2017, CASC.

[13]  H. Hong An efficient method for analyzing the topology of plane real algebraic curves , 1996 .

[14]  Anton Leykin,et al.  Robust Certified Numerical Homotopy Tracking , 2011, Foundations of Computational Mathematics.

[15]  Changbo Chen,et al.  Triangular decomposition of semi-algebraic systems , 2013, J. Symb. Comput..

[16]  Fabrice Rouillier,et al.  On the Topology of Real Algebraic Plane Curves , 2010, Math. Comput. Sci..

[17]  Fabrice Rouillier,et al.  Finding at Least One Point in Each Connected Component of a Real Algebraic Set Defined by a Single Equation , 2000, J. Complex..

[18]  Bernard Mourrain,et al.  On the computation of the topology of a non-reduced implicit space curve , 2008, ISSAC '08.

[19]  Changbo Chen,et al.  A Numerical Method for Computing Border Curves of Bi-parametric Real Polynomial Systems and Applications , 2016, CASC.

[20]  Yong Feng,et al.  Computing real witness points of positive dimensional polynomial systems , 2017, Theor. Comput. Sci..

[21]  Hélio Lopes,et al.  Robust adaptive polygonal approximation of implicit curves , 2002, Comput. Graph..

[22]  Wenyuan Wu,et al.  Finding points on real solution components and applications to differential polynomial systems , 2013, ISSAC '13.

[23]  Guillaume Moroz,et al.  A certified numerical algorithm for the topology of resultant and discriminant curves , 2017, J. Symb. Comput..

[24]  Lenore Blum,et al.  Complexity and Real Computation , 1997, Springer New York.

[25]  George E. Collins,et al.  Quantifier elimination for real closed fields by cylindrical algebraic decomposition , 1975 .

[26]  Fei Shen,et al.  Real Root Isolation of Polynomial Equations Based on Hybrid Computation , 2012, ASCM.

[27]  Changbo Chen,et al.  A Numerical Method for Analyzing the Stability of Bi-Parametric Biological Systems , 2016, 2016 18th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC).

[28]  Bo Dong,et al.  Robust continuation methods for tracing solution curves of parameterized systems , 2013, Numerical Algorithms.

[29]  Xiao-Shan Gao,et al.  On the Topology and Visualization of Plane Algebraic Curves , 2015, CASC.

[30]  Tsung-Lin Lee,et al.  HOM4PS-2.0: a software package for solving polynomial systems by the polyhedral homotopy continuation method , 2008, Computing.

[31]  Oliver Labs A List of Challenges for Real Algebraic Plane Curve Visualization Software , 2009 .

[32]  J. Hauenstein Numerically Computing Real Points on Algebraic Sets , 2011, Acta Applicandae Mathematicae.

[33]  Christian Stussak On reliable visualization algorithms for real algebraic curves and surfaces , 2013 .

[34]  Jack Bresenham,et al.  A linear algorithm for incremental digital display of circular arcs , 1977, CACM.

[35]  Raimund Seidel,et al.  On the exact computation of the topology of real algebraic curves , 2005, SCG.

[36]  Chee-Keng Yap,et al.  Complete subdivision algorithms, II: isotopic meshing of singular algebraic curves , 2008, ISSAC '08.

[37]  Anton Leykin,et al.  Newton's method with deflation for isolated singularities of polynomial systems , 2006, Theor. Comput. Sci..

[38]  Wenrui Hao,et al.  Algorithm 976 , 2017, ACM Trans. Math. Softw..