Geometric Constraint Solving in 2D

Similar documents
Strong Matching of Points with Geometric Shapes

GENERICITY OF INFINITE-ORDER ELEMENTS IN HYPERBOLIC GROUPS

Find the equation of a plane perpendicular to the line x = 2t + 1, y = 3t + 4, z = t 1 and passing through the point (2, 1, 3).

jf 00 (x)j ds (x) = [1 + (f 0 (x)) 2 ] 3=2 (t) = jjr0 (t) r 00 (t)jj jjr 0 (t)jj 3

MATH 6210: SOLUTIONS TO PROBLEM SET #3

Approximating min-max k-clustering

AP Calculus Testbank (Chapter 10) (Mr. Surowski)

GOOD MODELS FOR CUBIC SURFACES. 1. Introduction

Lecture 1.2 Pose in 2D and 3D. Thomas Opsahl

Topic 30 Notes Jeremy Orloff

Central Force Motion Challenge Problems

7. The olynomial P (x) =x, 9x 4 +9x, 4x +4x, 4 can be written in the form P (x) =(x, ) n Q(x), where n is a ositive integer and Q(x) is not divisible

Finite-State Verification or Model Checking. Finite State Verification (FSV) or Model Checking

Model checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]

A Note on Guaranteed Sparse Recovery via l 1 -Minimization

Multiplicative group law on the folium of Descartes

Uniform Law on the Unit Sphere of a Banach Space

A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands

Combinatorics of topmost discs of multi-peg Tower of Hanoi problem

CMSC 425: Lecture 4 Geometry and Geometric Programming

Introduction to Arithmetic Geometry Fall 2013 Lecture #10 10/8/2013

THE EXTREMUM OF A FUNCTION DEFINED ON THE EUCLIDEAN PLANE

Additive results for the generalized Drazin inverse in a Banach algebra

Matching points with geometric objects: Combinatorial results

HARMONIC EXTENSION ON NETWORKS

Circles. Exercise 9.1

Lecture 6. 2 Recurrence/transience, harmonic functions and martingales

Mollifiers and its applications in L p (Ω) space

19th Bay Area Mathematical Olympiad. Problems and Solutions. February 28, 2017

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

SECTION 12: HOMOTOPY EXTENSION AND LIFTING PROPERTY

Elementary theory of L p spaces

Solution sheet ξi ξ < ξ i+1 0 otherwise ξ ξ i N i,p 1 (ξ) + where 0 0

CMSC 425: Lecture 7 Geometric Programming: Sample Solutions

(Workshop on Harmonic Analysis on symmetric spaces I.S.I. Bangalore : 9th July 2004) B.Sury

Solvability and Number of Roots of Bi-Quadratic Equations over p adic Fields

Mathematics 2260H Geometry I: Euclidean geometry Trent University, Fall 2016 Solutions to the Quizzes

ISOSCELES TRIANGLES IN Q 3. Matt Noble Department of Mathematics, Middle Georgia State University, Macon, Georgia

Anisotropic elastoplasticity for cloth, knit and hair frictional contact supplementary technical document

The Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule

On Wrapping of Exponentiated Inverted Weibull Distribution

arxiv:math/ v4 [math.gn] 25 Nov 2006

E-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula

Extremal Polynomials with Varying Measures

0.1 Practical Guide - Surface Integrals. C (0,0,c) A (0,b,0) A (a,0,0)

Lecture 7: Introduction to syntax-based MT

MATH 361: NUMBER THEORY ELEVENTH LECTURE

Analysis of some entrance probabilities for killed birth-death processes

A Model for Randomly Correlated Deposition

Supplementary Materials for Robust Estimation of the False Discovery Rate

Principles of Computed Tomography (CT)

On the continuity property of L p balls and an application

FFTs in Graphics and Vision. Rotational and Reflective Symmetry Detection

Simplifications to Conservation Equations

Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems

216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are

CONGRUENCE PROPERTIES MODULO 5 AND 7 FOR THE POD FUNCTION

GRACEFUL NUMBERS. KIRAN R. BHUTANI and ALEXANDER B. LEVIN. Received 14 May 2001

PARTITIONS AND (2k + 1) CORES. 1. Introduction

Solutions of Benelux Mathematical Olympiad 2010

Improvement on the Decay of Crossing Numbers

UPPAAL tutorial What s inside UPPAAL The UPPAAL input languages

Mathematics. Class 12th. CBSE Examination Paper 2015 (All India Set) (Detailed Solutions)

