New below: Workshop preliminary program

How to come to IHP / Live stream

Registration

Special week

October 9 to 13, 2023

Lecture notes on the moment method

Location: Amphitheatre Darboux

Schedule Monday 10:00-12:00 and Tuesday 10:00-12:00

Lecture notes on resultants

Lecture notes on real algebraic geometry

Location: Amphitheatre Darboux

Schedule Tuesday 14:00-17:30 and Wednesday 10:00-12:00. Thursday 10:00-12:00 and new: 14:00-17:30 Friday 14:00-17:00

Location: Amphitheatre Darboux

Title: The quadratic equation revisited.

Abstract High school students learn how to express the solution of a quadratic equation in one unknown in terms of its three coefficients. What does this this formula matter? We offer an answer in terms of discriminants and data. This lecture invites the audience to a journey towards non-linear algebra.

Slides

Workshop: Geometry of Polynomial System Solving, Optimization and Topology

Organizers: C. D’Andrea, P. Lairez, M. Safey El Din, É. Schost, L. Zhi

October 16 to 20, 2023

Location: Amphitheatre Hermite

Polynomial systems encode a wide range of non-linear (but static) phenomena which arise in many applications. Non-linearity makes them non-trivial to handle, both from complexity and reliability viewpoints. Still, because of their importance for key applications e.g. in mechanism design and optimization amongst many others, various algorithmic approaches have been developed. During the last decades, tremendous achievements have been accomplished to design faster algorithms for polynomial system solving, extend their capabilities to tackle topological issues and understand their complexities. For instance, let us mention new families of algorithms to exploit algebraic and geometric properties of polynomial systems and their solution sets such as sparsity or weighted and multi-homogeneity, algorithms for understanding the topology of semi-algebraic sets (Betti numbers, connectivity queries), the raise of sums-of-squares certificates to certify emptiness over the reals of polynomial systems through symbolic-numeric approaches, and last but not least, qualitive advances in the line of Smale’s 17th problem.

Many challenges remain to be addressed to pave the way towards high performance polynomial system solvers tackling large scale applications. Topical issues lie in the combination of efficiency and certification, computing exact certificates of emptiness, understanding the geometry of polynomial systems and their solution sets to exploit better their properties algorithmically. This workshop will cover broadly all these topics.

Invited speakers

Christian Eder (Kaiserslautern, Germany) Recent advances in Gröbner basis algorithms and geometric applications

Abstract The tasks of designing innovative mathematical software and of solving complex research problems using computational methods are strongly mutually dependent. Developing a new generation of algorithms to considerably push the computational boundaries of nonlinear algebra, notably addressing polynomial system solving, is thus envitable. One important task of this process is to no longer use Gröbner bases only as a black box in higher level algorthms, but to optimize their computation with the geometric context in mind. In this talk, we will illustrate this idea by presenting a new algorithm for computing Gröbner bases of saturated polynomial ideals. Moreover, we introduce msolve, an open source software package build to provide this new generation of efficient and optimized algorithms for the community.
The contents of this talk are based on joint work with Jérémy Berthomieu and Mohab Safey El Din.

Elisenda Feliu (Copenhagen, Denmark) Positive solutions to polynomial systems and applications to reaction networks

Abstract The main object of study in the (algebraic) theory of reaction networks is the solution se of a system of parametric polynomial equations in the positive orthant. This system consists of polynomials with fixed support, the coefficients are linear in the parameters, but there might be some (proportionality) dependencies among the coefficients. The questions of interest concern properties of this system, and of its intersection with a family of parallel linear subspaces of complementary dimension. In this case, of relevance is to determine the possible number of zeros the system has when the parameters vary.
In this talk I will introduce the framework and the families of polynomial systems under study in full generality, and having the reaction networks as the main application. I will proceed to discuss recent results addressing the expected dimension of the solution sets and on how to decide whether the solution set admits a toric parametrization for all parameter values. The latter is relevant for the problem of counting solutions, and this connection will also be explained in the talk.

Jon Hauenstein (Notre Dame, USA) Some advances in numerical algebraic geometry for computing real solutions

