EN FR
EN FR

2023Activity reportProject-TeamAROMATH

RNSR: 201622154R
  • Research center Inria Centre at Université Côte d'Azur
  • In partnership with:National & Kapodistrian University of Athens
  • Team name: AlgebRa, geOmetry, Modeling and AlgoriTHms
  • Domain:Algorithmics, Programming, Software and Architecture
  • Theme:Algorithmics, Computer Algebra and Cryptology

Keywords

Computer Science and Digital Science

  • A5.5.1. Geometrical modeling
  • A6.1. Methods in mathematical modeling
  • A8.3. Geometry, Topology
  • A8.4. Computer Algebra

Other Research Topics and Application Domains

  • 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, Senior Researcher, HDR]
  • Evelyne Hubert [INRIA, Senior Researcher, HDR]
  • Angelos Mantzaflaris [INRIA, Researcher]

Faculty Members

  • Ioannis Emiris [UNIV NKUA, Greece, HDR]
  • André Galligo [UNIV. COTE D'AZUR, Emeritus]

Post-Doctoral Fellow

  • Christos Konaxis [UNIV. NKUA, Greece]

PhD Students

  • Emmanouil Christoforou [UNIV NKUA Greece]
  • Konstantinos Tertikas [ATHENA RC, UNIV. NKUA]
  • Ayoub Bel Hachmi [INRIA, SCHLUMBERGER]
  • Carles Checa [ATHENA RC, UNIV. NKUA, Greece]
  • Riccardo Di Dio [UNIV. COTE D'AZUR]
  • Pablo Gonzalez Mazon [INRIA]
  • Mehran Hatamzadeh [UNIV COTE D'AZUR]
  • Martin Jalard [INRIA]
  • Christos Karampatzakis [INRIA]
  • Thomas Laporte [UNIV. COTE D'AZUR, until Jan 2023]
  • Michelangelo Marsala [INRIA]

Technical Staff

  • Linh Tuy Nguyen [INRIA, until Jun 2023]

Administrative Assistant

  • Sophie Honnorat [INRIA]

Visiting Scientists

  • Sofia Imperatore [UNIV FLORENCE, from Apr 2023 until Jun 2023]
  • Hugo Verhelst [UNIV DELFT, from May 2023 until Jun 2023]

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 two-way. 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. CAD-CAM (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 daily-life 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 e-trade 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 CAD-CAM. 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 Computer Aided Geometric Design (CAGD) to describe complex geometry is based on parametric representations called NURBS (Non Uniform Rational B-Spline). The models are constructed by trimming and gluing together high order patches of algebraic surfaces. These models are built from the so-called B-Spline 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 so-called T-Splines, LR-Spline 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 skeleton-based 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 algebraic-geometric 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 input-output map is regular. In non-regular 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, multi-homogeneous, 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 semi-definite 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.

Computer-Aided Design (CAD) involves creating digital models defined by mathematical constructions, from geometric, functional or aesthetic considerations. Computer-aided 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.

CAD-CAM 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.

CAD-CAM 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, 5-axes 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 CAD-CAM 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 CAD-systems 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 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 time-consuming 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 B-Spline 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 so-called 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

5.1 Awards

Bernard Mourrain was awarded Solid Modeling Fellow of the Solid Modeling Association in 2023.

6 New software, platforms, open data

6.1 New software

6.1.1 G+Smo

  • Name:
    Geometry plus Simulation Modules
  • Keyword:
    Isogeometric analysis
  • Functional Description:

    G+Smo (pronounced gismo or gizmo) is a C++ library for isogeometric analysis (IGA).

    G+Smo (Geometry + Simulation Modules, pronounced "gismo") is an open-source 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 object-oriented, cross-platform, template C++ 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, open-source software to the forming isogeometric numerical simulation community and beyond. Geometry plus simulation modules aims at the seamless integration of Computer-aided Design (CAD) and high order Finite Element Analysis (FEA).

    The library and its documentation are available at https://gismo.github.io/

  • URL:
  • Contact:
    Angelos Mantzaflaris

6.1.2 MomentTools

  • Name:
    MomentTools
  • Keywords:
    Global optimization, Moment, Polynomial equations, Semi-algebraic set, Convex relaxation
  • Functional Description:

    The package provides efficient tools to build convex relaxations of moment sequences and their dual Sum-of-Squares relaxations, to optimize vectors of moment sequences that satisfy positivity constraints or mass constraints, to compute global minimizers of polynomial and moment optimization problems from moment sequences, polar ideals, approximate real radical. It also provides tools for computing minimum enclosing ellipsoids of basic semi-algebraic sets. It uses a connection with SDP solvers via the JuMP interface.

    The package is available at https://github.com/AlgebraicGeometricModeling/MomentTools.jl and its documentation at https://algebraicgeometricmodeling.github.io/MomentTools/

  • URL:
  • Contact:
    Bernard Mourrain
  • Participants:
    Lorenzo Baldi, Bernard Mourrain

6.1.3 TensorDec

  • Keywords:
    Tensor decomposition, Multivariate series, Low rank models, Hankel
  • Functional Description:

    TensorDec is a Julia package for the decomposition of tensors and polynomial-exponential series. It provides tools to compute rank decomposition or Waring decomposition of symmetric tensors or multivariate homogeneous, of multilinear tensors.

    It also allows computing low rank tensor approximations of given tensors, using Riemannian optimization techniques, with well-chosen initial start. It also provides tools to compute catalecticant or Hankel operators associated to tensors and their apolar ideal.

    The package is accessible at https://github.com/AlgebraicGeometricModeling/TensorDec.jl and its documentation at https://algebraicgeometricmodeling.github.io/TensorDec.jl/.

  • URL:
  • Contact:
    Bernard Mourrain
  • Participants:
    Rima Khouja, Bernard Mourrain

7 New results

7.1 Polynomial Optimization, Moments, and Applications

Participant: Michal Kocvara [Univ. Birmingham, UK], Bernard Mourrain, Cordian Riener [Univ. Tromso, Norway].

Polynomial optimization is a fascinating field of study that has revolutionized the way we approach nonlinear problems described by polynomial constraints. The applications of this field range from production planning processes to transportation, energy consumption, and resource control. This introductory book explores the latest research developments in polynomial optimization, presenting the results of cutting-edge interdisciplinary work conducted by the European network POEMA. For the past four years, experts from various fields, including algebraists, geometers, computer scientists, and industrial actors, have collaborated in this network to create new methods that go beyond traditional paradigms of mathematical optimization. By exploiting new advances in algebra and convex geometry, these innovative approaches have resulted in significant scientific and technological advancements. The book 31 aims to make these exciting developments accessible to a wider audience by gathering high-quality chapters on these hot topics. Aimed at both aspiring and established researchers, as well as industry professionals, this book will be an invaluable resource for anyone interested in polynomial optimization and its potential for real-world applications. The coeditors of the book are M. Kocvara and C. Riener from POEMA network.

7.2 Polynomial optimization in geometric modeling

Participant: Sodeeh Habibi [Univ. Birmingham, UK], Michal Kocvara [Univ. Birmingham, UK], Bernard Mourrain.

In the chapter 32 of the book 31, together with S. Habibi and M. Kocvara from POEMA network, we review applications of Polynomial Optimization techniques to Geometric Modeling problems. We present examples of topical problems in Geometric Modeling, illustrate their solution using Polynomial Optimization Tools, report some experimental results and analyse the behavior of the methods, showing what are their strengths and their limitations.

7.3 Algebraic Curves and Surfaces, A History of Shapes

Participants: Laurent Busé, Fabrizio Catanese [Bayreuth University], Elisa Postinghel [Trento University].

The aim of the book 30, co-authored with Fabrizio Catanese (Bayreuth University) and Elisa Postinghel (Trento University) is manifold, it intends to overview the wide topic of algebraic curves and surfaces (also with a view to higher dimensional varieties) from different aspects: the historical development that led to the theory of algebraic surfaces and the classification theorem of algebraic surfaces by Castelnuovo and Enriques; the use of such a classical geometric approach, as the one introduced by Castelnuovo, to study linear systems of hypersurfaces; and the algebraic methods used to find implicit equations of parametrized algebraic curves and surfaces, ranging from classical elimination theory to more modern tools involving syzygy theory and Castelnuovo-Mumford regularity. Since our subject has a long and venerable history, this book cannot cover all the details of this broad topic, theory and applications, but it is meant to serve as a guide for both young mathematicians to approach the subject from a classical and yet computational perspective, and for experienced researchers as a valuable source for recent applications.

7.4 G1 Spline Functions for Point Cloud Fitting

Participants: Michelangelo Marsala, Angelos Mantzaflaris, Bernard Mourrain.

In 20 we present a new construction of basis functions that generate the space of geometrically smooth splines on an unstructured quadrilateral mesh. The basis is represented in terms of biquintic Bézier polynomials on each quadrilateral face. The gluing along the face boundaries is achieved using quadratic gluing data functions, leading to globally G1-smooth spaces. We analyze the latter space and provide a combinatorial formula for its dimension as well as an explicit basis construction. Moreover, we assess the use of this basis in point cloud fitting problems. To apply G1 least squares fitting, a quadrilateral structure as well as parameters in each quadrilateral is required. Even though the general problem of segmenting and parametrizing point clouds is beyond the focus of the present work, we describe a procedure that produces such a structure as well as patch-local parameters. Our experiments demonstrate the accuracy and smoothness of the obtained reconstructed models in several challenging instances.

7.5 Variational Shape Reconstruction via Quadric Error Metrics

Participants: Laurent Busé, Tong Zhao [Titane], David Cohen-Steiner [Datashape], Tamy Boubekeur [Adobe Research Lab], Jean-Marc Thiery [Adobe Research Lab], Pierre Alliez [Titane].

Inspired by the strengths of quadric error metrics initially designed for mesh decimation, we proposed in 29 a concise mesh reconstruction approach for 3D point clouds. Our approach proceeds by clustering the input points enriched with quadric error metrics, where the generator of each cluster is the optimal 3D point for the sum of its quadric error metrics. This approach favors the placement of generators on sharp features, and tends to equidistribute the error among clusters. We reconstruct the output surface mesh from the adjacency between clusters and a constrained binary solver. We combine our clustering process with an adaptive refinement driven by the error. Compared to prior art, our method avoids dense reconstruction prior to simplification and produces immediately an optimized mesh.

7.6 A comparison of smooth basis constructions for isogeometric analysis

Participants: Angelos Mantzaflaris, Hugo M. Verhelst [TU Delft, The Netherlands], Pascal Weinmüller [MTU, Germany], Thomas Takacs [ RICAM - Johann Radon Institute for Computational and Applied Mathematics, Linz, Austria], Deepesh Toshniwal [TU Delft, The Netherlands].

In order to perform isogeometric analysis with increased smoothness on complex domains, trimming, variational coupling or unstructured spline methods can be used. The latter two classes of methods require a multi-patch segmentation of the domain, and provide continuous bases along patch interfaces. In the context of shell modeling, variational methods are widely used, whereas the application of unstructured spline methods on shell problems is rather scarce. In 22 we therefore provide a qualitative and a quantitative comparison of a selection of unstructured spline constructions, in particular the D-Patch, Almost-C1 , Analysis-Suitable C1 and the Approximate C1 constructions. Using this comparison, we aim to provide insight into the selection of methods for practical problems, as well as directions for future research. In the qualitative comparison, the properties of each method are evaluated and compared. In the quantitative comparison, a selection of numerical examples is used to highlight different advantages and disadvantages of each method. In the latter, comparison with weak coupling methods such as Nitsche's method or penalty methods is made as well. In brief, it is concluded that the Approximate C1 and Analysis-Suitable G1 converge optimally in the analysis of a bi-harmonic problem, without the need of special refinement procedures. Furthermore, these methods provide accurate stress fields. On the other hand, the Almost-C1 and D-Patch provide relatively easy construction on complex geometries. The Almost-C1 method does not have limitations on the valence of boundary vertices, unlike the D-Patch, but is only applicable to biquadratic local bases. Following from these conclusions, future research directions are proposed, for example towards making the Approximate C1 and Analysis-Suitable G1 applicable to more complex geometries.

7.7 Leveraging Moving Parameterization and Adaptive THB-Splines for CAD Surface Reconstruction of Aircraft Engine Components

Participants: Angelos Mantzaflaris, Carlotta Giannelli [DIMAI UniFI, Firenze, Italy], Sofia Imperatore [DIMAI UniFI, Firenze, Italy], Dominik Mokriš [MTU Aero Engines, Germany].

Reconstruction of highly accurate CAD models from point clouds is both paramount and challenging in industries such as aviation. Due to the acquisition process, this kind of data can be scattered and affected by noise, yet the reconstructed geometric models are required to be compact and smooth, while simultaneously capturing key geometric features of the engine parts. In 23, we present an iterative moving parameterization approach, which consists of alternating steps of surface fitting, parameter correction, and adaptive refinement using truncated hierarchical B-splines (THB-splines). We revisit two existing surface fitting methods, a global least squares approximation and a hierarchical quasi-interpolation scheme, both based on THB-splines. At each step of the adaptive loop, we update the parameter locations by solving a non-linear optimization problem to infer footpoints of the point cloud on the current fitted surface. We compare the behavior of different optimization settings for the critical task of distance minimization, by also relating the effectiveness of the correction step to the quality of the initial parameterization. In addition, we apply the proposed approach in the reconstruction of aircraft engine components from scanned point data. It turns out that the use of moving parameterization instead of fixed parameter values, when suitably combined with the adaptive spline loop, can significantly improve the resulting surfaces, thus outperforming state-of-the-art hierarchical spline model reconstruction schemes.

7.8 Learning meshless parameterization with (graph) convolutional neural networks

Participants: Angelos Mantzaflaris, Carlotta Giannelli [DIMAI UniFI, Firenze, Italy], Sofia Imperatore [DIMAI UniFI, Firenze, Italy], Dominik Mokriš [MTU Aero Engines, Germany], Felix Scholz [Institute of Applied Geometry, Linz, Austria].

The publication 24 proposes a deep learning approach for parameterizing an unorganized or scattered point cloud in R 3 with graph convolutional neural networks. It builds upon a graph convolutional neural network that predicts the weights (called parameterization weights) of certain convex combinations that lead to a mapping of the 3D points into a planar parameter domain. First, we compute a radius neighbors graph that yields proximity information to each 3D point in the cloud. This radius graph is then converted to its line graph, which encodes edge adjacencies, and is equipped with appropriate weights. The line graph is used as input to a graph convolutional neural network trained to predict optimal parameterizations. The proposed model outperforms closed-form choices of the parameterization weights and produces high quality parameterizations for surface reconstruction schemes. In a similar spirit, in 28 we propose a dimension independent method based on convolutional neural networks to assign parameter values to gridded point clouds of arbitrary size, without the need for additional data processing steps. We train the proposed networks by considering polynomial least squares approximations and demonstrate, both in the univariate and bivariate settings, that the accuracy of the final model properly scales when uniform and adaptive spline refinement is considered. A selection of numerical experiments on point clouds of different sizes highlights the performance of our parameterization scheme. Noisy data sets which simulate measurement errors are also considered.

7.9 Contact simulation of tooth flanks using Isogeometric Analysis

Participants: Angelos Mantzaflaris, Christos Karampatzakis [Aristotle University of Thessaloniki, Greece], Athanasios Mihailidis [Aristotle University of Thessaloniki, Greece], Christopher Provatidis [NTUA - National Technical University of Athens, Greece].

The Finite Element (FE) modeling of a tooth flank affects both the accuracy of the results as well as the computational time and resources required. In FE Analysis (FEA), the involute curve is modelled by linear segments, thus deviations from its actual geometry are introduced. In order to achieve accurate results, the mesh must be generated as dense as possible. This increases the number of degrees of freedom and the computational cost along with it. This inherent drawback of FEA unfavorably affects the calculation of the pressure distribution of the mating tooth flanks. Isogeometric Analysis (IGA) is a recent alternative to the FEA. It uses B-Splines, a technology used by Computer Aided Design (CAD) systems, to model the geometry as well as the solution field. Thus, no geometric error is introduced in the transition from CAD to analysis. An inherent characteristic of B-Splines is the continuity between adjacent elements. Furthermore, the smooth normal vector field of the surface is known at every point in the interior of elements. This is particularly advantageous for contact algorithms because this alleviates the need for special smoothing techniques such as those often used in FEA. In the study 25, a spur gear pair is simulated and the pressure at the contact area is calculated. The results are compared to those obtained by FEA in terms of both accuracy and computational cost.

7.10 Sampling – Practical volume approximation of high-dimensional convex bodies, applied to modeling portfolio dependencies and financial crises

Participants: Ioannis Emiris, Ludovic Calès [ JRC - European Commission - Joint Research Centre, Ispra, Italy], Apostolos Chalkis, Vissarion Fisikopoulos [ARC - Athena Research and Innovation Centre, Athens, Greece].

In 15 we applied basic algorithms for volume computation of general-dimensional polytopes and more general convex bodies, defined by the intersection of a simplex by a family of parallel hyperplanes, and another family of parallel hyperplanes or a family of concentric ellipsoids. Such convex bodies appear in modeling and predicting financial crises. The impact of crises on the economy (labor, income, etc.) makes its detection of prime interest for the public in general and for policy makers in particular. Certain features of dependencies in the markets clearly identify times of turmoil. We describe the relationship between asset characteristics by means of a copula; each characteristic is either a linear or quadratic form of the portfolio components, hence the copula can be estimated by computing volumes of convex bodies. We design and implement practical algorithms in the exact and approximate setting, and experimentally juxtapose them in order to study the trade-off of exactness and accuracy for speed. We also experimentally find an efficient parameter-tuning to achieve a sufficiently good estimation of the probability density of each copula. Our C++ software, based on Eigen and available on github, is shown to be very effective in up to 100 dimensions. Our results offer novel, effective means of computing portfolio dependencies and an indicator of financial crises, which is shown to correctly identify past crises.

7.11 Sampling – A Practical Algorithm for Volume Estimation based on Billiard Trajectories and Simulated Annealing

Participants: Ioannis Emiris, Apostolos Chalkis, Vissarion Fisikopoulos [ARC - Athena Research and Innovation Centre, Athens, Greece].

In 16 our most recent techniques tackle the problem of efficiently approximating the volume of convex polytopes, when these are given in 3 different representations: H-polytopes, which have been studied extensively, V-polytopes, and zonotopes (Z-polytopes). We design a novel practical Multiphase Monte Carlo algorithm that leverages random walks based on billiard trajectories, as well as a new empirical convergence test and a simulated annealing schedule of adaptive convex bodies. After tuning several parameters of our proposed method, we present a detailed experimental evaluation of our tuned algorithm using a rich dataset containing Birkhoff polytopes and polytopes from structural biology. Our open-source implementation tackles problems that have been intractable so far, offering the first software to scale up in thousands of dimensions for H-polytopes and in the hundreds for V- and Z-polytopes on moderate hardware. Last, we illustrate our software in evaluating Z-polytope approximations.

7.12 Sampling – Geometric algorithms for sampling the flux space of metabolic networks

Participants: Ioannis Emiris, Apostolos Chalkis, Vissarion Fisikopoulos [ARC - Athena Research and Innovation Centre, Athens, Greece], Elias Tsigaridas [Ouragan], Haris Zafeiropoulos [ UOC - University of Crete, Heraklion, Greece].

In 17 we apply random walk technology to systems biology. Metabolic networks and their reconstruction set a new era in the analysis of metabolic and growth functions in the various organisms. By modeling the reactions occurring inside an organism, metabolic networks provide the means to understand the underlying mechanisms that govern biological systems. Constraint-based approaches have been widely used for the analysis of such models and led to intriguing geometry-oriented challenges. In this setting, sampling uniformly points from polytopes derived from metabolic models (flux sampling) provides a representation of the solution space of the model under various conditions. However, the polytopes that result from such models are of high dimension (in the order of thousands) and usually considerably skinny. Therefore, to sample uniformly at random from such polytopes shouts for a novel algorithmic and computational framework specially tailored for the properties of metabolic models. We present a complete software framework to handle sampling in metabolic networks. Its backbone is a Multiphase Monte Carlo Sampling (MMCS) algorithm that unifies rounding and sampling in one pass, yielding both upon termination. It exploits an optimized variant of the Billiard Walk that enjoys faster arithmetic complexity per step than the original. We demonstrate the efficiency of our approach by performing extensive experiments on various metabolic networks. Notably, sampling on the most complicated human metabolic network accessible today, Recon3D, corresponding to a polytope of dimension 5335, took less than 30 hours. To the best of our knowledge, that is out of reach for existing software.

7.13 Deep Geometric learning – PartNeRF: Generating Part-Aware Editable 3D Shapes without 3D Supervision

Participants: Ioannis Emiris, Konstantinos Tertikas, Despoina Paschalidou [Stanford University, USA], Boxiao Pan [Stanford University, USA], Jeong Joon Park [Stanford University, USA], Mikaela Angelina Uy [Stanford University, USA], Yannis Avrithis [IARAI - Institute of Advanced Research in Artificial Intelligence, Vienna, Austria], Leonidas Guibas [Stanford University, USA].

Impressive progress in generative models and implicit representations gave rise to methods that can generate 3D shapes of high quality. However, being able to locally control and edit shapes is another essential property that can unlock several content creation applications. Local control can be achieved with part-aware models, but existing methods require 3D supervision and cannot produce textures. In 27, we devise PartNeRF, a novel part-aware generative model for editable 3D shape synthesis that does not require any explicit 3D supervision. Our model generates objects as a set of locally defined NeRFs, augmented with an affine transformation. This enables several editing operations such as applying transformations on parts, mixing parts from different objects etc. To ensure distinct, manipulable parts we enforce a hard assignment of rays to parts that makes sure that the color of each ray is only determined by a single NeRF. As a result, altering one part does not affect the appearance of the others. Evaluations on various ShapeNet categories demonstrate the ability of our model to generate editable 3D objects of improved fidelity, compared to previous part-based generative approaches that require 3D supervision or models relying on NeRFs.

7.14 Fairness in machine learning – Fairness Aware Counterfactuals for Subgroups

Participants: Ioannis Emiris, Loukas Kavouras [ARC - Athena Research and Innovation Centre, Athens, Greece], Konstantinos Tsopelas [ARC - Athena Research and Innovation Centre, Athens, Greece], Giorgos Giannopoulos [ARC - Athena Research and Innovation Centre, Athens, Greece], Dimitris Sacharidis [ULB - Université libre de Bruxelles, Belgique], Eleni Psaroudaki [NTUA - National Technical University of Athens, Greece], Nikolaos Theologitis [ARC - Athena Research and Innovation Centre, Athens, Greece], Dimitrios Rontogiannis [NKUA - National and Kapodistrian University of Athens, Greece], Dimitris Fotakis [].

In 26, we present Fairness Aware Counterfactuals for Subgroups (FACTS), a framework for auditing subgroup fairness through counterfactual explanations. We start with revisiting (and generalizing) existing notions and introducing new, more refined notions of subgroup fairness. We aim to (a) formulate different aspects of the difficulty of individuals in certain subgroups to achieve recourse, i.e. receive the desired outcome, either at the micro level, considering members of the subgroup individually, or at the macro level, considering the subgroup as a whole, and (b) introduce notions of subgroup fairness that are robust, if not totally oblivious, to the cost of achieving recourse. We accompany these notions with an efficient, model-agnostic, highly parameterizable, and explainable framework for evaluating subgroup fairness. We demonstrate the advantages, the wide applicability, and the efficiency of our approach through a thorough experimental evaluation of different benchmark datasets.

7.15 Discriminative sEMG-based features to assess damping ability and interpret activation patterns in lower-limb muscles of ACLR athletes

Participants: Mehran Hatamzadeh, Ali Sharifnezhad [ Department of Sport Biomechanics and Technology, Sport Sciences Research Institute, Tehran, Iran], Reza Hassannejad [ University of Tabriz, Iran], Raphael Zory [ LAMHESS - Laboratoire Motricité Humaine Expertise Sport Santé, Nice].

The main goal of the athletes who undergo anterior cruciate ligament reconstruction (ACLR) surgery is a successful return-to-sport. At this stage, identifying muscular deficits becomes important. Hence, in 18, we study three discriminative features based on surface electromyographic signals (sEMG) acquired in a dynamic protocol are introduced to assess the damping ability and interpret activation patterns in lower-limb muscles of ACLR athletes. Methods: The features include the median frequency of the power spectrum density (PSD), the relative percentage of the equivalent damping or equivalent stiffness derived from the median frequency, and the energy of the signals in the time-frequency plane of the pseudo-Wigner-Ville distribution (PWVD). To evaluate the features, 11 healthy and 11 ACLR athletes (6 months post-reconstruction surgery) were recruited to acquire the sEMG signals from the medial and the lateral parts of the hamstrings, quadriceps, and gastrocnemius muscles in pre- and post-fatigue single-leg landings. Results: A significant damping deficiency is observed in the hamstring muscles of ACLR athletes by evaluating the proposed features. This deficiency indicates that more attention should be paid to this muscle of ACLR athletes in pre-return-to-sport rehabilitations. Conclusion: The quality of electromyography-based pre-return-to-sport assessments on ACLR subjects depends on the sEMG acquisition protocol, as well as the type and nature of the extracted features. Hence, combinatorial application of both energy-based features (derived from the PWVD) and power-based features (derived from the PSD) could facilitate the assessment process by providing additional biomechanical information regarding the behavior of the muscles surrounding the knee.

8 Bilateral contracts and grants with industry

8.1 Bilateral grants with industry

Participants: Ioannis Emiris, Panagiotis Repouskos, Ioannis Psarros.

Geometric computing.

Ioannis Emiris coordinates a research contract with the industrial partner ANSYS Inc. (Greece), in collaboration with Athena Research Center. MSc students P. Repouskos, M. Dioletis, and T. Pappas, postdoc fellow I. Psarros, and Athena researcher George Ioannakis are 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 axis-aligned objects but the real challenge is the general case. The project, extended into 2024, 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 L-max metric, which is the primary focus of the contract, (2) random sampling of the free space among objects, (3) data-driven algorithmic design for problems concerning data-structures and their construction and initialization. The implementation of prototypes has led into the development of a software library including implementations of parallel algorithms (Cuda).

It has continued into 2023 along with a tripartite grant from the Greek ministry of Development (Athena RC, ANSYS, and University of Patras, Greece).

Participants: Ayoub Belhachmi, Bernard Mourrain.

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.

Participants: Laurent Busé, Angelos Mantzaflaris.

Spline surfaces for optics applications.OpenReality is a startup company that is developing curved photonic optics to enable high-fidelity Augmented Reality solutions. In the frame of a consulting contract we collaborated on the development of 3D curved photonic surfaces represented by splines.

9 Partnerships and cooperations

9.1 International initiatives

9.1.1 Participation in other International Programs

PHC Galileo (Hubert Curien Partnership, 2022–2023)

Participants: Angelos Mantzaflaris, Michelangelo Marsala.

Collaboration with the Universiry of Florence. Efficient high order learning for deep geometric design networks. The research program focuses on the interaction between the computational side of geometric models, and their application-oriented side devoted to the design and analysis of efficient adaptive spline approximation schemes. The work plan addresses important challenges in the area of geometric modeling and processing, creating a connection with suitable machine learning applications.

9.2 European initiatives

9.2.1 H2020 projects

POEMA

POEMA project on cordis.europa.eu

  • Title:
    Polynomial Optimization, Efficiency through Moments and Algebra
  • Duration:
    From January 1, 2019 to June 30, 2023
  • Partners:
    • INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET AUTOMATIQUE (INRIA), France
    • UNIVERSITETET I TROMSOE - NORGES ARKTISKE UNIVERSITET (UiT), Norway
    • FRIEDRICH-ALEXANDER-UNIVERSITAET ERLANGEN-NUERNBERG (FAU), Germany
    • UNIVERSITA DEGLI STUDI DI FIRENZE (UNIFI), Italy
    • THE UNIVERSITY OF BIRMINGHAM (UoB), United Kingdom
    • UNIVERSITAT KONSTANZ (UKON), Germany
    • TILBURG UNIVERSITY- UNIVERSITEIT VAN TILBURG (TILBURG UNIVERSITY), Netherlands
    • ARTELYS, France
    • CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE CNRS (CNRS), France
    • STICHTING NEDERLANDSE WETENSCHAPPELIJK ONDERZOEK INSTITUTEN (NWO-I), Netherlands
    • SORBONNE UNIVERSITE, France
  • Inria contact:
    Bernard Mourrain
  • Coordinator:
    Bernard Mourrain
  • Summary:

    Non-linear optimization problems are present in many real-life 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 Semi-Definite Programming problems. This approach provides a powerful way to compute global solutions of non-linear 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. The network will train a new squad of 15 young researchers to master high-level mathematics, algorithm design, scientific computation and software development, and to solve optimization problems for real-world applications. It will advance the research on algebraic methods for moment approaches, tackle mixed integer non-linear optimization problems and enhance the efficiency and robustness of moment relaxation methods. Specific applications of these approaches to optimization problems are related to smarter cities challenges, such as water distribution network management, energy flow in power systems, urban traffic management, as well as to oceanography and environmental monitoring and finance.

GRAPES

GRAPES project on cordis.europa.eu

  • Title:
    learninG, pRocessing, And oPtimising shapES
  • Duration:
    From December 1, 2019 to May 31, 2024
  • Partners:
    • INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET AUTOMATIQUE (INRIA), France
    • UNIVERSITA DEGLI STUDI DI ROMA TOR VERGATA (UNITOV), Italy
    • RHEINISCH-WESTFAELISCHE TECHNISCHE HOCHSCHULE AACHEN (RWTH AACHEN), Germany
    • ATHINA-EREVNITIKO KENTRO KAINOTOMIAS STIS TECHNOLOGIES TIS PLIROFORIAS, TON EPIKOINONION KAI TIS GNOSIS (ATHENA - RESEARCH AND INNOVATION CENTER), Greece
    • UNIVERSITAT LINZ (JOHANNES KEPLER UNIVERSITAT LINZ UNIVERSITY OF LINZ JOHANNES KEPLER UNIVERSITY OF LINZ JKU), Austria
    • SINTEF AS (SINTEF), Norway
    • VILNIAUS UNIVERSITETAS (Vilniaus universitetas), Lithuania
    • UNIVERSITA DELLA SVIZZERA ITALIANA (USI), Switzerland
    • UNIVERSITAT DE BARCELONA (UB), Spain
    • GEOMETRY FACTORY SARL, France
    • UNIVERSITY OF STRATHCLYDE, United Kingdom
  • Inria contact:
    Laurent Buse
  • Coordinator:
  • Summary:

    GRAPES aims at considerably advancing the state of the art in Mathematics, Computer-Aided Design, and Machine Learning in order to promote game changing approaches for generating, optimizing, and learning 3D shapes, along with a multisectoral training for young researchers. Recent advances in the above domains have solved numerous tasks concerning multimedia and 2D data. However, automation of 3D geometry processing and analysis lags severely behind, despite their importance in science, technology and everyday life, and the well-understood underlying mathematical principles. The CAD industry, although well established for more than 20 years, urgently requires advanced methods and tools for addressing new challenges.

    The scientific goal of GRAPES is to bridge this gap based on a multidisciplinary consortium composed of leaders in their respective fields. Top-notch 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 Network-wide workshops.

    Horizontally, our results lead to open-source, 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 visualisation, urban planning and autonomous driving.

10 Dissemination

Participants: Laurent Busé, Carles Checa, Ioannis Emiris, Evelyne Hubert, Pablo Gonzalez Mazon, Angelos Mantzaflaris, Michelangelo Marsala, Bernard Mourrain.

10.1 Promoting scientific activities

10.1.1 Scientific events: organisation

  • Evelyne Hubert was a co-organizer of the conference Symmetry, Stability and Computation held in CIRM (Marseille)
  • Laurent Busé, Carlès Checa and Pablo Gonzalez-Mazon co-organized a mini-symposium on «Castelnuovo-Mumford Regularity and Applications», at the SIAM Conference on Applied Algebraic Geometry, Eindhoven, The Netherlands, July 10-14, 2023.
  • Angelos Mantzaflaris co-organized with Felix Scholz the mini-symposium “Deep learning for geometric design” at the SIAM Conference on Computational Geometric Design (GD23), Genova, Italy, July 3-7, 2023.

    Angelos Mantzaflaris co-organized with Victor Calo, Pablo Antolin and Mattia Tani the mini-symposium “Fast formation and solution techniques for large-scale IGA” 11th International Conference on Isogeometric Analysis (IGA 2023) in Lyon, France, on 18-21 June 2023.

    Angelos Mantzaflaris co-organized with Mattia Tani, John Evans and Stefan Takacs the mini-symposium “Efficient methods for Isogeometric Analysis” at the 10th International Congress on Industrial and Applied Mathematics (ICIAM), Tokyo, Japan, August 20-24, 2023.

  • Ioannis Emiris has organized the GRAPES ITN Software Workshop in February 2023 in Athens, Greece.
General chair, scientific chair
  • Evelyne Hubert is chair of the society Foundation Computational Mathematics for the period 2023-2026.

10.1.2 Scientific events: selection

Member of the conference program committees
  • Bernard Mourrain, Laurent Busé and Angelos Mantzaflaris were members of the program committee of GMP23 and SPM 2023.
  • Angelos Mantzaflaris was a member of the program committees of the Isogeometric Analysis (IGA) 2023 in Lyin, France.
  • Evelyne Hubert was on the selection committees for the plenary speakers of the FoCM 2023 conference and the Steven Smale prize.
Reviewer
  • Bernard Mourrain and Angelos Mantzafalris reviewed submissions for the conferences GMP23, SPM 23.
  • Laurent Busé reviewed submissions for the conference GMP23, SPM 23 and SIGGRAPH ASIA 2023.

10.1.3 Journal

Member of the editorial boards
  • Laurent Busé is an associated editor of the Journal of Pure and Applied Algebra, and of Maple Transactions.
  • Evelyne Hubert is associated editor for the journal Foundations of Computational Mathematics and the Journal of Symbolic Computation. She is a co-editor of the Special Issue on Computational Geometry and Algebra in honor of the memory of A. Szanto.
  • Bernard Mourrain is associated editor for the Journal of Symbolic Computation.
  • Ioannis Emiris is associate editor of Journal of Symbolic Computation, Elsevier, Mathematics for Computer Science, Springer-Birkhauser, and Journal of Applicable Algebra in Engineering, Communication & Computing (Springer).
Reviewer - reviewing activities
  • Bernard Mourrain reviewed submissions for the following journals: Journal of Algebra and Applications, Journal of Symbolic Computation, Mathematics of Computation, Mathematical Programming, Journal of Computational and Applied Mathematics, SIAM Journal on Applied Algebra and Geometry, Applied Mathematics and Computation, Foundations of Computational Mathematics, Journal of Computer-Aided Design, Journal of Computer Aided Geometric Design, International Journal of Computer Vision.
  • Evelyne Hubert reviewed submissions for the following journals: Foundations of COmputational Mathematics, Journal of Symbolic Computatiom, SIAM Journal on Applied Mathematics, Computational and Applied Mathematics, Advances in Geometry, Computer Aided Geometric Design.
  • Laurent Busé reviewed submissions for the following journals: Journal of the London Mathematical Society, Foundations of COmputational Mathematics, Mathematics of Computation, Journal of Computational Algebra, Examples and Counterexamples, Journal of Algebraic Combinatorics, Journal of Algebra, Journal of Symbolic Computation, Computer Aided Geometric Design.
  • Angelos Mantzaflaris reviewed submissions for the following journals: Computer Methods in Applied Mechanics and Engineering, Computer Aided Geometric Design, Applied Numerical Mathematics, Journal of Symbolic Computation, Applied Mathematics and Computation, Engineering with Computers.

10.1.4 Invited talks

  • Bernard Mourrain was an invited speaker at the Workshop Real Algebraic Geometry with a View toward Koopman Operator Methods (13-17 March, Oberwolfach, Germany), Foundation of Computational Mathematics conference, Computational Algebraic Geometry (15-17 June 2023, Paris), Huawei Lavender Summit 2023 - Towards 2030 Wireless Communication Workshop (5-7 Jul. 2023, Cannes), SIAM conference on Applied Algebra and Geometry, Computational Real Algebraic Geometry (10-14 Jul. 2023, Eindhoven, The Netherlands), Workshop Computational and Applied Mathematics 2023 (28 Aug. - 1 Sep. 2023, Selva di Fasano, Italy), Recent Trends in Computer Algebra, Geometry of Polynomial System Solving, Optimization and Topology (16-20 Oct. 2023, IHP, Paris), POP23 - Future Trends in Polynomial OPtimization (13-17 Nov. 2023, Toulouse).
  • Evelyne Hubert was an invited speaker at the conferences New Directions in Real Algebraic Geometry (19-25 March, Oberwolfach), Random Algebraic Geometry (16-21 April BIRS Banff), Computational Algebraic Geometry (15-17 June Sorbonne Université) Symbolic Analysis (FoCM) (19-22 June Sorbonne Université), Equations Fonctionnelles et Interactions (3-5 July Université de Rennes), Symmetry, Stability and Computation (13-17 November CIRM Marseille), Géométrie Differentielle et Mécanique (22-24 novembre Sorbonne Université).
  • Laurent Busé gave an invited talk the Algebra & Geometry seminar at the University of Genova, March 15th.
  • Angelos Mantzaflaris gave an invited talk at the 20 Years Johann Radon Institute Workshop, Linz, Austria (March 2023).
  • Ioannis Emiris was an invited speaker at the Graph rigidity Workshop (Lancaster University, UK April 2023), the Laszlo Endrenyi memorial meeting on Bioavailability and Bioequivalence (Univ. of Athens, Greece, May 2023), the Dagstuhl seminar on Computational geometry (Germany, May 2023), the SIAM Conf. on Algebraic Geometry (Eindhoven, Holland, July 2023), the Workshop on Online Algorithms, Learning and Games (Koutsoupias Fest, Athens, July 2023), Recent Trends in Computer Algebra, Institut Henri Poincaré (IHP, Paris, France, October 2023), The Cyprus Institute Colloquim (Nicosia, Cyprus, November 2023).
  • Michelangelo Marsala was an invited speaker in the minisymposia at the GIMC SIMAI YOUNG Workshop (29-30 Sep. 2023, Pavia, Italy), SIAM Conference on Computational Science and Engineering (27 Feb - 3 Mar. 2023, Amsterdam, The Netherlands), International Conference on Approximation Theory and Beyond (15-19 May 2023, Nashville, TN, USA), 11th International Conference on Isogeometric Analysis (18-21 Jun. 2023, Lyon, France), SIAM Conference on Computational Geometric Design (3-7 Jul, Genova, Italy), Advances in Computational Mechanics (22-25 Oct. 2023, Austin, TX, USA)

10.1.5 Leadership within the scientific community

  • Ioannis Emiris is the coordinator of the European MCSA network GRAPES (learninG, pRocessing, And oPtimising shapES).
  • Bernard Mourrain is the coordinator of the European MCSA network POEMA (Polynomial Optimization, Efficiency through Moments and Algebra).

10.1.6 Scientific expertise

Bernard Mourrain was

  • a member of the committee of evaluation of LIP6 for HCERES (November 2023),
  • reviewer of a proposal for the Swiss National Science Foundation,
  • reviewer for a promotion at the Academy of Mathematics and Systems Science (Chinese Academy of Sciences),
  • reviewer for a tenured track promotion in a US university.

Evelyne Hubert was requested for

  • evaluation of a grant proposal for the Austrian Science Fund,
  • evaluation for a grant proposal for the Natural Sciences and Engineering Research Council of Canada,
  • recommendation for a professorship promotion in a US university.

Laurent Busé was

  • a member of the selection committee for CRCN positions at Inria Saclay (June),
  • chair of the selection committee for the Chair Professeur Junior (CPJ) on computational geometry at the laboratory of mathematics J.-A. Dieudonné of Université Côte d'Azur (UniCA).

10.1.7 Research administration

Angelos Mantzaflaris is a member of the Bureau of AMIES (Agence pour les Mathématiques en Interaction avec l'Entreprise et la Société) and member of the Comité du Centre of Inria d'Université Côte d'Azur in 2023.

Bernard Mourrain is member of the Bureau du Comité des Equipes Projets (BCEP).

Laurent Busé is co-chair, with Clément Pernet, of the french computer algebra research group of CNRS (GT Calcul Formel du GDR Informatique-Mathématiques)

10.2 Teaching - Supervision - Juries

10.2.1 Teaching

  • License : Martin Jalard, Analyse, 64h (L1), Polytech Nice Sophia - Université Côte d'Azur
  • Licence : Ioannis Emiris, Algorithms and complexity, 52 h (L2), NKU Athens
  • Licence : Ioannis Emiris, Software development, 26 h (L3), NKU Athens
  • Licence : Angelos Mantzaflaris, Mathematics 1, 64 h (L1), ESDHEM, SKEMA Business School in Sophia-Antipolis
  • 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 : Angelos Mantzaflaris, Geometric Modeling 6 h (M2) and Numerical Interpolation 20 h (M2), Polytech Nice Sophia - Univ Côte d'Azur.

10.2.2 Supervision

  • 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, co-supervised by Bernard Mourrain.
  • PhD: Pablo Gonzalez Mazon; Generation of valid high-order curved meshes. GRAPES Marie Skłodowska-Curie ITN, started in December 2020, defended December 2023, Inria, supervised by Laurent Busé.
  • PhD in progress: Mehran Hatamzadeh; An innovative gait analysis technology. PhD grant from the EU CoFUND BoostUrCareer program of UniCA, co-supervised by Laurent Busé and Raphaël Zory (LAMHESS, UniCA).
  • PhD in progress: Martin Jalard. Stratification of orbit space by orbit type : a constructive approach through equivariants. Funded for the most part by Ecole Normale Supérieure de Rennes through a contrat doctoral Inria. Started October 1st 2022, supervised by Evelyne Hubert.
  • PhD: Thomas Laporte, Towards a 4D model of the respiratory system. Fellowship from ED SFA/UniCA. Started on October 2019, defended in January 2023 33, co-supervised by Benjamin Mauroy (UniCA) and Angelos Mantzaflaris.
  • PhD: Riccardo DiDio, Building a digital twin of the respiratory system; study of the effects of inertial and constrictions on the airflow at mouth. Started on November 2019, defended in December 2023, co-supervised by Benjamin Mauroy (UniCA) and Angelos Mantzaflaris.
  • PhD: Michelangelo Marsala, Modelling and simulation using analysis-suitable subdivision surfaces and solids. GRAPES Marie Skłodowska-Curie ITN, started in November 2020, Defended in December 2023, Inria, supervised by Angelos Mantzaflaris and Bernard Mourrain.
  • PhD: Tong Zhao, Learning priors and metrics for 3D reconstruction of large-scale scenes. Started in October 2019, defended in March 2023; supervised by Pierre Alliez (Titane) and Laurent Busé.
  • PhD: Emmanouil Christoforou, since November 2021. Unsupervised learning algorithms in geometry. Supervised by Ioannis Emiris.
  • PhD: Carles Checa, since November 2020. Algebraic computing for geometric predicates. GRAPES MSCA ITN. Co-supervised by Ioannis Emiris and Bernard Mourrain.
  • PhD: Konstantinos Tertikas, since November 2020. Deep learning for representing, searching, generating 3D shapes. GRAPES MSCA ITN. Supervised by Ioannis Emiris.
  • PhD: Efi Malesiou, since March 2022. Machine learning for robust personalized long-term therapeutical treatment. Greek Ministy of education fellowship. Supervised by Ioannis Emiris.
  • PhD: Petros Stavropoulos, since April 2023. Reproducibility in science through information extraction from the scientific literature. Supervised by Ioannis Emiris.

10.2.3 Juries

Bernard Mourrain was

  • a referee and member of the jury for the Habilitation à Diriger des Recherches of Jeremy Berthomieu (LIP6) entitled Contributions to polynomial system solving : Recurrences and Gröbner bases.
  • a referee and member of the jury of the PhD thesis of Christof Vermeersch (Univ. Leuven) entitled The (Block) Macaulay Matrix: Solving Systems of Multivariate Polynomial Equations and Multiparameter Eigenvalue Problems

Evelyne Hubert

  • was a referee and member of the jury for the PhD of Christina Katsamaki, Sorbonne Université, entitled Exact Algebraic and Geometric Computations for Parametric Curves
  • chaired the jury for the PhD of Subhayan Saha, Ecole Normale de Lyon, Algebraic and Numerical Algorithms for Symmetric Tensor Decompositions

Laurent Busé

  • was a member of the jury for the PhD of Jorge García Fonán, Sorbonne Université, entitled Singularity and stability analysis of vision-based controllers and defended on January 24, 2023.
  • was a referee and member of the jury for the PhD of Rémi Prébet, Sorbonne Université, entitled Connexité dans les ensemble algébriques réels algorithmes et applications and defended on December 20th, 2023.

11 Scientific production

11.1 Major publications

  • 1 articleL.Lorenzo Baldi and B.Bernard Mourrain. On the Effective Putinar’s Positivstellensatz and Moment Approximation.Mathematical Programming, Series ASeptember 2022HALDOI
  • 2 articleE.Evangelos Bartzos, I. Z.Ioannis Z. Emiris and R.Raimundas Vidunas. New upper bounds for the number of embeddings of minimally rigid graphs.Discrete and Computational Geometry6832022, 796HALDOI
  • 3 articleL.Laurent Busé, Y.Yairon Cid-Ruiz and C.Carlos D'Andrea. Degree and birationality of multi-graded rational maps.Proceedings of the London Mathematical Society12142020, 743-787HALDOI
  • 4 articleL.Laurent Busé and A.Anna Karasoulou. Resultant of an equivariant polynomial system with respect to the symmetric group.Journal of Symbolic Computation762016, 142-157HALDOI
  • 5 articleA.Apostolos Chalkis, I. Z.Ioannis Z. Emiris and V.Vissarion Fisikopoulos. A Practical Algorithm for Volume Estimation based on Billiard Trajectories and Simulated Annealing.ACM Journal of Experimental Algorithmics28May 2023, 1-34HALDOI
  • 6 articleA. J.Alvaro Javier Fuentes Suárez and E.Evelyne Hubert. Scaffolding skeletons using spherical Voronoi diagrams: feasibility, regularity and symmetry.Computer-Aided Design102May 2018, 83 - 93HALDOI
  • 7 articleA.Alessandro Giust, B.Bert Jüttler and A.Angelos Mantzaflaris. Local (T)HB-spline projectors via restricted hierarchical spline fitting.Computer Aided Geometric Design80June 2020, 101865HALDOI
  • 8 articleE.Evelyne Hubert and E.Erick Rodriguez Bazan. Algorithms for fundamental invariants and equivariants: (of finite groups).Mathematics of Computation913372022, 2459-2488HALDOI
  • 9 articleZ.Zbigniew Jelonek and A.André Galligo. Elimination ideals and Bezout relations.Journal of Algebra5622020, 621-626HALDOI
  • 10 articleA.Angelos Mantzaflaris, B.Bert Jüttler, B.Boris Khoromskij and U.Ulrich Langer. Low Rank Tensor Methods in Galerkin-based Isogeometric Analysis.Computer Methods in Applied Mechanics and Engineering316April 2017, 1062-1085HALDOI
  • 11 articleB.Bernard Mourrain. Polynomial-Exponential Decomposition from Moments.Foundations of Computational Mathematics186December 2018, 1435--1492HALDOI
  • 12 inproceedingsK.Konstantinos Tertikas, D.Despoina Paschalidou, B.Boxiao Pan, J. J.Jeong Joon Park, M. A.Mikaela Angelina Uy, I. Z.Ioannis Z. Emiris, Y.Yannis Avrithis and L.Leonidas Guibas. PartNeRF: Generating Part-Aware Editable 3D Shapes without 3D Supervision.CVPR 2023 - IEEE/CVF Conference on Computer Vision and Pattern RecognitionProceedings of the 2023 IEEE/CVF Conference on Computer Vision and Pattern RecognitionVancouver, CanadaIEEEMarch 2023, 4466-4478HALDOI

11.2 Publications of the year

International journals

International peer-reviewed conferences

  • 23 inproceedingsC.Carlotta Giannelli, S.Sofia Imperatore, A.Angelos Mantzaflaris and D.Dominik Mokriš. Leveraging Moving Parameterization and Adaptive THB-Splines for CAD Surface Reconstruction of Aircraft Engine Components.Italian Chapter Conference 2023 - Smart Tools and Apps in GraphicsSTAG 2023 - annual international conference Smart Tools and Applications in Graphics - Eurographics Italian Chapter Conference 2023Matera, ItalyThe Eurographics Association2023, 125-134HALDOIback to text
  • 24 inproceedingsC.Carlotta Giannelli, S.Sofia Imperatore, A.Angelos Mantzaflaris and F.Felix Scholz. Learning meshless parameterization with graph convolutional neural networks.Lectures Notes in Networks and Systems2023 World Conference on Smart Trends in Systems, Security and SustainabilityLondon, United Kingdom2023HALback to text
  • 25 inproceedingsC.Christos Karampatzakis, A.Athanasios Mihailidis, A.Angelos Mantzaflaris and C.Christopher Provatidis. Contact simulation of tooth flanks using Isogeometric Analysis: Leveraging the power of high order, increased smoothness and curved geometries.VDI-BerichteInternational Conference on Gears 20232422Munich, GermanyVDI Verlag2023, 413-428HALDOIback to text
  • 26 inproceedingsL.Loukas Kavouras, K.Konstantinos Tsopelas, G.Giorgos Giannopoulos, D.Dimitris Sacharidis, E.Eleni Psaroudaki, N.Nikolaos Theologitis, D.Dimitrios Rontogiannis, D.Dimitris Fotakis and I. Z.Ioannis Z. Emiris. Fairness Aware Counterfactuals for Subgroups.NeurIPS 2023 - 37th Conference on Neural Information Processing SystemsProceedings Thirty-seventh Conference on Neural Information Processing SystemsNew-Orleans, Lousiane, United StatesJune 2023HALback to text
  • 27 inproceedingsK.Konstantinos Tertikas, D.Despoina Paschalidou, B.Boxiao Pan, J. J.Jeong Joon Park, M. A.Mikaela Angelina Uy, I. Z.Ioannis Z. Emiris, Y.Yannis Avrithis and L.Leonidas Guibas. PartNeRF: Generating Part-Aware Editable 3D Shapes without 3D Supervision.CVPR 2023 - IEEE/CVF Conference on Computer Vision and Pattern RecognitionProceedings of the 2023 IEEE/CVF Conference on Computer Vision and Pattern RecognitionVancouver, CanadaIEEEMarch 2023, 4466-4478HALDOIback to text
  • 28 inproceedingsM.Michele de Vita, C.Carlotta Giannelli, S.Sofia Imperatore and A.Angelos Mantzaflaris. Parameterization learning with convolutional neural networks for gridded data fitting.Lecture Notes in Networks and SystemsFICC 2024 - Future of Information and Communication ConferenceBerlin, Germany2023HALback to text
  • 29 inproceedingsT.Tong Zhao, L.Laurent Busé, D.David Cohen-Steiner, T.Tamy Boubekeur, J.-M.Jean-Marc Thiery and P.Pierre Alliez. Variational Shape Reconstruction via Quadric Error Metrics.SIGGRAPH 2023 - The 50th International Conference & Exhibition On Computer Graphics & Interactive TechniquesLos Angeles, United StatesAugust 2023HALDOIback to text

Scientific books

Scientific book chapters

  • 32 inbookS.Soodeh Habibi, M.Michal Kočvara and B.Bernard Mourrain. Polynomial optimization in geometric modeling.Polynomial Optimization, Moments, and ApplicationsSpringer2023, 163-186HALback to text

Doctoral dissertations and habilitation theses

  • 33 thesisT.Thomas Laporte. From CT-Scan to numerical simulations : building of a personalized model of the lung envelope and of the whole bronchial tree, application to freediving.Université Côte d'AzurJanuary 2023HALback to text

Reports & preprints

11.3 Cited publications

  • 43 articleM.Mehran Hatamzadeh, A.Ali Sharifnezhad, R.Reza Hassannejad and R.Raphael Zory. Discriminative sEMG-based features to assess damping ability and interpret activation patterns in lower-limb muscles of ACLR athletes.Biomedical Signal Processing and Control83May 2023, 104665HALDOIback to text