Opposite-quadrant depth in the plane

Dependence on Initial Conditions of Attainable Sets of Control Systems with p-integrable Controls

Proof Nets and Boolean Circuits

Optimal Approximate Polytope Membership

LECTURE 1: ADVERTISEMENT LECTURE. (0) x i + 1 2! x i x j

( 1 ) Show that P ( a, b + c ), Q ( b, c + a ) and R ( c, a + b ) are collinear.

Robustness of classifiers to uniform l p and Gaussian noise Supplementary material

#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS

Complex Analysis Homework 1

Topic 2 [312 marks] The rectangle ABCD is inscribed in a circle. Sides [AD] and [AB] have lengths

Recognise the Equation of a Circle. Solve Problems about Circles Centred at O. Co-Ordinate Geometry of the Circle - Outcomes

SUBORBITAL GRAPHS FOR A SPECIAL SUBGROUP OF THE NORMALIZER OF. 2p, p is a prime and p 1 mod4

Gaps in Semigroups. Université Pierre et Marie Curie, Paris 6, Equipe Combinatoire - Case 189, 4 Place Jussieu Paris Cedex 05, France.

Additional Mathematics Lines and circles

SECTION 5: FIBRATIONS AND HOMOTOPY FIBERS

THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN

On the minimax inequality for a special class of functionals

Numerical Linear Algebra

2.6 Primitive equations and vertical coordinates

Pretest (Optional) Use as an additional pacing tool to guide instruction. August 21

MAZUR S CONSTRUCTION OF THE KUBOTA LEPOLDT p-adic L-FUNCTION. n s = p. (1 p s ) 1.

Pythagorean triples and sums of squares

Class IX Chapter 5 Introduction to Euclid's Geometry Maths

d) Find the equation of the circle whose extremities of a diameter are (1,2) and (4,5).

7. Two Random Variables

Answers Investigation 2

Computational Solid Mechanics Computational Plasticity

INJECTIVITY RADIUS AND DIAMETER OF THE MANIFOLDS OF FLAGS IN THE PROJECTIVE PLANES

On the relationship between sound intensity and wave impedance

CTL, the branching-time temporal logic

CHAPTER 2: SMOOTH MAPS. 1. Introduction In this chapter we introduce smooth maps between manifolds, and some important

Inference for Empirical Wasserstein Distances on Finite Spaces: Supplementary Material

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

CO-ORDINATE GEOMETRY. 1. Find the points on the y axis whose distances from the points (6, 7) and (4,-3) are in the. ratio 1:2.

On the Field of a Stationary Charged Spherical Source

2. (i) Find the equation of the circle which passes through ( 7, 1) and has centre ( 4, 3).

Convex Analysis and Economic Theory Winter 2018

Transcription:

Universitat Politècnica de Catalunya Deartament de Llenguatges i Sistemes Informàtics PhD Thesis Geometric Constraint Solving in 2D November, 1998 Antoni Soto i Riera Advisor: Robert Joan i Arinyo

Contents 1. Motivation. 2. Geometric Constraint Problem (GCP). 3. Aroaches to Geometric Constraint Solving. 4. Constructive technique. 5. Hybrid technique. 6. Conclusions and future work. 1

Motivation d7 d2 a4 t5 d6 rd5 d5 t4 Feature-based CAD systems + Editable reresentation (Ere) 2-dimensional constraint-based editor 2

Geometric Constraint Problem y C x D L BC L AC h α A L AB B d y = x v v = 0.5 A geometric constraint roblem (GCP) consists of A set of geometric elements {A,, L AB, }. A set of values {d, α, h} A set of dimensional variables {x, y}. A set of external variables {v}. A set of valuated and symbolic constraints. A set of equations. 3

Reresentation: Geometric Constraint Grah The geometric constraints of a GCP can be reresented by a constraint grah G = (E, V ). y C D L AC x h L BC α L AB A d B The vertices in V are two-dimensional geometric elements with two degrees of freedom. The edges in E are constraints that reduce by one the degrees of freedom. 4

Reresentation: First-order Logic y C x D L BC L AC h α A L AB B d y = x v v = 0.5 A geometric constraint roblem can be reresented by a formula in first-order logic. ϕ(a, B, C, D, L AB, L AC, L BC, x, y, v) d(a, B) = d on(a, L AB ) on(b, L AB ) on(a, L AC ) on(c, L AC ) on(d, L AC ) on(b, L BC ) on(c, L BC ) h(c, L AB ) = h a(l AB, L BC ) = α d(a, C) = x d(c, D) = y y = x v v = 0.5 5

