Keywords
Computer Science and Digital Science
 A1.2.4. QoS, performance evaluation
 A6.1.4. Multiscale modeling
 A6.2.3. Probabilistic methods
 A8.1. Discrete mathematics, combinatorics
 A8.2. Optimization
 A8.3. Geometry, Topology
 A8.6. Information theory
 A8.7. Graph theory
 A8.8. Network science
 A8.9. Performance evaluation
 A9.2. Machine learning
 A9.7. AI algorithmics
Other Research Topics and Application Domains
 B4.3. Renewable energy production
 B6.2.2. Radio technology
 B6.3.4. Social Networks
1 Team members, visitors, external collaborators
Research Scientists
 Bartlomiej Blaszczyszyn [Team leader, INRIA, Senior Researcher, HDR]
 François Baccelli [INRIA, Senior Researcher, HDR]
 Antoine Baker [Inria, Starting Research Position]
 Ana Busic [INRIA, Researcher]
 Christine Fricker [INRIA, Researcher, HDR]
 Ali Khezeli [INRIA, Starting Research Position]
 Marc Lelarge [INRIA, Senior Researcher, HDR]
 Laurent Massoulié [INRIA, Senior Researcher, HDR]
 Sean Meyn [ECE FLORIDA, from Jun 2022 until Aug 2022, Inria International Chair]
 Kevin Scaman [INRIA, Researcher]
PostDoctoral Fellow
 Ke Feng [INRIA]
PhD Students
 Claire Bizon Monroc [INRIA]
 Matthieu Blanke [INRIA]
 Romain Cosson [Inria, from Sep 2022]
 Michel Davydov [INRIA]
 Bastien Dubail [ENS PARIS]
 Mathieu Even [INRIA]
 Roman Gambelin [INRIA]
 Sayeh Khaniha [INRIA]
 Thomas Le Corre [PSL University]
 Maxime Leiber [SAFRAN]
 Jakob Maier [Inria, from Oct 2022]
 Pierre Popineau [INRIA]
 David Robin [INRIA]
 Bharath Roy Choudhury [INRIA]
 Ilia Shilov [INRIA]
 Guodong Sun [NOKIA]
 Amaury Triboulin [INRIA]
 Lucas Weber [DGA]
Technical Staff
 Thibault La BatideAlanore [INRIA, Engineer, from Nov 2022, Prethèse]
Interns and Apprentices
 Alessia Rigonat [Inria, from Apr 2022 until Oct 2022]
 Solène DelannoyPavy [Inria, from Mar 2022 until Jul 2022]
 Thibault La BatideAlanore [Inria, from Feb 2022 until Jul 2022]
 Jakob Maier [Inria, from Apr 2022 until Sep 2022]
Administrative Assistants
 Helene Bessin Rousseau [INRIA]
 Julien Guieu [Inria, for ERC Nemo]
 Helene Milome [INRIA]
Visiting Scientist
 Mir Omid Haji Mirsadeghi [SUT SHARIF]
External Collaborators
 Charles Bordenave [CNRS, from Jul 2022]
 Pierre Bremaud [Emeritus EPFL, HDR]
 MarcOlivier Buob [Nokia Bell Labs]
 Serguei Foss [UNIV HERIOTWATT, from May 2022]
 Alexandre Gaudilliere [CNRS, from Mar 2022]
 Subhroshekhar Ghosh [UNIV SINGAPOUR, from Jun 2022]
 Thomas Hutchcroft [UNIV CAMBRIDGE]
 Sanket Kalamkar [Qualcomm]
 Panogiotis Konstantopoulos [UNIV LIVERPOOL, from Feb 2022]
 Fabien Mathieu [SWAPCARD, HDR]
 Samuel Mellick [ENS DE LYON, from Feb 2022]
 Elizabeth O'Reilly [CALTECH]
2 Overall objectives
The general scientific focus of DYOGENE is on the development of network mathematics. The following theories lie within our research interest: dynamical systems, queuing theory, optimization and control, information theory, stochastic processes, random graphs, stochastic geometry.
Our theoretical developments are motivated by and applied in the context of communication networks (Internet, wireless, mobile, cellular, peertopeer), social and economic networks, power grids, and, recently, infectious diseases.
We collaborate with many industrial partners. Our current industrial relations involve EDF, Huawei, Microsoft, Nokia, Orange, Safran.
More specifically, the scientific focus of DYOGENE defined in 2013 was on geometric network dynamics arising in communications. By geometric networks we understand networks with a nontrivial, discrete or continuous, geometric definition of the existence of links between the nodes. In stochastic geometric networks, this definition leads to random graphs or stochastic geometric models.
A first type of geometric network dynamics is the one where the nodes or the links change over time according to an exogeneous dynamics (e.g. node motion and geometric definition of the links). We will refer to this as dynamics of geometric networks below. A second type is that where links and/or nodes are fixed but harbor local dynamical systems (in our case, stemming from e.g. information theory, queuing theory, social and economic sciences). This will be called dynamics on geometric networks. A third type is that where the dynamics of the network geometry and the local dynamics interplay. Our motivations for studying these systems stem from many fields of communications where they play a central role, and in particular: message passing algorithms; epidemic algorithms; wireless networks and information theory; device to device networking; distributed content delivery; social and economic networks, neural networks, and power grids.
3 Research program
Here are our main research accesses.
 Models of infectious diseases 8.1
 Distributed network control and smartgrids 8.2
 Stochastic matching and queueing systems 8.3
 Reinforcement learning 8.4
 Mathematics of wireless cellular networks 8.5
 Highdimensional statistical inference and distributed learning 8.6
 Distributed optimization for machine learning 8.7
 Random Geometry 8.8
 Mathematics of stochastic networks via meanfield analysis 8.9
4 Application domains
4.1 Physical communication networks
Internet, wireless, mobile, cellular networks, transportation networks, distributed systems (cloud, call centers). In collaboration with Nokia Bell Labs and Orange Labs.
4.2 Abstract networks
Social interactions, human communities, economic networks.
4.3 Power grids
Energy networks. In collaboration with EDF and Vito (Belgium).
5 Social and environmental responsibility
6 Highlights of the year
6.1 Awards
 Solène DelannoyPavy received a research internship award (section Applied mathematics) from Ecole Polytechnique for her internship entitled "Designing call centers with two levels of service".
 Michel Davydov received the MNA award for the best junior researcher presentation during the workshop "Mathematical modeling and Statistical Analysis in Neuroscience” at the Institut Henri Poincaré, Paris France;
7 New software and platforms
7.1 New software
7.1.1 pp_syn

Name:
Particle gradient descent model for point process generation

Keywords:
2D, Point process generation, Gradient descent, Wavelet phase harmonics

Scientific Description:
This paper introduces a generative model for planar point processes in a square window, built upon a single realization of a stationary, ergodic point process observed in this window. Inspired by recent advances in gradient descent methods for maximum entropy models, we propose a method to generate similar point patterns by jointly moving particles of an initial Poisson configuration towards a target counting measure. The target measure is generated via a deterministic gradient descent algorithm, so as to match a set of statistics of the given, observed realization. Our statistics are estimators of the multiscale wavelet phase harmonic covariance, recently proposed in image modeling. They allow one to capture geometric structures through multiscale interactions between wavelet coefficients. Both our statistics and the gradient descent algorithm scale better with the number of observed points than the classical knearest neighbour distances previously used in generative models for point processes, based on the rejection sampling or simulatedannealing. The overall quality of our model is evaluated on point processes with various geometric structures through spectral and topological data analysis, compared in particular to [Tscheschel, Stoyan (2006) Statistical reconstruction of random point patterns].

Functional Description:
Code for the paper "Particle gradient descent model for point process generation" in Statistics and Computing 32 (3) June 2022 Statistics and Computing. Using Kymatio/phaseharmonics software.
 URL:
 Publication:

Contact:
Sixin Zhang

Participants:
Antoine Brochard, Sixin Zhang, Bartlomiej Blaszczyszyn
7.1.2 wavelettexturesynthesis

Name:
Generalized Rectifier Wavelet Covariance Model For texture Synthesis

Keywords:
Texture Synthesis, Wavelet phase harmonics, Generative Models

Scientific Description:
Stateoftheart maximum entropy models for texture synthesis are built from statistics relying on image representations defined by convolutional neural networks (CNN). Such representations capture rich structures in texture images, outperforming waveletbased representations in this regard. However, conversely to neural networks, wavelets offer meaningful representations, as they are known to detect structures at multiple scales (e.g. edges) in images. In this work, we propose a family of statistics built upon nonlinear wavelet based representations, that can be viewed as a particular instance of a onelayer CNN, using a generalized rectifier nonlinearity. These statistics significantly improve the visual quality of previous classical waveletbased models, and allow one to produce syntheses of similar quality to stateoftheart models, on both grayscale and color textures. We further provide insights on memorization effects in these models.

Functional Description:
Code for the paper: "Generalized Rectifier Wavelet Covariance Model For texture Synthesis" In ICLR 2022  10th International Conference on Learning Representation, France, April 2022 ICLR 2022
 URL:
 Publication:

Contact:
Sixin Zhang

