Download
1 / 45
450 likes | 562 Views
Computing the Genus of a Curve Numerically. Andrew Sommese University of Notre Dame In collaboration with Daniel Bates, IMA Christopher Peterson, Colorado State University Charles Wampler, General Motors R&D Center. Reference on Numerical Algebraic Geometry up to 2005:
E N D
Computing the Genus of a Curve Numerically Andrew Sommese University of Notre Dame In collaboration with Daniel Bates, IMA Christopher Peterson, Colorado State University Charles Wampler, General Motors R&D Center
Reference on Numerical Algebraic Geometry up to 2005: • A.J. Sommese and C.W. Wampler, Numerical solution of systems of polynomials arising in engineering and science, (2005), World Scientific Press. • Website with more information and some articles: • www.nd.edu/~sommese
Overview • What is an algebraic curve? • What is the genus of a curve? • An example from applications • Four-bar-planar linkage coupler curves • Background on Numerical Algebraic Geometry • How to represent Positive Dimensional Solution Sets • Numerical approach by using Hurwitz’s formula
What is an algebraic curve? • Curves are ancient, beautiful, and useful. • They are pervasive in mathematics: in number theory, algebraic, analytic, and differential geometry, complex analysis, numerical analysis, topology,... • The arise everywhere in applications.
What is the genus of a curve? • Let’s start with x2 + y2 = 1:
Interlude between slides • Including blackboard discussion of • the examples; • the Euler characteristic; and • Hurwitz’s formula.
Some examples from applications • One dimensional component of a polynomial system on :
Planar four-bar coupler curves • A four-bar planar linkage is a planar quadrilateral with a rotational joint at each vertex. • They are useful for converting one type of motion to another. • They occur everywhere.
P0 δ P y x C D D′ µ λ u v b C′
P0 δ P θ y yeiθ C b b-δ µ v v = y – b veiμ = yeiθ- (b - δ) = yeiθ+ δ- b C′
We use complex numbers (as is standard in this area) Summing over vectors we have two equations plus their two conjugates
This gives four equations: in the variables
Multiplying each side by its complex conjugate and letting we get three equations in the four variables
We can solve for in the first two equations using Cramer’s rule, and substitute into the last Equation. This gives us an equation of degree (3,3) in We replace by
Background on Numerical Algebraic Geometry • Find all isolated solutions in of a system on n polynomials:
Solving a system • Homotopy continuation is our main tool: • Start with known solutions of a known start system and then track those solutions as we deform the start system into the system that we wish to solve.
Path Tracking This method takes a system g(x) = 0, whose solutions we know, and makes use of a homotopy, e.g., Hopefully, H(x,t) defines “nice paths” x(t) as t runs from 1 to 0. They start at known solutions of g(x) = 0 and end at the solutions of f(x) at t = 0.
The paths satisfy the Davidenko equation • To compute the paths: use ODE methods to predict and Newton’s method to correct.
x1(t) • x2(t) • x3(t) • x4(t) • Known solutions of g(x)=0 • Solutions of • f(x)=0 • t=0 • H(x,t) = (1-t) f(x) + t g(x) • t=1
x* { • xj(t) • 0 • Dt • prediction • Newton correction 1
Continuation’s Core Computation • Given a system f(x) = 0 of n polynomials in n unknowns, continuation computes a finite set S of solutions such that: • any isolated root of f(x) = 0 is contained in S; • any isolated root “occurs” a number of times equal to its multiplicity as a solution of f(x) = 0; • S is often larger than the set of isolated solutions.
Positive Dimensional Solution Sets We now turn to finding the positive dimensional solution sets of a system
How to represent positive dimensional components? • S. + Wampler in ’95: • Use the intersection of a component with generic linear space of complementary dimension.
Use a generic flag of affine linear spaces to get witness point supersets This approach has 19th century roots in algebraic geometry, e.g., adjunction theory and the use of hyperplane sections to study varieties.
Numerical approach by using Hurwitz’s formula • Given an irreducible curve X that is a component of V(f) of a polynomial system, • X is represented by a generic hyperplane L = 0 plus W, the set of degree X points where L meets X. • L = 0 is the fiber of a general projection p from CN to C.
g = -deg X + 1 + r/2 • Steps: • Find a finite set B of points of C containing all singular points of X plus all branch points of X. • Do monodromy of the fiber of p around sufficiently small loops around the images R of B under p. • Compute the number of sets gr which the monodromy breaks the fiber over r in R into. • The sum of the (deg X - gr) as r runs over R equals r for the induced map from the desingularization of X down to C.
Fermat curves: xd + yd = 1 • we use a general linear projection. So we have • a degree d projection; • d(d-1) branchpoints in different fibers all contributing 1 to r: the genus g is d(d-3)/2 + 1. • So g = 6 for d = 5; g = 10 for d = 6; and g = 15 for d = 7.
genus 0: 10 contributions of one each tor zero contributions 7 double points 1 sextuple point g = -6 + 10/2 + 1 = 0
A general planar four-bar coupler curve C • C is a (3,3) curve in P1 x P1. Arithmetic genus at most 4 with at least 3 singularities: so geometric genus is at most 1. • We treat it as a degree 6 curve in P2. We compute • 30 potential branch points counting multiplicities; • 2 potential branch points with multiplicity 6 each make a local contribution 0 to r; • 3 potential branch points with multiplicity 2 each make a contribution of 0 to r; • 12 potential branch points contribute one each to r.
Summary • The numerical approach gives a direct method of computing the geometric genus of one-dimensional irreducible components of the solutions set of a system of polynomials.
