- By
**mahon** - Follow User

- 177 Views
- Uploaded on

Download Presentation
## PowerPoint Slideshow about 'Introduction to Gröbner Bases for Geometric Modeling' - mahon

**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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

### Introduction to Gröbner Bases for Geometric Modeling

Geometric & Solid Modeling

1989

Christoph M. Hoffmann

Algebraic Geometry

- Branch of mathematics.
- Express geometric facts in algebraic terms in order to interpret algebraic theorems geometrically.
- Computations for geometric objects using symbolic manipulation.
- Surface intersection, finding singularities, and more…
- Historically, methods have been computationally intensive, so they have been used with discretion.

source: Hoffmann

Goal

- Operate on geometric object(s) by solving systems of algebraic equations.
- “Ideal”: (informal partial definition) Set of polynomials describing a geometric object symbolically.
- Considering algebraic combinations of algebraic equations (without changing solution) can facilitate solution.
- Ideal is the set of algebraic combinations (to be defined more rigorously later).
- Gröbner basis of an ideal: special set of polynomials defining the ideal.
- Many algorithmic problems can be solved easily with this basis.
- One example (focus of our lecture): abstract ideal membership problem:
- Is a given polynomial g in a given ideal I ?
- Equivalently: can g be expressed as an algebraic combination of the fjfor some polynomials hj?
- Answer this using Gröbner basis of the ideal.
- Rough geometric interpretation: g can be expressed this way when surface g = 0 contains all points that are common intersection of surfaces fj = 0.

source: Hoffmann

Overview

- Algebraic Concepts
- Fields, rings, polynomials
- Field extension
- Multivariate polynomials and ideals
- Algebraic sets and varieties
- Gröbner Bases
- Lexicographic term ordering and leading terms
- Rewriting and normal-form algorithms
- Membership test for ideals
- Buchberger’s theorem and construction of Gröbner bases
- For discussion of geometric modeling applications of Gröbner bases, see Hoffmann’s book.
- e.g. Solving simultaneous algebraic expressions to find:
- surface intersections
- singularities

source: Hoffmann

Algebraic Concepts:Fields, Rings, and Polynomials

- Consider single algebraic equation:
- Values of xi’s are from a field.(Recall from earlier in semester.)
- Elements can be added, subtracted, multiplied, divided*.
- Ground fieldk is the choice of field .
- Univariate polynomial over k is of form:
- Coefficients are numbers in k.
- k[x] = all univariate polynomials using x’s.
- It is a ring (recall from earlier in semester): addition, subtraction, multiplication, but not necessarily division.
- Can a given polynomial be factored?
- Depends on ground field
- e.g. x2+1 factors over complex numbers but not real numbers.
- Reducible: polynomial can be factored over ground field.
- Irreducible: polynomial cannot be factored over ground field.

* for non-0 elements

source: Hoffmann

Algebraic Concepts:Field Extension

- Field extension: enlarging a field by adjoining (adding) new element(s) to it.
- Algebraic Extension:
- Adjoin an element u that is a root of a polynomial (of degree m) in k[x].
- Resulting elements in extended field k(u)are of form:
- e.g. extending real numbers to complex numbers by adjoining i
- i is root of x2+1, so m=2 and extended field elements are of form a + bi
- e.g. extending rational numbers to algebraic numbers by adjoining roots of all univariate polynomials (with rational coefficients)
- Transcendental Extension:
- Adjoin an element (such as p) that is not the root of any polynomial in k[x].

source: Hoffmann

Algebraic Concepts:Multivariate Polynomials

- Multivariate polynomial over k is of form:
- Coefficients are numbers in k.
- Exponents are nonnegative integers.
- k[x1,…,xn] = all multivariate polynomials using x’s.
- It is a ring: addition, subtraction, multiplication, but not necessarily division.
- Can a given polynomial be factored?
- Depends on ground field (as in univariate case)
- Reducible: polynomial can be factored over ground field.
- Irreducible: polynomial cannot be factored over ground field.
- Absolutely Irreducible: polynomial cannot be factored over any ground field.
- e.g.