Abstract. Numerical algebraic geometry provides a collection of algorithms for computing and analyzing solution sets of polynomial systems. This talk will discuss new techniques that have been developed in numerical algebraic geometry for focusing on real solution sets of polynomial systems. Several applications of these techniques will be presented such as computing smooth points on algebraic sets, approximate synthesis of mechanisms, and path planning for output mode switching.

Martin Helmer (Raleigh, USA) Conormal Spaces and Whitney Stratifications

Abstract. We describe a new algorithm for computing Whitney stratifications of complex projective varieties. The main ingredients are (a) an algebraic criterion, due to Lê and Teissier, which reformulates Whitney regularity in terms of conormal spaces and maps, and (b) a new interpretation of this conormal criterion via ideal saturations, which can be practically implemented on a computer. We show that this algorithm improves upon the existing state of the art by several orders of magnitude, even for relatively small input varieties. This is joint work with Vidit Nanda (Oxford).

Teresa Krick (Buenos Aires, Argentina) Non-negativity and rational sums of squares over zero-dimensional varieties

Abstract In this work in progress with Lorenzo Baldi and Bernard Mourrain, we extend previous results on univariate rational sums of squares, obtained with Bernard and Agnes Szanto, to the case of a non-negative rational polynomial on a basic zero-dimensional semi-algebraic set defined by rational polynomials.

Monique Laurent (Centrum Wiskunde & Informatica (CWI) Amsterdam, and Tilburg University, The Netherlands) Sums of squares approximations in polynomial optimization: performance analysis and degree bounds

Abstract Polynomial optimization deals with optimizing a polynomial function over a feasible region defined by polynomial inequalities, thus modeling a broad range of hard nonlinear nonconvex optimization problems. Hierarchies of tractable semidefinite relaxations have been introduced that are based on using sums of squares of polynomials as a ``proxy” for global nonnegativity. These hierarchies give bounds on the global minimum of the original problem with asymptotic convergence (under a minor compactness assumption). In this lecture we discuss recent results on the performance analysis of these hierarchies and related effective degree bounds for dedicated sums of squares representations of positive polynomials on some classes of compact semi-algebraic sets (including the hypercube, the sphere or the ball).

Anton Leykin (Atlanta, USA) u-generation: solving systems of polynomials equation-by-equation

Abstract We develop a new method that improves the efficiency of equation-by-equation homotopy continuation methods for solving polynomial systems. Our method is based on a novel geometric construction and reduces the total number of homotopy paths that must be numerically continued. These improvements may be applied to the basic algorithms of numerical algebraic geometry in the settings of both projective and multiprojective varieties. (This is joint work with T. Duff and J. I. Rodriguez.)

Zijia Li (Beijing, China) Polynomial Systems Arising in Paradoxical 6R Linkages

Abstract In this talk, we first provide a comprehensive definition of closed n-linkages and explain their mobility, typically denoted as n-6. We then focus on the intriguing subset of closed n-linkages with a mobility higher than n-6, known as paradoxical linkages. Based on the powerful tools of Bond Theory and the freezing technique, we present a thorough classification of n-linkages with a mobility of n-4 or higher, incorporating revolute, prismatic, or helical joints. Additionally, we explicitly derive strong necessary conditions for nR-linkages with a mobility of n-5. Utilizing these necessary conditions, we explore and discuss possible polynomial systems that arise in paradoxical 6R linkages.

Fatemeh Mohammadi (Leuven, Belgium) Polynomial systems arising in the formal verification of programs

Abstract. Multiple classical problems in the formal verification of programs such as reachability, termination, and template-based synthesis can be reduced to solving polynomial systems of equations. In this talk, I will describe the primary objects and these connections. In particular, I will show how the algebraic and geometric techniques can be applied, enhancing the scalability and completeness for such problems.

