Solving nonlinear polynomial systems in the barycentric Bernstein basis
Author(s)Reuter, Martin; Mikkelsen, Tarjei Sigurd; Sherbrooke, Evan C.; Maekawa, Takashi; Patrikalakis, Nicholas M.
MetadataShow full item record
We present a method for solving arbitrary systems of N nonlinear polynomials in n variables over an n-dimensional simplicial domain based on polynomial representation in the barycentric Bernstein basis and subdivision. The roots are approximated to arbitrary precision by iteratively constructing a series of smaller bounding simplices. We use geometric subdivision to isolate multiple roots within a simplex. An algorithm implementing this method in rounded interval arithmetic is described and analyzed. We find that when the total order of polynomials is close to the maximum order of each variable, an iteration of this solver algorithm is asymptotically more efficient than the corresponding step in a similar algorithm which relies on polynomial representation in the tensor product Bernstein basis. We also discuss various implementation issues and identify topics for further study.
DepartmentMassachusetts Institute of Technology. Department of Mechanical Engineering
Reuter, Martin et al. “Solving Nonlinear Polynomial Systems in the Barycentric Bernstein Basis.” The Visual Computer 24.3 (2007) : 187-200. © 2007 Springer-Verlag
Author's final manuscript