Algorithms for polynomials in Bernstein form
Computer Aided Geometric Design
What every computer scientist should know about floating-point arithmetic
ACM Computing Surveys (CSUR)
A highly parallel algorithm for approximating all zeros of a polynomial with only real zeros
Communications of the ACM
Recursive de Casteljau bisection and rounding errors
Computer Aided Geometric Design
The Bernstein polynomial basis: A centennial retrospective
Computer Aided Geometric Design
Hi-index | 0.00 |
The surface/curve intersection problem, through the resultants process results in a high degree (n=100) polynomial equation on [0,1] in the Bernstein basis. The knowledge of multiplicities of the roots is critical for the topological coherence of the results. In this aim, we propose an original two-steps algorithm based on successive differentiations which separates any root (even multiple) and guarantees that the assumptions of Newton global convergence theorem are satisfied. The complexity is @q(n^4) but the algorithm can easily be parallelized. Experimental results show its efficiency when facing ill-conditioned polynomials.