1 / 43

Solving Multivariate Polynomial Systems using Hyperplane Arithmetic and Linear Programming

Solving Multivariate Polynomial Systems using Hyperplane Arithmetic and Linear Programming. Iddo Hanniel. The Problem. Consider the following set of n polynomial equations: in R n . We seek the simultaneous solution, ,

bglass
Download Presentation

Solving Multivariate Polynomial Systems using Hyperplane Arithmetic and Linear Programming

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Solving Multivariate Polynomial Systemsusing Hyperplane Arithmetic and Linear Programming IddoHanniel

  2. The Problem Consider the following set ofnpolynomial equations: inRn. We seek the simultaneous solution, , such that for alli =1,…,n.

  3. Example Application – Apollonius 10th Problem 3

  4. Apollonius Problem Applications Voronoi diagram of circles Artillery sound range locator 4

  5. Example Voronoi Diagram Vertices p 2 distance constraint + 3 tangency constraint = 5 equations in 5 unknowns (can we do in less?) 5

  6. Example application – Kinematics of Parallel Mechanisms 6

  7. Previous Work on Multivariate Polynomial Solvers Algebraic Solvers: Grobner bases [Cox et al., 1992]. Multivariate Sturm sequences [Milne, 1992]. Homotopy methods [Wampler 05]. Subdivision Solvers: Interval arithmetic solvers [Merlet 2000] Subdivision solvers based on the Bernstein/Bezier/B-spline properties [Sherbrooke and Patrikalakis, 1993], [Elber and Kim, 2001]. 7

  8. Interval Arithmetic For two intervals [a, b] and [c, d] that are subsets of (-∞,+ ∞): [a, b] + [c, d] = [a + c, b + d] [a, b] − [c, d] = [a − d, b − c] [a, b] × [c, d] = [min (a × c, a × d, b × c, b × d), max (a × c, a × d, b × c, b × d)] [a, b] ÷ [c, d] = [min (a ÷ c, a ÷ d, b ÷ c, b ÷ d), max (a ÷ c, a ÷ d, b ÷ c, b ÷ d)] when 0 is not in [c, d]. 8

  9. Interval of Powers Odd n: Even n: 9

  10. Interval Based Multivariate Bisection Algorithm root_isolation_in_box(Multivariate F(x1,..,xn), Box B[a1,b1]×.. ×[an,bn], Output {Boxes}) If (max(bi-ai) < ε) append B to output boxes and return. Evaluate [Fmin,Fmax] = F([a1,b1], ..,[an,bn]). If [Fmin,Fmax] does not contain 0, return. Split B into two sub-boxes, B1,B2 at 0.5(bm+am), where m = {i: (bi-ai) is maximal}. root_isolation_in_box(F,B1,{Boxes}) root_isolation_in_box(F,B2,{Boxes}) 10

  11. Problem with Intervals 11

  12. The Basic Bezier/B-Spline Solver A subdivision step in a multidimensional space, using B-spline/Bezier subdivision. A multivariate Newton-Raphson (NR)numeric step. [Elber and Kim, 2001] 12

  13. Subdivision Illustration – Univariate Case Subdivision. CH containment. 13

  14. Multivariate Bezier/B-spline Tensor Representation 14

  15. Basic Bezier Based Multivariate Subdivision Algorithm root_isolation_in_box(Multivariate F(x1,..,xn), Box B[a1,b1]×.. ×[an,bn], Output {Boxes}) • If (max(bi-ai) < ε) append B to output boxes and return. • If CH does not contain 0, return. • Split B into two sub-boxes, B1,B2 at 0.5(bm+am), where m = {i: (bi-ai) is maximal}. • root_isolation_in_box(F,B1,{Boxes}) • root_isolation_in_box(F,B2,{Boxes}) 15

  16. Research Questions When do we stop the subdivision? Single solution test [Hanniel & Elber 07, Wei et al. 11]. • Can we handle univariate solutions? [Barton, Hanniel & Elber 11] • What about bivariate solutions? [Mizrahi & Elber 15] • Can preconditioning improve the algorithm? [Mourrain & Pavone 09] • Can we reduce the domain more efficiently? • Projected Polyhedron [Sherbrooke & Patrikalakis 93]. • Projected curves [Mourrain & Pavone 09]. • Linear Programming on Bernstein Polyhedra [Funfzig et al. 09] 16

  17. Domain Reduction using Bounding Hyperplanes The basic idea: Bound the function Fi=0 by a pair of parallel hyper-planes in Rd. The root can only be at the intersection of the bounding hyperplanes. F1 F2 17

  18. Linear Programming Under a set of linear constraints (inequalities) of type: Maximize the linear objective function of type: 18

  19. Linear Programming Example Constraints: • Maximize f(x1,x2)=x2 • Maximize f(x1,x2)=x1 19

  20. Bounding Hyperplane Construction • Promote Fi from a scalar function to a hyper-surface in Rn+1, . • Evaluate the gradient (i.e., the normal) at the sub-domain mid-point and project all control points of the hyper-surface onto it. • The two hyper-planes orthogonal to the gradient and passing through the extreme projection points, bound Fi. 20

  21. Bounding Hyperplane Construction Illustration (d=1) [ ] xn+1 = 0 21

  22. Domain Reduction Procedure using LP Bound Fi from above and below by two hyperplanes: HPmin ≤Fi ≤HPmax . Perform 2n LP queries (for objective functions min/maxxi ). Reduce the domain to box with min/maxxi. 22

  23. Why is Tensor Representation a Problem? [Bamberger and Shoham, 2007] 23

  24. Why is Tensor Representation a Problem? Grows exponentially with number of variables n (O(l+1)n). Dense representation: Even for sparse systems, all coefficients are required for Bezier properties (e.g., CH). Standard power form requires only O(nl) and in many cases less (e.g., linear, constant). 24

  25. Research Question Can we use a sparse representation and still benefit the features of Bezier representation? Previous work: • Elber & Grandine 09: Expression Trees • Funfzig et al. 09: Bernstein polyhedra and LP 25

  26. From Interval Arithmetic to Hyperplane Arithmetic Interval arithmetic: • Generalization to bounding functions:

  27. Hyperplane Arithmetic Example • Bounding sum: • Bounding first monom: • Bounding second monom:

  28. Hyperplane Arithmetic Example x1 - 0.25 ≤ x12≤x1 x2 - 0.25 ≤ x22≤ x2 x1 + x2 – 1.5 ≤x12 + x22 - 1≤x1+x2– 1

  29. LP Algorithm using Hyperplane Arithmetic Bound Fi polynomials by two parallel hyperplanes, using hyperplane arithmetic of bounds on its monomial terms. Perform 2n LP computations to compute new [ximin,ximax] values for each xi. Continue recursively on the new reduced domain.

  30. Experimental Results Tested Algorithms: • IRIT Naïve - basic with mid parameter subdivision – no reduction. • IRIT PP - variant of [Sherbrooke & Patrikalakis 93]. • Bernst LP - variant of [Funfzig et al. 09]. • Naïve LP - our algorithm with dense Bezier representation. • HPA LP - our algorithm with hyper-plane arithmetic. Test cases: • Easily scalable by number of variables • Easily computable roots. 30

  31. Experimental Results – Hyper-cylinders What does it look like? How many roots are there? 31

  32. Results Hyper-cylinders (I) ms 32 n

  33. Results Hyper-cylinders (II) ms 33 n

  34. Extension to Degree 3 34

  35. Results – Degree 3 (I) ms 35 n

  36. Results – Degree 3 (Opt vs. No Opt) ms n Tighter (optimized) bounding hyperplanes specialized for x3. 36

  37. Broyden Tridiagonal System • Classic function in optimization community (where not all root are necessarily required). • Easily scaled by number of variables. • Two roots in [-2,2]n . 37

  38. Broyden Results (I) ms 38 n

  39. Broyden Results (II) ms 39 n

  40. A Dense System • Taken from [Tsidon, Hanniel & Kesslassi 2012]. • Dense system easily scaled by number of variables. • One root in (0,1]n . 40

  41. Dense System (cont.) Solution formula: What does the geometry of these surfaces look like?

  42. Results – Dense System ms 42 n

More Related