Keywords
 A5.5.1. Geometrical modeling
 A6.1. Methods in mathematical modeling
 A8.3. Geometry, Topology
 A8.4. Computer Algebra
 B9.5.1. Computer science
 B9.5.2. Mathematics
1 Team members, visitors, external collaborators
Research Scientists
 Bernard Mourrain [Team leader, Inria, Senior Researcher, HDR]
 Laurent Busé [Inria, Researcher, HDR]
 Evelyne Hubert [Inria, Senior Researcher, HDR]
 Angelos Mantzaflaris [Inria, Researcher]
Faculty Members
 Ioannis Emiris [Université nationale et capodistrienne d'Athènes, HDR]
 André Galligo [Univ Côte d'Azur, Professor]
 Adam Parusinski [Univ Côte d'Azur, Professor, until Feb 2020, HDR]
PostDoctoral Fellows
 Navid Nemati [Inria, from Mar 2020]
 Xiao Xiao [Inria, until Oct 2020]
PhD Students
 Lorenzo Baldi [Inria]
 Ayoub Belhachmi [Schlumberger Limited, from Aug 2020]
 Riccardo Di Dio [Univ Côte d'Azur]
 Edgar Fuentes Figueroa [Artelys]
 Pablo Gonzalez Mazon [Inria, from Oct 2020]
 Mehran Hatamzadeh [Univ Côte d'Azur, from Nov 2020]
 Rima Khouja [Univ de Lorraine]
 Thomas Laporte [Univ Côte d'Azur]
 Michelangelo Marsala [Inria, from Nov 2020]
 Tobias Metzlaff [Inria]
 Erick David Rodriguez Bazan [Inria]
 Fatma Yildirim [Univ de Nice  Sophia Antipolis, until Sep 2020]
Technical Staff
 Thuy Linh Nguyen [ERCIM, Engineer, until Jan 2020]
Administrative Assistant
 Sophie Honnorat [Inria]
External Collaborator
 Adam Parusinski [Univ Côte d'Azur, from Mar 2020, HDR]
2 Overall objectives
Our daily life environment is increasingly interacting with digital information. An important amount of this information is of geometric nature. It concerns the representation of our environment, the analysis and understanding of “real” phenomena, the control of physical mechanisms or processes. The interaction between physical and digital worlds is twoway. Sensors are producing digital data related to measurements or observations of our environment. Digital models are also used to “act” on the physical world. Objects that we use at home, at work, to travel, such as furniture, cars, planes, ... are nowadays produced by industrial processes which are based on digital representation of shapes. CADCAM (Computer Aided Design – Computer Aided Manufacturing) software is used to represent the geometry of these objects and to control the manufacturing processes which create them. The construction capabilities themselves are also expanding, with the development of 3D printers and the possibility to create dailylife objects “at home” from digital models.
The impact of geometry is also important in the analysis and understanding of phenomena. The 3D conformation of a molecule explains its biological interaction with other molecules. The profile of a wing determines its aeronautic behavior, while the shape of a bulbous bow can decrease significantly the wave resistance of a ship. Understanding such a behavior or analyzing a physical phenomenon can nowadays be achieved for many problems by numerical simulation. The precise representation of the geometry and the link between the geometric models and the numerical computation tools are closely related to the quality of these simulations. This also plays an important role in optimisation loops where the numerical simulation results are used to improve the “performance” of a model.
Geometry deals with structured and efficient representations of information and with methods to treat it. Its impact in animation, games and VAMR (Virtual, Augmented and Mixed Reality) is important. It also has a growing influence in etrade where a consumer can evaluate, test and buy a product from its digital description. Geometric data produced for instance by 3D scanners and reconstructed models are nowadays used to memorize old works in cultural or industrial domains.
Geometry is involved in many domains (manufacturing, simulation, communication, virtual world...), raising many challenging questions related to the representations of shapes, to the analysis of their properties and to the computation with these models. The stakes are multiple: the accuracy in numerical engineering, in simulation, in optimization, the quality in design and manufacturing processes, the capacity of modeling and analysis of physical problems.
3 Research program
3.1 High order geometric modeling
The accurate description of shapes is a long standing problem in mathematics, with an important impact in many domains, inducing strong interactions between geometry and computation. Developing precise geometric modeling techniques is a critical issue in CADCAM. Constructing accurate models, that can be exploited in geometric applications, from digital data produced by cameras, laser scanners, observations or simulations is also a major issue in geometry processing. A main challenge is to construct models that can capture the geometry of complex shapes, using few parameters while being precise.
Our first objective is to develop methods, which are able to describe accurately and in an efficient way, objects or phenomena of geometric nature, using algebraic representations.
The approach followed in CAGD, to describe complex geometry is based on parametric representations called NURBS (Non Uniform Rational BSpline). The models are constructed by trimming and gluing together high order patches of algebraic surfaces. These models are built from the socalled BSpline functions that encode a piecewise algebraic function with a prescribed regularity at knots. Although these models have many advantages and have become the standard for designing nowadays CAD models, they also have important drawbacks. Among them, the difficulty to locally refine a NURBS surface and also the topological rigidity of NURBS patches that imposes to use many such patches with trims for designing complex models, with the consequence of the appearing of cracks at the seams. To overcome these difficulties, an active area of research is to look for new blending functions for the representation of CAD models. Some examples are the socalled TSplines, LRSpline blending functions, or hierarchical splines, that have been recently devised in order to perform efficiently local refinement. An important problem is to analyze spline spaces associated to general subdivisions, which is of particular interest in higher order Finite Element Methods. Another challenge in geometric modeling is the efficient representation and/or reconstruction of complex objects, and the description of computational domains in numerical simulation. To construct models that can represent efficiently the geometry of complex shapes, we are interested in developing modeling methods, based on alternative constructions such as skeletonbased representations. The change of representation, in particular between parametric and implicit representations, is of particular interest in geometric computations and in its applications in CAGD.
We also plan to investigate adaptive hierarchical techniques, which can locally improve the approximation of a shape or a function. They shall be exploited to transform digital data produced by cameras, laser scanners, observations or simulations into accurate and structured algebraic models.
The precise and efficient representation of shapes also leads to the problem of extracting and exploiting characteristic properties of shapes such as symmetry, which is very frequent in geometry. Reflecting the symmetry of the intended shape in the representation appears as a natural requirement for visual quality, but also as a possible source of sparsity of the representation. Recognizing, encoding and exploiting symmetry requires new paradigms of representation and further algebraic developments. Algebraic foundations for the exploitation of symmetry in the context of non linear differential and polynomial equations are addressed. The intent is to bring this expertise with symmetry to the geometric models and computations developed by aromath.
3.2 Robust algebraicgeometric computation
In many problems, digital data are approximated and cannot just be used as if they were exact. In the context of geometric modeling, polynomial equations appear naturally as a way to describe constraints between the unknown variables of a problem. An important challenge is to take into account the input error in order to develop robust methods for solving these algebraic constraints. Robustness means that a small perturbation of the input should produce a controlled variation of the output, that is forward stability, when the inputoutput map is regular. In nonregular cases, robustness also means that the output is an exact solution, or the most coherent solution, of a problem with input data in a given neighborhood, that is backward stability.
Our second long term objective is to develop methods to robustly and efficiently solve algebraic problems that occur in geometric modeling.
Robustness is a major issue in geometric modeling and algebraic computation. Classical methods in computer algebra, based on the paradigm of exact computation, cannot be applied directly in this context. They are not designed for stability against input perturbations. New investigations are needed to develop methods which integrate this additional dimension of the problem. Several approaches are investigated to tackle these difficulties.
One relies on linearization of algebraic problems based on “elimination of variables” or projection into a space of smaller dimension. Resultant theory provides a strong foundation for these methods, connecting the geometric properties of the solutions with explicit linear algebra on polynomial vector spaces, for families of polynomial systems (e.g., homogeneous, multihomogeneous, sparse). Important progress has been made in the last two decades to extend this theory to new families of problems with specific geometric properties. Additional advances have been achieved more recently to exploit the syzygies between the input equations. This approach provides matrix based representations, which are particularly powerful for approximate geometric computation on parametrized curves and surfaces. They are tuned to certain classes of problems and an important issue is to detect and analyze degeneracies and to adapt them to these cases.
A more adaptive approach involves linear algebra computation in a hierarchy of polynomial vector spaces. It produces a description of quotient algebra structures, from which the solutions of polynomial systems can be recovered. This family of methods includes Gröbner Basis , which provides general tools for solving polynomial equations. Border Basis is an alternative approach, offering numerically stable methods for solving polynomial equations with approximate coefficients . An important issue is to understand and control the numerical behavior of these methods as well as their complexity and to exploit the structure of the input system.
In order to compute “only” the (real) solutions of a polynomial system in a given domain, duality techniques can also be employed. They consist in analyzing and adding constraints on the space of linear forms which vanish on the polynomial equations. Combined with semidefinite programming techniques, they provide efficient methods to compute the real solutions of algebraic equations or to solve polynomial optimization problems. The main issues are the completness of the approach, their scalability with the degree and dimension and the certification of bounds.
Singular solutions of polynomial systems can be analyzed by computing differentials, which vanish at these points. This leads to efficient deflation techniques, which transform a singular solution of a given problem into a regular solution of the transformed problem. These local methods need to be combined with more global root localisation methods.
Subdivision methods are another type of methods which are interesting for robust geometric computation. They are based on exclusion tests which certify that no solution exists in a domain and inclusion tests, which certify the uniqueness of a solution in a domain. They have shown their strength in addressing many algebraic problems, such as isolating real roots of polynomial equations or computing the topology of algebraic curves and surfaces. The main issues in these approaches is to deal with singularities and degenerate solutions.
4 Application domains
4.1 Geometric modeling for Design and Manufacturing.
The main domain of applications that we consider for the methods we develop is Computer Aided Design and Manufacturing.
ComputerAided Design (CAD) involves creating digital models defined by mathematical constructions, from geometric, functional or aesthetic considerations. Computeraided manufacturing (CAM) uses the geometrical design data to control the tools and processes, which lead to the production of real objects from their numerical descriptions.
CADCAM systems provide tools for visualizing, understanding, manipulating, and editing virtual shapes. They are extensively used in many applications, including automotive, shipbuilding, aerospace industries, industrial and architectural design, prosthetics, and many more. They are also widely used to produce computer animation for special effects in movies, advertising and technical manuals, or for digital content creation. Their economic importance is enormous. Their importance in education is also growing, as they are more and more used in schools and educational purposes.
CADCAM has been a major driving force for research developments in geometric modeling, which leads to very large software, produced and sold by big companies, capable of assisting engineers in all the steps from design to manufacturing.
Nevertheless, many challenges still need to be addressed. Many problems remain open, related to the use of efficient shape representations, of geometric models specific to some application domains, such as in architecture, naval engineering, mechanical constructions, manufacturing ...Important questions on the robustness and the certification of geometric computation are not yet answered. The complexity of the models which are used nowadays also appeals for the development of new approaches. The manufacturing environment is also increasingly complex, with new type of machine tools including: turning, 5axes machining and wire EDM (Electrical Discharge Machining), 3D printer. It cannot be properly used without computer assistance, which raises methodological and algorithmic questions. There is an increasing need to combine design and simulation, for analyzing the physical behavior of a model and for optimal design.
The field has deeply changed over the last decades, with the emergence of new geometric modeling tools built on dedicated packages, which are mixing different scientific areas to address specific applications. It is providing new opportunities to apply new geometric modeling methods, output from research activities.
4.2 Geometric modeling for Numerical Simulation and Optimization
A major bottleneck in the CADCAM developments is the lack of interoperability of modeling systems and simulation systems. This is strongly influenced by their development history, as they have been following different paths.
The geometric tools have evolved from supporting a limited number of tasks at separate stages in product development and manufacturing, to being essential in all phases from initial design through manufacturing.
Current Finite Element Analysis (FEA) technology was already well established 40 years ago, when CADsystems just started to appear, and its success stems from using approximations of both the geometry and the analysis model with low order finite elements (most often of degree $\le 2$).
There has been no requirement between CAD and numerical simulation, based on Finite Element Analysis, leading to incompatible mathematical representations in CAD and FEA. This incompatibility makes interoperability of CAD/CAM and FEA very challenging. In the general case today this challenge is addressed by expensive and timeconsuming human intervention and software developments.
Improving this interaction by using adequate geometric and functional descriptions should boost the interaction between numerical analysis and geometric modeling, with important implications in shape optimization. In particular, it could provide a better feedback of numerical simulations on the geometric model in a design optimization loop, which incorporates iterative analysis steps.
The situation is evolving. In the past decade, a new paradigm has emerged to replace the traditional Finite Elements by BSpline basis element of any polynomial degree, thus in principle enabling exact representation of all shapes that can be modeled in CAD. It has been demonstrated that the socalled isogeometric analysis approach can be far more accurate than traditional FEA.
It opens new perspectives for the interoperability between geometric modeling and numerical simulation. The development of numerical methods of high order using a precise description of the shapes raises questions on piecewise polynomial elements, on the description of computational domains and of their interfaces, on the construction of good function spaces to approximate physical solutions. All these problems involve geometric considerations and are closely related to the theory of splines and to the geometric methods we are investigating. We plan to apply our work to the development of new interactions between geometric modeling and numerical solvers.
5 Highlights of the year
Fatmanur Yildirim defended her PhD thesis 40 titled "Finite fibers of rational maps by means of matrix representations with applications to distance problem" at University Côte d'Azur on February 3, 2020.
Ahmed Blidia defended his PhD thesis 38 titled "Geometric continuity for surfaces and scalar fields" at University Côte d'Azur on April 8, 2020. He was supervized by Bernard Mourrain and the committee included Ioannis Emiris.
Clément Laroche defended his PhD thesis 39 titled "Compact and efficient implicit representations" ("Représentations implicites compactes et efficaces") at University of Athens on April 30, 2020. He was supervized by Ioannis Emiris and the examination committee included Bernard Mourrain.
Ioannis Emiris was cogeneral chair of International Symposium on Symbolic and Algebraic Computation (ISSAC 2020), held online on July 2023, 2020.
6 New results
6.1 Elimination ideals and Bézout relations
Participant: André Galligo.
Let $k$ be an infinite field and $I\subset k[{x}_{1},...,{x}_{n}]$ be a nonzero ideal such that $dimV\left(I\right)=q\ge 0$. Denote by $({f}_{1},...,{f}_{s})$ a set of generators of $I$. One can see that in the set $I\cap k[{x}_{1},...,{x}_{q+1}]$ there exist nonzero polynomials, depending only on these $q+1$ variables. We aim to bound the minimal degree of the polynomials of this type, and of a Bézout (i.e. membership) relation expressing such a polynomial as a combination of the ${f}_{i}$. In particular we show in 27 that if $deg{f}_{i}={d}_{i}$ where ${d}_{1}\ge {d}_{2}\ge \cdots \ge {d}_{s}$, then there exist a nonzero polynomial $\varphi \left(x\right)\in I\cap k[{x}_{1},...,{x}_{q+1}]$, such that $deg\left(\varphi \right)\le {d}_{s}{\prod}_{i=1}^{nq1}{d}_{i}$. We also give a relative version of this theorem.This is a joint work with Zbigniew Jelonek, Instytut Matematyczny PAN, Warszawa, Poland.
6.2 Separation bounds for polynomial systems
Participant: Ioannis Emiris, Bernard Mourrain, Elias Tsigaridas (OURAGAN).
We rely on aggregate separation bounds for univariate polynomials to introduce, in 24, novel worstcase separation bounds for the isolated roots of zerodimensional, positivedimensional, and overdetermined polynomial systems. We exploit the structure of the given system, as well as bounds on the height of the sparse (or toric) resultant, by means of mixed volume, thus establishing adaptive bounds. Our bounds improve upon Canny’s Gap theorem. Moreover, they exploit sparseness and they apply without any assumptions on the input polynomial system. To evaluate the quality of the bounds, we present polynomial systems whose root separation is asymptotically not far from our bounds. We apply our bounds to three problems. First, we use them to estimate the bitsize of the eigenvalues and eigenvectors of an integer matrix; thus we provide a new proof that the problem has polynomial bit complexity. Second, we bound the value of a positive polynomial over the simplex: we improve by at least one order of magnitude upon all existing bounds. Finally, we asymptotically bound the number of steps of any purely subdivisionbased algorithm that isolates all real roots of a polynomial system.
6.3 Matrix formulae for Resultants and Discriminants of Bivariate Tensorproduct Polynomials
Participant: Laurent Busé, Angelos Mantzaflaris, Elias Tsigaridas (OURAGAN).
The construction of optimal resultant formulae for polynomial systems is one of the main areas of research in computational algebraic geometry. However, most of the constructions are restricted to formulae for unmixed polynomial systems, that is, systems of polynomials which all have the same support. Such a condition is restrictive, since mixed systems of equations arise frequently in many problems. Nevertheless, resultant formulae for mixed polynomial systems is a very challenging problem. In 21 we present a square, Koszultype, matrix, the determinant of which is the resultant of an arbitrary (mixed) bivariate tensorproduct polynomial system. The formula generalizes the classical Sylvester matrix of two univariate polynomials, since it expresses a map of degree one, that is, the elements of the corresponding matrix are up to sign the coefficients of the input polynomials. Interestingly, the matrix expresses a primaldual multiplication map, that is, the tensor product of a univariate multiplication map with a map expressing derivation in a dual space. In addition we prove an impossibility result which states that for tensorproduct systems with more than two (affine) variables there are no universal degreeone formulae, unless the system is unmixed. Last but not least, we present applications of the new construction in the efficient computation of discriminants and mixed discriminants.
6.4 Fibers of multigraded rational maps and orthogonal projection onto rational surfaces
Participant: Laurent Busé, Fatmanur Yildirim.
In 16 we contribute a new algebraic method for computing the orthogonal projections of a point onto a rational algebraic surface embedded in the three dimensional projective space. This problem is first turned into the computation of the finite fibers of a generically finite dominant rational map: a congruence of normal lines to the rational surface. Then, an indepth study of certain syzygy modules associated to such a congruence is presented and applied to build elimination matrices that provide universal representations of its finite fibers, under some genericity assumptions. These matrices depend linearly in the variables of the three dimensional space. They can be precomputed so that the orthogonal projections of points are approximately computed by means of fast and robust numerical linear algebra calculations.
This is a joint work with Nicolas Notbol, the University of Buenos Aires, Argentina and Marc Chardin, CNRS and Institut de Mathématiques de Jussieu, Sorbonne Universités, Paris, France.
6.5 Degree and birationality of multigraded rational maps
Participant: Laurent Busé.
In 18 we give formulas and effective sharp bounds for the degree of multigraded rational maps and provide some effective and computable criteria for birationality in terms of their algebraic and geometric properties. We also extend the Jacobian dual criterion to the multigraded setting. Our approach is based on the study of blowup algebras, including syzygies, of the ideal generated by the defining polynomials of the rational map. A key ingredient is a new algebra that we call the saturated special fiber ring, which turns out to be a fundamental tool to analyze the degree of a rational map. We also provide a very effective birationality criterion and a complete description of the equations of the associated Rees algebra of a particular class of plane rational maps.
This is a joint work Carlos D'Andrea and Yairon Cid Ruiz from the University of Barcelona, Spain.
6.6 Freeness and invariants of rational plane curves
Participant: Laurent Busé.
Given a parameterization φ of a rational plane curve C, in 20 we study some invariants of C via φ. We first focus on the characterization of rational cuspidal curves, in particular we establish a relation between the discriminant of the pullback of a line via φ, the dual curve of C and its singular points. Then, by analyzing the pullbacks of the global differential forms via φ, we prove that the (nearly) freeness of a rational curve can be tested by inspecting the Hilbert function of the kernel of a canonical map. As a by product, we also show that the global Tjurina number of a rational curve can be computed directly from one of its parameterization, without relying on the computation of an equation of C.
This is a joint work Alexandru Dimca from the University Côté d'Azur, France, and Gabriel Sticlaru from the Faculty of Mathematics and Informatics of Constanta, Romania.
6.7 The geometry of the flex locus of a hypersurface
Participant: Laurent Busé.
In 19 we give a formula in terms of multidimensional resultants for an equation for the flex locus of a projective hypersurface, generalizing a classical result of Salmon for surfaces in P3. Using this formula, we compute the dimension of this flex locus, and an upper bound for the degree of its defining equations. We also show that, when the hypersurface is generic, this bound is reached, and that the generic flex line is unique and has the expected order of contact with the hypersurface.
This is a joint work Carlos D'Andrea and Martin Sombra from the University of Barcelona, Spain, and Martin Weimann from the University of Caen, France.
6.8 Regularity of bicyclic graphs and their powers
Participant: Navid Nemati.
In 22 we characterise the CastelnuovoMumford regularity of edge ideal $I\left(G\right)$ of a bicyclic graph $G$ in terms of the induced matching numbers of $G$. For the base case of this family of graphs, i.e. dumbbell graphs, we explicitly compute the induced matching number. Moreover, we prove that $\mathrm{reg}I{\left(G\right)}^{q}=2q+\mathrm{reg}I\left(G\right)2$, for all $q\ge 1$, when $G$ is a dumbbell graph with a connecting path having no more than two vertices.
This is a joint work with Yairon CidRuiz from Universitat de Barcelona, Spain, Sepehr Jafari from University Degli Studi di Genova, and Beatrice Picone from University of Catania, Italy.
6.9 Lefschetz properties of monomial algebras with almost linear resolution
Participant: Navid Nemati.
In 22 we study the WLP and SLP of artinian monomial ideals in $S=k[{x}_{1},\cdots ,{x}_{n}]$ via studying their minimal free resolutions. We study the Lefschetz properties of such ideals where the minimal free resolution of $S/I$ is linear for at least $n2$ steps. We give an affirmative answer to a conjecture of Eisenbud, Huneke and Ulrich for artinian monomial ideals with almost linear resolutions.
This is a joint work with Nasrin Altafi from KTH Royal Institute of Technology, Sweden.
6.10 Punctual Hilbert Schemes and Certified Approximate Singularities
Participant: Angelos Mantzaflaris, Bernard Mourrain.
In 35 we provide a new method to certify that a nearby polynomial system has a singular isolated root with a prescribed multiplicity structure. More precisely, given a polynomial system $f=({f}_{1},\cdots ,{f}_{N})\in \u2102{[{x}_{1},...,{x}_{n}]}^{N}$, we present a Newton iteration on an extended deflated system that locally converges, under regularity conditions, to a small deformation of $f$ such that this deformed system has an exact singular root. The iteration simultaneously converges to the coordinates of the singular root and the coefficients of the so called inverse system that describes the multiplicity structure at the root. We use $\alpha $theory test to certify the quadratic convergence, and to give bounds on the size of the deformation and on the approximation error. The approach relies on an analysis of the punctual Hilbert scheme, for which we provide a new description. We show in particular that some of its strata can be rationally parametrized and exploit these parametrizations in the certification. We show in numerical experimentation how the approximate inverse system can be computed as a starting point of the Newton iterations and the fast numerical convergence to the singular root with its multiplicity structure, certified by our criteria.
This is a joint work with Agnes Szanto, Department of mathematics, North Carolina State University, USA.
6.11 Ideal Interpolation, HBases and Symmetry
Participant: Evelyne Hubert, Erick Rodriguez Bazan.
Multivariate Lagrange and Hermite interpolation are examples of ideal interpolation. The interpolation problem is defined by a set of linear forms, on the polynomial ring, whose kernels intersect into an ideal. Another class of examples thus arises when the linear forms are defined by the normal forms modulo an ideal. For an ideal interpolation problem with symmetry, we address in 36 the simultaneous computation of a symmetry adapted basis of the least interpolation space and the related symmetry adapted Hbasis. The algorithm exploits symmetry at all stages, in addition to preserving it in the output. Defined by de Boor and Ron (1990), the least interpolation space is a canonically defined representation of the quotient of the polynomial ring by the ideal. As a result it has great properties that cannot be achieved by a space spanned by monomials. In particular it is invariant under a group action when the ideal is; Equivalently, when the set of linear forms is. Its natural complement is a Hbasis of the ideal, as introduced by Macaulay (1916). As a contrast with Groebner bases, which are tied to a monomial term order and are notorious for not preserving symmetry, we show that symmetry can be naturally preserved and represented in Hbases, as well as exploited in the course of their computation.
6.12 Computing Free Noncommutative Groebner Bases over Z with Singular:Letterplace
Participant: Tobias Metzlaff.
The extension of Groebner bases concept from polynomial algebras over fields to polynomial rings over rings allows to tackle numerous applications, both of theoretical and of practical importance. Groebner and GroebnerShirshov bases can be defined for various noncommutative and even nonassociative algebraic structures. In 34 we study the case of associative rings and aim at free algebras over principal ideal rings. We concentrate ourselves on the case of commutative coefficient ring without zero divisors (i. e. a domain). Even working over Z allows one to do computations, which can be treated as universal for fields of arbitrary characteristic. By using the systematic approach, we revisit the theory and present the algorithms in the implementable form. We show drastic differences in the behavior of Gröbner bases between free algebras and algebras, close to commutative. Even the formation of critical pairs has to be reengineered, together with the criteria for their quick discarding. We present an implementation of algorithms in the Singular subsystem called Letterplace, which internally uses Letterplace techniques (and Letterplace Gröbner bases), due to La Scala and Levandovskyy. Interesting examples accompany our presentation.This is joint work with Viktor Levandovskyy and Karim Zeid, RWTH Aachen University, Germany.
6.13 On minimal decompositions of low rank symmetric tensors
Participant: Bernard Mourrain.
In 30, we use an algebraic approach to construct minimal decompositions of symmetric tensors with low rank. This is done by using Apolarity Theory and by studying minimal sets of reduced points apolar to a given symmetric tensor, namely, whose ideal is contained in the apolar ideal associated to the tensor. In particular, we focus on the structure of the Hilbert function of these ideals of points. We give a procedure which produces a minimal set of points apolar to any symmetric tensor of rank at most 5. This procedure is also implemented in the algebra software Macaulay2.This is a joint work with Alessandro Oneto, Univ. de Trento, Italy.
6.14 Computational geometry in general dimension
Participant: Ioannis Emiris, Ioannis Psarros.
In 25 we examine general products of Euclidean metrics and construct a general framework in order to study such metrics and design nearest neighbor (NN) query algorithms for them. Then, we apply our approach to the Frechet distance metric and the Dynamic time warping distance function and show they are both special cases of our framework. Our methods are applied to proximity problems and approximate NN queries among generaldimensional polygonal curves, thus improving upon the existing asymptotic complexity of such operations.Participant: Ioannis Emiris, Ioannis Psarros.
In 13 we design efficient algorithms for highdimensional approximate $r$nets that improve upon the asymptotic complexity of existing approaches. The problem asks to find a subset of a given pointset that shall define the centers of balls of radius $r$ which shall cover the input points. We expect the simplicity of our algorithms to allow for a practical implementation of our techniques.This is joint work with Zeta Avarikioti (ETH Zurich), and Loukas Kavouras (Athena RC).
Participant: Ioannis Emiris, Evangelos Bartzos.
In 14 we employ the multihomogeneous Bézout (mBezout) number in order to derive new asymptotic bounds on the number of Euclidean embeddings of generically minimally rigid graphs. These are distance graphs that admit a finite number of realizations and generalize the class of Laman graphs that are rigid in the plane. Despite a number of studies on this problem, our bounds are the first to improve upon the trivial Bézout bound for dimension larger than 4 as well as 3D graphs that are graphtheoretically planar. The crucial observation is that the mBezout number is given by matrix permanents that have a strong combinatorial structure which can be revealed and exploited.This is joint work with Josef Schicho (U Linz).
6.15 Complete Classification and Efficient Determination of Arrangements Formed by Two Ellipsoids
Participant: Bernard Mourrain.
Arrangements of geometric objects refer to the spatial partitions formed by the objects and they serve as an underlining structure of motion design, analysis and planning in CAD/CAM, robotics, molecular modeling, manufacturing and computerassisted radiosurgery. Arrangements are especially useful to collision detection, which is a key task in various applications such as computer animation , virtual reality, computer games, robotics, CAD/CAM and computational physics. Ellipsoids are commonly used as bounding volumes in approximating complex geometric objects in collision detection. In 28, we present an indepth study on the arrangements formed by two ellipsoids. Specifically, we present a classification of these arrangements and propose an efficient algorithm for determining the arrangement formed by any particular pair of ellipsoids. A stratification diagram is also established to show the connections among all the arrangements formed by two ellipsoids. Our results for the first time elucidate all possible relative positions between two arbitrary ellipsoids and provides an efficient and robust algorithm for determining the relative position of any two given ellipsoids, therefore providing the necessary foundation for developing practical and trustworthy methods for processing ellipsoids for collision analysis or simulation in various applications.
This is a joint work with Xiaohong Jia, KLMM  Key Laboratory of Mathematics Mechanization, Beijing, China; Changhe Tu, School of Computer Science and Technology, Qingdao, China; Wenping Wang, HKU  The University of Hong Kong, China.
6.16 Interpolatory CatmullClark volumetric subdivision over unstructured hexahedral meshes for modeling and simulation applications
Participant: Bernard Mourrain.
Volumetric modeling is an important topic for material modeling and isogeometric simulation. In 33, two kinds of interpolatory CatmullClark volumetric subdivision approaches over unstructured hexahedral meshes are proposed based on the limit point formula of CatmullClark subdivision volume. The basic idea of the first method is to construct a new control lattice, whose limit volume by the CatmullClark subdivision scheme interpolates vertices of the original hexahedral mesh. The new control lattice is derived by the local pushback operation from one CatmullClark subdivision step with modified geometric rules. This interpolating method is simple and efficient, and several shape parameters are involved in adjusting the shape of the limit volume. The second method is based on progressiveiterative approximation using limit point formula. At each iteration step, we progressively modify vertices of an original hexahedral mesh to generate a new control lattice whose limit volume interpolates all vertices in the original hexahedral mesh. The convergence proof of the iterative process is also given. The interpolatory subdivision volume has C${}^{2}$smoothness at the regular region except around extraordinary vertices and edges. Furthermore, the proposed interpolatory volumetric subdivision methods can be used not only for geometry interpolation, but also for material attribute interpolation in the field of volumetric material modeling. The application of the proposed volumetric subdivision approaches on isogeometric analysis is also given with several examples.
This is a joint work with Jin Xie, Jinlan Xu, Zhenyu Dong, Gang Xu, Chongyang Deng, HDU  College of Computer Science, Hangzhou, China; Yongjie Jessica Zhang, Department of Mechanics, Carnegie Mellon University, USA.
6.17 Geometrically smooth spline bases for data fitting and simulation
Participant: Ahmed Blidia, Bernard Mourrain.
Given a topological complex $M$ with glueing data along edges shared by adjacent faces, we study the associated space of geometrically smooth spline functions that satisfy differentiability properties across shared edges. In 15, we present new and efficient constructions of basis functions of the space of G${}^{1}$spline functions on quadrangular meshes, which are tensor product bspline functions on each quadrangle and with bspline transition maps across the shared edges. This new strategy for constructing basis functions is based on a local analysis of the edge functions, and does not depend on the global topology of $M$. We show that the separability of the space of G${}^{1}$ splines across an edge allows to determine the dimension and a basis of the space of G${}^{1}$ splines on $M$. This leads to explicit and effective constructions of basis functions attached to the vertices, edges and faces of $M$. This basis construction has important applications in geometric modeling and simulation. We illustrate it by the fitting of point clouds by G${}^{1}$ splines on quadrangular meshes of complex topology and in Isogeometric Analysis methods for the solution of diffusion equations. The ingredients are detailed and experimentation results showing the behavior of the method are presented.This is joint work with Gang Xu, HDU  College of Computer Science, Hangzhou, China.
6.18 An overview of Geometry plus Simulation Modules
Participant: Angelos Mantzaflaris.
In 37 we give an overview of the opensource library "G+Smo". G+Smo is a C++ library that brings together mathematical tools for geometric design and numerical simulation. It implements the relatively new paradigm of isogeometric analysis, which suggests the use of a unified framework in the design and analysis pipeline. G+Smo is an objectoriented, crossplatform, fully templated library and follows the generic programming principle, with a focus on both efficiency and ease of use. The library aims at providing access to high quality, opensource software to the community of numerical simulation and beyond.6.19 Local (T)HBspline projectors via restricted hierarchical spline fitting
Participant: Angelos Mantzaflaris.
The paper 26 is devoted to techniques for adaptive spline projection via quasiinterpolation, enabling the efficient approximation of given functions. We employ local leastsquares fitting in restricted hierarchical spline spaces to establish novel projection operators for hierarchical splines of degree p. This leads to efficient spline projectors that require $O\left({p}^{d}\right)$ floating point operations and $O\left(1\right)$ evaluations of the given function per degree of freedom, while providing essentially the same accuracy as global approximation. Our spline projectors are based on a unifying framework for quasiinterpolation in hierarchical spline spaces. We present a detailed comparison with the scheme of Speleers and Manni (2016).This is a joint work with Alessandro Giust and Bert Jüttler for the Johannes Kepler University of Linz, Austria.
6.20 Efficient Matrix Assembly in Isogeometric Analysis with Hierarchical Bsplines
Participant: Angelos Mantzaflaris.
Hierarchical Bsplines that allow local refinement have become a promising tool for developing adaptive isogeometric methods. Unfortunately, similar to tensorproduct Bsplines, the computational cost required for assembling the system matrices in isogeometric analysis with hierarchical Bsplines is also high, particularly if the spline degree is increased. To address this issue, we propose in 31 an efficient matrix assembly approach for bivariate hierarchical Bsplines based the previously proposed IIL approach. The new algorithm consists of three stages: approximating the integrals by quasiinterpolation, building three compact lookup tables and assembling the matrices via sumfactorization. A detailed analysis shows that the complexity of our method has the order $O\left(N{p}^{3}\right)$ under a mild assumption about mesh admissibility, where $N$ and $p$ denote the number of degrees of freedom and spline degree respectively. Finally, several experimental results are demonstrated to verify the theoretical results and to show the performance of the proposed method.This is a joint work with Maodong Pan and Bert Jüttler for the Johannes Kepler University of Linz, Austria.
7 Bilateral contracts and grants with industry
7.1 Bilateral contracts with industry
$\u2022$ NURBSFIX: Repairing the topology of a NURBS model in view of its approximation. We have a research contract with the industrial partner GeometryFactory, in collaboration with the projectteam Titane (Pierre Alliez). The postdoc of Xiao Xiao is funded by this research contract together with a PEPS from the labex AMIES. We continue to develop a robust algorithm which offers the capability to repair 155 while meshing CAD models, without requiring a valid surface mesh as input, and 156 to generate volume meshes which are valid by design and 157 independent of the input NURBS patch layout, without resorting 158 to postprocessing steps such as mesh quilting or remeshing.
$\u2022$ Geometric computing. Ioannis Emiris coordinates a research contract with the industrial partner ANSYS (Greece), in collaboration with Athena Research Center. PhD candidate A. Chalkis is partially funded.
Electronic design automation (EDA) and simulating Integrated Circuits requires robust geometric operations on thousands of electronic elements (capacitors, resistors, coils etc) represented by polyhedral objects in 2.5 dimensions, not necessarily convex. A special case may concern axisaligned objects but the real challenge is the general case. The project, extended into 2021, focuses on 3 axes: (1) efficient data structures and prototype implementations for storing the aforementioned polyhedral objects so that nearest neighbor queries are fast in the Lmax metric, which is the primary focus of the contract, (2) random sampling of the free space among objects, (3) datadriven algorithmic design for problems concerning datastructures and their construction and initialization.
7.2 Bilateral grants with industry
$\u2022$ Interactive construction of 3D models  Application to the modeling of complex geological structures.
CIFRE collaboration between Schlumberger Montpellier (A. Azzedine) and Inria Sophia Antipolis (B. Mourrain). The PhD candidate is A. Belhachmi. The objective of the work is the development of a new spline based high quality geomodeler for reconstructing the stratigraphy of geological layers from the adaptive and efficient processing of large terrain information.
8 Partnerships and cooperations
8.1 International initiatives
High order methods for computational engineering and data analysis
 Title: High order methods for computational engineering and data analysis
 Program: Partenariats Hubert Curien (PHC) Alliance
 Duration: January 2020  December 2022
 Other partner: U. Swansea (UK)
 Inria contact: Angelos Mantzaflaris
 Summary: The aim of this project is to develop a mathematical framework for the integration of geometric modeling and simulation using splinebased finite elements of high degree of smoothness. Highorder methods are known to provide a robust and efficient methodology to tackle complex challenges in multiphysics simulations, shape optimization, and the analysis of largescale datasets arising in datadriven engineering and design. However, the analysis and design of highorder methods is a daunting task requiring a concurrent effort from diverse fields such as applied algebraic geometry, approximation theory and splines, topological data analysis, and computational mathematics. Our strategic vision is to create a research team combining uniquely broad research expertise in these areas by establishing a link between the AROMATH and Swansea University.
8.2 European initiatives
8.2.1 FP7 & H2020 Projects
ARCADES
 Title: Algebraic Representations for ComputerAided Design of Complex Shapes
 Project acronym: ARCADES
 Program: Marie SkłodowskaCurie ETN
 Duration: January 2016  December 2019
 Coordinator: I.Z. Emiris, ATHENA Research Innovation Center

Other partners:
 U. Barcelona (Spain),
 J. Kepler University, Linz (Austria),
 SINTEF Institute, Oslo (Norway),
 U. Strathclyde, Glascow (UK),
 Technische U. Wien (Austria),
 Evolute GmBH, Vienna (Austria).
 Inria contact: Laurent Busé

Webpage:
http://
arcadesnetwork. eu/  Summary: ARCADES aims at disrupting the traditional paradigm in ComputerAided Design (CAD) by exploiting cuttingedge research in mathematics and algorithm design. Geometry is now a critical tool in a large number of key applications; somewhat surprisingly, however, several approaches of the CAD industry are outdated, and 3D geometry processing is becoming increasingly the weak link. This is alarming in sectors where CAD faces new challenges arising from fast point acquisition, big data, and mobile computing, but also in robotics, simulation, animation, fabrication and manufacturing, where CAD strives to address crucial societal and market needs. The challenge taken up by ARCADES is to invert the trend of CAD industry lagging behind mathematical breakthroughs and to build the next generation of CAD software based on strong foundations from algebraic geometry, differential geometry, scientific computing, and algorithm design. Our gamechanging methods lead to realtime modelers for architectural geometry and visualisation, to isogeometric and designthroughanalysis software for shape optimisation, and marine design and hydrodynamics, and to tools for motion design, robot kinematics, path planning, and control of machining tools.
POEMA
 Title: Polynomial Optimization, Efficiency through Moments and Algebra
 Project acronym: POEMA
 Program: Marie SkłodowskaCurie ITN
 Duration: January 2019  December 2022 (48 months)
 Coordinator: Bernard Mourrain, Aromath, Inria Sophia Antipolis

Other partners:
 LAAS  CNRS, Toulouse (France),
 Sorbonne Université, Paris (France),
 Centrum Wiskunde & Informatica, Amsterdam (The Netherlands),
 Stichting Katholieke Universiteit Brabant, Tilburd (The Netherlands),
 Universität Konstanz (Germany),
 Università degli Studi di Firenze (Italy),
 University of Birmingham (United Kingdom),
 Friedrich Alexander University ErlangenNuremberg (Germany),
 Universitet I Tromsø (Norway),
 ARTELYS SAS, Paris (France).
 Inria contact: Bernard Mourrain

Webpage:
http://
poemanetwork. eu/ 
Summary:
Nonlinear optimization problems are present in many reallife applications and in scientific areas such as operations research, control engineering, physics, information processing, economy, biology, etc. However, efficient computational procedures, that can provide the guaranteed global optimum, are lacking for them. The project will develop new polynomial optimization methods, combining moment relaxation procedures with computational algebraic tools to address this type of problems. Recent advances in mathematical programming have shown that the polynomial optimization problems can be approximated by sequences of SemiDefinite Programming problems. This approach provides a powerful way to compute global solutions of nonlinear optimization problems and to guarantee the quality of computational results. On the other hand, advanced algebraic algorithms to compute all the solutions of polynomial systems, with efficient implementations for exact and approximate solutions, were developed in the past twenty years. The network combines the expertise of active European teams working in these two domains to address important challenges in polynomial optimization and to show the impact of this research on practical applications.
POEMA aims to train scientists at the interplay of algebra, geometry and computer science for polynomial optimization problems and to foster scientific and technological advances, stimulating interdisciplinary and intersectoriality knowledge exchange between algebraists, geometers, computer scientists and industrial actors facing reallife optimization problems.
GRAPES
 Title: learninG, Representing, And oPtimizing shapES.
 Duration: 2019  2023
 Coordinator: ATHENA

Partners:
 ATHENA, ATHINAEREVNITIKO KENTRO KAINOTOMIAS STIS TECHNOLOGIES TIS PLIROFORIAS, TON EPIKOINONION KAI TIS GNOSIS (Greece)
 AROMATH, INRIA SOPHIAANTIPOLIS (France)
 VISUAL COMPUTING, RHEINISCHWESTFAELISCHE TECHNISCHE HOCHSCHULE AACHEN (Germany)
 OSLO, SINTEF AS (Norway)
 MATH DEPARTMENT, UNIVERSITA DEGLI STUDI DI ROMA TOR VERGATA (Italy)
 INFORMATICS DEPARTMENT), UNIVERSITA DELLA SVIZZERA ITALIANA (Switzerland)
 MATH DEPARTMENT, UNIVERSITAT DE BARCELONA (Spain)
 MATH Dept, UNIVERSITAT LINZ (Austria)
 Naval engineering, University of Strathclyde (UK)
 Math DEPARTMENT, VILNIAUS UNIVERSITETAS (Lithuania)
 GeometryFactory SARL (France)
 Inria contact: Laurent Busé

Summary:
GRAPES aims at considerably advancing the state of the art in Mathematics, CAD, and Machine Learning in order to promote game changing approaches for generating, optimising, and learning 3D shapes, along with a multisectoral training for young researchers. The scientific goals of GRAPES rely on a multidisciplinary consortium composed of leaders in their respective fields. Topnotch research is also instrumental in forming the new generation of European scientists and engineers. Their disciplines span the spectrum from Computational Mathematics, Numerical Analysis, and Algorithm Design, up to Geometric Modelling, Shape Optimisation, and Deep Learning. This allows the 15 PhD candidates to follow either a theoretical or an applied track and to gain knowledge from both research and innovation through a nexus of intersectoral secondments and Networkwide workshops.
Horizontally, our results lead to opensource, prototype implementations, software integrated into commercial libraries as well as open benchmark datasets. These are indispensable for dissemination and training but also to promote innovation and technology transfer. Innovation relies on the active participation of SMEs, either as a beneficiary hosting an ESR or as associate partners hosting secondments. Concrete applications include simulation and fabrication, hydrodynamics and marine design, manufacturing and 3D printing, retrieval and mining, reconstruction and urban planning.
8.3 National initiatives
GdR EFI and GDM
Evelyne Hubert is part of the Scientific Committee of the GdR Equations Fonctionnelles et Interactions (http://
9 Dissemination
9.1 Promoting scientific activities
9.1.1 Scientific events: organisation
General chair, scientific chair
Ioannis Emiris was general cochair of the International Symposium on Symbolic and Algebraic Computation (ISSAC) 2020.
Bernard Mourrain was cochair of the 2020 International conference on Geometric Modeling and Processing (GMP).
Member of the organizing committees
Laurent Busé was a member of the organizing committee of the Journées Nationales de Calcul Formel (JNCF), Luminy, France, March 26.
Angelos Mantzaflaris was proceedings editor of the International Symposium on Symbolic and Algebraic Computation (ISSAC) July 2020, and a cooorganizer of the 4th G+Smo Developer days Workshop, December 2020.
Member of the conference program committees
Evelyne Hubert was part of the program committee for the selection of submissions to the International Symposium on Symbolic and Algebraic Computation (ISSAC) 2020.
Evelyne Hubert was part of the committee to select the plenary speakers for the triennal conference Foundation of Computational Mathematics. The conference was canceled but the Board of Directors is organizing a monthly webinar based on the selected pleanry speakers.
Laurent Busé, Angelos Mantzaflaris and Bernard Mourrain were members of the program committee of the 2020 Symposium on Physical and Solid Modeling (SPM), Strasbourg, France, June 24.
Angelos Mantzaflaris was a member of the program committee of the 2020 Geometric Modeling and Processing (GMP) Conference, Okinawa, Japan, May 1315, 2020.
Reviewer
Laurent Busé reviewed for Symposium on Physical and Solid Modeling conference and SIGGRAPH conference.
Bernard Mourrain reviewed for SPM conference, GMP conference, ISSAC conference.
9.1.2 Journal
Member of the editorial boards
Ioannis Emiris is editor of Journal of Symbolic Computation, Elsevier; Mathematics for Computer Science, SpringerBirkhauser; Journal of Applicable Algebra in Engineering, Communication & Computing, Springer.
Evelyne Hubert is associate editor for the journal of Foundations of Computational Mathematics (Springer) and Journal of Symbolic Computation (Elsevier)
Bernard Mourrain is associate editor of the Journal of Symbolic Computation (since 2007) and of the SIAM Journal on Applied Algebra and Geometry (since 2016).
Reviewer  reviewing activities
Evelyne Hubert reviewed for the journal Foundations of Computational Mathematics, the Journal of Symbolic Computation, Advances in Applied Mathematics, Mathematics and Computers in Simulation, Symmetry, Computers and Graphics. She also write reviews for MathSciNet.
Laurent Busé reviewed for MS Mathematics of Computations journal, Journal of Symbolic Computation, Israel Journal of Mathematics, Computer AidedDesign, Computer Aided Geometric Design, Journal of Computational and Applied Mathematics, Portugaliae Mathematica, Mathematica Scandinavica and Boletín de la Sociedad Matemática Mexicana.
Bernard Mourrain reviewed for the Journal of Symbolic Computation, for the Journal of Computational and Applied Mathematics, the journal ComputerAided Desig, the journal Computer Aided Geometric Design, the Journal of Algebra , the journal Foundations of Computational Mathematics, the FWF Austrian Science Fund.
Angelos Mantzaflaris reviewed for Computer Graphics Forum, Journal of Numerical Methods in Engineering (Wiley), Journal on Applied Algebra and Geometry (SIAM), Computer Aided Geometric Design, Computer Methods in Applied Mechanics and Engineering (Elsevier).
9.1.3 Invited talks
Evelyne Hubert was invited to give a talk at the workshop Symmetry, Randomness, and Computations in Real Algebraic Geometry hosted by the Institute for Computational and Experimental Research in Mathematics (Providence, USA) and at the annual meeting of the GdR Géométrie Différentielle et Mécanique held online.
Angelos Mantzaflaris gave an invited talk at the INdAM Workshop on Geometric Challenges in Isogeometric Analysis January 2731, 2020, Rome, Italy.
9.1.4 Leadership within the scientific community
Bernard Mourain is vice chair of the SIAM Algebraic Geometry group.
9.1.5 Scientific expertise
Evelyne Hubert was part of the selection committee for a Professor position in Mathematics at Université Versaille Saint Quentin.
Evelyne Hubert was part of the working group around the creation of the Inria Project team Calisto (Mireille Bossy).
9.1.6 Research administration
Bernard Mourrain is member of the BCEP (Bureau du Comité des Equipes Projet) of the center Inria Sophia Antipolis.
Laurent Busé is a member of the admnistrative and scientific committee of the labex AMIES and a member of the executive committee of the MSI of UCA. He is also member of the CDT at Inria SophiaAntipolis and the CPRH of the Mathematics Laboratory JeanAlexandre Dieudonné of the University of Nice.
9.2 Teaching  Supervision  Juries
9.2.1 Teaching
 Licence : Ioannis Emiris, Algorithms and complexity, 52 h (L2), NKU Athens
 Licence : Ioannis Emiris, Software development, 26 h (L3), NKU Athens
 Master : Ioannis Emiris, Geometric data science, 52 h (M2), NKU Athens
 Master : Ioannis Emiris, Structural bioinformatics, 39 h (M2), NKU Athens
 Master : Laurent Busé, Geometric Modeling, 18h (M2), Polytech Nice Sophia  Univ Côte d'Azur.
 Master : Laurent Busé, Computational Algebraic Geometry, 30h (M2), Department of Mathematics of University Côte d'Azur.
9.2.2 Supervision
 PhD in progress: Lorenzo Baldi, Structure of moment problems and applications to polynomial optimization. POEMA Marie SkłodowskaCurie ITN, started in October 2019, supervised by Bernard Mourrain.
 PhD in progress: Evangelos Bartzos, Algebraic elimination and Distance graphs. ARCADES Marie SkłodowskaCurie ITN, started in June 2016, NKUA, supervised by Ioannis Emiris.
 PhD in progress: Ayoub Belhachmi, Interactive construction of 3D models  Application to the modeling of complex geological structures. CIFRE, started in August 2020, Inria/Schlumberger, cosupervised by Bernard Mourrain.
 PhD defended on April 8, 2020: Ahmed Blidia, New geometric models for the design and computation of complex shapes. ARCADES Marie SkłodowskaCurie ITN, started in September 2016, supervised by Bernard Mourrain.
 PhD in progress: Apostolos Chalkis, Sampling in highdimensional convex regions, Google Summer of Code and Pegasus national Project, started in June 2018, NKUA, supervised by Ioannis Emiris.
 PhD in progress: Emmanouil Christoforou, Geometric approximation algorithms for clustering, Structural Bioinformatics national infrastructure "Inspire" and Pegasus national Project, started in Jan. 2018, NKUA, supervised by Ioannis Emiris.
 PhD in progress: Carles Checa, Algebraic computing for geometric predicate operations. GRAPES Marie SkłodowskaCurie ITN, started in December 2020, NKUA, supervised by Ioannis Emiris.
 PhD in progress: Riccardo Di Dio, Building a diagnosis tool to detect bronchoconstrictions, BoostUrCAreer Marie SkłodowskaCurie COFUND fellowship. Started on November 2019, cosupervised by Benjamin Mauroy (UCA) and Angelos Mantzaflaris.
 PhD in progress: Edgar Fuentes Figueroa, Polynomial Optimization Techniques for Energy Network Operation and Design, Artelys, POEMA Marie SkłodowskaCurie ITN, started in December 2019, cosupervised by Michael Gabay (Artelys) and Bernard Mourrain.
 PhD in progress: Pablo Gonzalez Mazon; Generation of valid highorder curved meshes. GRAPES Marie SkłodowskaCurie ITN, started in December 2020, Inria, supervised by Laurent Busé.
 PhD in progress: Mehran Hatamzadeh; An innovative gait analysis technology. PhD grant from the EU CoFUND BoostUrCareer program of UCA, cosupervised by Laurent Busé and Raphaël Zory (LAMHESS, UCA).
 PhD in progress: Rima Khouja, Tensor decomposition, best approximations, algorithms and applications. Cotutelle Univ. Liban, started in November 2018, cosupervised by Houssam Khalil and Bernard Mourrain.
 PhD in progress: Thomas Laporte, Towards a 4D model of the respiratory system. Fellowship from ED SFA/UCA. Started on October 2019, cosupervised by Benjamin Mauroy (UCA) and Angelos Mantzaflaris.
 PhD defended on April, 30, 2020: Clément Laroche, Algebraic representations of geometric objects. ARCADES Marie SkłodowskaCurie ITN, started in Nov. 2016, NKUA; supervised by Ioannis Emiris.
 PhD in progress: Tobias Metzlaff. Multivariate orthogonal polynomials and applications to global optimization. POEMA Marie SkłodowskaCurie ITN, started in December 2019, supervised by Evelyne Hubert.
 PhD in progress: Michelangelo Marsala, Modelling and simulation using analysissuitable subdivision surfaces and solids. GRAPES Marie SkłodowskaCurie ITN, started in November 2020, Inria, supervised by Angelos Mantzaflaris and Bernard Mourrain.
 PhD: Erick Rodriguez Bazan defended on December 2020. Symmetry in Algebraic Computations. Inria. Started in November 2017, supervised by Evelyne Hubert.
 PhD in progress: Konstantinos Tertikas, Machine learning for geometric shapes. GRAPES Marie SkłodowskaCurie ITN, started in December 2020, NKUA, supervised by Ioannis Emiris.
 PhD in progress: Tong Zhao; Learning priors and metrics for 3D reconstruction of largescale scenes. PhD grant from the 3IA Institut of UCA, cosupervised by Pierre Alliez (TITANE) and Laurent Busé.
9.2.3 Juries
Bernard Mourrain was a member of the PhD committee of Clement Laroche, Algebraic representations of geometric objects. NKUA, April, 30, 2020.
Evelyne Hubert was a reviewer for the Habilitation thesis of Georg Regensbuger, Johannes Kepler University (Austria): Algebraic and algorithmic Approached to Analysis: Integrodifferential equations, positive steady states, and wavelets.
Laurent Busé was a reviewer and member of the committee of the PhD thesis of This Xuan Vu, Homotopy algorithms for solving structured determinantal systems, University of Waterloo, Waterloo, Ontario, Canada, December 9th.
Angelos Mantzaflaris was a reviewer of the PhD thesis of Gabriele Loli, Efficient Solvers for Isogeometric Analysis, University of Pavia, Italy.
10 Scientific production
10.1 Major publications
 1 articleOn the multihomogeneous Bézout bound on the number of embeddings of minimally rigid graphsApplicable Algebra in Engineering, Communication and Computing31562020, 325357
 2 articleDegree and birationality of multigraded rational mapsProceedings of the London Mathematical Society12142020, 743787
 3 articleResultant of an equivariant polynomial system with respect to the symmetric groupJournal of Symbolic Computation762016, 142157
 4 articleProducts of Euclidean Metrics, Applied to Proximity Problems among CurvesACM Transactions on Spatial Algorithms and Systems64August 2020, 120
 5 articleScaffolding skeletons using spherical Voronoi diagrams: feasibility, regularity and symmetryComputerAided Design102May 2018, 83  93
 6 articleLocal (T)HBspline projectors via restricted hierarchical spline fittingComputer Aided Geometric Design80June 2020, 101865
 7 articleRational invariants of even ternary forms under the orthogonal groupFoundations of Computational Mathematics192019, 13151361
 8 articleElimination ideals and Bezout relationsJournal of Algebra5622020, 621626
 9 articleLow Rank Tensor Methods in Galerkinbased Isogeometric AnalysisComputer Methods in Applied Mechanics and Engineering316April 2017, 10621085
 10 articlePolynomialExponential Decomposition from MomentsFoundations of Computational Mathematics186December 2018, 14351492
 11 articleSolving Polynomial Systems via a Stabilized Representation of Quotient AlgebrasSIAM Journal on Matrix Analysis and Applications393October 2018, 14211447
10.2 Publications of the year
International journals
 12 articleLefschetz properties of monomial algebras with almost linear resolutionCommunications in Algebra484April 2020, 14991509
 13 articleHighDimensional Approximate rNetsAlgorithmica826June 2020, 16751702
 14 articleOn the multihomogeneous Bézout bound on the number of embeddings of minimally rigid graphsApplicable Algebra in Engineering, Communication and Computing31562020, 325357
 15 articleGeometrically smooth spline bases for data fitting and simulationComputer Aided Geometric Design78March 2020, 101814
 16 articleFibers of multigraded rational maps and orthogonal projection onto rational surfacesSIAM Journal on Applied Algebra and Geometry422020, 322353
 17 article Determinantal tensor product surfaces and the method of moving quadrics Transactions of the American Mathematical Society 2020
 18 articleDegree and birationality of multigraded rational mapsProceedings of the London Mathematical Society12142020, 743787
 19 articleThe geometry of the flex locus of a hypersurfacePacific Journal of Mathematics30422020, 419437
 20 articleFreeness and invariants of rational plane curvesMathematics of Computation892020, 15251546
 21 articleMatrix formulae for Resultants and Discriminants of Bivariate Tensorproduct PolynomialsJournal of Symbolic Computation98June 2020, 6583
 22 articleRegularity of bicyclic graphs and their powersJournal of Algebra and Its Applications1903March 2020, 2050057
 23 articleMultilinear Polynomial Systems: Root Isolation and Bit ComplexityJournal of Symbolic Computation1052021, 145164
 24 articleSeparation bounds for polynomial systemsJournal of Symbolic Computation1012020, 128151
 25 articleProducts of Euclidean Metrics, Applied to Proximity Problems among CurvesACM Transactions on Spatial Algorithms and Systems64August 2020, 120
 26 articleLocal (T)HBspline projectors via restricted hierarchical spline fittingComputer Aided Geometric Design80June 2020, 101865
 27 articleElimination ideals and Bezout relationsJournal of Algebra5622020, 621626
 28 articleComplete Classification and Efficient Determination of Arrangements Formed by Two EllipsoidsACM Transactions on Graphics393June 2020, 112
 29 article Practical isogeometric shape optimization: Parameterization by means of regularization Journal of computational design and engineering January 2021
 30 articleOn minimal decompositions of low rank symmetric tensorsLinear Algebra and its Applications607December 2020, 347377
 31 articleEfficient Matrix Assembly in Isogeometric Analysis with Hierarchical BsplinesJournal of Computational and Applied Mathematics3902020, 113278
 32 article Multivariate Interpolation: Preserving and Exploiting Symmetry Journal of Symbolic Computation 2021
 33 articleInterpolatory CatmullClark volumetric subdivision over unstructured hexahedral meshes for modeling and simulation applicationsComputer Aided Geometric Design80June 2020, 101867
International peerreviewed conferences
 34 inproceedings Computation of free noncommutative Gröbner Bases over Z with Singular:Letterplace ISSAC 2020  International Symposium on Symbolic and Algebraic Computation Kalamata, Greece July 2020
 35 inproceedings Punctual Hilbert Schemes and Certified Approximate Singularities ISSAC 2020  International Symposium on Symbolic and Algebraic Computation Kalamata, Greece July 2020
 36 inproceedings Ideal Interpolation, HBases and Symmetry ISSAC 2020  International Symposium on Symbolic and Algebraic Computation Kalamata, Greece July 2020
Scientific book chapters
 37 inbookAn overview of Geometry plus Simulation ModulesMathematical Aspects of Computer and Information Sciences. MACIS 2019. Lecture Notes in Computer Science, vol 11989.March 2020, 453456
Doctoral dissertations and habilitation theses
 38 thesis Geometric continuity for surfaces and scalar fields Université Côte d'Azur July 2020
 39 thesis Compact and efficient implicit representations Université d'Athènes April 2020
 40 thesis Finite fibers of rational maps by means of matrix representations with applications to distance problem Université Côte d'Azur February 2020
Reports & preprints
 41 misc Computing real radicals by moment optimization February 2021
 42 misc Exact Moment Representation in Polynomial Optimization February 2021

43
misc
The Hessian polynomial and the Jacobian ideal of a reduced surface in
${}^{3}$ ' 2020  44 misc Formulas for the eigendiscriminants of ternary and quaternary forms June 2020
 45 misc Efficient Sampling from Feasible Sets of SDPs and Volume Approximation October 2020
 46 misc Curve Valuations and Mixed Volumes in the Implicitization of Rational Varieties October 2020
 47 misc Sparse Interpolation in Terms of Multivariate Chebyshev Polynomials January 2020
 48 misc A Riemannian Newton Optimization Framework for the Symmetric Tensor Rank Approximation Problem July 2020
 49 misc Computing Free Noncommutative Groebner Bases over Z with Singular:Letterplace December 2020
 50 misc A certified iterative method for isolated singular roots December 2020
 51 misc Symmetry in Multivariate Ideal Interpolation February 2021
 52 misc StretchBased Hyperelastic Material Formulations for Isogeometric KirchhoffLove Shells July 2020
 53 misc Completeness characterization of TypeI box splines November 2020