Computations in Algebraic Geometry with Macaulay 2
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
The Sagemanifolds Project
Tensor calculus with free softwares: the SageManifolds project Eric´ Gourgoulhon1, Micha l Bejger2 1Laboratoire Univers et Th´eories (LUTH) CNRS / Observatoire de Paris / Universit´eParis Diderot 92190 Meudon, France http://luth.obspm.fr/~luthier/gourgoulhon/ 2Centrum Astronomiczne im. M. Kopernika (CAMK) Warsaw, Poland http://users.camk.edu.pl/bejger/ Encuentros Relativistas Espa~noles2014 Valencia 1-5 September 2014 Eric´ Gourgoulhon, Micha l Bejger (LUTH, CAMK) SageManifolds ERE2014, Valencia, 2 Sept. 2014 1 / 44 Outline 1 Differential geometry and tensor calculus on a computer 2 Sage: a free mathematics software 3 The SageManifolds project 4 SageManifolds at work: the Mars-Simon tensor example 5 Conclusion and perspectives Eric´ Gourgoulhon, Micha l Bejger (LUTH, CAMK) SageManifolds ERE2014, Valencia, 2 Sept. 2014 2 / 44 Differential geometry and tensor calculus on a computer Outline 1 Differential geometry and tensor calculus on a computer 2 Sage: a free mathematics software 3 The SageManifolds project 4 SageManifolds at work: the Mars-Simon tensor example 5 Conclusion and perspectives Eric´ Gourgoulhon, Micha l Bejger (LUTH, CAMK) SageManifolds ERE2014, Valencia, 2 Sept. 2014 3 / 44 In 1969, during his PhD under Pirani supervision at King's College, Ray d'Inverno wrote ALAM (Atlas Lisp Algebraic Manipulator) and used it to compute the Riemann tensor of Bondi metric. The original calculations took Bondi and his collaborators 6 months to go. The computation with ALAM took 4 minutes and yield to the discovery of 6 errors in the original paper [J.E.F. Skea, Applications of SHEEP (1994)] In the early 1970's, ALAM was rewritten in the LISP programming language, thereby becoming machine independent and renamed LAM The descendant of LAM, called SHEEP (!), was initiated in 1977 by Inge Frick Since then, many softwares for tensor calculus have been developed.. -
Using Macaulay2 Effectively in Practice
Using Macaulay2 effectively in practice Mike Stillman ([email protected]) Department of Mathematics Cornell 22 July 2019 / IMA Sage/M2 Macaulay2: at a glance Project started in 1993, Dan Grayson and Mike Stillman. Open source. Key computations: Gr¨obnerbases, free resolutions, Hilbert functions and applications of these. Rings, Modules and Chain Complexes are first class objects. Language which is comfortable for mathematicians, yet powerful, expressive, and fun to program in. Now a community project Journal of Software for Algebra and Geometry (started in 2009. Now we handle: Macaulay2, Singular, Gap, Cocoa) (original editors: Greg Smith, Amelia Taylor). Strong community: including about 2 workshops per year. User contributed packages (about 200 so far). Each has doc and tests, is tested every night, and is distributed with M2. Lots of activity Over 2000 math papers refer to Macaulay2. History: 1976-1978 (My undergrad years at Urbana) E. Graham Evans: asked me to write a program to compute syzygies, from Hilbert's algorithm from 1890. Really didn't work on computers of the day (probably might still be an issue!). Instead: Did computation degree by degree, no finishing condition. Used Buchsbaum-Eisenbud \What makes a complex exact" (by hand!) to see if the resulting complex was exact. Winfried Bruns was there too. Very exciting time. History: 1978-1983 (My grad years, with Dave Bayer, at Harvard) History: 1978-1983 (My grad years, with Dave Bayer, at Harvard) I tried to do \real mathematics" but Dave Bayer (basically) rediscovered Groebner bases, and saw that they gave an algorithm for computing all syzygies. I got excited, dropped what I was doing, and we programmed (in Pascal), in less than one week, the first version of what would be Macaulay. -
A NOTE on COMPLETE RESOLUTIONS 1. Introduction Let
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 138, Number 11, November 2010, Pages 3815–3820 S 0002-9939(2010)10422-7 Article electronically published on May 20, 2010 A NOTE ON COMPLETE RESOLUTIONS FOTINI DEMBEGIOTI AND OLYMPIA TALELLI (Communicated by Birge Huisgen-Zimmermann) Abstract. It is shown that the Eckmann-Shapiro Lemma holds for complete cohomology if and only if complete cohomology can be calculated using com- plete resolutions. It is also shown that for an LHF-group G the kernels in a complete resolution of a ZG-module coincide with Benson’s class of cofibrant modules. 1. Introduction Let G be a group and ZG its integral group ring. A ZG-module M is said to admit a complete resolution (F, P,n) of coincidence index n if there is an acyclic complex F = {(Fi,ϑi)| i ∈ Z} of projective modules and a projective resolution P = {(Pi,di)| i ∈ Z,i≥ 0} of M such that F and P coincide in dimensions greater than n;thatis, ϑn F : ···→Fn+1 → Fn −→ Fn−1 → ··· →F0 → F−1 → F−2 →··· dn P : ···→Pn+1 → Pn −→ Pn−1 → ··· →P0 → M → 0 A ZG-module M is said to admit a complete resolution in the strong sense if there is a complete resolution (F, P,n)withHomZG(F,Q) acyclic for every ZG-projective module Q. It was shown by Cornick and Kropholler in [7] that if M admits a complete resolution (F, P,n) in the strong sense, then ∗ ∗ F ExtZG(M,B) H (HomZG( ,B)) ∗ ∗ where ExtZG(M, ) is the P-completion of ExtZG(M, ), defined by Mislin for any group G [13] as k k−r r ExtZG(M,B) = lim S ExtZG(M,B) r>k where S−mT is the m-th left satellite of a functor T . -
The Geometry of Syzygies
The Geometry of Syzygies A second course in Commutative Algebra and Algebraic Geometry David Eisenbud University of California, Berkeley with the collaboration of Freddy Bonnin, Clement´ Caubel and Hel´ ene` Maugendre For a current version of this manuscript-in-progress, see www.msri.org/people/staff/de/ready.pdf Copyright David Eisenbud, 2002 ii Contents 0 Preface: Algebra and Geometry xi 0A What are syzygies? . xii 0B The Geometric Content of Syzygies . xiii 0C What does it mean to solve linear equations? . xiv 0D Experiment and Computation . xvi 0E What’s In This Book? . xvii 0F Prerequisites . xix 0G How did this book come about? . xix 0H Other Books . 1 0I Thanks . 1 0J Notation . 1 1 Free resolutions and Hilbert functions 3 1A Hilbert’s contributions . 3 1A.1 The generation of invariants . 3 1A.2 The study of syzygies . 5 1A.3 The Hilbert function becomes polynomial . 7 iii iv CONTENTS 1B Minimal free resolutions . 8 1B.1 Describing resolutions: Betti diagrams . 11 1B.2 Properties of the graded Betti numbers . 12 1B.3 The information in the Hilbert function . 13 1C Exercises . 14 2 First Examples of Free Resolutions 19 2A Monomial ideals and simplicial complexes . 19 2A.1 Syzygies of monomial ideals . 23 2A.2 Examples . 25 2A.3 Bounds on Betti numbers and proof of Hilbert’s Syzygy Theorem . 26 2B Geometry from syzygies: seven points in P3 .......... 29 2B.1 The Hilbert polynomial and function. 29 2B.2 . and other information in the resolution . 31 2C Exercises . 34 3 Points in P2 39 3A The ideal of a finite set of points . -
How to Make Free Resolutions with Macaulay2
How to make free resolutions with Macaulay2 Chris Peterson and Hirotachi Abo 1. What are syzygies? Let k be a field, let R = k[x0, . , xn] be the homogeneous coordinate ring of Pn and let X be a projective variety in Pn. Consider the ideal I(X) of X. Assume that {f0, . , ft} is a generating set of I(X) and that each polynomial fi has degree di. We can express this by saying that we have a surjective homogenous map of graded S-modules: t M R(−di) → I(X), i=0 where R(−di) is a graded R-module with grading shifted by −di, that is, R(−di)k = Rk−di . In other words, we have an exact sequence of graded R-modules: Lt F i=0 R(−di) / R / Γ(X) / 0, M |> MM || MMM || MMM || M& || I(X) 7 D ppp DD pp DD ppp DD ppp DD 0 pp " 0 where F = (f0, . , ft). Example 1. Let R = k[x0, x1, x2]. Consider the union P of three points [0 : 0 : 1], [1 : 0 : 0] and [0 : 1 : 0]. The corresponding ideals are (x0, x1), (x1, x2) and (x2, x0). The intersection of these ideals are (x1x2, x0x2, x0x1). Let I(P ) be the ideal of P . In Macaulay2, the generating set {x1x2, x0x2, x0x1} for I(P ) is described as a 1 × 3 matrix with gens: i1 : KK=QQ o1 = QQ o1 : Ring 1 -- the class of all rational numbers i2 : ringP2=KK[x_0,x_1,x_2] o2 = ringP2 o2 : PolynomialRing i3 : P1=ideal(x_0,x_1); P2=ideal(x_1,x_2); P3=ideal(x_2,x_0); o3 : Ideal of ringP2 o4 : Ideal of ringP2 o5 : Ideal of ringP2 i6 : P=intersect(P1,P2,P3) o6 = ideal (x x , x x , x x ) 1 2 0 2 0 1 o6 : Ideal of ringP2 i7 : gens P o7 = | x_1x_2 x_0x_2 x_0x_1 | 1 3 o7 : Matrix ringP2 <--- ringP2 Definition. -
Luis David Garcıa Puente
Luis David Garc´ıa Puente Department of Mathematics and Statistics (936) 294-1581 Sam Houston State University [email protected] Huntsville, TX 77341–2206 http://www.shsu.edu/ldg005/ Professional Preparation Universidad Nacional Autonoma´ de Mexico´ (UNAM) Mexico City, Mexico´ B.S. Mathematics (with Honors) 1999 Virginia Polytechnic Institute and State University Blacksburg, VA Ph.D. Mathematics 2004 – Advisor: Reinhard Laubenbacher – Dissertation: Algebraic Geometry of Bayesian Networks University of California, Berkeley Berkeley, CA Postdoctoral Fellow Summer 2004 – Mentor: Lior Pachter Mathematical Sciences Research Institute (MSRI) Berkeley, CA Postdoctoral Fellow Fall 2004 – Mentor: Bernd Sturmfels Texas A&M University College Station, TX Visiting Assistant Professor 2005 – 2007 – Mentor: Frank Sottile Appointments Colorado College Colorado Springs, CO Professor of Mathematics and Computer Science 2021 – Sam Houston State University Huntsville, TX Professor of Mathematics 2019 – 2021 Sam Houston State University Huntsville, TX Associate Department Chair Fall 2017 – 2021 Sam Houston State University Huntsville, TX Associate Professor of Mathematics 2013 – 2019 Statistical and Applied Mathematical Sciences Institute Research Triangle Park, NC SAMSI New Researcher fellowship Spring 2009 Sam Houston State University Huntsville, TX Assistant Professor of Mathematics 2007 – 2013 Virginia Bioinformatics Institute (Virginia Tech) Blacksburg, VA Graduate Research Assistant Spring 2004 Virginia Polytechnic Institute and State University Blacksburg, -
Persistence in Algebraic Specifications
Persistence in Algebraic Specifications Freek Wiedijk Persistence in Algebraic Specifications Persistence in Algebraic Specifications Academisch Proefschrift ter verkrijging van de graad van doctor aan de Universiteit van Amsterdam, op gezag van de Rector Magnificus prof. dr. P.W.M. de Meijer, in het openbaar te verdedigen in de Aula der Universiteit (Oude Lutherse Kerk, ingang Singel 411, hoek Spui), op donderdag 12 december 1991 te 12.00 uur door Frederik Wiedijk geboren te Haarlem Promotores: prof. dr. P. Klint & prof. dr. J.A. Bergstra Faculteit Wiskunde en Informatica Het hier beschreven onderzoek werd mede mogelijk gemaakt door steun van de Ne- derlandse organisatie voor Wetenschappelijk Onderzoek (voorheen de nederlandse organisatie voor Zuiver-Wetenschappelijk Onderzoek) binnen project nummer 612- 317-013 getiteld Executeerbaar maken van algebraïsche specificaties. voor Jan Truijens Contents Contents Acknowledgements 1 Introduction 1 1.1 Motivation and general overview 2 1.2 Examples of erroneous specifications 13 1.3 Background and related work 22 2. Theory 25 2.1 Basic notions 27 2.1.1. Equational specifications 28 2.1.2 Term rewriting systems 29 2.2 Termination 30 2.2.1 Weak and strong termination 33 2.2.2 Path orderings 34 2.2.3 Decidability 40 2.3 Confluence 41 2.3.1 Weak and strong confluence 42 2.3.2 Overlapping 44 2.3.3 Decidability 45 2.4 Reduction strategies 46 2.5 Persistence 48 2.5.1 Weak and strong persistence 52 2.5.2 Term approximations and bases 54 2.5.3 Normal form analysis 56 2.5.4 Decidability 60 2.6 Primitive recursive -
Matchgates Revisited
THEORY OF COMPUTING, Volume 10 (7), 2014, pp. 167–197 www.theoryofcomputing.org RESEARCH SURVEY Matchgates Revisited Jin-Yi Cai∗ Aaron Gorenstein Received May 17, 2013; Revised December 17, 2013; Published August 12, 2014 Abstract: We study a collection of concepts and theorems that laid the foundation of matchgate computation. This includes the signature theory of planar matchgates, and the parallel theory of characters of not necessarily planar matchgates. Our aim is to present a unified and, whenever possible, simplified account of this challenging theory. Our results include: (1) A direct proof that the Matchgate Identities (MGI) are necessary and sufficient conditions for matchgate signatures. This proof is self-contained and does not go through the character theory. (2) A proof that the MGI already imply the Parity Condition. (3) A simplified construction of a crossover gadget. This is used in the proof of sufficiency of the MGI for matchgate signatures. This is also used to give a proof of equivalence between the signature theory and the character theory which permits omittable nodes. (4) A direct construction of matchgates realizing all matchgate-realizable symmetric signatures. ACM Classification: F.1.3, F.2.2, G.2.1, G.2.2 AMS Classification: 03D15, 05C70, 68R10 Key words and phrases: complexity theory, matchgates, Pfaffian orientation 1 Introduction Leslie Valiant introduced matchgates in a seminal paper [24]. In that paper he presented a way to encode computation via the Pfaffian and Pfaffian Sum, and showed that a non-trivial, though restricted, fragment of quantum computation can be simulated in classical polynomial time. Underlying this magic is a way to encode certain quantum states by a classical computation of perfect matchings, and to simulate certain ∗Supported by NSF CCF-0914969 and NSF CCF-1217549. -
Research in Algebraic Geometry with Macaulay2
Research in Algebraic Geometry with Macaulay2 1 Macaulay2 a software system for research in algebraic geometry resultants the community, packages, the journal, and workshops 2 Tutorial the ideal of a monomial curve elimination: projecting the twisted cubic quotients and saturation 3 the package PHCpack.m2 solving polynomial systems numerically MCS 507 Lecture 39 Mathematical, Statistical and Scientific Software Jan Verschelde, 25 November 2019 Scientific Software (MCS 507) Macaulay2 L-39 25November2019 1/32 Research in Algebraic Geometry with Macaulay2 1 Macaulay2 a software system for research in algebraic geometry resultants the community, packages, the journal, and workshops 2 Tutorial the ideal of a monomial curve elimination: projecting the twisted cubic quotients and saturation 3 the package PHCpack.m2 solving polynomial systems numerically Scientific Software (MCS 507) Macaulay2 L-39 25November2019 2/32 Macaulay2 D. R. Grayson and M. E. Stillman: Macaulay2, a software system for research in algebraic geometry, available at https://faculty.math.illinois.edu/Macaulay2. Funded by the National Science Foundation since 1992. Its source is at https://github.com/Macaulay2/M2; licence: GNU GPL version 2 or 3; try it online at http://habanero.math.cornell.edu:3690. Several workshops are held each year. Packages extend the functionality. The Journal of Software for Algebra and Geometry (at https://msp.org/jsag) started its first volume in 2009. The SageMath interface to Macaulay2 requires its binary M2 to be installed on your computer. Scientific -
An Alternative Algorithm for Computing the Betti Table of a Monomial Ideal 3
AN ALTERNATIVE ALGORITHM FOR COMPUTING THE BETTI TABLE OF A MONOMIAL IDEAL MARIA-LAURA TORRENTE AND MATTEO VARBARO Abstract. In this paper we develop a new technique to compute the Betti table of a monomial ideal. We present a prototype implementation of the resulting algorithm and we perform numerical experiments suggesting a very promising efficiency. On the way of describing the method, we also prove new constraints on the shape of the possible Betti tables of a monomial ideal. 1. Introduction Since many years syzygies, and more generally free resolutions, are central in purely theo- retical aspects of algebraic geometry; more recently, after the connection between algebra and statistics have been initiated by Diaconis and Sturmfels in [DS98], free resolutions have also become an important tool in statistics (for instance, see [D11, SW09]). As a consequence, it is fundamental to have efficient algorithms to compute them. The usual approach uses Gr¨obner bases and exploits a result of Schreyer (for more details see [Sc80, Sc91] or [Ei95, Chapter 15, Section 5]). The packages for free resolutions of the most used computer algebra systems, like [Macaulay2, Singular, CoCoA], are based on these techniques. In this paper, we introduce a new algorithmic method to compute the minimal graded free resolution of any finitely generated graded module over a polynomial ring such that some (possibly non- minimal) graded free resolution is known a priori. We describe this method and we present the resulting algorithm in the case of monomial ideals in a polynomial ring, in which situ- ation we always have a starting nonminimal graded free resolution. -
Prizes and Awards Session
PRIZES AND AWARDS SESSION Wednesday, July 12, 2021 9:00 AM EDT 2021 SIAM Annual Meeting July 19 – 23, 2021 Held in Virtual Format 1 Table of Contents AWM-SIAM Sonia Kovalevsky Lecture ................................................................................................... 3 George B. Dantzig Prize ............................................................................................................................. 5 George Pólya Prize for Mathematical Exposition .................................................................................... 7 George Pólya Prize in Applied Combinatorics ......................................................................................... 8 I.E. Block Community Lecture .................................................................................................................. 9 John von Neumann Prize ......................................................................................................................... 11 Lagrange Prize in Continuous Optimization .......................................................................................... 13 Ralph E. Kleinman Prize .......................................................................................................................... 15 SIAM Prize for Distinguished Service to the Profession ....................................................................... 17 SIAM Student Paper Prizes .................................................................................................................... -
PLS Toolbox 7.0 Pltternf - Plots a 3D Ternary Diagram with Frequency of Occurrence
plotgui - Interactive data viewer. logdecay - Mean centers and variance scales a matrix using the log plttern - Plots a 2D ternary diagram. decay of the variable axis. PLS_Toolbox 7.0 pltternf - Plots a 3D ternary diagram with frequency of occurrence. lsq2top - Fits a polynomial to the top/(bottom) of data. querydb - Executes a query on a database defined by connection mdcheck - Missing Data Checker and infiller. Quick-Reference Card string. med2top - Fits a constant to top/(bottom) of data. Copyright © Eigenvector Research, 2012 reportwriter - Write a summary of the analysis including medcn - Median center scales matrix to median zero. associated figures to html/word/powerpoint. mncn - Scale matrix to mean zero. Help and Information rwb - Red white and blue color map. mscorr - Multiplicative scatter/signal correction (MSC). setpath - Modifies and saves current directory to the MATLAB normaliz - Normalize rows of matrix. helppls - Context related help on the PLS_Toolbox. search path. npreprocess - Preprocessing of multi-way arrays. readme - Release notes for Version 4.1 of PLS_Toolbox. snabsreadr - Reads Stellarnet ABS XY files. oscapp - Applies OSC model to new data. demos - Demo list for the PLS_Toolbox. spcreadr - Reads a Galactic SPC file. osccalc - Calculates orthogonal signal correction (OSC). evricompatibility - Tests for inter-product compatibility of trendtool - Univariate trend analysis tool. poissonscale - Perform Poisson scaling with scaling offset. Eigenvector toolboxes. vline - Adds vertical lines to figure at specified locations. polyinterp - Polynomial interpolation, smoothing, and evridebug - Checks the PLS_Toolbox installation for problems. writeasf - Writes AIT ASF files from a dataset object. differentiation. evridir - Locate and or create EVRI home directory. writecsv - Export a DataSet object to a comma-separated values preprocess - Selection and application of standard preprocessing evriinstall - Install Eigenvector Research Product.