Participants:
Antoine Brochard, Sixin Zhang
8 New results
Participants: All Dyogene.
8.1 Models of infectious diseases
1. A computational framework for evaluating the role of mobility on the propagation of epidemics on point processes 4 This paper is focused on SIS (SusceptibleInfectedSusceptible) epidemic dynamics (also known as the contact process) on populations modeled by homogeneous Poisson point processes of the Euclidean plane, where the infection rate of a susceptible individual is proportional to the number of infected individuals in a disc around it. The main focus of the paper is a model where points are also subject to some random motion. Conservation equations for moment measures are leveraged to analyze the stationary regime of the point processes of infected and susceptible individuals. A heuristic factorization of the third moment measure is then proposed to obtain simple polynomial equations allowing one to derive closed form approximations for the fraction of infected individuals in the steady state. These polynomial equations also lead to a phase diagram which tentatively delineates the regions of the space of parameters (population density, infection radius, infection and recovery rate, and motion rate) where the epidemic survives and those where there is extinction. A key takeaway from this phase diagram is that the extinction of the epidemic is not always aided by a decrease in the motion rate. These results are substantiated by simulations on large two dimensional tori. These simulations show that the polynomial equations accurately predict the fraction of infected individuals when the epidemic survives. The simulations also show that the proposed phase diagram accurately predicts the parameter regions where the mean survival time of the epidemic increases (resp. decreases) with motion rate.
2. Disparities in accessibility to oncology care centers in France 42 Background Cancer caused nearly 10 million deaths in 2020. While most of the ongoing research focuses on finding new treatments, accessibility to oncology care receives less attention. However, access to health services plays a key role in cancer survival. Spatial accessibility methods have been successfully applied to measure accessibility to primary care. Yet, little research to date focused on oncology care specifically. Methods We focused on all care centers with medicine, surgery, or obstetric activity in metropolitan France. We propose a clustering algorithm to automatically label the hospitals in terms of oncology specialization. Then, we computed an accessibility score to these hospitals for every municipality in metropolitan France. Finally, we proposed an optimization algorithm to increase the oncology accessibility by identifying centers which should increase their capacity. Results We labelled 1,662 care centers into 8 clusters. Half of them were eligible for oncology care and 118 centers were identified as experts. We computed the oncology accessibility score for 34,877 municipalities in metropolitan France. Half of the population lived in the top 20% accessibility areas, and 6.3% in the bottom 20% zones. Accessibility was higher near dense cities, where the experts care centers were located. By combining the care centers clusters and the accessibility distributions, our optimization algorithm could identify hospitals to grow, to reduce accessibility disparities. Conclusion Our method made it possible to quantify oncology care accessibility across all metropolitan France, as well as to make suggestions on where to increase hospital capacity to improve accessibility, especially in more populated suburban areas. Our approach was deliberately nonspecific to cancer type nor to the kind of stays, but it could be adapted to more specific scenarios. We packaged our method into a web application allowing the users to run the algorithms with various parameters and visualize the results. Highlights We computed the oncology accessibility score for 34,877 municipalities and highlighted disparities. Our optimization algorithm can identify hospitals to grow, to reduce accessibility disparities. We packaged our algorithms and results into a web application, opened to healthcare professionals.
3. Geographic and sociodemographic disparities in oncology care pathways 34 In France, during year 2018, there was 382,000 new cancer cases and 157,400 deaths. The 20142019 Cancer Plan sets new objectives for cancer care in France. In particular, objectives 2 and 7 emphasize the quality of the care pathways: they aim respectively to guarantee the qual ity and safety of care and ensure comprehensive and personalized care. In order to standardize the care pathways while personalizing care, care trajectories have been introduced. The definition of these optimal care trajectories is based on national and international good practice recommendations. We propose to study in detail the geographical and sociodemographic disparities in the care pathways of cancer patients in France. First, we will try to characterize the health care institutions in France based on their oncology activity. Then, we will study the distribution of these centers on the territory, in order to highlight possible disparities in access to them. Finally, we will try to propose a care center recommendation algorithm, using the French social security database (SNDS). This algorithm will aim at guiding patients towards the optimal center, in order to maximize the quality of the care pathways.
8.2 Distributed network control and smartgrids
4. An Equilibrium Analysis of RiskHedging Strategies in Decentralized Electricity Markets 50 We investigate equilibrium problems arising in various decentralized designs of the electricity market involving riskaverse prosumers. The prosumers have the possibility to hedge their risks through financial contracts that they can trade with peers or purchase from an insurance company. We build several market designs of increasing complexity, from a onestage market design with interagent financial contract trading to a Stackelberg game where an insurance company acts as a leader and prosumers are followers. We derive riskhedging pricing scheme for each model and show that the Stackelberg game pessimistic formulation might have no solution. We propose an equivalent reformulation as a parametrizated generalized Nash equilibrium problem, and characterize the set of equilibria. We prove that the insurance company can design price incentives that guarantee the existence of a solution of the pessimistic formulation, which is ε close to the optimistic one. We then derive economic properties of the Stackelberg equilibria such as fairness, equity, and efficiency. We also quantify the impact of the insurance company incomplete information on the prosumers' riskaversion levels on its cost and social cost. Finally, we evaluate numerically the proposed riskhedging market models, using residential data provided by Pecan Street.
5. Feature Projection for Optimal Transport 41 Optimal transport is now a standard tool for solving many problems in statistics and machine learning. The optimal "transport of probability measures" is also a recurring theme in stochastic control and distributed control, where in the latter application the probability measure corresponds to an empirical distribution associated with a large collection of distributed agents, subject to local and global control. The goal of this paper is to make precise these connections, which inspires new relaxations of optimal transport for application in new and traditional domains. The proposed relaxation replaces a target measure with a "moment class": a set of probability measures defined by generalized moment constraints. This is motivated by applications to control, outlier detection, and to address computational complexity. The main conclusions are (i) A characterization of the solution is obtained, similar to Kantorovich duality, in which one of the dual functions in the classical theory is replaced by a linear combination of the features defining the generalized moments. Hence the dimension of the optimization problem coincides with the number of constraints, even with an uncountable state space; (ii) By introducing regularization in the form of relative entropy, the solution can be interpreted as replacing a maximum with a softmax in the dual; (iii) In applications such as control for which it is not known apriori if the moment class is nonempty, a relaxation is proposed whose solution admits a similar characterization; (iv) The gradient of the dual function can be expressed in terms of the expectation of the features under a tilted probability measure, which motivates MonteCarlo techniques for computation.
6. Can locational disparity of prosumer energy optimization due to inverter rules be limited? 15 To mitigate issues related to the growth of variable smart loads and distributed generation, distribution system operators (DSO) now make it binding for prosumers with inverters to operate under preset rules. In particular, the maximum active and reactive power set points for prosumers are based on local voltage measurements to ensure that inverter output does not cause voltage violations. However, such actions, as observed in this work, restrict the range available for local energy management, with more adverse losses on arbitrage profits for prosumers located farther away from the substation. The goal of the paper is threefold: (a) to develop an optimal local energy optimization algorithm for activation of load flexibility and inverterinterfaced solar PV and energy storage under timevarying electricity prices; (b) to quantify the locational impact on prosumer arbitrage gains due to inverter injection rules prevalent in different energy markets; (c) to propose a computationally efficient hybrid inverter control policy which provides voltage regulation while substantially reducing locational disparity. Using numerical simulations on three identical prosumers located at different parts of a radial feeder, we show that our control policy is able to minimize locational disparity in arbitrage gains between customers at the beginning and end of the feeder to 1.4%, while PV curtailment is reduced by 91.7% compared to the base case with restrictive voltVar and voltwatt policy.
7. A Generalized Nash Equilibrium analysis of the interaction between a peertopeer financial market and the distribution grid 32 A Generalized Nash Equilibrium analysis of the interaction between a peertopeer financial market and the distribution grid.
8.3 Stochastic matching and queueing systems
8. Analysis of an optimal policy in dynamic bipartite matching models 9 A dynamic bipartite matching model is given by a bipartite matching graph which determines the possible matchings between the various types of supply and demand items. Both supply and demand items arrive to the system according to a stochastic process. Matched pairs leave the system and the others wait in the queues, which induces a holding cost. We model this problem as a Markov Decision Process and study the discounted cost and the average cost problem. We assume that the cost function is linear on the queue sizes. We show that for the shaped matching graph, an optimal matching control prioritizes the matchings in the pendant edges and is of threshold type for the diagonal edge. In addition, for the average cost problem, we compute the optimal threshold value. We then show how the obtained results can be used to characterize the structure of an optimal matching control for a quasicomplete graph with an arbitrary number of nodes. For arbitrary bipartite graphs, we show that, when the cost of the pendant edges is larger than in the neighbors, an optimal matching policy prioritizes the items in the pendant edges. We also study the shaped matching graph and, when the cost of the pendant edges is larger than the cost of the middle edge, we conjecture that an optimal matching policy is also of threshold type with priority to the pendant edges; however, when the cost of the middle edge is larger, we present simulations that show that it is not optimal to prioritize items in the pendant edges.
9. Flexibility can hurt dynamic matching system performance 10 We study the performance of general dynamic matching models. This model is defined by a connected graph, where nodes represent the class of items and the edges the compatibilities between items. Items of different classes arrive one by one to the system according to a given probability distribution. Upon arrival, an item is matched with a compatible item according to the First Come First Served discipline and leave the system immediately, whereas it is enqueued with other items of the same class, if any. We show that such a model may exhibit a non intuitive behavior: increasing the services ability by adding new edges in the matching graph may lead to a larger average population. This is similar to a Braess paradox. We first consider a quasicomplete graph with four nodes and we provide values of the probability distribution of the arrivals such that when we add an edge the mean number of items is larger. Then, we consider an arbitrary matching graph and we show sufficient conditions for the existence or nonexistence of this paradox. We conclude that the analog to the Braess paradox in matching models is given when specific independent sets are in saturation, i.e., the system is close to the stability condition.
10. Product form solution for the steadystate distribution of a Markov chain associated with a general matching model with selfloops 20 We extend the general matching graph model to deal with matching graph where every node has a self loop. Thus the states on the Markov chain are associated with the independent sets of the matching graph. We prove that under i.i.d. arrivals assumptions the steadystate distribution of the Markov chain has a product form solution.
11. Appariements, polytopes et dons d’organes 21 Les programmes de dons croisés de reins donnent lieu à un problème d'appariement complexe pour lequel une solution optimale est inconnue. Ce problème peut être décrit par un modèle d'appariement dynamique et aléatoire où des éléments de différentes classes arrivent selon des processus de Poisson indépendants. Chaque élément représente un couple donneurreceveur où le donneur, prêt à donner un rein pour le receveur, est incompatible avec lui ; et la classe encode des propriétés du couple, telles que les groupes sanguins du donneur et du receveur, qui influencent sa compatibilité. Les contraintes d'appariement sont décrites par un graphe nonorienté dont les sommets sont les classes : deux couples ont des classes voisines dans le graphe s'ils peuvent être appariés dans le sens où chaque donneur peut donner son rein au receveur de l'autre couple. Les éléments en attente sont stockés dans une file d'attente. Nous étudions l'efficacité des politiques d'appariement, non seulement en termes de stabilité du système, mais aussi de contrôle des taux d'appariement entre différentes classes. Ce dernier point est important en vue de privilégier les appariements maximisant d'autres mesures de performance, comme la qualité de vie après la transplantation ou le taux de survie du greffon. Cet article présente les principales contributions sur les appariements dynamiques prépubliés et implantés par les auteurs [CMB22, Mat]. Par souci de concision, les preuves sont omises dans la présente version. La partie 1 formalise le problème et introduit l'équation de conservation qui caractérise les taux d'appariement de toute politique stable. La partie 2 présente une relation simple reliant la structure du graphe, l'équation de conservation, et la stabilité du système. Cette relation nous permet de caractériser dans la partie 3 l'ensemble des solutions de l'équation de conservation comme un polytope convexe. Finalement, la partie 4 montre que tout point du polytope est atteignable ou approchable par des politiques simples qui éliminent certaines arêtes.
12. Stochastic dynamic matching: A mixed graphtheory and linearalgebra approach 40 The stochastic dynamic matching problem has recently drawn attention in the stochasticmodeling community due to its numerous applications, ranging from supplychain management to kidney exchange programs. In this paper, we consider a matching problem in which items of different classes arrive according to independent Poisson processes. Unmatched items are stored in a queue, and compatibility constraints are described by a simple graph on the classes, so that two items can be matched if their classes are neighbors in the graph. We analyze the efficiency of matching policies, not only in terms of system stability, but also in terms of matching rates between different classes. Our results rely on the observation that, under any stable policy, the matching rates satisfy a conservation equation that equates the arrival and departure rates of each item class. Our main contributions are threefold. We first introduce a mapping between the dimension of the solution set of this conservation equation, the structure of the compatibility graph, and the existence of a stable policy. In particular, this allows us to derive a necessary and sufficient stability condition that is verifiable in polynomial time. Secondly, we describe the convex polytope of nonnegative solutions of the conservation equation. When this polytope is reduced to a single point, we give a closedform expression of the solution; in general, we characterize the vertices of this polytope using again the graph structure. Lastly, we show that greedy policies cannot, in general, achieve every point in the polytope. In contrast, nongreedy policies can reach any point of the interior of this polytope, and we give a condition for these policies to also reach the boundary of the polytope.
13. Stationary queues in space and time 5
8.4 Reinforcement learning
14. DelayAware Decentralized Qlearning for Wind Farm Control 27 Wind farms are subject to the socalled "wake effect", where upstream turbines facing the wind create suboptimal wind conditions for turbines located downstream. One strategy to address this issue is to use yaw actuators to misalign the wind turbines with regard to the incoming wind direction, thus deflecting wakes away from downstream turbines. Tractable models for yaw optimization are however subject to inaccuracies, ignore wake dynamics and lack adaptability. This incentivizes the use of modelfree methods. In this paper, we propose a delayaware decentralized Qlearning algorithm for yaw control on wind farms. We introduce a strategy to handle delayed cost collection, and show that our method significantly increases power production in simulations with realistic wake dynamics. We validate our results for two farm layouts on midfidelity wind farm simulator FAST.Farm.
15. Particle gradient descent model for point process generation 8 This paper presents a statistical model for stationary ergodic point processes, estimated from a single realization observed in a square window. With existing approaches in stochastic geometry, it is very difficult to model processes with complex geometries formed by a large number of particles. Inspired by recent works on gradient descent algorithms for sampling maximumentropy models, we describe a model that allows for fast sampling of new configurations reproducing the statistics of the given observation. Starting from an initial random configuration, its particles are moved according to the gradient of an energy, in order to match a set of prescribed moments (functionals). Our moments are defined via a phase harmonic operator on the wavelet transform of point patterns. They allow one to capture multiscale interactions between the particles, while controlling explicitly the number of moments by the scales of the structures to model. We present numerical experiments on point processes with various geometric structures, and assess the quality of the model by spectral and topological data analysis.
8.5 Mathematics of wireless cellular networks
16. Energy and Delay TradeOffs of EndtoEnd Vehicular Communications using a Hyperfractal Urban Modelling 16, 37 We characterize tradeoffs between the endtoend communication delay and the energy in urban vehicular communications with infrastructure assistance. Our study exploits the selfsimilarity of the location of communication entities in cities by modeling them with an innovative model called "hyperfractal". We show that the hyperfractal model can be extended to incorporate roadside infrastructure and provide stochastic geometry tools to allow a rigorous analysis. We compute theoretical bounds for the endtoend communication hop count considering two different energyminimizing goals: either total accumulated energy or maximum energy per node. We prove that the hop count for an endtoend transmission is bounded by $O\left({n}^{1\alpha /({d}_{F}1)}\right)$ where $\alpha <1$ and ${d}_{F}>2$ is the fractal dimension of the mobile nodes process. This proves that for both constraints the energy decreases as we allow choosing routing paths of higher length. The asymptotic limit of the energy becomes significantly small when the number of nodes becomes asymptotically large. A lower bound on the network throughput capacity with constraints on path energy is also given. We show that our model fits real deployments where open data sets are available. The results are confirmed through simulations using different fractal dimensions in a Matlab simulator.
17. Modeling and Analysis of 2Tier Heterogeneous Vehicular Networks Leveraging Roadside Units and Vehicle Relays 39 While roadside units (RSUs) play an essential role in vehicletoeverything (V2X) by communicating with users, some users in congestion areas may not be wellserved due to data traffic, signal attenuation, and interference. In these cases, vehicle relays can be employed to enhance the network topology to better serve those users. This paper leverages stochastic geometry to propose a novel framework for the performance analysis of heterogeneous vehicular networks with RSUs, vehicle relays, and vehicle users. We present a twodimensional analytical model where the spatial dependence between RSUs, vehicle relays, vehicle users, and roads is accurately taken into account through a Cox point process structure. Assuming relays are backhauled to RSUs over a reserved wireless resource and users are associated with the closest RSU or relay, we derive the probability that the typical user is associated with either an RSU or a relay. Then, we derive the signaltointerference ratio (SIR) coverage probability of the typical user. Finally, using the derived formulas, we evaluate the average effective rate of the typical user in the network. This allows us to determine the gain of the average effective rate of users that results from the deployment of relays in the network.
18. A Stochastic Geometry Model for Spatially Correlated Blockage in Vehicular Networks 11 This paper presents a novel spatially consistent approach for modeling lineofsight (LOS) paths in vehicular networks. We use stochastic geometry to model transmitters, obstacles, and receivers located in three parallel lines, respec tively. Their geometric interactions are leveraged to characterize the existence of LOS paths. Specifically, the proposed approach focuses on the role of obstacles in blocking one or more LOS paths, which has been overlooked in most statistical models for blockage. Under the proposed framework, we derive the probability that a typical vehicle is in LOS with respect to transmitters with received signaltonoise ratios greater than a threshold. The proposed framework and LOS coverage analysis are instrumental to the analysis of LOScritical applications such as positioning or mmWave communications in vehicular networks.
19. A User Centric Blockage Model for Wireless Networks 3 This paper proposes a cascade blockage model for analyzing the vision that a user has of a wireless network. This model, inspired by the classical multiplicative cascade models, has a radial structure meant to analyze blockages seen by the receiver at the origin in different angular sectors. The main novelty is that it is based on the geometry of obstacles and takes the joint blockage phenomenon into account. We show on a couple of simple instances that the Laplace transforms of total interference satisfies a functional equation that can be solved efficiently by an iterative scheme. This is used to analyze the coverage probability of the receiver and the effect of blockage correlation and penetration loss in both dense and sparse blockage environments. Furthermore, this model is used to investigate the effect of blockage correlation on user beamforming techniques. Another functional equation and its associated iterative algorithm are proposed to derive the coverage performance of the best beam selection in this context. In addition, the conditional coverage probability is also derived to evaluate the effect of beam switching. The results not only show that beam selection is quite efficient for multibeam terminals, but also show how the correlation brought by blockages can be leveraged to accelerate beam sweeping and pairing.
8.6 Highdimensional statistical inference and distributed learning
20. Partial Recovery in the Graph Alignment Problem 14 In this paper, we consider the graph alignment problem, which is the problem of recovering, given two graphs, a onetoone mapping between nodes that maximizes edge overlap. This problem can be viewed as a noisy version of the wellknown graph isomorphism problem and appears in many applications, including social network deanonymization and cellular biology. Our focus here is on partial recovery, i.e., we look for a onetoone mapping which is correct on a fraction of the nodes of the graph rather than on all of them, and we assume that the two input graphs to the problem are correlated ErdősRényi graphs of parameters $(n,q,s)$. Our main contribution is then to give necessary and sufficient conditions on $(n,q,s)$ under which partial recovery is possible with high probability as the number of nodes $n$ goes to infinity. In particular, we show that it is possible to achieve partial recovery in the $nqs=\Theta \left(1\right)$ regime under certain additional assumptions.
21. Online greedy identification of linear dynamical systems 18 This work addresses the problem of exploration in an unknown environment. For multiinput multioutput, linear timeinvariant dynamical systems, we use an experimental design framework and introduce an online greedy policy where the control maximizes the information of the next step. We evaluate our approach experimentally and compare it with more elaborate gradientbased methods. In a setting with a limited number of observations, our algorithm has low complexity and shows competitive performances.
22. Spectral alignment of correlated Gaussian random matrices 13 In this paper we analyze a simple spectral method (EIG1) for the problem of matrix alignment, consisting in aligning their leading eigenvectors: given two matrices $A$ and $B$, we compute ${v}_{1}$ and ${v}_{1}^{\text{'}}$ two corresponding leading eigenvectors. The algorithm returns the permutation $\widehat{\pi}$ such that the rank of coordinate $\widehat{\pi}\left(i\right)$ in ${v}_{1}$ and that of coordinate $i$ in ${v}_{1}^{\text{'}}$ (up to the sign of ${v}_{1}^{\text{'}}$) are the same. We consider a model of weighted graphs where the adjacency matrix $A$ belongs to the Gaussian Orthogonal Ensemble (GOE) of size $N\times N$, and $B$ is a noisy version of $A$ where all nodes have been relabeled according to some planted permutation $\pi $, namely $B={\Pi}^{T}(A+\sigma H)\Pi $, where $\Pi $ is the permutation matrix associated with $\pi $ and $H$ is an independent copy of $A$. We show the following zeroone law: with high probability, under the condition $\sigma {N}^{7/6+\u03f5}\to 0$ for some $\u03f5>0$, EIG1 recovers all but a vanishing part of the underlying permutation $\pi $, whereas if $\sigma {N}^{7/6\u03f5}\to \infty $, this method cannot recover more than $o\left(N\right)$ correct matches. This result gives an understanding of the simplest and fastest spectral method for matrix alignment (or complete weighted graph alignment), and involves proof methods and techniques which could be of independent interest.
23. Accelerating abelian random walks with hyperbolic dynamics 12 Given integers $d\ge 2,n\ge 1$, we consider affine random walks on torii ${(\mathbb{Z}/n\mathbb{Z})}^{d}$ defined as ${X}_{t+1}=A{X}_{t}+{B}_{t}$ mod $n$, where $A\in {\mathrm{GL}}_{d}\left(\mathbb{Z}\right)$ is an invertible matrix with integer entries and ${\left({B}_{t}\right)}_{t\ge 0}$ is a sequence of iid random increments on ${\mathbb{Z}}^{d}$. We show that when $A$ has no eigenvalues of modulus 1, this random walk mixes in $O(lognloglogn)$ steps as $n\to \infty $, and mixes actually in $O(logn)$ steps only for almost all $n$. These results generalize those on the socalled ChungDiaconisGraham process, which corresponds to the case $d=1$. Our proof is based on the initial arguments of Chung, Diaconis and Graham, and relies extensively on the properties of the dynamical system $x\mapsto {A}^{\top}x$ on the continuous torus ${\mathbb{R}}^{d}/{\mathbb{Z}}^{d}$. Having no eigenvalue of modulus one makes this dynamical system a hyperbolic toral automorphism, a typical example of a chaotic system known to have a rich behaviour. As such our proof sheds new light on the speedup gained by applying a deterministic map to a Markov chain.
26. The graph alignment problem: fundamental limits and efficient algorithms 35 This thesis focuses on statistical inference in graphs (or matrices) in high dimension and studies the graph alignment problem which aims to recover a hidden underlying matching between the nodes of two correlated random graphs. Similarly to many other inference problems in planted models, we are interested in understanding the fundamental informationtheoretical limits as well as the computational hardness of graph alignment. First, we study the Gaussian setting, when the graphs are complete and the signal lies on correlated Gaussian edges weights. We prove that the exact recovery task exhibits a sharp informationtheoretic threshold, characterize it, and study a simple and natural spectral method for recovery, EIG1, which consists in aligning the leading eigenvectors of the adjacency matrices of the two graphs. While most of the recent work on the subject was dedicated to recovering the hidden signal in dense graphs, we next explore graph alignment in the sparse regime, where the mean degrees are constant, not scaling with the graph size. In this particularly challenging setting, for sparse ErdosRényi graphs, only a fraction of the nodes can be correctly matched by any algorithm. Our second contribution is an informationtheoretical result which characterizes a regime where even this partial alignment is impossible, and gives upper bounds on the reachable overlap between any estimator and the true planted matching. We next propose an algorithm that performs partial alignment, NTMA, which is based on a measure of similarity – called the tree matching weight – between treelike neighborhoods of the nodes in the graphs. Under this local approach in the sparse regime, we are brought to study a related problem: correlation detection in random unlabeled trees. This hypothesis testing problem consists in testing whether two trees are correlated or independent. The tree matching weight yields a first method for this question as well; another contribution is to study an optimal test based on the likelihood ratio. In a correlated GaltonWatson model, which is wellknown to be the local approximation of sparse ErdosRényi graphs, we characterize the regimes of performance of this test. Finally, we come back to graph alignment and propose a messagepassing algorithm, MPAlign, naturally inspired by the study of the related problem on trees. This messagepassing algorithm is analyzed and provably recovers a fraction of the planted signal in some regimes of parameters.
27. Statistical limits of correlation detection in trees 46 In this paper we address the problem of testing whether two observed trees $(t,{t}^{\text{'}})$ are sampled either independently or from a joint distribution under which they are correlated. This problem, which we refer to as correlation detection in trees, plays a key role in the study of graph alignment for two correlated random graphs. Motivated by graph alignment, we investigate the conditions of existence of onesided tests, i.e. tests which have vanishing type I error and nonvanishing power in the limit of large tree depth. For the correlated GaltonWatson model with Poisson offspring of mean $\lambda >0$ and correlation parameter $s\in (0,1)$, we identify a phase transition in the limit of large degrees at $s=\sqrt{\alpha}$, where $\alpha \sim 0.3383$ is Otter's constant. Namely, we prove that no such test exists for $s\le \sqrt{\alpha}$, and that such a test exists whenever $s>\sqrt{\alpha}$, for $\lambda $ large enough. This result sheds new light on the graph alignment problem in the sparse regime (with $O\left(1\right)$ average node degrees) and on the performance of the MPAlign method studied in Ganassali et al. (2021), Piccioli et al. (2021), proving in particular the conjecture of Piccioli et al. (2021) that MPAlign succeeds in the partial recovery task for correlation parameter $s>\sqrt{\alpha}$ provided the average node degree $\lambda $ is large enough.
8.7 Distributed optimization for machine learning
28. Muffliato: PeertoPeer Privacy Amplification for Decentralized Optimization and Averaging 22 Decentralized optimization is increasingly popular in machine learning for its scalability and efficiency. Intuitively, it should also provide better privacy guarantees, as nodes only observe the messages sent by their neighbors in the network graph. But formalizing and quantifying this gain is challenging: existing results are typically limited to Local Differential Privacy (LDP) guarantees that overlook the advantages of decentralization. In this work, we introduce pairwise network differential privacy, a relaxation of LDP that captures the fact that the privacy leakage from a node $u$ to a node $v$ may depend on their relative position in the graph. We then analyze the combination of local noise injection with (simple or randomized) gossip averaging protocols on fixed and random communication graphs. We also derive a differentially private decentralized optimization algorithm that alternates between local gradient descent steps and gossip averaging. Our results show that our algorithms amplify privacy guarantees as a function of the distance between nodes in the graph, matching the privacyutility tradeoff of the trusted curator, up to factors that explicitly depend on the graph topology. Finally, we illustrate our privacy gains with experiments on synthetic and realworld datasets.
29. On Sample Optimality in Personalized Collaborative and Federated Learning 23 In personalized federated learning, each member of a potentially large set of agents aims to train a model minimizing its loss function averaged over its local data distribution. We study this problem under the lens of stochastic optimization, focusing on a scenario with a large number of agents, that each possess very few data samples from their local data distribution. Specifically, we prove novel matching lower and upper bounds on the number of samples required from all agents to approximately minimize the generalization error of a fixed agent. We provide strategies matching these lower bounds, based on a gradient filtering approach: given prior knowledge on some notion of distance between local data distributions, agents filter and aggregate stochastic gradients received from other agents, in order to achieve an optimal biasvariance tradeoff. Finally, we quantify the impact of using rough estimations of the distances between local distributions of agents, based on a very small number of local samples.
30. Correlation detection in trees for planted graph alignment 24 Motivated by alignment of correlated sparse random graphs, we introduce a hypothesis testing problem of deciding whether or not two random trees are correlated. We obtain sufficient conditions under which this testing is impossible or feasible. We propose MPAlign, a messagepassing algorithm for graph alignment inspired by the tree correlation detection problem. We prove MPAlign to succeed in polynomial time at partial alignment whenever tree detection is feasible. As a result our analysis of tree detection reveals new ranges of parameters for which partial alignment of sparse random graphs is feasible in polynomial time. We then conjecture that graph alignment is not feasible in polynomial time when the associated tree detection problem is impossible. If true, this conjecture together with our sufficient conditions on tree detection impossibility would imply the existence of a hard phase for graph alignment, i.e. a parameter range where alignment cannot be done in polynomial time even though it is known to be feasible in nonpolynomial time.
31. Robustness in MultiObjective Submodular Optimization: a Quantile Approach 25 The optimization of multiobjective submodular systems appears in a wide variety of applications. However, there are currently very few techniques which are able to provide a robust allocation to such systems. In this work, we propose to design and analyse novel algorithms for the robust allocation of submodular systems through lens of quantile maximization. We start by observing that identifying an exact solution for this problem is computationally intractable. To tackle this issue, we propose a proxy for the quantile function using a softmax formulation, and show that this proxy is well suited to submodular optimization. Based on this relaxation, we propose a novel and simple algorithm called SOFTSAT. Theoretical properties are provided for this algorithm as well as novel approximation guarantees. Finally, we provide numerical experiments showing the efficiency of our algorithm with regards to stateoftheart methods in a test bed of realworld applications, and show that SOFTSAT is particularly robust and wellsuited to online scenarios.
32. Asynchronous SGD Beats Minibatch SGD Under Arbitrary Delays 26 The existing analysis of asynchronous stochastic gradient descent (SGD) degrades dramatically when any delay is large, giving the impression that performance depends primarily on the delay. On the contrary, we prove much better guarantees for the same asynchronous SGD algorithm regardless of the delays in the gradients, depending instead just on the number of parallel devices used to implement the algorithm. Our guarantees are strictly better than the existing analyses, and we also argue that asynchronous SGD outperforms synchronous minibatch SGD in the settings we consider. For our analysis, we introduce a novel recursion based on "virtual iterates" and delayadaptive stepsizes, which allow us to derive stateoftheart guarantees for both convex and nonconvex objectives.
33. Convergence beyond the overparameterized regime using Rayleigh quotients 29 In this paper, we present a new strategy to prove the convergence of deep learning architectures to a zero training (or even testing) loss by gradient flow. Our analysis is centered on the notion of Rayleigh quotients in order to prove KurdykaŁojasiewicz inequalities for a broader set of neural network architectures and loss functions. We show that Rayleigh quotients provide a unified view for several convergence analysis techniques in the literature. Our strategy produces a proof of convergence for various examples of parametric learning. In particular, our analysis does not require the number of parameters to tend to infinity, nor the number of samples to be finite, thus extending to test loss minimization and beyond the overparameterized regime.
34. Periodic Signal Recovery with Regularized Sine Neural Networks 30 We consider the problem of learning a periodic onedimensional signal with neural networks, and designing models that are able to extrapolate the signal well beyond the training window. First, we show that multilayer perceptrons with ReLU activations are provably unable to perform this extrapolation task, and lead to poor performance in practice even close to the training window. Then, we propose a modified training procedure for twolayer architectures with sine activations with a more diverse feature initialization and wellchosen nonconvex regularization, that is able to extrapolate the signal with low error well beyond the training window. This procedure yields results several orders of magnitude better than its competitors for distant extrapolation (beyond 100 periods of the signal), while being able to accurately recover the frequency spectrum of the signal in a multitone setting.
35. Convergence Rates of NonConvex Stochastic Gradient Descent Under a Generic Łojasiewicz Condition and Local Smoothness 31 Training overparameterized neural networks involves the empirical minimization of highly nonconvex objective functions. Recently, a large body of works provided theoretical evidence that, despite this nonconvexity, properly initialized overparameterized networks can converge to a zero training loss through the introduction of the PolyakŁojasiewicz condition. However, these analyses are restricted to quadratic losses such as mean square error, and tend to indicate fast exponential convergence rates that are seldom observed in practice. In this work, we propose to extend these results by analyzing stochastic gradient descent under more generic Łojasiewicz conditions that are applicable to any convex loss function, thus extending the current theory to a larger panel of losses commonly used in practice such as crossentropy. Moreover, our analysis provides highprobability bounds on the approximation error under subGaussian gradient noise and only requires the local smoothness of the objective function, thus making it applicable to deep neural networks in realistic settings.
36. Sample Optimality and Allforall Strategies in Personalized Federated and Collaborative Learning 44 In personalized Federated Learning, each member of a potentially large set of agents aims to train a model minimizing its loss function averaged over its local data distribution. We study this problem under the lens of stochastic optimization. Specifically, we introduce information theoretic lower bounds on the number of samples required from all agents to approximately minimize the generalization error of a fixed agent. We then provide strategies matching these lower bounds, in the allforone and allforall settings where respectively one or all agents desire to minimize their own local function. Our strategies are based on a gradient filtering approach: provided prior knowledge on some notions of distances or discrepancies between local data distributions or functions, a given agent filters and aggregates stochastic gradients received from other agents, in order to achieve an optimal biasvariance tradeoff.
37. On NonLinear operators for Geometric Deep Learning 49
8.8 Random Geometry
38. A Unified Framework for Generalizing the GromovHausdorff Metric 47 In this paper, a general approach is presented for generalizing the GromovHausdorff metric to consider metric spaces equipped with some additional structure. A special case is the GromovHausdorffProkhorov metric which considers measured metric spaces. This abstract framework also unifies several existing generalizations which consider metric spaces equipped with a measure, a point, a closed subset, a curve or a tuple of such structures. It can also be useful for studying new examples of additional structures. The framework is provided both for compact metric spaces and for boundedlycompact pointed metric spaces. In addition, completeness and separability of the metric is proved under some conditions. This enables one to study random metric spaces equipped with additional structures, which is the main motivation of this work.
39. An Improved Lower Bound on the Largest Common Subtree of Random LeafLabeled Binary Trees 48 It is known that the size of the largest common subtree (i.e., the maximum agreement subtree) of two independent random binary trees with $n$ given labeled leaves is of order between ${n}^{0.366}$ and ${n}^{1/2}$. We improve the lower bound to order ${n}^{0.4464}$ by constructing a common subtree recursively and by proving a lower bound for its asymptotic growth. The construction is a modification of an algorithm proposed by D. Aldous by splitting the tree at the centroid and by proceeding recursively.
40. Waveletbased representations of point processes for modelling and statistical learning 33 This dissertation presents a class of representations of spatial point processes. Inspired from the success of wavelet methods in signal processing, these descriptors rely on the convolution of a point process with a family of wavelet filters. From these convolutions are built sets of statistical descriptors of stationary point process, by applying nonlinear operators, followed by a spatial averaging. Much like classical summary characteristics for point processes, these statistics are designed to extract information about the process with a relatively small number of numerical values, by describing its geometry. Their goal is to describe whether the atoms of the process tend to repel each other, or cluster together, and by doing so, form possibly complex geometric shapes. By construction, these descriptors enjoy several properties that make them suitable for statistical analysis and learning tasks. To illustrate the quality of these representations as statistical descriptors, we study several problems involving statistical analysis of point processes. In a first experiment, we seek to estimate an unknown function that takes as input a point pattern, and returns a marked version of this pattern, where a numerical value is associated to each atom of the pattern. We use a waveletbased representation of point patterns to estimate the relation between their nonmarked and marked version. We then study, in a second experiment, the ability of such representations to model the distribution of a point process, by defining a maximum entropy model defined by a set of waveletbased statistics, computed on a single observation. For these two problems, we observe that our representations lead to better performance than summary statistics commonly used in the literature on point processes. Finally, to study to what extent such representations can capture geometric structures of texture images, we define a maximum entropy model relying on similar wavelet statistics, yielding syntheses of similar visual quality to stateoftheart models based on deep convolutional neural networks representations.
41. Counter examples to invariant circle packing 17 In this work, a unimodular random planar triangulation is constructed that has no invariant circle packing. This disputes a problem asked in [9]. A natural weaker problem is the existence of pointstationary circle packings for a graph, which are circle packings that satisfy a certain mass transport principle. It is shown that the answer to this weaker problem is also false. Two examples are provided with two diﬀerent approaches: Using indistinguishability and ﬁnite approximations.
42. Generalized Rectifier Wavelet Covariance Models For Texture Synthesis 19 Stateoftheart maximum entropy models for texture synthesis are built from statistics relying on image representations defined by convolutional neural networks (CNN). Such representations capture rich structures in texture images, outperforming waveletbased representations in this regard. However, conversely to neural networks, wavelets offer meaningful representations, as they are known to detect structures at multiple scales (e.g. edges) in images. In this work, we propose a family of statistics built upon nonlinear wavelet based representations, that can be viewed as a particular instance of a onelayer CNN, using a generalized rectifier nonlinearity. These statistics significantly improve the visual quality of previous classical waveletbased models, and allow one to produce syntheses of similar quality to stateoftheart models, on both grayscale and color textures.
43. Coupling from the Past for the Null Recurrent Markov Chain 36 The Doeblin Graph of a countable state space Markov chain describes the joint pathwise evolutions of the Markov dynamics starting from all possible initial conditions, with two paths coalescing when they reach the same point of the state space at the same time. Its Bridge Doeblin subgraph only contains the paths starting from a tagged point of the state space at all possible times. In the irreducible, aperiodic, and positive recurrent case, the following results are known: the Bridge Doeblin Graph is an infinite tree that is unimodularizable. Moreover, it contains a single biinfinite path, which allows one to build a perfect sample of the stationary state of the Markov chain. The present paper is focused on the null recurrent case. It is shown that when assuming irreducibility and aperiodicity again, the Bridge Doeblin Graph is either an infinite tree or a forest made of a countable collection of infinite trees. In the first case, the infinite tree in question has a single end, is not unimodularizable in general, but is always locally unimodular. These key properties are used to study the stationary regime of several measurevalued random dynamics on this Bridge Doeblin Tree. The most important ones are the taboo random dynamics, which admits as steady state a random measure with mean measure equal to the invariant measure of the Markov chain, and the potential random dynamics which is a random extension of the classical potential measure, with a mean measure equal to infinity at every point of the state space.
44. Nash equilibrium structure of Cox process Hotelling games 1 We study an Nplayer game where a pure action of each player is to select a nonnegative function on a Polish space supporting a finite diffuse measure, subject to a finite constraint on the integral of the function. This function is used to define the intensity of a Poisson point process on the Polish space. The processes are independent over the players, and the value to a player is the measure of the union of its open Voronoi cells in the superposition point process. Under randomized strategies, the process of points of a player is thus a Cox process, and the nature of competition between the players is akin to that in Hotelling competition games. We characterize when such a game admits Nash equilibria and prove that when a Nash equilibrium exists, it is unique and comprised of pure strategies that are proportional in the same proportions as the total intensities. We give examples of such games where Nash equilibria do not exist. A better understanding of the criterion for the existence of Nash equilibria remains an intriguing open problem.
8.9 Mathematics of stochastic networks via meanfield analysis
45. A new stochastic model for carsharing suited to freefloating 45 Freefloating has an increasing popularity in carsharing but imbalance problem shows the need of stochastic modelling and analysis. In this paper, a new stochastic model suited to freefloating is proposed, taking into account the sharing of public space between private and freefloating cars.The capacity seen by freefloating cars in such a model turns out to be random. We show that unlike stationbased carsharing systems, it is not bounded. Moreover, a stochastic averaging principle governs the freefloating car behavior. We exhibit a phase transition between a nonsaturated regime where freefloating cars can always be parked and a saturated one where freefloating cars do not find an available parking space with positive probability. This probability is completely determined by the environment which implies that the operator cannot act on the proportion of saturated areas, i.e. without parking parking spaces. We solve the dimensioning problem; the more freefloating cars in the system, the more satisfied users are.
46. ReplicaMeanField Limits of FragmentationInteractionAggregation Processes 2 Network dynamics with pointprocessbased interactions are of paramount modeling interest. Unfortunately , most relevant dynamics involve complex graphs of interactions for which an exact computational treatment is impossible. To circumvent this difficulty, the replicameanfield approach focuses on randomly interacting replicas of the networks of interest. In the limit of an infinite number of replicas , these networks become analytically tractable under the socalled "Poisson Hypothesis". However, in most applications, this hypothesis is only conjectured. Here, we establish the Poisson Hypothesis for a general class of discretetime, pointprocessbased dynamics, that we propose to call fragmentationinteractionaggregation processes, and which are introduced in the present paper. These processes feature a network of nodes, each endowed with a state governing their random activation. Each activation triggers the fragmentation of the activated node state and the transmission of interaction signals to downstream nodes. In turn, the signals received by nodes are aggregated to their state. Our main contribution is a proof of the Poisson Hypothesis for the replicameanfield version of any network in this class. The proof is obtained by establishing the propagation of asymptotic independence for state variables in the limit of an infinite number of replicas. Discrete time GalvesLöcherbach neural networks are used as a basic instance and illustration of our analysis.
47. Propagation of chaos and Poisson Hypothesis for replica meanfield models of intensitybased neural networks 43 Neural computations arising from myriads of interactions between spiking neurons can be modeled as network dynamics with punctuate interactions. However, most relevant dynamics do not allow for computational tractability. To circumvent this difficulty, the Poisson Hypothesis regime replaces interaction times between neurons by Poisson processes. We prove that the Poisson Hypothesis holds at the limit of an infinite number of replicas in the replicameanfield model, which consists of randomly interacting copies of the network of interest. The proof is obtained through a novel application of the ChenStein method to the case of a random sum of Bernoulli random variables and a fixed point approach to prove a law of large numbers for exchangeable random variables.
48. Mean field analysis of stochastic networks with reservation 38 The problem of reservation in a large distributed system is analyzed via a new mathematical model. A typical application is a stationbased carsharing system which can be described as a closed stochastic network where the nodes are the stations and the customers are the cars. The user can reserve the car and the parking space. In the paper, we study the evolution of the system when the reservation of parking spaces and cars is effective for all users. The asymptotic behavior of the underlying stochastic network is given when the number N of stations and the fleet size increase at the same rate. The analysis involves a Markov process on a state space with dimension of order N 2. It is quite remarkable that the state process describing the evolution of the stations, whose dimension is of order N , converges in distribution, although not Markov, to an nonhomogeneous Markov process. We prove this meanfield convergence. We also prove, using combinatorial arguments, that the meanfield limit has a unique equilibrium measure when the time between reserving and picking up the car is sufficiently small. This result extends the case where only the parking space can be reserved.
49. Mean field analysis of an incentive algorithm for a closed stochastic network 28 The paper deals with a loadbalancing algorithm for a closed stochastic network with two zones with different demands. The algorithm is motivated by an incentive algorithm for redistribution of cars in a largescale carsharing system. The service area is divided into two zones. When cars stay too much long in the lowdemand zone, users are encouraged to pick up them and return them in the highdemand zone. The zones are divided in cells called stations. The cars are the network customers. The meanfield limit solution of an ODE gives the large scale distribution of the station state in both clusters for this incentive policy in a discrete Markovian framework. An equilibrium point of this ODE is characterized via the invariant measure of a random walk in the quarterplane. The proportion of empty and saturated stations measures how the system is balanced. Numerical experiments illustrate the impact of the incentive policy. Our study shows that the incentive policy helps when the highdemand zone observes a lack of cars but a saturation must be prevented especially when the highdemand zone is small.
9 Bilateral contracts and grants with industry
9.1 Bilateral contracts with industry
Participants: François Baccelli, Bartłomiej Błaszczyszyn, Ana Bušić, Guodong Sun.
9.1.1 Contract with EDF
Collaborative research in the area of demand dispatch of flexible loads. PI: A. Bušić.
9.1.2 CRE with Orange Labs
Two year contract titled Guaranteed throughput and millimeter waves in the dimensioning of 5G cellular networks between Inria and Orange Labs started 2022 and attached to the joint InriaOrange lab "IOLab". It is a part of a longterm collaboration between TREC/DYOGENE, represented by B. Błaszczyszyn and Orange Labs, represmnted by M. K. Karray on the development of mathematical methods and engineering analytical tools for the dimensioning of wireless cellular networks enabling operators to solve critical technical and economic issues related to the main business related to the permanent evolution of radio technology. They capture the macroscopic relation between antenna deployment, frequency allocation, the volume of traffic carried on the network, and QoS parameters, such as the average and variation in bandwidth available to end users. The math solutions developed in cooperation with Orange Labs are implemented by Orange Labs in their internal toolbox and used for dimensioning studies by Direction des Affaires Réglementaires of Orange, Orange France in relation to ARCEP (Autorité de Régulation des Communications Électroniques, des Postes et de la Distribution de la Presse) and Orange affiliates.
9.1.3 CIFRE with Nokia Bell Labs
Contract with Nokia Bell Labs France started in 2021 for the coadvising by F. Baccelli of a PhD student Guodong Sun.
10 Partnerships and cooperations
10.1 International initiatives
10.1.1 IFCAM poject
Participants: Bartłomiej Błaszczyszyn, Antoine Brochard, Romain Gambelin.
IndoFrench Centre for Applied Mathematics (IFCAM) project Geometric statistics of stationary point processes started in 2018 and extended to March 2023 with coordinators B. Błaszczyszyn (Inria), Yogeshwaran D. (ISI Bangalore). The project is investigating the asymptotic properties of the geometric statistics of point processes, which have recently gained importance in stochastic geometry, in particular for more general point processes, for which theory still lags far behind Poisson point processes. Our main contribution and results are as follows. (1) In the theory of point processes, we introduce the concept of decorrelation to general marked point processes and show that this "asymptotic independence" of the structure — an essential element of limit theory — appears when points and marks have this property separately. (2) This allows us to study limit results for more general models, in particular spatially dependent interacting particle systems in Euclidean space, permitting for dependent sites and initial states, unbounded interaction range and general (and nonnecessarily Markovian) updates. (3) We prove or mention the results of the central limit theorem for several examples of these modes (as sequential adsorption, ballistic deposition, majority dynamics, epidemic models) allowing for nonPoisson locations of particles and their dependent initial states. We shall also investigate some applications to the theory of signal processing. (4) In general, limits of variance or asymptotics of covariance are commonly accepted in central limit theorems. In two complementary projects, we show that volumeorder variance asymptotics are preserved under quasilocal perturbations of a stationary random measure satisfying some mixing condition. We also show possible covariance asymptotics for different test statistics of general stationary random measures. (5) Finally, in a related project, we are also interested in optimization of large particle systems considered in the framework of stationary marked point process to formalize novel optimization concepts: intensityoptimal and locally optimal marking. Our main results provide general conditions when the intensityoptimal solution exists and when both optimization solutions are equivalent. We discuss several applications ranging from theoretical probability to applications in statistical physics, combinatorial optimization and communications.10.1.2 Alliance Communauto Montreal
Participants: Christine Fricker.
Dyogene participates in 4year research project (20232027) cofunded by Communauto, NSERC and MItacs in Montreal.10.2 International research visitors
Scientists speaking at the DYOGENE's seminar in 2022
 Céline Compte (Eindhoven University of Technology),
 Senya Shlosman (Centre de Physique Theorique; Skoltech Center of Advance Studies),
 Nicolas Curien (Université ParisSaclay),
 Samuel Mellick (ENS de Lyon),
 Umberto de Ambroggio (University of Bath),
 Takis Konstantopoulos (University of Liverpool),
 Alexandre Gaudillière (CNRS),
 JeanFrançois Le Gall (Université ParisSaclay),
 David Aldous (UC Berkeley),
 Eva Löcherbach ( Université Paris 1),
 Vsevolod Shneer (HeriotWatt University),
 Sergey Foss (HeriotWatt University),
 Hermann Thorisson (University of Iceland),
 Sergey Foss (HeriotWatt University),
 Pierre Calka (University of Rouen),
 Matteo D'Achille (Université ParisEst  Créteil ValdeMarne),
 Eliza O'Reilly (California Institute of Technology),
 Charles Bordenave (Institut de Mathématiques de Marseille),
 Milad Barzegar (Sharif University of Technology),
 Anita Winter (University of DuisburgEssen),
 Payam Delgosha (University of Illinois),
 Payam Delgosha (University of Illinois),
 Nicolas Gast (Inria).
Inria International Chair
Participants: Ana Busic, Sean Meyn.
Prof Sean Meyn (University of Florida) is holding Inria International Chair from 2019 to 2023 in association with Dyogene and Sierra teams.10.3 European initiatives
10.3.1 ERC NEMO
Participants: François Baccelli, Bartłomiej Błaszczyszyn, Bharath Roy Choudhury, Ke Feng, Sayeh Khaniha, Ali Khezeli, Pierre Popineau.
NEMO NEtwork MOtion, Inria/NEMO is an ERC Advanced Grant (2019 – 2024, PI François Baccelli). It is an interdisciplinary proposal centered on network dynamics. The interdisciplinarity spans from communication engineering to mathematics, with an innovative interplay between the two. NEMO’s aim is to introduce dynamics in stochastic geometry. General mathematical tools combining stochastic geometry, random graph theory, and the theory of dynamical systems will be developed. NEMO will leverage interactions of Inria with Ecole Normale Supérieure on the mathematical side, and with Nokia Bell Labs and Orange on the engineering side.
10.3.2 Collaboration with VITO (Belgium)
Participants: Ana Busic, Ilia Shilov.
Collaboration with VITO Belgium. Coadvising of PhD student I. Shilov. Started: Nov 2019. Topic: “Algorithmic Games and Distributed Learning for PeertoPeer Energy Trading”. PhD scholarship by VITO.10.4 National initiatives
10.4.1 PRAIRIE Institute
Participants: Marc Lelarge, Laurent Massoulié.
The Prairie Institute (PaRis AI Research InstitutE) is one of the four French Institutes for Interdisciplinary Artificial Intelligence Research (3IA), which were created as part of the national French initiative on AI announced by President Emmanuel Macron on May 29, 2018. It brings together five academic partners (CNRS, Inria, Institut Pasteur, PSL University, and University of Paris) as well as 14 industrial partners.10.4.2 GdR GeoSto
Participants: François Baccelli, Bartłomiej Błaszczyszyn, Bharath Roy Choudhury, Michel Davydov, Romain Gambelin, Sayeh Khaniha, Ali Khezeli.
Members of Dyogene participate in Research Group GeoSto (Groupement de recherche, GdR 3477) gdrgeostoch.math.cnrs.fr on Stochastic Geometry currenly led by Hermine Biermé [Université de Tours]. This is a collaboration framework for all French research teams working in the domain of spatial stochastic modeling, both on theory development and in applications.
10.4.3 GdR RO
Participants: Ana Busic, Thomas Le Corre, Lucas Weber.
Members of Dyogene participate in the group Recherche Opérationelle GdRRO working on Stochastic optimization and control, modeling and simulation COSMOSlead by A. Busic and E. Hyon (LIP 6).10.4.4 Défi Inria FedMalin
Participants: Kevin Scaman, Laurent Massoulié.
Members of Dyogene participate in the FedMalin Inria défi on Federated Learning.10.5 Regional initiatives
10.5.1 LINCS
Participants: All Dyogene.
Laboratory of Information, Networking and Communication Sciences Dyogene participates in LINCS, a research centre cofounded by Inria, Institut Mines Télécom, UPMC and AlcatelLucent Bell Labs (currently Nokia Bell Labs) dedicated to research and innovation in the domains of future information and communication networks, systems and services.11 Dissemination
All Dyogene
11.1 Promoting scientific activities
11.1.1 LINCS+NEMO workshop on vehicular networks
ERC Nemo coorganized an international Workshop on vehicular networks, hosted a LINCS in Paris; March 1011. The 14 guest speakers covered a wide range of questions relating to the scientific and industrial challenges of this type of communication network:
 Analysis of the reliability of the various radio links (safety messages, exchanged between vehicles, pedestrians or cyclists). Presentations by C.S. Choi (Hongik Univ.), H. Dhillon (Virginia Tech.), M. Haenggi (ETH), V. Mancuso (IMDEA) and O. Shagdar (VEDECOM).
 Integrating communications and measurements: N. GonzalezPrélcic (NCSU), V. Petrov (Nokia) and H. Wymeersh (Chalmers).
 Distributed systems aspects, including distributed usage of AI at the edge of the network. C. Chiasserini (PoliTO), F. Dressler (TU Berlin) and R. Heath (NCSU).
 Industrial issues: O. Altintas (Toyota) and T.V. Nguyen (Qualcomm).
 The standardization: T. Ernst (YoGoKo).
The Presentations and recorded talks as well as the discussions are available on LINCSNEMO coference website.
11.1.2 Other scientific events
 ERC Nemo organized a Workshop on random measures, unimodularity and Palm Theory at Inria Paris; May 910.
 F. Baccelli participated in the organization of the Inria national day on communications networks on September 13, 2023 in Rungis.
 A. Khezeli organized the weekly DYOGENE seminar. See see Section 10.2 for the scientists speaking at the DYOGENEs seminar in 2022.
 M. Even et M. Blanke orgnaised the weekly ARGO seminar.
 B. RoyChowdhury organized the NEMO reading group
Reviewer
 B. Blaszczyszyn: for Electronic Journal of Statistics, Statistics and Computing, ACM Symposium on Theory of Computing, NWO  The Dutch Research Council.
 K. Scaman: for Journal of Machine Learning Research (2 papers) and Agence Nationale de la Recherche (1 proposal)
11.1.3 Invited talks
 Invited lecture on Particle gradient descent model for point process generation at Journées communes GéoSto MIA, SaintÉtienneduRouvray 2223 sept. 2022, B. Blaszczyszyn
 Invited lecture on Limit theory for asymptotically decorrelated dynamic spatial random models at Recent Trends in Spatial Stochastic Processes, Eindhoven, Netherlands Oct 0307, 2022, B. Blaszczyszyn
 Invited lecture on replica meanfield limits for neural networks at the Mathematics Colloquium, University of Houston, Houston, USA, January 31, 2022, F. Baccelli.
 Keynote lecture on vehicular networks at the COSTINTERACT 22 meeting, Lyon, France, June 1315, 2022, F.Baccelli.
 Invited lecture on contagionmigration processes at "Atelier d'Évaluation de Performance 22", Grenoble, France, July 0405, 2022, F. Baccelli.
 Keynote lecture on wireless stochastic geometry at IEEE Wiopt, Torino, Italy, September 1923, 2022, F. Baccelli.
 Invited lecture on unimodular random graphs at the Workshop "Recent Trends in Spatial Stochastic Processes", Eindhoven, Netherlands Oct 0307, 2022, F. Baccelli.
 Online seminar on stationarity and ergodicity, IPM Geometry and Topology seminars, Feb 2022, A. Khezeli.
 Lecture on unimodular continuum Spaces at the Workshop "Recent Trends in Spatial Stochastic Processes", Eindhoven, Netherlands Oct 0307, 2022, A. Khezeli.
 Lecture entitled "An Improved Lower Bound on the Largest Common Subtree of Random LeafLabeled Binary Trees", IPM Combinatorics, Oct 2022, A. Khezeli.
11.1.4 Scientific expertise
 A. Busic: member of Comission d'Emplois Scientifiques (CES) et Commission de Développement Technologique (CDT)
11.2 Teaching  Supervision  Juries
11.2.1 Teaching
 Licence: B. Błaszczyszyn (Cours) Théorie de l'information et du codage 24 heqTD, L3, ENS Paris. moodle.psl.eu
 Licence: A. Busic (Cours) and Romain Cosson (Cours + TD) Structures et algorithmes aléatoires 60 heqTD, L3, ENS Paris.
 Licence: M. Davydov (Cours) Introduction aux probabilités (L2 Mathématiques) 38h(eqTD) and Bases des mathématiques : semaine de prérentrée (L1 Mathématiques & Informatique) 15h(eqTD), University ParisCité.
 Licence: R. Gambelin, CM de "Probabilités approfondies" en seconde année de la CPES de PSL (filière mathématiques et physique).
 Licence: L. Massoulié (Cours) Social and Communication networks 60 heqTD, L3, l'X.
 Master: B. Błaszczyszyn (Cours) "Random Geometric Models", jointly at M2 “Probabilities and Random Models”, Sorbonne University and M2 “Applied and Theoretical Mathematics”, University ParisDauphinePSL (39heqTD). moodle.psl.eu
 Master: A. Busic (Cours) and L. Weber (TD) Modèles et algorithmes de réseaux 60 heqTD, M1, ENS Paris.
 Master: A. Busic and L. Massoulié (Cours) Fondements de la modélisation des réseaux 37.5 heqTD, M2 MPRI.
 Master: M. Lelarge (Cours) Deep Learning Do it Yourself, X, XHEC. mlelarge.github.io/dataflowrweb/
 Master: L. Massoulié (Cours) Inference in large random graphs, M2 Université d'Orsay.
 Master: K. Scaman with the participation of M. Lelarge (Cours) Deep Learning, 24 heqTD, M1, ENS Paris.
 Master: K. Scaman (Cours) and M. Even (TD) Mathematics of Deep Learning, 30 heqTD, M2 MASH, Université Dauphine.
 Master: K. Scaman (TD) Machine Learning 2, 21 heqTD, M2 XHEC, Institut Polytechnique de Paris.
 Master: R. Gambelin, Chargé de TD "Méthodes quantitatives", L2 Eco & Gestion, l'université Paris Descartes.
11.2.2 Supervision
PhD defended:
 Antoine Brochard, "Waveletbased representations of point processes for modelling and statistical learning", Ecole Normale Supérieure (Paris),March 2022, 33; supervised by B. Błaszczyszyn
 Eric Daoud, "Geographic and sociodemographic disparities in oncology care pathways", Université ParisSaclay, December 2022, 34; supervised by M. Lelarge
 Luca Ganassali, "The graph alignment problem: fundamental limits and efficient algorithms", PSL Research University; Ecole normale supérieure, September 2022, 35; supervised by Laurent Massoulié and Marc Lelarge
PhD in progress:
 Claire Bizon Monroc, since Nov 2021, supervised by A. Bušić
 Matthieu Blanke since Sep 2021, supervised by M. Lelarge
 Romain Cosson, since Sep 2022, supervised by Laurent Massoulié
 Bharath Roy since Sep 2019, supervised by F. Baccelli
 Michel Davydov, since 2020, supervised by F. Baccelli
 Bastien Dubail since Sep 2020, supervised by Laurent Massoulié
 Mathieu Even since Sep 2021, supervised by Laurent Massoulié
 Roman Gambelin since Sep 2020, supervised by B. Blaszczyszyn
 Cedric Gerbelot since Sep 2021, supervised by M. Lelarge
 Sayeh Khaniha, since Feb 2019, supervised by F. Baccelli
 Thomas Le Corre since Nov 2021, supervised by Bušić
 Maxime Leiber since Feb 2020, supervised by Laurent Massoulié
 Jakob Maier, since Oct 2022, supervised by Laurent Massoulié
 Pierre Popineau, since Sep 2019, supervised by F. Baccelli
 David Robin, since Oct 2021, supervised by Laurent Massoulié
 Ilia Shilov since Sep 2019, supervised by A. Bušić
 Guodong Sun, since 2021, supervised by F. Baccelli
 Amaury Triboulin since Sep 2021, supervised by M. Lelarge
 Lucas Weber since Oct 2021, supervised by A. Bušić
11.2.3 Juries
 F. Baccelli: member of the thesis committee of Quentin Gontier (ULB, Bruxelles).
 B. Błaszczyszyn: member of the thesis committee of Iman Hmedoush (Université Paris Saclay) and Safa Ladgham (Université de Paris Cité).
12 Scientific production
12.1 Publications of the year
International journals
 1 articleNash equilibrium structure of Cox process Hotelling games.Advances in Applied Probability542June 2022, 570598
 2 articleReplicaMeanField Limits of FragmentationInteractionAggregation Processes.Journal of Applied Probability591March 2022, 3859
 3 articleA User Centric Blockage Model for Wireless Networks.IEEE Transactions on Wireless Communications21102022, 10 p.
 4 articleA computational framework for evaluating the role of mobility on the propagation of epidemics on point processes.Journal of Mathematical Biology8412January 2022, 4
 5 articleStationary queues in space and time.Queueing Systems10034April 2022, 501503
 8 articleParticle gradient descent model for point process generation.Statistics and Computing323June 2022
 9 articleAnalysis of an optimal policy in dynamic bipartite matching models.Performance Evaluation154April 2022, 102286
 10 articleFlexibility can hurt dynamic matching system performance.ACM SIGMETRICS Performance Evaluation Review493March 2022
 11 articleA Stochastic Geometry Model for Spatially Correlated Blockage in Vehicular Networks.IEEE Internet of Things Journal920October 2022, 1988119889
 12 articleAccelerating abelian random walks with hyperbolic dynamics.Probability Theory and Related Fields18434December 2022, 939968
 13 articleSpectral alignment of correlated Gaussian random matrices.Advances in Applied Probability54Issue 1March 2022, 279  310
 14 articlePartial Recovery in the Graph Alignment Problem.Operations ResearchJanuary 2022
 15 articleCan locational disparity of prosumer energy optimization due to inverter rules be limited?IEEE Transactions on Power Systems2022, 115
 16 articleEnergy and Delay TradeOffs of EndtoEnd Vehicular Communications using a Hyperfractal Urban Modelling.Annals of Telecommunications  annales des télécommunications2023
 17 articleCounter examples to invariant circle packing.Annales de l'Institut Henri Poincaré (B) Probabilités et Statistiques584November 2022
International peerreviewed conferences
 18 inproceedingsOnline greedy identification of linear dynamical systems.2022 IEEE 61st Conference on Decision and Control (CDC)CDC 2022  61st Conference on Decision and Control61Cancun, MexicoIEEE2022, 53635368
 19 inproceedingsGeneralized Rectifier Wavelet Covariance Models For Texture Synthesis.ICLR 2022  10th International Conference on Learning RepresentationsVirtual, FranceApril 2022
 20 inproceedingsProduct form solution for the steadystate distribution of a Markov chain associated with a general matching model with selfloops.EPEW 2022: 18th European Performance Engineering WorkshopSanta Pola, Alicante, Spain2022
 21 inproceedingsAppariements, polytopes et dons d'organes.AlgoTel 2022  24èmes Rencontres Francophones sur les Aspects Algorithmiques des TélécommunicationsSaintRémyLèsChevreuse, FranceMay 2022
 22 inproceedingsMuffliato: PeertoPeer Privacy Amplification for Decentralized Optimization and Averaging.Advances in Neural Information Processing Systems 35 (NeurIPS)New Orleans, United States2022
 23 inproceedingsOn Sample Optimality in Personalized Collaborative and Federated Learning.NeurIPS 2022  36th Conference on Neural Information Processing SystemNew Orleans, United StatesNovember 2022
 24 inproceedingsCorrelation detection in trees for planted graph alignment.ITCS 2022  13th Innovations in Theoretical Computer Science ConferenceBerkeley, United StatesJanuary 2022
 25 inproceedingsRobustness in MultiObjective Submodular Optimization: a Quantile Approach.Proceedings of the 39 th International Conference on Machine LearningICML 2022  39th International Conference on Machine LearningBaltimore, United StatesJuly 2022
 26 inproceedingsAsynchronous SGD Beats Minibatch SGD Under Arbitrary Delays.NeurIPS 2022  Thirtysixth Conference on Neural Information Processing SystemsAdvances in Neural Information Processing SystemsNew Orleans, United StatesNovember 2022
 27 inproceedingsDelayAware Decentralized Qlearning for Wind Farm Control.CDC 2022  IEEE 61st Conference on Decision and ControlCancun, MexicoIEEEDecember 2022, 807813
 28 inproceedingsMean field analysis of an incentive algorithm for a closed stochastic network.33rd International Conference on Probabilistic, Combinatorial and Asymptotic Methods for the Analysis of Algorithms (AofA 2022)AofA 2022225Leibniz International Proceedings in Informatics (LIPIcs)Philadelphia, PA, United StatesSchloss Dagstuhl – LeibnizZentrum für InformatikJune 2022, 13:113:17
 29 inproceedingsConvergence beyond the overparameterized regime using Rayleigh quotients.NeurIPS 2022  36th Conference on Neural Information Processing SystemNew Orleans, United StatesNovember 2022
 30 inproceedingsPeriodic Signal Recovery with Regularized Sine Neural Networks.Proceedings Track 2022 NeurIPS Workshop on Symmetry and Geometry in Neural RepresentationsNeurIPS 2022 Workshop on Symmetry and Geometry in Neural RepresentationsNew Orleans, United StatesDecember 2022
 31 inproceedingsConvergence Rates of NonConvex Stochastic Gradient Descent Under a Generic Łojasiewicz Condition and Local Smoothness.Proceedings of Machine Learning ResearchICML 2022  39th International Conference on Machine Learning162Baltimore, United StatesJuly 2022
 32 inproceedingsA Generalized Nash Equilibrium analysis of the interaction between a peertopeer financial market and the distribution grid.ROADEF 2022  23ème congrès annuel de la Société Française de Recherche Opérationnelle et d'Aide à la DécisionVilleurbanne  Lyon, FranceFebruary 2022
Doctoral dissertations and habilitation theses
 33 thesisWaveletbased representations of point processes for modelling and statistical learning.Ecole Normale Supérieure (Paris)March 2022
 34 thesisGeographic and sociodemographic disparities in oncology care pathways.Université ParisSaclayDecember 2022
 35 thesisThe graph alignment problem: fundamental limits and efficient algorithms.PSL Research University; Ecole normale supérieureSeptember 2022
Reports & preprints
 36 miscCoupling from the Past for the Null Recurrent Markov Chain.April 2022
 37 miscEnergy and Delay TradeOffs of EndtoEnd Vehicular Communications using a Hyperfractal Urban Modelling.January 2022
 38 miscMean field analysis of stochastic networks with reservation.January 2022
 39 miscModeling and Analysis of 2Tier Heterogeneous Vehicular Networks Leveraging Roadside Units and Vehicle Relays.April 2022
 40 miscStochastic dynamic matching: A mixed graphtheory and linearalgebra approach.February 2022
 41 miscFeature Projection for Optimal Transport.2022
 42 miscDisparities in accessibility to oncology care centers in France.January 2022
 43 miscPropagation of chaos and Poisson Hypothesis for replica meanfield models of intensitybased neural networks.November 2022
 44 miscSample Optimality and Allforall Strategies in Personalized Federated and Collaborative Learning.February 2022
 45 miscA new stochastic model for carsharing suited to freefloating.January 2023
 46 miscStatistical limits of correlation detection in trees.January 2023
 47 miscA Unified Framework for Generalizing the GromovHausdorff Metric.2022
 48 miscAn Improved Lower Bound on the Largest Common Subtree of Random LeafLabeled Binary Trees.2022
 49 miscOn NonLinear operators for Geometric Deep Learning.July 2022
 50 miscA Stackelberg Game Analysis of RiskHedging Strategies in Decentralized Electricity Markets.January 2023