Computations in Algebraic Geometry with Macaulay 2 PDF
Preview Computations in Algebraic Geometry with Macaulay 2
Algorithms and Computation in Mathematics • Volume 8 Editors E. Becker M. Bronstein H. Cohen D. Eisenbud R. Gilman Springer-Verlag Berlin Heidelberg GmbH David Eisenbud Daniel R. Grayson Michael Stillman Bernd Sturmfels (Eds.) Computations in Algebraic Geometry with Macaulay 2 , Springer Editors David Eisenbud Daniel R. Grayson Mathematical Sciences University of Illinois Research Institute at Urbana-Champaign 1000 Centennial Drive Department of Mathematics Berkeley, CA 94720, USA Urbana, IL 61801, USA e-mail: [email protected] e-mail: [email protected] Michael Stillman Bernd Sturmfels Cornell University University of California Department of Mathematics Department of Mathematics Ithaca, NY 14853, USA Berkeley, CA 94720, USA e-mail: [email protected] e-mail: [email protected] Mathematics Subject Classification (2000): 13-04, 13P, 14-04, 14Q, 16Z05, 68W30, nY, 12Y05, 14P, 12DI0, 14M15, 14N15 Library of Congress Cataloging-in-Publication Data applied for Die Deutsche Bibliothek -CIP-Einheitsaufnahme Computations in algebraic geometry with macaulay 21 David Eisenbud ... (ed.). -Berlin; Heidelberg; New York; Barcelona; Hong Kong; London; Milan; Paris; Tokyo: Springer, 2002 (Algorithms and computation in mathematics: Vol. 8) ISBN 978-3-642-07592-6 ISBN 978-3-662-04851-1 (eBook) DOl 10.1007/978-3-662-04851-1 This work is subject to copyright. All rights are reserved. whether the whole or part of the material is concerned. specifically the rights of translation. reprinting. reuse of illustrations. recitation. broadcasting. reproduction on microfilm or in any other way. and storage in data banks. Duplication of this publication or parts thereofis permitted only under the provisions of the German Copyright Law of September 9. 1965. in its current version. and permission for use must always be obtained from Springer-Verlag. Violations are liable for prosecution under the German Copyright Law. @Springer-Verlag Berlin Heidelberg 2002 Originally published by Springer-Verlag Berlin Heidelberg New York in 2002. Softcover reprint of the hardcover 1st edition 2002 The use of general descriptive names, registered names. trademarks. etc. in this publication does not intply. even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations and therefore free for general use. Typeset by the authors using a Springer M'sX-package Cover design: design 0-production GmbH. Heidelberg SPIN: 10842488 46/3142LK -5 43210 - Printed on acid-free paper Preface Systems of polynomial equations arise throughout mathematics, science, and engineering. Algebraic geometry provides powerful theoretical techniques for studying the qualitative and quantitative features of their solution sets. Re cently developed algorithms have made theoretical aspects of the subject accessible to a broad range of mathematicians and scientists. The algorith mic approach to the subject has two principal aims: developing new tools for research within mathematics, and providing new tools for modeling and solv ing problems that arise in the sciences and engineering. A healthy synergy emerges, as new theorems yield new algorithms and emerging applications lead to new theoretical questions. This book presents algorithmic tools for algebraic geometry and experi mental applications of them. It also introduces a software system in which the tools have been implemented and with which the experiments can be carried out. Macaulay 2 is a computer algebra system devoted to supporting research in algebraic geometry, commutative algebra, and their applications. The reader of this book will encounter Macaulay 2 in the context of concrete applications and practical computations in algebraic geometry. The expositions of the algorithmic tools presented here are designed to serve as a useful guide for those wishing to bring such tools to bear on their own problems. A wide range of mathematical scientists should find these expositions valuable. This includes both the users of other programs similar to Macaulay 2 (for example, Singular and CoCoA) and those who are not interested in explicit machine computations at all. The chapters are ordered roughly by increasing mathematical difficulty. The first part of the book is meant to be accessible to graduate students and computer algebra users from across the mathematical sciences and is pri marily concerned with introducing Macaulay 2. The second part emphasizes the mathematics: each chapter exposes some domain of mathematics at an accessible level, presents the relevant algorithms, sometimes with proofs, and illustrates the use of the program. In both parts, each chapter comes with its own abstract and its own bibliography; the index at the back of the book covers all of them. One of the first computer algebra packages aimed at algebraic geometry was Macaulay, the predecessor of Macaulay 2, written during the years 1983- 1993 by Dave Bayer and Mike Stillman. Worst-case estimates suggested that trying to compute Grabner bases might be a hopeless approach to solving problems. But from the first prototype, Macaulay was successful surprisingly often, perhaps because of the geometrical origin of the problems attacked. Macaulay improved steadily during its first decade. It helped transform the theoretical notion of a projective resolution into an exciting new practical vi Preface research tool, and became widely used for research and teaching in com mutative algebra and algebraic geometry. It was possible to write routines in the top-level language, and many important algorithms were added by David Eisenbud and other users, enhancing the system and broadening its usefulness. There were certain practical drawbacks for the researcher who wanted to use Macaulay effectively. A minor annoyance was that only finite prime fields were available as coefficient rings. The major problem was that the language made available to users was primitive and barely supported high level development of new algorithms; it had few basic data types and didn't support the addition of new ones. Macaulay 2 is based on experience gained from writing and using its pre decessor Macaulay, but is otherwise a fresh start. It was written by Dan Grayson and Mike Stillman with the generous financial support of the U.S. National Science Foundation, with the work starting in 19931. It also incor porates some code from other authors: the package SINGULAR-FACTORy2 provides for factorization of polynomials; SINGULAR-LIBFAC3 uses FAC TORY to enable the computation of characteristic sets and thus the decom position of subvarieties into their irreducible components; and GNU Mp4 by Torbjorn Granlund and others provides for multiple precision arithmetic. Macaulay 2 aims to support efficient computation associated with a wide variety of high level mathematical objects, including Galois fields, number fields, polynomial rings, exterior algebras, Weyl algebras, quotient rings, ide als, modules, homomorphisms of rings and modules, graded modules, maps between graded modules, chain complexes, maps between chain complexes, free resolutions, algebraic varieties, and coherent sheaves. To make the system easily accessible, standard mathematical notation is followed closely. As with Macaulay, it was hoped that users would join in the further development of new algorithms for Macaulay 2, so the developers tried to make the language available to the users as powerful as possible, yet easy to use. Indeed, much of the high-level part of the system is written in the same language available to the user. This ensures that the user will find it just as 1 NSF grants DMS 92-10805, 92-10807, 96-23232, 96-22608, 99-70085, and 99- 70348. 2 SINGULAR-FACTORY, a subroutine library for factorization, by G.-M. Greuel, R. Stobbe, G. Pfister, H. Schoenemann, and J. Schmidt; available at ftp:/ /helios.mathematik. uni-kl.de/pub/Math/Singular /Factory /. 3 SINGULAR-LIBFAC, a subroutine library for characteristic sets and irreducible decomposition, by M. Messollen; available at ftp:/ /helios.mathematik.uni kl.de/pub/Math/Singular/Libfac/ . 4 GMP, a library for arbitrary precision arithmetic, by Torbjorn Granlund, John Amanatides, Paul Zimmermann, Ken Weber, Bennet Yee, Andreas Schwab, Robert Harley, Linus Nordberg, Kent Boortz, Kevin Ryde, and Guillaume Han rot; available at ftp://ftp.gnu.org/gnu/gmp/. Preface vii easy as the developers did to implement a new type of mathematical object or to modify the high-level aspects of the current algorithms. The language available to the user is interpreted. The interpreter itself is written in a convenient language designed to be mostly type-safe and to handle memory allocation and initialization automatically. For maximum effi ciency, the core mathematical algorithms are written in C++ and compiled, not interpreted. This includes the arithmetic operations of rings, modules, and matrices, the Grabner basis algorithm (in several enhanced versions, tai lored for various situations), several algorithms for computing free resolutions of modules, the algorithm for computing the Hilbert series of a graded ring or module, the algorithms for computing determinants and Pfaffians, the basis reduction algorithm, factoring, etc. In one way Macaulay 2 is like a standard computer algebra system, such as Mathematica or Maple: the user enters mathematical expressions at the keyboard, and the program computes the value of the expression and displays the answer. Here is the first input prompt offered to the user. i1 : In response to the prompt, the user may enter, for example, a simple arith metic expression. it : 3/5 + 7/11 68 01 = -- 55 01 : QQ The answer itself is displayed to the right of the output label 01 = and its type (or class) is displayed to the right of the following label. 01 : The symbol QQ appearing in this example denotes the class of all rational numbers, and is meant to be reminiscent of the notation Q. Macaulay 2 often finds itself being run in a window with horizontal scroll bars, so by default it does not wrap output lines, but instead lets them grow without bound. This book was produced by an automated mechanism that submits code provided by the authors to Macaulay 2 and incorporates the result into the text. Output lines that exceed the width of the pages of this book are indicated with ellipses, as in the following example. i2 : 100! 02 = 93326215443944152681699238856266700490715968264381621468592963895 ... viii Preface Next we describe an important difference between general computer al gebra systems (such as Maple and Mathematica) and Macaulay 2. Before entering an expression involving variables (such as x+y) into Macaulay 2 the user must first create a ring containing those variables. Rings are important objects of study in algebraic geometry; quotient rings of polynomial rings, for example, encapsulate the essential information about a system of polyno mial equations, including, for example, the field from which the coefficients are drawn. Often one has several rings under consideration at once, along with ring homomorphisms between them, so it is important to treat them as first-class objects in the computer, capable of being named and manipu lated the same way numbers and characters can be manipulated in simpler programming languages. Let's give a hint of the breadth of types of mathematical objects available in Macaulay 2 with some examples. In Macaulay 2 one defines a quotient ring of a polynomial ring R over the rational numbers by entering a command such as the one below. i3 : R = QQ[x,y,z]/(x"3-y"3-z"3) 03 = R 03 : QuotientRing Having done that, we can compute in the ring. i4 : (x+y+z)"3 2 2 3 2 2 2 2 3 04 = 3x Y + 3x*y + 2y + 3x z + 6x*y*z + 3y z + 3x*z + 3y*z + 2z 04 : R We can make matrices over the ring. i5 : b = vars R 05 = I x y z I 3 05 Matrix R <--- R i6 c = matrix {{x"2,y"2,z"2}} 06 = I x2 y2 z2 I 1 3 06 : Matrix R <--- R We can make modules over the ring. i7 : M = coker b 07 = cokernel I x y z 1 07 R-module, quotient of R i8 N = ker c 08 = image {2} x 0 -y2 -z2 {2} -y -z2 x2 0 {2} -z y2 0 x2 Preface ix 3 08 : R-module, submodule of R We can make projective resolutions of modules. i9 : res M 1 3 4 4 4 09 = R <-- R <-- R <-- R <-- R o 1 2 3 4 09 : ChainComplex We can make projective varieties. il0 : X = Proj R 010 = X 010 : ProjectiveVariety We can make coherent sheaves and compute their cohomology. ill : HH~l cotangent Sheaf X 1 011 = QQ all : QQ-module, free At this writing, Macaulay 2 is available for GNU /Linux and other flavors of Unix, and also for Microsoft Windows and the Macintosh operating system. Although it can be used as a free-standing program, it is most convenient to use it in an editor's buffer; Emacs (on Unix or Windows systems) or MPW on Macintosh systems are currently the editors of choice. To obtain Macaulay 2, download it from the website5 and unpack the file. Among the resulting files will be a file called Macaulay2/README. txt, which you should read. It will tell you how to run the setup script, and how to install a few lines of code in your emacs init file to enable you to run M2 in an emacs buffer and to edit Macaulay 2 code. A system administrator of a Unix system may optionally arrange for those lines of code to be available to every emacs user. The editors thank the authors of the chapters for their valuable contri butions and hard work, and the National Science Foundation for funding the development of Macaulay 2 and for partial funding of the authors who have contributed to this volume. May, 2001 David Eisenbud Daniel R. Grayson Michael E. Stillman Bernd Sturmfels 5 Macaulay 2, a software system for research in algebraic geometry, by Daniel R. Grayson and Michael E. Stillman; available online in source code form and compiled for various architectures at http://www . math. uiue. edu/Maeaulay2/. Table of Contents Preface .......................... " .. .. . . . . .. . . .. . . . . . . . . . . . . . . v List of Contributors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. xv Part I Introducing Macaulay 2 Ideals, Varieties and Macaulay 2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 Bernd Sturmfels 1 A Curve in Affine Three-Space . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 2 Intersecting Our Curve With a Surface. . . . . . . . . . . . . . . . . . . . . . . . . 4 3 Changing the Ambient Polynomial Ring. . . . . . . . . . . . . . . . . . . . . . . . 6 4 Monomials Under the Staircase. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 5 Pennies, Nickels, Dimes and Quarters . . . . . . . . . . . . . . . . . . . . . . . . .. 12 References .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 15 Projective Geometry and Homological Algebra . . . . . . . . . . . . . .. 17 David Eisenbud 1 The Twisted Cubic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 18 2 The Cotangent Bundle of jp» . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 20 3 The Cotangent Bundle of a Projective Variety .................. 24 4 Intersections by Serre's Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 26 5 A Mystery Variety in jp» ..................................... 28 Appendix A. How the "Mystery Variety" was Made . . . . . . . . . . . . . . . .. 37 References .... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 40 Data Types, Functions, and Programming . . . . . . . . . . . . . . . . . . .. 41 Daniel R. Grayson and Michael E. Stillman 1 Basic Data Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 41 2 Control Structures. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 44 3 Input and Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 46 4 Hash Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 48 5 Methods. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 52 6 Pointers to the Source Code . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 53 References ..................................................... 53 Teaching the Geometry of Schemes . . . . . . . . . . . . . . . . . . . . . . . . . .. 55 Gregory G. Smith and Bernd Sturmfels 1 Distinguished Open Sets ..................................... 55 2 Irreducibility. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 56 3 Singular Points ............................................. 58
The list of books you might like