Geometric Constraint Solving Geometric constraint solving (GCS) consists in roving the truth of the formula A B C D L AB L AC L BC x y v ϕ(a, B, C, D, L AB, L AC, L BC, x, y, v) by finding the osition of the geometric elements and the values of tags and external variables that satisfy the constraints. 6

Over-constrained Geometric Constraint Problem Theorem 1 (Laman, 1970) Let G = (P, D) a geometric constraint grah where the vertices in P are oints in the two dimensional Euclidean sace and the edges D P P are distance constraints. G is generically well constrained if and only if for all G = (P, D ), subgrah of G induced by P P, 1. D 2 P 3, and 2. D = 2 P 3. r r r 7

Structurally over-constrained Geometric Constraint Problem Definition 1 A geometric constraint grah is structurally over-constrained if and only if exists an induced subgrah with n vertices and m edges such that m > 2 n 3. Definition 2 A geometric constraint grah is structurally well-constrained if and only if it is not structurally over-constrained and E = 2 V 3. Definition 3 A geometric constraint grah is structurally under-constrained if and only if it is not structurally over-constrained and E 2 V 3. 8

Aroaches to Geometric Constraint Solving Solving systems of equations Numerical Constraint Solvers Symbolic Constraint Solvers Proagation Methods Structural analysis Constructive Constraint Solvers Grah based Rule based Degrees of freedom analysis Geometric theorem roving 9

Constructive technique: Ruler-and-comass constructibility A oint P is constructible if there exists a finite sequence P 0, P 1,..., P n = P of oints in the lane with the following roerty. Let S j = {P 0, P 1,..., P j }, for 1 j n. For each 2 j n is either 1. the intersection of two distinct straight lines, each joining two oints of S j 1, or 2. a oint of intersection of a straight line joining two oints of S j 1 and a circle with centre a oint of S j 1 and radius the distance between two oints of S j 1, or 3. a oint of intersection of two distinct circles, each with centre a oint of S j 1 and radius the distance between two oints of S j 1. 10

Constructive technique: Constraints sets A CA set is a air of oriented segments which are mutually constrained in angle. 4 3 1 2 A CD set is a set of oints with mutually constrained distances. 1 3 5 2 4 A CH set is a oint and a segment constrained by the erendicular distance from the oint to the segment. 1 2 3 11

Constructive technique: Geometric locus cd 2 1 cd 1 CR ca 1 ca 1 1 2 1 2 3 cd 1 CC cd 1 RA ch 1 ch 1 cd 1 1 2 cd 1 1 2 RP RT 12

cd 1 Constructive technique: Set of rules cd 2 cd 3 1 2 cd 1 CR-CR cd 2 ca 1 cd 2 ca 1 1 1 cd 1 CR-CC CR-RA ch 1 cd 2 cd 2 ch 1 1 1 cd 1 cd 1 CR-RP CR-RT 13

Constructive technique: Analysis and Construction hases ϕ( 1, 2, 3 ) d( 1, 2 ) = d 1 d( 2, 3 ) = d 2 d( 3, 1 ) = d 3 Geometric Constraint Problem 3 analysis d 3 d 2 1 d 1 2 ψ( 1, 2, 3 ) 1 = (0, 0) 2 = (d 1, 0) 3 = inter(circle( 1, d 3 ), circle( 2, d 2 )) Constructive Formula construction υ d 1 = 200 d 2 = 180 d 3 = 220 Values ψ υ ( 1, 2, 3 ) 1 = (0, 0) 2 = (200, 0) 3 = (140, 169.706) Valuated Formula 14

Constructive technique: Structure of the Analysis hase ϕ( 1, 2, 3 ) d( 1, 2 ) = d 1 d( 2, 3 ) = d 2 d( 3, 1 ) = d 3 translation C2(cd 1, 1, 2, d 1 ) C2(cd 2, 2, 3, d 2 ) C2(cd 3, 1, 3, d 3 ) E 0 cd 1 = { 1, 2 } cd 2 = { 2, 3 } cd 3 = { 1, 3 } analysis DDD(cd 4, 1, 2, 3, cd 1, cd 2, cd 3 ) E f cd 4 = { 1, 2, 3 } generation ψ( 1, 2, 3 ) 1 = (0, 0) 2 = (d 1, 0) 3 = inter(circle( 1, d 3 ), circle( 2, d 2 )) 15