source: Hoffmann

Algebraic Concepts:Ideals

- For ground field k, let:
- kn be the n-dimensional affine space over k.
- mathematical physicist John Baez: "An affine space is a vector space that's forgotten its origin”.
- Points in kn are n-tuples (x1,…,xn), with xi’s having values in k.
- f be an irreducible multivariate polynomial in k[x1,…,xn]
- gbe a multivariate polynomial in k[x1,…,xn]
- f = 0 be the hypersurface in kn defined by f
- Since hypersurface gf = 0 includes f = 0, view f as intersection of all surfaces of form gf = 0
- is an ideal*
- g varies over k[x1,…,xn]
- Consider the ideal as the description of the surface f.
- Ideal is closed under addition and subtraction.
- Product of an element of k[x1,…,xn]with a polynomial in the ideal is in the ideal.

source: Hoffmann and others

*Ideals are defined more generally in algebra.

Algebraic Concepts:Ideals (continued)

- Let F be a finite set of polynomials f1, f2,…, fr in k[x1,…,xn]
- Algebraic combinations of the fi form an ideal generated by F (a generating set*):
- generators: { f, g }
- Goal: find generating sets, with special properties, that are useful for solving geometric problems.

* Not necessarily unique.

source: Hoffmann

Algebraic Concepts:Algebraic Sets

- Let be the ideal generated by the finite set of polynomials F = { f1, f2,…, fr }.
- Let p = (a1,…, an) be a point in kn such that g(p) = 0 for every g in I.
- Set of all such points p is the algebraic setV(I) of I.
- It is sufficient that fi(p) = 0 for every generator fi in F.
- In 3D, the algebraic surface f = 0 is the algebraic set of the ideal .

source: Hoffmann

Algebraic Concepts:Algebraic Sets (cont.)

- Intersection of two algebraic surfaces f, g in 3D is an algebraic space curve.
- The curve is the algebraic set of the ideal.
- But, not every algebraic space curve can be defined as the intersection of 2 surfaces.
- Example where 3 are needed*: twisted cubic (in parametric form):
- Can define twisted cubic using 3 surfaces: paraboloid with two cubic surfaces:
- Motivation for considering ideals with generating sets containing > 2 polynomials.

*see Hoffman’s Section 7.2.6 for subtleties related to this statement.

source: Hoffmann

Algebraic Concepts:Algebraic Sets and Varieties (cont.)

- Given generatorsF = { f1, f2,…, fr }, the algebraic set defined by F in kn has dimension n-r
- If equations fi = 0 are algebraically independent.
- Complication: some of ideal’s components may have different dimensions.

source: Hoffmann

Algebraic Concepts:Algebraic Sets and Varieties (cont.)

- Consider algebraic set V(I) for ideal I in kn.
- V(I) is reducible when V(I) is union of > 2 point sets, each defined separately by an ideal.
- Analogous to polynomial factorization:
- Multivariate polynomial f that factors describes surface consisting of several components
- Each component is an irreducible factor of f.
- V(I) is irreducible implies V(I) is a variety.

source: Hoffmann

Algebraic Concepts:Algebraic Sets and Varieties (cont.)

- Example: Intersection curve of 2 cylinders:
- Intersection lies in 2 planes:

and

- Irreducible ellipse in plane is is algebraic set in ideal generated by { f1,g1 }.
- Irreducible ellipse in plane is is algebraic set in ideal generated by { f1,g2 }.
- Ideal is reducible.
- Decomposes into and
- Algebraic set
- Varieties: V(I2) and V(I3)

source: Hoffmann

Algebraic Concepts:Algebraic Sets and Varieties (cont.)