Marc Moreno Maza (London, Canada) Modular algorithms for computing triangular decompositions of polynomial systems [abstract](https://rtca2023.github.io/pages_Paris/files_m5/abstract_moreno-maza.pdf)
Chenqi Mou (Beijing, China) Chordal Graphs in Triangular Decomposition in Top-Down Style

Abstract In this talk, I will present the connections between chordal graphs from graph theory and triangular decomposition in top-down style from symbolic computation, including the underlying theories, algorithms, and applications in biology. Viewing triangular decomposition in top-down style as polynomial generalization of Gaussian elimination, we show that all the polynomial sets, including all the computed triangular sets, appearing in several typical top-down algorithms for triangular decomposition have associated graphs which are subgraphs of the chordal graph of the input polynomial set. These theoretical results can be interpreted as “triangular decomposition in top-down style preserves chordality” and are further used to design sparse triangular decomposition for polynomial sets which are sparse with respect to their variables. Sparse triangular decomposition is shown to be more efficient than ordinary one experimentally, and its application on computing equilibria of biological dynamic systems will also be reported.
This talk is based on the joint work with Yang Bai, Jiahua Lai, and Wenwen Ju.

Bernard Mourrain (Inria, Sophia-Antipolis) Solving by duality

Abstract. Finding the common roots of a set of polynomial equations is a problem that appears in many contexts and applications. Standard approaches for solving this difficult question, such as Grobner bases, border basis, triangular sets, etc. are based on polynomial reductions but their instability against numerical approximations can be critical. In this talk, we will describe a dual approach which focuses on linear functionals vanishing at the roots. We will review the properties of Truncated Normal Forms, the connexion with classical computer algebra approaches and resultants. We will also detail the dual approach in the context of optimisation problems and for analysing isolated singularities. Examples from geometric modeling, robotics and tensor decomposition will illustrate the numerical behavior of these dual methods.

Cordian Riener (Tromso, Norway) The geometry of the Vandermonde map at infinity

Abstract The Vandermonde map is the polynomial map given by the power-sum polynomials. We study the geometry of the image of the nonnegative orthant under under this map and focus on the limit as the number of variables approaches infinity. We will show, the geometry of this limit is the key to new undecidability results in nonnegativity of symmetric polynomials and deciding validity of trace inequalities in linear algebra.

Pierre-Jean Spaenlehauer (Nancy, France) Dimension results for sparse systems homogenized via rational polytopes

Abstract. A classical method to compute with sparse polynomials is to homogenize them with respect to Newton polytopes, regarding them as homogeneous elements of Cartier degrees in the Cox ring of a projective toric variety. In this talk, we investigate subvarieties defined by generic polynomial systems in the Cox ring when the degrees are non-necessarily Cartier, with a view towards identifying alternative toric homogenizations that are suitable for practical computations. Joint work with Matías Bender.

Topical day: Mechanism Design and Computer Algebra

October 24, 2023

Organizer: J. Schicho

Location: Amphitheatre Darboux

Speakers. Maria Alberich-Carraminana, Didier Henrion, Hans-Peter Schröcker, Josef Schicho

Schedule:

Abstract. Let $G$ be a graph with $n$ vertices and $e$ edges. The computation of the position of $n$ points in the plane such that for any two vertices in the graph connected by an edge, the distance between the two corresponding points is given, is equivalent to the inverse kinematic problem for a (highly parallel) planar mechanism with revolute joints. If the graph is a Laman graph, then the solution set is generically a finite set of orbits under the group of Euclidean displacements, and can be assigned a Galois group (which is associated to the field extension needed to express the solutions exactly). We explain some geometric ideas for analyzing the Galois group. Using these ideas, we determine the number of components of the solution set for graphs that have the property that the above position problem is generically solvable.
Abstract. Robotic manipulation of cloth presents a complex challenge due to the infinite-dimensional shape-state space of textiles. This complexity makes accurate state estimation a daunting task. To address this issue, we introduce the concept of dGLI Cloth Coordinates—a finite, low-dimensional representation of cloth states. This novel approach enables us to effectively distinguish among a wide range of folded states, paving the way for efficient learning methods in cloth manipulation planning and control.


Preliminary Program of the workshop - All talks in Amphitheater Hermite

Live stream

    Monday October 16th, 2023
8:45 Welcome coffee
9:15-9:30 Opening
Chair: Carlos D'Andrea
9:30-10:30 Pierre-Jean Spaenlehauer
Dimension results for sparse systems homogenized via rational polytopes
Abstract. A classical method to compute with sparse polynomials is to homogenize them with respect to Newton polytopes, regarding them as homogeneous elements of Cartier degrees in the Cox ring of a projective toric variety. In this talk, we investigate subvarieties defined by generic polynomial systems in the Cox ring when the degrees are non-necessarily Cartier, with a view towards identifying alternative toric homogenizations that are suitable for practical computations. Joint work with Matías Bender.
10:30-11:00 Coffee break
11:00-12:00 Chenqi Mou
Chordal Graphs in Triangular Decomposition in Top-Down Style
Abstract. In this talk, I will present the connections between chordal graphs from graph theory and triangular decomposition in top-down style from symbolic computation, including the underlying theories, algorithms, and applications in biology. Viewing triangular decomposition in top-down style as polynomial generalization of Gaussian elimination, we show that all the polynomial sets, including all the computed triangular sets, appearing in several typical top-down algorithms for triangular decomposition have associated graphs which are subgraphs of the chordal graph of the input polynomial set. These theoretical results can be interpreted as “triangular decomposition in top-down style preserves chordality” and are further used to design sparse triangular decomposition for polynomial sets which are sparse with respect to their variables. Sparse triangular decomposition is shown to be more efficient than ordinary one experimentally, and its application on computing equilibria of biological dynamic systems will also be reported. This talk is based on the joint work with Yang Bai, Jiahua Lai, and Wenwen Ju.
12:00-14:30 Lunch break
Chair: Teresa Krick
14:30-15:30 Marc Moreno Maza
Modular algorithms for computing triangular decompositions of polynomial systems
abstract_moreno-maza.pdf
15:30-16:00 Coffee break
    Tuesday October 17th, 2023
Chair: Jérémy Berthomieu
9:30-10:30 Fatemeh Mohammadi
Polynomial systems arising in the formal verification of programs
Abstract. Multiple classical problems in the formal verification of programs such as reachability, termination, and template-based synthesis can be reduced to solving polynomial systems of equations. In this talk, I will describe the primary objects and these connections. In particular, I will show how the algebraic and geometric techniques can be applied, enhancing the scalability and completeness for such problems.
10:30-11:00 Coffee break
11:00-12:00 Elisenda Feliu
Positive solutions to polynomial systems and applications to reaction networks
Abstract. The main object of study in the (algebraic) theory of reaction networks is the solution se of a system of parametric polynomial equations in the positive orthant. This system consists of polynomials with fixed support, the coefficients are linear in the parameters, but there might be some (proportionality) dependencies among the coefficients. The questions of interest concern properties of this system, and of its intersection with a family of parallel linear subspaces of complementary dimension. In this case, of relevance is to determine the possible number of zeros the system has when the parameters vary. In this talk I will introduce the framework and the families of polynomial systems under study in full generality, and having the reaction networks as the main application. I will proceed to discuss recent results addressing the expected dimension of the solution sets and on how to decide whether the solution set admits a toric parametrization for all parameter values. The latter is relevant for the problem of counting solutions, and this connection will also be explained in the talk.
12:00-14:30 Lunch break
Chair: Mohab Safey El Din
14:30-15:30 Zijia Li
Polynomial Systems Arising in Paradoxical 6R Linkages
Abstract. In this talk, we first provide a comprehensive definition of closed n-linkages and explain their mobility, typically denoted as n-6. We then focus on the intriguing subset of closed n-linkages with a mobility higher than n-6, known as paradoxical linkages. Based on the powerful tools of Bond Theory and the freezing technique, we present a thorough classification of n-linkages with a mobility of n-4 or higher, incorporating revolute, prismatic, or helical joints. Additionally, we explicitly derive strong necessary conditions for nR-linkages with a mobility of n-5. Utilizing these necessary conditions, we explore and discuss possible polynomial systems that arise in paradoxical 6R linkages.
15:30-16:00 Coffee
    Wednesday October 18th, 2023
Chair: Pierre Lairez
9:30-10:30 Jon Hauenstein
Some advances in numerical algebraic geometry for computing real solutions
Abstract. Numerical algebraic geometry provides a collection of algorithms for computing and analyzing solution sets of polynomial systems. This talk will discuss new techniques that have been developed in numerical algebraic geometry for focusing on real solution sets of polynomial systems. Several applications of these techniques will be presented such as computing smooth points on algebraic sets, approximate synthesis of mechanisms, and path planning for output mode switching.
10:30-11:00 Coffee break
11:00-12:00 Anton Leykin
u-generation: solving systems of polynomials equation-by-equation
Abstract. We develop a new method that improves the efficiency of equation-by-equation homotopy continuation methods for solving polynomial systems. Our method is based on a novel geometric construction and reduces the total number of homotopy paths that must be numerically continued. These improvements may be applied to the basic algorithms of numerical algebraic geometry in the settings of both projective and multiprojective varieties. (This is joint work with T. Duff and J. I. Rodriguez.)
Free afternoon
18:30 Reception.
    Thursday October 19th, 2023
Chair: Lihong Zhi
9:30-10:30 Monique Laurent
Sums of squares approximations in polynomial optimization: performance analysis and degree bounds
Abstract. Polynomial optimization deals with optimizing a polynomial function over a feasible region defined by polynomial inequalities, thus modeling a broad range of hard nonlinear nonconvex optimization problems. Hierarchies of tractable semidefinite relaxations have been introduced that are based on using sums of squares of polynomials as a ``proxy” for global nonnegativity. These hierarchies give bounds on the global minimum of the original problem with asymptotic convergence (under a minor compactness assumption). In this lecture we discuss recent results on the performance analysis of these hierarchies and related effective degree bounds for dedicated sums of squares representations of positive polynomials on some classes of compact semi-algebraic sets (including the hypercube, the sphere or the ball).
10:30-11:00 Coffee break
11:00-12:00 Bernard Mourrain
Solving by duality
Abstract. Finding the common roots of a set of polynomial equations is a problem that appears in many contexts and applications. Standard approaches for solving this difficult question, such as Grobner bases, border basis, triangular sets, etc. are based on polynomial reductions but their instability against numerical approximations can be critical. In this talk, we will describe a dual approach which focuses on linear functionals vanishing at the roots. We will review the properties of Truncated Normal Forms, the connexion with classical computer algebra approaches and resultants. We will also detail the dual approach in the context of optimisation problems and for analysing isolated singularities. Examples from geometric modeling, robotics and tensor decomposition will illustrate the numerical behavior of these dual methods.
12:00-14:30 Lunch break
Chair: Elisenda Feliu
14:30-15:30 Teresa Krick
Non-negativity and rational sums of squares over zero-dimensional varieties
Abstract. In this work in progress with Lorenzo Baldi and Bernard Mourrain, we extend previous results on univariate rational sums of squares, obtained with Bernard and Agnes Szanto, to the case of a non-negative rational polynomial on a basic zero-dimensional semi-algebraic set defined by rational polynomials.
15:30-16:00 Coffee break
16:00-17:00 Martin Helmer
Conormal Spaces and Whitney Stratifications
Abstract. We describe a new algorithm for computing Whitney stratifications of complex projective varieties. The main ingredients are (a) an algebraic criterion, due to Lê and Teissier, which reformulates Whitney regularity in terms of conormal spaces and maps, and (b) a new interpretation of this conormal criterion via ideal saturations, which can be practically implemented on a computer. We show that this algorithm improves upon the existing state of the art by several orders of magnitude, even for relatively small input varieties. This is joint work with Vidit Nanda (Oxford).
    Friday October 20th, 2023
9:30-10:30 Cordian Riener
The geometry of the Vandermonde map at infinity
Abstract. The Vandermonde map is the polynomial map given by the power-sum polynomials. We study the geometry of the image of the nonnegative orthant under under this map and focus on the limit as the number of variables approaches infinity. We will show, the geometry of this limit is the key to new undecidability results in nonnegativity of symmetric polynomials and deciding validity of trace inequalities in linear algebra.
10:30-11:00 Coffee break
11:00-12:00 Christian Eder
Recent advances in Gröbner basis algorithms and geometric applications
Abstract. The tasks of designing innovative mathematical software and of solving complex research problems using computational methods are strongly mutually dependent. Developing a new generation of algorithms to considerably push the computational boundaries of nonlinear algebra, notably addressing polynomial system solving, is thus envitable. One important task of this process is to no longer use Gröbner bases only as a black box in higher level algorthms, but to optimize their computation with the geometric context in mind. In this talk, we will illustrate this idea by presenting a new algorithm for computing Gröbner bases of saturated polynomial ideals. Moreover, we introduce msolve, an open source software package build to provide this new generation of efficient and optimized algorithms for the community. The contents of this talk are based on joint work with Jérémy Berthomieu and Mohab Safey El Din.
End of the workshop

Code of Conduct