Constructive technique: Correctness (1) Let R be a set of tules (D, X) where D is a set of CD sets, and X is a set of CA sets and CH sets. Let ρ = DDD DDX DXX be a reduction relation. We define the abstract reduction system R = R, ρ. We roof termination and confluence that imlies canonicity and unique normal form roerty. 16

Constructive technique: Correctness (2) d 4 q s r s d 3 d 2 q a l h m r DDD q DAH s q r d 1 DAH s r DDD d 5 d 4 17

Constructive technique: Reduction rules (1) (D, X) DDD ((D {d 1, d 2, d 3 }) {d 1 d 2 d 3 }, X) if {d 1, d 2, d 3 } D d 1 d 2 = {1} d 2 d 3 = {2} d 1 d 3 = {3} 1 2 3 18

Constructive technique: Reduction rules (2) (D, X) DDX ((D {d 1, d 2 }) {d 1 d 2 unts(x 1 )}, X {x 1 }) {d 1, d 2 } D d 1 d 2 = { 1 } d if 2 x 1 X unts(x 1 ) d 1 = {} 1 19

Constructive technique: Reduction rules (3) (D, X) DXX ((D {d 1 }) {d 1 unts(x 1 ) unts(x 2 )}, X {x 1, x 2 }) d 1 D {x if 1, x 2 } X unts(x 1 ) d 1 = {} unts(x 2 ) d 1 = {} 20

Hybrid technique: Introduction Goal Solve symbolic constraints keeing the two hases of the constructive technique for valuated constraints. Idea Federate a constructive solver and an equational solver. Required A technique of analysis of systems of equations. 21

Hybrid technique: Analysis of systems of equations 1. Reresent the structure of the systems of equations by a biartite grah (bigrah). x = 5 x = 2y x = y 2y = z yt = u 3t = uu = w v x y z t u v w 2. Comute the Dulmage-Mendelsohn decomosition of the bigrah. V 0 is the over-determined art, V is the underdetermined art. V 1,, V n are the consistent art. x = 5 x = 2y x = y 2y = z yt = u 3t = uu = w v x y z t u v w V 0 V 1 V2 V 22

Hybrid technique: Motivation (1) A geometric constraint roblem that can not be solved without considering geometric variables. d 1 b a o 1, t 1 l d 2 n α 1 g α 2 m d o 2, t 2 c ϕ(a, b, c, d, g, l, m, n, r) d(a, b) = d 1 d(b, c) = d 2 on(a, l) on(b, l) on(b, n) on(c, n) on(c, m) on(d, m) a(l, n) = α 1 a(l, m) = α 2 d(g, c) = r h(g, l) = r h(g, m) = r d(g, d) = r 23

Hybrid technique: Motivation (2) Final state of the geometric analyzer. C 1 d 1 b a o 1, t 1 l d 2 n α 1 g α 2 m d o 2, t 2 c Final state of the equational analyzer. d(g, a) = r h(g, l) = r h(g, m) = r d(g, d) = r r 24

Hybrid technique: Technique (1) 1. Reresent geometric variables in the bigrah (B). d(g, a) = r h(g, l) = r h(g, m) = r d(g, d) = r g x a x a y g y l x l y r n x n y d x d x 2. Comute R(B, C 1 ), the restriction of bigrah B by CD set C 1. The equations are analyzed with resect to a coordinate system (CD set). d(g, a) = r h(g, l) = r h(g, m) = r d(g, d) = r g x g y r d x d x V 1 V 25

Hybrid technique: Technique (2) 3. For each solved dimensional variable, add a new constraint set to the state of the geometric analyzer. 4. For each air of solved geometric variables (v x, v y ), add the geometric element v to the rojection CD set C 1. 5. Remove solved variables and equations from the bigrah B. C 1 d 1 b a o 1, t 1 l d 2 n α 1 g α 2 m d o 2, t 2 c C 2 26

Hybrid technique: Correctness (1) Let S be a set of tules (D, X, B) where D is a set of CD sets, X is a set of CA sets and CH sets, and B is a bigraf reresenting symbolic geometric constraints and equations. Let ρ relation. be the constructive reduction Let κ be the equational analysis reduction relation. We define the abstract reduction system S = S, ρ κ. We roof termination and confluence that imlies canonicity and unique normal form roerty. 27

Conclusions A correct ruler-and-comass constructive method. A clean hase structure. A correct hybrid method combining a constructive method and an equation analysis method. A rototye imlementation. 28

Future work Study the domain of constructive methods. Extending the domain of our constructive method. Selection of the solution. Determine the range of values of a constraint. 29