- Example: Intersection curve of 2 cylinders:
- Intersection curve is not reducible
- These 2 component curves cannot be defined separately by polynomials.
- Rationale: Bezout’s Theorem implies intersection curve has degree 4. Furthermore:
- Union of 2 curves of degree m and n is a reducible curve of degree m + n.
- If intersection curve were reducible, then consider degree combinations for component curves (total = 4):
- 1 + 3: illegal since neither has degree 1.
- 2 + 2: illegal since neither is planar.
- Conclusion: intersection curve irreducible.
- Bezout’s Theorem also implies that twisted cubic cannot be defined algebraically as intersection of 2 surfaces:
- Twisted cubic has degree 3.
- Bezout’s Theorem would imply it is intersection of plane and cubic surface.
- But twisted cubic is not planar; hence contradiction.

Bezout’s Theorem*: 2 irreducible surfaces of degree m and n intersect in a curve of degree mn. *allowing complex coordinates, points at infinity

source: Hoffmann

Gröbner Bases:Formulating Ideal Membership Problem

- Can help to solve geometric modeling problems such as intersection of implicit surfaces (see Hoffmann Sections 7.4-7.8).
- Here we only treat the ideal membership problem for illustrative purposes:
- “Given a finite set of polynomials F = { f1, f2,…, fr }, and a polynomial g, decide whether g is in the ideal generated by F; that is, whether g can be written in the form where the hi are polynomials.”
- Strategy: rewrite g until original question can be easily answered.

source: Hoffmann

Gröbner Bases:Lexicographic Term Ordering and Leading Terms

- Need to judge if “this polynomial is simpler than that one.”
- Power Product:
- Lexicographic ordering of power products:
- x
- If then for all power products w.
- If u and v are not yet ordered by rules 1 and 2, then order them lexicographically as strings.

Example for n=2 on board...

source: Hoffmann

Gröbner Bases:Lexicographic Term Ordering and Leading Terms

- Each term in a polynomial g is a coefficient combined with a power product.
- Leading term lt(g) of g: term whose power product is largest with respect to ordering
- lcf (g) =leading coefficient of lt(g)
- lpp (g) =leading power product of lt(g)
- Definition: Polynomial f is simpler than polynomial g if:

Example 7.1 on board...

source: Hoffmann

Gröbner Bases:Rewriting and Normal-Form Algorithms

- Given polynomial g and set of polynomials F = { f1, f2,…, fr }
- Rewrite/simplify g using polynomials in F.
- gis in normal formNF(g, F) if it cannot be reduced further. Note: normal form need not be unique.

source: Hoffmann

Example 7.2 on board...

Gröbner Bases:Rewriting and Normal-Form Algorithms

- If normal form from rewriting algorithm is unique
- then g is in ideal when NF(g, F) = 0.
- This motivates search for generating sets that produce unique normal forms.

source: Hoffmann

Gröbner Bases:A Membership Test for Ideals

- Goal: Rewrite g to decide whether g is in the ideal generated by F.
- Gröbner basisG of ideal
- Set of polynomials generating F.
- Rewriting algorithm using G produces unique normal forms.
- Ideal membership algorithm usingG:

source: Hoffmann

Example 7.3 on board...

Gröbner Bases:Buchberger’s Theorem & Construction

- Algorithm will consist of 2 operations:
- Consider a polynomial, and bring it into normal form with respect to some set of generators G.
- From certain generator pairs, compute S-polynomials (see definition on next slide) and add their normal forms to the generator set.
- G starts as input set F of polynomials
- G is transformed into a Gröbner basis.
- Some Implementation Issues:
- Coefficient arithmetic must be exact.
- Rational arithmetic can be used.
- Size of generator set can be large.
- Reduced Gröbner bases can be developed.

source: Hoffmann

Gröbner Bases:Buchberger’s Theorem & Construction (continued)

Example 7.4 on board...

source: Hoffmann

Gröbner Bases:Buchberger’s Theorem & Construction (continued)

Buchberger’s Theorem:

foundation of algorithm

Gröbner basis construction algorithm

Example 7.5 on board...

source: Hoffmann

Download Presentation

Connecting to Server..