1 / 34

Subdivision Termination Criteria in Subdivision Multivariate Solvers

Subdivision Termination Criteria in Subdivision Multivariate Solvers. Iddo Hanniel, Gershon Elber CGGC, CS, Technion. The Problem. Consider the following set of d polynomial equations: In R d . We seek the simultaneous solution, ,

rhea
Download Presentation

Subdivision Termination Criteria in Subdivision Multivariate Solvers

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. Subdivision Termination Criteria in Subdivision Multivariate Solvers Iddo Hanniel, Gershon Elber CGGC, CS, Technion

  2. The Problem Consider the following set ofdpolynomial equations: InRd. We seek the simultaneous solution, , such that for alli =1,…,d.

  3. Motivation for the Solver • Many geometric problems are reduced to the • simultaneous solution of a set of constraints. • Ray-surface, curve-curve and surface-surface intersections. • Voronoi diagram of curves. • Curve-curve bi-tangents, convex hull of curves. • Minimum enclosing circle/sphere of a set of curves/surfaces.

  4. Example – Intersection of 3 Explicit Surfaces in R3 Here, one can find four solutions that satisfy the equations. We are interested in finding all solutions.

  5. Previous Work on Multivariate Polynomial Solvers Algebraic Solvers: Grobner bases [Cox et al., 1992]. Multivariate Sturm sequences [Milne, 1992]. Geometric Solvers: Based on the Bernstein/Bezier/Bspline properties. Bezier clipping methods [Sherbrooke and Patrikalakis, 1993], [Mourrain and Pavone, 2005]. Subdivision step and numeric improvement [Elber and Kim, 2001].

  6. The Bezier/B-Spline Solver A subdivision stage in a multidimensional space, using B-spline/Bezier subdivision. A multivariate Newton-Raphson (NR)numeric step. Elber and Kim, 2001: Geometric Constraint Solver using Multivariate Rational Spline Functions

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

  8. The Termination Question Question: When do we stop the subdivision? We can identify domain cells with no root by the CH property. Can we do any better? • If we know there is at most a single root in some domain cell, we can move to the NR step. How can we know that?

  9. The Hodograph – Univariate Case

  10. Single Solution Guarantee

  11. Extension to Higher Dimensions Theorem: Givendimplicit hyper-surfacesFi(u) = 0, i = 1,….,d, inRd, there is at most one common solution if where0is the origin. Elber and Kim, 2001: Geometric Constraint Solver using Multivariate Rational Spline Functions

  12. Implementing the Single Solution Test How can this condition be tested efficiently inRd? Clearly the gradient field of Fi(u) provides a bound on the normal space of Fi(u). Optimal bounding cones over vectors in Rd could be found in average linear time.

  13. Implementing the SST ( Cont.) The complementary cone can easily be derived, given the normal cone. and share the same axis but with complementary angles.

  14. Implementing the SST (Cont.) But now we have to test for the intersection of d complementary cones in Rd. A difficult task, even in R3 ! d = 3

  15. Implementing the SST (Cont.) Reexamine the univariate case (two curves): …intersect iff the parallel lines domain’s Intersection… ….is not confined to the unit circle Two cones of two planar curves…

  16. Two infinite parallel hyper-planes in Rd. Implementing the SST (Cont.) The efficient solution is found by using a different representation for the complementary cones.

  17. Implementing the SST - Main Lemma Lemma: contains a vector other than {0}, if and only if the intersection of the unit hyper-sphere Sd-1 with the regions between the bounding hyper-planes is non-empty. Implication: In order to test for a single solution, we do not have to intersect d-dimensional cones.

  18. Implementing the SST (Cont.) Two options: The entire intersection is inside the unit sphere. Then SST holds. Otherwise, SST failed. Hence, only test if all the vertices of the intersection are inside the unit sphere.

  19. Implementing the SST (Cont.) SST holds SST failed

  20. Implementing the SST (Cont.) The set of 2d hyper-plane intersection (HPI) equations to solve is: Or in matrix (d d) form:

  21. Implementing the SST (Cont.) Due to symmetry, only 2d-1hyper-plane intersection (HPI) tests are required. 4HPI tests for d = 3 (x, y, z) constraints in R3.

  22. Implementing the SST (Cont.) An efficient way to solve the 2d-1 systems of equations is due to the inherent symmetry of the problem: Therefore, given one solution (e.g., {-,-,…,-,-}), computing a second solution that differs by a single sign (e.g., {-,-,…,-,+}), takes only O(d) operations and not O(d2 ).

  23. Implementing the SST – Example Without SST With SST

  24. Implementing the SST – Example (Zoom) Without SST With SST

  25. Purging Away Zero-solution Domains Problem: SST does not guarantee that a solution exists in the sub-domain. We might start the numerical iterations only to find that no root exists in this sub-domain.

  26. Purging Away Zero-solution Domains (Cont.) Therefore, we seek to purge away, as much as possible, sub-domains that contain no solution. We present a second criterion, in addition to the CH criterion, for identifying no-solution sub-domains.

  27. Purging Away Zero-solution Domains (Cont.) The basic idea: Bound the function Fi=0 by a pair of parallel hyper-planes in Rd. If the intersection of the hyper-planes is entirely outside the sub-domain, then the solution is outside the sub-domain. F1 F2

  28. Purging Away Zero-solution Domains (Cont.) Problem: In order to bound Fi=0 inRd we need to have at least a sample point uonFi=0. Solution: Bound FiinRd+1 and intersect the resulting (d+1)-dimensional hyper-planes with ud+1=0, resulting bounding hyper-planes in Rd.

  29. Purging Away Zero-solution Domains (Cont.) • Promote Fi from a scalar function to a hyper-surface in Rd+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.

  30. Purging Away Zero-solution Domains Illustration (d=1) [ ] ud+1 = 0 Center for Graphics and Geometric Computing, Technion

  31. Purging Away Zero-solution Domains (Cont.) How to determine whether the intersection is entirely outside the sub-domain? Possible solutions: Linear programming techniques. Similarly to the SST for finding all vertices of the intersection.

  32. Purging Away Zero-solution sub-domains – Example With hyper-plane test Without hyper-plane test

  33. Future Work • Compare the SST with algebraic/numeric termination criteria (Kantorovich). • Use the vertices, which were computed for purging away zero-solutions, for clipping the sub-domain. • Extend the presented ideas for under- and over-constrained systems. • Search for tighter bounding volumes on the solutions (not just cones).

  34. End

More Related