The Software and Industrial Workshop II (SIW2) is “GRAPES developer days”, allowing for senior researchers and ESRs to partake in coding sessions, completing small coding tasks and learning from each other. These sessions aim to facilitate knowledge transfer towards first solutions of real-world problems posed by our industrial partners.
Confirmed invited speakers include: RISC-Software, Tutorial: Machine Learning for metrology (XMANAI team, Athena RC), Leonidas Guibas (tentative) (Stanford University), André Lieutier (Geometry Engineering, Dassault Systèmes).
The ESR Days is an important aspect of networking and training, collocated with SIW2. The ESRs, through their organisation, coordinated by their representatives and assisted by the Coordinator and the organisers of SIW2, shall have to make choices and assign roles to each other. They shall choose GRAPES and external speakers, shall organise oral or poster presentations by themselves, undertake public engagement activities, and shall interact so as to foster ties useful for their future careers.
The meeting will take place physically at the Fresh Hotel in downtown Athens, see the map below. The hotel is close to the Athens City Hall and within walking distance form the Omonia square and Monastiraki square and the corresponding metro stations. Lunches and coffee breaks will be served on site.
Schedule of the event
The event will start on Tuesday morning and will close on Friday early afternoon. A preliminary program can be found below.
Parallel WP sessions
Educational Committee meeting
Supervisory Board meeting
Invited talks and tutorials
André Lieutier (Geometry Engineering) – From topological inference to meshing algorithms
In this tutorial we consider the general problem of retrieving a topologically correct model from partial geometric information. Applications areas span from automated 3D scene faithful reconstruction from scanned points clouds to abstract submanifolds inference from high dimensional data in the context of machine learning.
In a first part, the tutorial will recall several topological notions and consider their possible inference from data. In a second part, we focus on the computation of simplicial complexes homeomorphic to manifolds known through point samples, in other words, manifolds triangulation algorithms with topological guarantees.
Directions from/To the Athens International Airport “E. Venizelos”
The best option is to use Metro Line 3 (subway, blue line), and get off at “Syntagma” stop. Timetable: 05:30-24:00. Frequency: 36′
Taxis: the taxi stand extends from Door 4 to Door 1 at the Arrivals Level. Fee: 38 or 54 Euro at daytime or evenings, resp.
Express Airport Bus X95 (to Syntagma sq.): Timetable. Runs on a 24-hour basis. The duration of the route from and to the Eleftherios Venizelos Airport depends on traffic conditions and on the time of use of the Express bus line. Εstimated duration of route for X95 bus line is approximately 60’.
Directions to the venue (“Fresh Hotel”):
From the airport take Metro Line 3 (blue line) or Express Airport Bus X95 or a taxi and at Syntagma Metro Station change to Metro Line 2 (red line). Get off at “Omonia” stop. From there it is a 6′ walk (500m) to the venue.
ALternatively, from the airport take Metro Line 3 (blue line) and get off at “Monastiraki” stop. From there it is a 8′ walk (700m) to the venue.
Schedule of the metro and tram can be found here and of the buses here.
We have reserved a few rooms at the Fresh hotel where the workshop will take place. Here are some other suggestions for lodging in Athens in the vicinity of the workshop’s venue:
The GRAPES Network has four (4) SHORT TERM (4-9 months) open research training positions at Athena RC (Greece), Inria Research Center -Sophia Antipolis Mediterranee (France), University of Strathclyde (UK) and Università della Svizzera italiana (Switzerland). For details see http://grapes-network.eu/phd-positions/, or contact directly the corresponding advisors.
Tor Dokken will present a poster showcasing GRAPES and the results obtained so far by our ESRs at the I4MS event ‘European leadership and digital game-changers of the manufacturing industry’ in Budapest, on 19th of October 2022. Many thanks to the fellows that contributed their results and to Georg Muntingh for compiling them to a poster!
The aims of the event are to bring researchers closer to the general public and to increase public awareness of research and innovation activities with a view to supporting the public recognition of research and innovation activities, to creating an understanding of the impact of researchers’ work on citizens’ daily lives and to encouraging young people to embark on research careers.
The event is therefore particularly aimed at young people and their families.
The organisers of these events always seek the involvement of researchers funded by H2020 and Horizon Europe, including the Marie Skłodowska-Curie actions (MSCA).
The participation of the GRAPES fellows in the event is therefore strongly encouraged. The link European Researchers’ Night includes the websites of each organisers that the fellow can contact directly and offer his/her participation.
We are very hapy to welcome to GRAPES the Cognition and Computation in Design (Co-Design) Lab (https://codesign.berkeley.edu/) of the Department of Mechanical Engineering, University California Berkeley! The Co-Design Lab, led by Prof. Kosa Goucher-Lambert has just joined our Project as a Partner Organization.
The second Doctoral School will take place *physically* from Monday, June 13 to Friday, June 17 at the Università della Svizzera italiana (USI) in Lugano. It will start with lunch on Monday and end after lunch on Friday.
The event will focus on *CAD* in its broadest sense, covering the historical development of the field, latest research results, applications, and the relevance for industry. It features keynote talks by GRAPES members and invited high-profile speakers, short presentation by the ESRs, work package meetings, supervisory board and educational committee meetings, and social activities.
Confirmed speakers include Fehmi Cirak (U. Cambridge), Adi Levin (Meta), and Ulrich Reif (TU Darmstadt).
You can download a pdf version of the programme here and the book of abstracts here.
To register for the workshop, just send an e-mail with your name, affiliation, arrival and departure date, and dietary restrictions to firstname.lastname@example.org before *May 20, 2022*.
The main event will happen at room D1.14, on the first floor (one above the ground floor) of the D sector of the East Campus building. The registration desk will be on the ground floor of the D sector.
To travel to Lugano, you will fly into Zurich or Milano Malpensa airport and then take a train (approximately 2 hours) to Lugano. Take bus #6 (direction Cornaredo) to get from the train station to the university (get off at stop “Campus Universitario”).
Lugano has numerous hotels, but we recommend the following hotels:
The first Software and Industrial Workshop will offer intersectoral training in technological aspects with links to the industry. It will include lectures and demo sessions so that the ESRs can acquire hands-on experience with software developed within the Network on computational geometry and geometry processing (CGAL), isogeometric analysis (G+Smo) and machine learning (PyTorch). Industrial partners will also be invited to give talks and discuss some of their open problems.
The meeting will take place physically at the Inria research center at Sophia Antipolis on December 6-10, 2021. The last two days consist of the annual G+Smo Developer Days workshop. Address : 2004 Route des Lucioles, 06902 Valbonne, France
The social dinner will take place on Wednesday evening in Antibes, at 19:30, at restaurant Alcyon (13 Av. du 11 Novembre). See also this link.
For people going back to Nice note that the last train to Nice is at 22:00.
The agenda of the workshop includes: Tutorials on Software, Industrial presentations, Contributed presentations, and ESR presentations.
Dominik Mokris (MTU) – Geometric modelling in the CAE pipeline of MTU Aero Engines AG
Together, we will have a look at how geometry is handled in the CAE pipeline at MTU. I will show you what our geometry generator does, its scope, inputs and outputs as well as why and how it has become indispensable for the design phase of aircraft engines. After discussing some architectural decisions and implementation specifics you will hopefully be better equipped to see whether developing such a tool in your organisation would make sense or not. In the second half, we will turn our attention to the problem of using 3D scans of physical parts as inputs into the process. Based on the examples of some of the recent [1, 2] and less recent  scientific results you will get an idea of how we work and how long does it take to integrate a scientific idea into day-to-day engineering practice. Joint work with David Großmann (MTU Aero Engines AG, Germany).
 Lisa Groiss, Bert Jüttler and Dominik Mokriš. 27 variants of Tutte’s theorem
for plane near-triangulations and an application to periodic spline surface
fitting. Comput. Aided Geom. Design, 85, 101975, 2021.
 Cesare Bracco, Carlotta Giannelli, David Großmann, Sofia Imperatore, Do-
minik Mokriš and Alessandra Sestini. THB-spline approximations for turbine blade design with local B-spline approximations. Accepted to SEMA-SIMAI granted to MACMAS 2019. Springer.
 Gábor Kiss, Carlotta Giannelli, Urška Zore, Bert Jüttler, David Großmann
and Johannes Barner. Adaptive CAD model (re-) construction with THB-
splines. Graphical models, 76(5):273–288, 2014.
Mark Gammon (ITI) – B-rep, Subdivision and Facet – Towards a Hybrid Geometry Engine
The 3D digital geometry landscape has changed rapidly in recent years with novel forms of geometry reaching industrial maturity. High fidelity 3D scans are now commonplace. Super-flexible subdivision surface geometry powers intuitive free-form design tools. Lattice-based implicit surface geometry driven by simulation algorithms create designs unimaginable with traditional MCAD tools. This rapid expansion of the industrial geometry frontier is exciting and promising, but geometry exchange weaknesses are limiting opportunities. In this talk I will describe progress ITI is making towards developing a multi-representational geometry engine designed to support alternative, but concurrent, geometrical forms of the same master design. An overview of the complex transformations and intelligent associativity between the various representations will be described. Some industrial examples of moving a 3D design between different forms to enable the most optimal representation to be used for a particular operation will be shown.
Matthias Fey (TU Dortmund) – Graph Neural Networks within PyTorch Geometric: Applications, Implementation and Scalability
Graph Neural Networks (GNNs) recently emerged to a powerful approach for representation learning on relational data such as social networks, molecular graphs or geometry. Similar to the concepts of convolutional and pooling layers on regular domains, GNNs are able to (hierarchically) extract localized embeddings by passing, transforming, and aggregating information between nodes. In this lecture, I will provide a broad overview of this highly active research field, and will cover relevant topics such as scalability and applications based on GNNs. In a hands-on session, you will learn to implement and train GNNs from scratch using the PyTorch Geometric library [https://github.com/pyg-team/pytorch_geometric]. You will learn how to apply GNNs to your own problems and how PyTorch Geometric enables high GPU throughput on highly sparse and irregular data of varying size.
PyG was started under the name Pytorch Geometric by Matthias Fey and Jan Eric Lenssen during their PhDs at TU Dortmund University. After its release in November 2017, the state-of-the-art in Graph Neural Networks has been quickly growing in leaps and bounds over the past few years, including significant contributions from Stanford University such as GraphSAGE, the Open Graph Benchmark and GraphGym. To support the continuous appetite of new features and to ensure longevity, TU Dortmund and Stanford University have joined forces and are now committed to the maintenance and enhancement of the package in accordance with the latest trends in academic research. GNN tools from both parties are unified under the umbrella PyG, which now involves a team of core developers extended by a community of over 170 contributors across the world.
Pierre Alliez (Inria), Andreas Fabri (GeometryFactory) – CGAL Tutorial
In this tutorial on CGAL library (https://www.cgal.org/) we will present various CGAL packages from the user perspective, that is we will not present the theory behind the algorithms and data structures. By “CGAL packages”, we mean the different chapters in the online manual as you can find them in the CGAL Package Overview page: https://doc.cgal.org/latest/Manual/packages.html. The goal is to be able to write code during the workshop where you will use CGAL, whether you have never used CGAL before or are an advanced user.
Jean Feydy (Inria Paris) – Fast geometric learning with symbolic matrices
Sparse representations such as 3D point clouds have a key position in the computer vision toolbox. They complement bitmap images effectively, enabling fast geometric computations for e.g. shape registration. In this talk, I will present extensions for PyTorch, NumPy, Matlab and R that speed up fundamental computations on (generalized) point clouds by several orders of magnitude, compared with PyTorch, TF and JAX GPU baselines. These software tools allow researchers to break through major computational bottlenecks in the field and have been downloaded more than 100k times over the last few years. The presentation will be of interest to all researchers who deal with point clouds, time series and segmentation maps, with a special focus on:
1. Fast and scalable computations with (generalized) distance matrices.
2. Efficient and robust solvers for the optimal transport (= “Earth Mover’s”) problem.
3. Applications to shape analysis and geometric deep learning, with a case study on the “pixel-perfect” registration of lung vessel trees.
The tutorial will be a hands-on introduction to fast geometric computations with the KeOps library. We will use Google Colab as a simple GPU sandbox.
– “Geometric data analysis, beyond convolutions”: https://www.jeanfeydy.com/geometric_data_analysis.pdf
– “Fast geometric learning with symbolic matrices”: http://jeanfeydy.com/Papers/KeOps_NeurIPS_2020.pdf
– “Accurate point cloud registration with robust optimal transport”: https://www.jeanfeydy.com/Papers/RobOT_NeurIPS_2021.pdf
– KeOps library (geometric computations): http://kernel-operations.io/keops/index.html
– GeomLoss library (optimal transport): https://www.kernel-operations.io/geomloss/
The tutorials on Thursday will give an overview of the open-source library “G+Smo” and its capabilities, as well as hands-on session in the afternoon. 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. The library aims at providing access to high quality, open-source software to the community of numerical simulation and beyond.
Contributed talks for G+Smo Developer Days
Andrea Farahat (RICAM, Austria)
Isogeometric analysis with C1-smooth functions over multi-patch surfaces.
A framework for the construction of a C1-smooth isogeometric spline spaces over particular G1 multi-patch surfaces, called analysis-suitable G1 (AS-G1) surfaces , will be presented. The class of AS-G1 multi-patch geometries is of importance since it includes exactly those G1-smooth multi-patch geometries which allow the design of C1-smooth isogeometric spline spaces with optimal approximation properties. The method extends the construction for AS-G1 planar multi-patch parametrizations to the AS-G1 multi-patch surface case. We also generate for the C1-smooth isogeometric spline space a local basis which is used to solve, by a standard Galerkin approach, the biharmonic equation — a particular fourth order partial differential equation — over several AS-G1 multi-patch surfaces. The obtained numerical results exhibit optimal convergence orders in the L2, H1 and H2 norm, and demonstrate the potential of our C1-smooth isogeometric spline functions for solving fourth order partial differential equations over multi-patch surfaces.
Pascal Weinmüller (JKU Linz) – Solving fourth order equations on multpatches with isogeometric analysis in G+Smo
Isogeometric analysis (IGA) is a numerical method that uses spline-based geometry parameterizations to solve partial differential equations (PDEs). Since IGA is based on B-splines, it is simple to achieve high order smoothness within a single patch. However, to represent more complex geometries one usually uses a multi-patch construction. In this case, the global continuity for the basis functions is in general only C^0. Therefore, for C^1-smooth isogeometric functions, a special construction for the basis is needed. Such spaces are of interest when solving numerically fourth-order PDE problems, such as the biharmonic equation or Kirchhoff-Love plate/shell formulations, using an isogeometric Galerkin method.
There are several different methods to overcome the problem of reduced smoothness. One way is to use isogeometric spaces that are globally C^1 over multi-patch domains where different constructions are studied in [1,2,3].
However, each approaches need a different construction, but end with the same problem: solving the PDEs. In this talk, a method in G+Smo is shown which handles the different spaces by using the class “gsMSpline” for constructing the spaces. The main idea is that it uses basis transformation for describing the space which makes it in general simple for any bases. The handling of the class is kept straightforward and we show an example with different basis constructions which can be used for solving the biharmonic equation on multi-patch domains.
Joint work with Hugo Verhelst, TU Delft, Netherlands, Andrea Farahat, RICAM Linz, Austria and Angelos Mantzaflaris, Inria, France
 Hughes, T. JR and Sangalli, G. and Takacs, T. and Toshniwal, D., Smooth multi-patch discretizations in Isogeometric Analysis, Handbook of Numerical Analysis, (2021) Vol 22: 467–543.
 Kapl, M., Sangalli, G., Takacs, T. An isogeometric C^1 subspace on unstructured multi-patch planar domains. Computer Aided Geometric Design, (2019) Vol 69: 55–75.
 Weinmüller, P. and Takacs T., Construction of approximate C^1 bases for isogeometric analysis on two-patch domains, Computer Methods in Applied Mechanics and Engineering, (2021) Vol 385: 114017.
Sofia Imperatore (U. Florence) – Data-driven spline parametrization and fitting problems
Free-form curve fitting plays an important role in Computer Aided Geometric Design (CAGD) and Geometric Modelling. In particular, in many applications geometric models are built on measured data and subsequently employed for further design. A necessary pre-processing step to fitting a given sequence of points with spline curves is to first compute suitable parameterizations. The talk will investigate different neural network architectures to address this problem for different kinds of point sequences.
Physics-informed neural networks (PINNs) are emerging technologies that strive to revolutionize the field of computer-aided analysis of scientific and engineering problems by directly exploiting physical laws to drive network optimization. In this short talk we propose a novel approach to embed the PINN paradigm into the framework of Isogeometric Analysis. In contrast to classical PINNs which predict point-wise solution values to (initial-)boundary-value problems directly, our IGA-PINNs learn solutions in terms of their expansion coefficients relative to a given B-Spline basis. This approach is furthermore used to encode the geometry and other problem parameters such as boundary conditions and feed them into the network as inputs, which allows the user to analyze different problem configurations effectively.
Felix Weber (RWTH Aachen) – Isogeometric Analysis for Micromechanical Analysis of Cast Iron
Material properties such as fatigue strength highly depend on the local microstructure. With increasing complexity of microstructure models finite element based approaches become more costly. Therefore, IGA-suitable 2D-modelling of microstructures becomes relevant and allows for a decrease in required degrees of freedom making further analysis feasible.
Hugo Verhelst (TU Delft) – A Parallel Adaptive Arc-Length Method
Parallel computing is omnipresent in today’s scientific computing landscape starting at multicore processors in desktop computers up to massively parallel clusters. While domain decomposition methods have a long tradition in computational mechanics to decompose spatial problem into many subproblems that can be solved in parallel, advancing solution schemes for dynamics or quasi-statics are inherently serial processes. Methods like Parareal and Multigrid-reduction-in-time (MGRIT) are based on a multigrid approach over the temporal domain, enabling parallelization in time. These techniques have been applied successfully in many fields of study, including power network analysis , cardiac fluid-structure interaction  or the training of neural networks .
As MGRIT is based on dynamic, their temporal multigrid scheme is based on a time discretization error that evolves over time. For quasi-static simulations, solutions at different time steps are independent so that the time discretization error and thus the motivation for MGRIT vanishes.
In this talk, we will present an MGRIT-inspired approach to parallelize quasi-static computations. Due to the parametrization of the arc-length instead of time, the multi-level error for the arc-length parametrization is formed by the load parameter and the solution norm; enabling to measure errors like MGRIT. By applying local refinements in the arc-length parameter, this multi-level adaptive arc-length method refines solutions where the non-linearity in the load-response space is maximal, i.e. where the path is curved. To reduce the number of levels needed for reference errors, a higher-dimensional spline is adaptively fitted through the load-response solution space. Furthermore, knot insertion algorithms will assure that bifurcations can nicely be embedded in the presented method.
 S. Günther, R. D. Falgout, P. Top, C. S. Woodward and J. B. Schroder, “Parallel-in-Time Solution of Power Systems with Unscheduled Events,” 2020 IEEE Power & Energy Society General Meeting (PESGM), 2020, pp. 1-5, doi: 10.1109/PESGM41954.2020.9281595.
 Hessenthaler, A., Falgout, R. D., Schroder, J. B., de Vecchi, A., Nordsletten, D., & Röhrle, O. (2021). Time-periodic steady-state solution of fluid-structure interaction and cardiac flow problems through multigrid-reduction-in-time. arXiv preprint arXiv:2105.00305.
 Cyr, E. C., Günther, S., & Schroder, J. B. (2019). Multilevel initialization for layer-parallel deep neural network training. arXiv preprint arXiv:1912.08974.
The meeting will take place at the Inria research center at Sophia Antipolis: Inria Sophia Antipolis Méditerranée 2004, route des Lucioles BP 93 06902 Sophia Antipolis Cedex. You can find travel directions to Inria here.
Here are some suggestions for lodging in Antibes with an easy access by bus to Sophia Antipolis by bus.