The Departement of Mathematical Engineering organizes a series of seminars. The seminars are usually held on Tuesday from 2:00pm to 3:00pm in the Euler lecture room, Building EULER, av. Georges Lemaître 4-6, Louvain-la-Neuve (Parking 13). Be mindful that exceptions may occur; see the talk annoucements.
If you wish to receive the seminar announcements by email, please send an email to Etienne Huens.
Master students can take this seminar for credit in either of the two semesters; see LINMA2120 for more information.
Seminars to come
|21/06/2022 (14h) [Location: Euler building (room A.002)]|
Abdou Sene (Université virtuelle du Sénégal - UVS)
COVID-19 propagation mathematical modeling: the case of Senegal
The outburst of the COVID-19 pandemic has raised several questions leading to a complex system in terms of modeling. Indeed, the modeling of the epidemic, at the level of a country, needs considering each of the different sources of contamination as well as the public health authorities strategy, in a specific way. With this in mind, we have developed a mathematical model of the COVID-19 epidemic in Senegal. In the model, the population is subdivided into five compartments: susceptible, infected but asymptomatic, symptomatic, quarantined, and recovered immune people. In addition, due to its important impact in the propagation of the disease, we have added one more variable: the number of infected objects. Based on the senegalese territory COVID-19 data, we simulate various scenarios as for the evolution of the epidemic in the country, in order to predict the peak and its magnitude with regard to the application of barrier measures.
|07/06/2022 (14h) [Location: Euler building (room a.002)]|
Nelson Maculan (Federal University of Rio de Janeiro)
Mixed Integer Nonlinear Programming (MINLP) Models for the Euclidean Steiner Tree Problem in R^n
Abstract:New Mixed Integer Nonlinear Programming (MINLP) models for the Euclidean Steiner Tree Problem in R^n will be presented. The novelty of these models is the introduction of constraints that represent second-order cones, avoiding the problem of non-differentiability of continuous relaxation, which appears in other models. Computational results provided a more effective way to solve the Euclidean Steiner Tree Problem in R^n.
|10/05/2022 (14h) [Location: Euler building (room a.002)] (UCLouvain (EPL)) |
Presentation from Master's Students
Abstract:This seminar will feature the presentations by the students who followed the course LINMA2120 on their elected topics. The presentation will be organised in co-modal, as two students are currently abroad and will present remotely. The channel to follow the session on Teams is as usual (O365G-cours-linma2120-seminars).
|26/04/2022 (14h) [Location: Euler building (room a.002)]|
Tijl De Bie (UGent)
Automating data exploration
Abstract: Machine learning is increasingly becoming a commodity, available as a building block to programmers with minimal required understanding about detailed modelling aspects. In large part, this is thanks to advances in automating the labor-intensive aspects of machine learning. In contrast, the more exploratory tasks in data science still remain highly labor-intensive, often requiring a deeper understanding of the intricacies of the methods used. In this talk, I will survey some past results on the automation of data exploration tasks, and outline what I view as some of the more important challenges for the future.
|19/04/2022 (14h) [Location: Euler building (room a.002)]|
Raphaël Berthier (École Polytechnique Fédérale de Lausanne)
A Continuized View on Nesterov Acceleration for Stochastic Gradient Descent and Randomized Gossip
Abstract: We introduce the continuized Nesterov acceleration, a close variant of Nesterov acceleration whose variables are indexed by a continuous time parameter. The two variables continuously mix following a linear ordinary differential equation and take gradient steps at random times. This continuized variant benefits from the best of the continuous and the discrete frameworks: as a continuous process, one can use differential calculus to analyze convergence and obtain analytical expressions for the parameters; and a discretization of the continuized process can be computed exactly with convergence rates similar to those of Nesterov original acceleration. We show that the discretization has the same structure as Nesterov acceleration, but with random parameters. We provide continuized Nesterov acceleration under deterministic as well as stochastic gradients, with either additive or multiplicative noise. Finally, using our continuized framework and expressing the gossip averaging problem as the stochastic minimization of a certain energy function, we provide the first rigorous acceleration of asynchronous gossip algorithms.
|12/04/2022 (14h) [Location: Euler building (room a.002)]|
Peyman Esfahani (Delft Centre for Systems and Control (TU Delft))
Data-driven Decision-Making in Dynamic Environments
Abstract: In this seminar, we start with a broad class of anomaly detection for large-scale nonlinear dynamical systems. Noting a connection between the diagnosis filter and the so-called behavioral sets of dynamical systems, we leverage tools from the traditional model-based approaches and modern data-driven analytics to address the inherent complexity of the problem. We then shift our attention to the performance guarantees of our proposed solution. In this part, we study this topic in a general context of data-driven decision-making with a particular focus on the distributionally robust optimization framework. We will discuss the role of convexity from the different viewpoints of computational, statistical, and real-time implementation.
|29/03/2022 (14h) [Location: Euler building (room a.002)]|
Hugues Goosse (UCLouvain (ELI))
Reconstructing the climate of the past millennium by combining indirect climate observations and simulation results
Abstract: It is essential to characterize well past natural climate variations to detect the potential anthropogenic contribution in recent changes. Past variations also allow studying the dynamics of the climate system for a range of conditions broader than the one observed recently. The past millennium appears particularly interesting in this framework as the conditions are similar to those observed currently, except that the human impact was much lower. The reconstruction of past climate changes is based on indirect observations as the instrumental data cover less than 150 years in most regions. Different techniques exists. The goal here is to illustrate how climate models results and indirect records derived from natural archives can be combined to obtain estimates of past states of the climate system and of the mechanisms responsible for the past climate changes.
|22/03/2022 (14h) [Location: Euler building (room a.002)]|
Bálint Daróczy (UCLouvain (INMA))
Learning from pairwise comparisons and gradient representations of neural networks
Abstract: During the seminar we will consider two separate problems. First, we study of learning problems in which we would like to learn intrinsic values of objects based on pairwise comparisons. We suggest an algorithm and deterime a minimax rate and show that both the upper and the lower error bounds are connected to the trace of the Moore-Penrose inverse of the weighted Laplacian of the comparison graph. In the second part of the seminar, we consider feed-forward neural networks and investigate how the tangent space of the network can be exploited to refine the decision in case of ReLU (Rectified Linear Unit) activations by forming a similarity function with a Riemannian metric.
|15/03/2022 (14h) [Location: Euler building (room a.002)]|
Alessandro Abate (University of Oxford)
Certified learning, or learning for verification?
Abstract: We are witnessing an increased, inter-disciplinary convergence between areas underpinned by model-based reasoning and by data-driven learning. Work across these areas is not only scientifically justified, but also motivated by industrial applications where access to information-rich data has to be traded off with a demand for safety criticality: cyber-physical systems are exemplar applications. In this talk, I will report on ongoing initiatives in this cross-disciplinary domain. According to the dual perspective in the title of this talk, I will sketch, on the one hand, results where formal methods can provide certificates to learning algorithms, and on the other hand, results where learning can bolster formal verification and strategy synthesis objectives.
|08/03/2022 (14h) [Location: Euler building (room a.002)] |
Title: Towards tight convergence rates for optimization methods on hypoconvex functions
Speaker: Teodor Rotaru (UCLouvain/KULeuven)
Abstract: Using the framework of performance estimation (PEP), we established the first tight convergence rates of the gradient method for smooth hypoconvex (or weakly-convex) functions. These functions’ curvature (i.e., maximum Hessian eigenvalue) belongs to the interval [µ, L], where µ is negative and L is positive. With the help of PEP, we obtained mathematical proofs for a large range of step sizes. As a direct application, we recommend the optimal step size that minimizes the convergence rate.
Title: Accelerating large-scale Kernel Support Vector Machines
Speaker: Sofiane Tanji (UCLouvain)
Abstract:Kernel methods provide an elegant extension to well-known linear statistical learning. Due to their poor scalability in time and memory, they have limited applications in large scale learning. We propose Snacks : a kernel SVM solver which can tackle large-scale datasets. Our approach consists in using an accelerated version of the stochastic subgradient method to solve the primal optimization problem on smaller random subspaces. The computational savings do not lead to any degradation in terms of learning performance and we demonstrate the effectiveness of the proposed algorithm on benchmark datasets.
|1/03/2022 (14h) [Location: Euler building (room a.002)]|
Estelle Massart (UCLouvain,ICTEAM)
Exploring geometry for deep learning
Abstract: We propose to use stochastic Riemannian coordinate descent on the orthogonal group for recurrent neural network training. The algorithm rotates successively two columns of the recurrent matrix, an operation that can be efficiently implemented as a multiplication by a Givens matrix. In the case when the coordinate is selected uniformly at random at each iteration, we prove the convergence of the proposed algorithm under standard assumptions on the loss function, stepsize and minibatch noise. In addition, we numerically demonstrate that the Riemannian gradient in recurrent neural network training has an approximately sparse structure. Leveraging this observation, we propose a faster variant of the proposed algorithm that relies on the Gauss-Southwell rule. Experiments on a benchmark recurrent neural network training problem are presented to demonstrate the effectiveness of the proposed algorithm. Time permitting, I will also present another application of geometry in deep learning, namely, exploiting invariant curves in the loss landscape when training Wasserstein generative adversarial neural networks.
|22/02/2022 (14h) [Location: Euler building (room a.002)]|
Renaud Ronsse (UCLouvain,IMMC)
Adaptive oscillator as a template for the assistance of cyclical movements.
Abstract:In this lecture, I will review the theory of adaptive oscillators and how this tool has been initially developed by generating primitive locomotion movements in robotics. Then I will explain how this framework can be extended to provide different types of assistance in cyclical movements (mainly for the lower-limb, but not exclusively). I will further illustrate the coding of a simple adaptive oscillator on a haptic device that can be manipulated by the participants.
|15/02/2022 (14h) [Location: Euler building (room a.207)]|
Anne-Katrin Schmuck (Max Plank Institute for Software Systems, Kaiserslautern, Germany)
Let's play! - Solving Controller Synthesis Games for Cyber-Physical System Design
Abstract:Cyber-Physical Systems (CPS) are technical systems where a large software stack orchestrates the interaction of physical and digital components. Such systems are omnipresent in our daily life and their correct behavior is crucial. However, developing safe, reliable and performant CPS is challenging. A promising research direction towards this goal is the combination of formal methods from computer science and controller synthesis techniques from automation. In my talk, I will discuss how infinite two-player games over finite graphs, originating from the formal methods community, can be utilized and enhanced for higher layer control of CPS. In particular, I will discuss how the use of environment assumptions - used to model particularities of the system under control within these games - has to be rethought in order to effectively solve controller synthesis tasks for CPS.
|08/02/2022 (14h) [Location: Euler building (room a.002)] |
Title: Extension of the Performance Estimation framework through constraint interpolation.
Speaker: Anne Rubbens(PhD UCLouvain/INMA)
Abstract: Selecting parameters involved in an optimization algorithm is a major challenge often without satisfactory solution, and can have a dramatic impact on the algorithm performance. For instance, selecting parameters by optimizing performance bounds that are not exact may lead to misguided choices.Recently, an approach called Performance Estimation Problem (PEP) has been developed to automatically compute exact worst-case bounds on the performance of a wide class of optimization algorithms whose input belongs to a given class of functions, for instance convex functions.This project consists in extending the field of application of this approach to classes of functions for which no satisfactory PEP formulation exists yet. PEP formulations rely on necessary and sufficient interpolation conditions, ensuring the existence of a function of a given class interpolating a finite set of data. Hence, we will provide a novel approach to derive interpolation conditions for various classes of functions.
Title: Optimization of compositional functions
Speaker: Nizar Bousselmi(PhD UCLouvain/INMA)
Abstract:We consider functions with compositional structure, i.e., functions that can be expressed using one or several compositions of simple elementary component functions. These functions naturally appear in many practical problems. We would like to analyze the impact of the compositional structure of objective functions on the performance of the existing methods and to develop new methods exploiting this structure.
Title: Automatic quality control of weather and climate time series
Speaker: Benoît Loucheur(PhD UCLouvain/INMA)
Abstract:In Belgium, the Royal Meteorological Institute (RMI) is the national meteorological service that provide weather and climate services based on observations and scientific research. The RMI collects and archives meteorological observations in Belgium since the 19th century. Currently, air temperature is monitored in Belgium in about 30 synoptic automatic weather stations as well as in 110 manual climatological stations. All observations are routinely checked for errors, inconsistencies and missing values by the RMI staff. Misleading data are corrected and gaps are filled by estimations. This quality control tasks require a lot of human intervention. With the forthcoming deployment of low-cost weather stations and the subsequent increase in the volume of data to verify, the process of data quality control and completion should become as automated as much as possible. The aim of this project is to develop algorithmic tools to address this quality control problem in a fully automatic way.
|14/12/2021 (14h) [Location: Euler building (room a.002)]|
Nina Miolane (UC Santa Barbara)
Geomstats: a Python Package for Riemannian Geometry in Statistics and Machine Learning
Abstract: We introduce Geomstats, an open-source Python package for computations and statistics on nonlinear manifolds that appear in machine learning applications, such as: hyperbolic spaces, spaces of symmetric positive definite matrices, Lie groups of transformations, and many more. We provide object-oriented and extensively unit-tested implementations. Manifolds come equipped with families of Riemannian metrics with associated exponential and logarithmic maps, geodesics, and parallel transport. Statistics and learning algorithms provide methods for estimation, clustering, and dimension reduction on manifolds. All associated operations provide support for different execution backends --- namely NumPy, Autograd, PyTorch, and TensorFlow. This talk presents the package, compares it with related libraries, and provides relevant examples. We show that Geomstats provides reliable building blocks to both foster research in differential geometry and statistics and democratize the use of Riemannian geometry in statistics and machine learning. The source code is freely available under the MIT license at https://github.com/geomstats/geomstats.
|07/12/2021 (14h) [Location: Euler building (room a.002)]|
Pierre Ablin (CNRS, Université Paris-Dauphine)
The symbiotic relationship between optimization and deep learning
Abstract: Optimization is one of the cornerstones of deep learning: most deep neural networks are trained by optimizing a cost function. The purpose of this talk is to cover other fruitful interactions between the fields of optimization and deep learning that are less obvious. First, I will discuss how deep neural networks can be designed to quickly and efficiently solve classical optimization problems, with a focus on the Lasso. I will then argue that the tools of optimization allow us to describe accurately what weights are learned by these networks during the training process. In a second part, I will discuss how classical optimization ideas like momentum acceleration can be translated to deep learning, and allow us to develop novel architectures that perform well and are memory-efficient.
|30/11/2021 (14h) [Location: Euler building (room a.002)]|
Geovani Nunes Grapiglia (UCLouvain/INMA)
A Generalized Worst-Case Complexity Analysis for Non-Monotone Line Searches
Abstract: In this talk we discuss the worst-case complexity of a wide class of non-monotone line search methods for non-convex unconstrained minimization problems. For the algorithms in this class, the non monotonicity is controlled by a sequence of nonnegative parameters. We prove complexity bounds to achieve approximate first-order optimality even when this sequence is not summable. As a by-product, we obtain a unified global convergence result. Our generalized results allow more freedom for the development of new non-monotone line search algorithms. As an example, we design a non-monotone scheme related to the Metropolis rule. Preliminary numerical experiments suggest that the new method is suitable to nonconvex problems with many non-global local minimizers.
|23/11/2021 (14h) [Location: Euler building (room a.002)]|
Flavio Abreu Araujo (UCLouvain)
Implementing spintronic based neuromorphic computing hardware under the reservoir computing approach
Abstract: The brain displays many signatures of non-linear dynamical behavior, including synchronization and complex transient behavior. These observations have inspired a whole class of neuromorphic concepts based on complex networks of interconnected non-linear nodes. Such non-linearity has been identified as one of the main ingredients to achieve excellent performance in cognitive tasks such as spoken digit recognition. In this work, we have quantified in detail the contribution of the acoustic filtering and the neural network, respectively, for spoken digit recognition task using three different frequency decomposition methods: Cochlear, MFCC and Spectrogram. In a first step, we have demonstrated that Cochlear and MFCC are powerful stand-alone features extractors, and they can achieve for themselves very high recognition level: up to 95.8% and 77.2% for cochlear and MFCC, respectively. We have found that such high recognition level is mainly due to the non-linear character of these frequency decomposition methods. First, we have investigated the non- linear dependence of the Spectrogram showing a huge increase of recognition rate from 10% (linear) to 85.6%. In a second step, we have evaluated the gain of the recognition rate provided by the neural network. For simplicity, we have modeled a neural network based on the non- linear dynamics of oscillators in the framework of the reservoir computing approach. The reservoir is generally composed by a large number of fixed and random interconnected non- linear oscillators which generates very complex non-linear dynamics. The key insight behind reservoir computing is that only the external connections (between the reservoir and the output layer) should be trained to obtain the desirable target. Reservoir computing has been identified to be very suitable for different hardware implementations: optical, photonic and spintronic devices. We have found that the contribution of the neural network is dominant for the linear Spectrogram filter but not in the other two cases, i.e., cochlear and MFCC. Finally, we have carried out experiments using non-linear and tunable spin-torque nano-oscillators exhibiting an excellent agreement with our simulations.
|16/11/2021 (14h) [Location: Euler building (room a.002)]|
Radu Dragomir (UCL/INMA)
Gradient methods with Bregman distances
Abstract: Large-scale optimization problems from signal processing and machine learning are typically solved with gradient methods, because of their low cost per iteration and their simplicity. In this talk, we study a generalization of the standard gradient descent, which consists in replacing the Euclidean distance by a more general Bregman divergence induced by some simple convex reference function. This function is chosen to be adapted to the geometry of the problem at hand through the so-called relative smoothness condition. We present some advances in this recent line of work, including the study of worst-case complexity through performance estimation problems (PEPs), as well as applications to low-rank matrix optimization.
|09/11/2021 (14h) [Location: Euler building (room a.002)]|
Renato Vizuete (CentraleSupélec and GIPSA-lab, France )
Graphons for the Estimation of Performance Indices in Large Networks
Abstract: In recent years, the analysis of large networks has received increasing attention in many scientific fields due to the continuous evolution of the world towards a networked environment with a large number of connections. In this type of networks, uncertainties are almost present, and even the graph topology may not be available, which makes the analysis more complicated. One of the most promising tools to address these problems are graphons, defined as the limits of convergent sequences of dense graphs. In this talk, we will present the fundamentals of graphons and their applications in the estimation of performance indices associated with large networks sampled from graphons. In the first part, we will analyze the stability of a SIS epidemic over a network and its robustness to noise using the properties of the graphon operator. For the second application, we will present a study of the spectrum of the Laplacian matrix of a graph using the degree function of the graphon and the estimation of the well known average effective resistance of sampled networks.
|02/11/2021 (14h) [Location: Euler building (room a.002)] |
Title: Using photogrammetry for the objective study of ancient bowed instruments: a machine learning approach
Speaker: Philémon Beghin(UCLouvain/INMA)
Abstract: The morphology of today’s violin differs greatly from that of the first instruments of the late 16th century. Indeed, in order to meet the standards suggested by famous orchestras and conservatories, many ancient violins have been recut. It is important for musicologists, violin makers and museum curators to analyse the alterations they have undergone. Specialists agree that the instruments have been reshaped, but have difficulty to prove it rigorously. Moreover, the historical testimonies about this process are imprecise. It is therefore necessary to find an objective way to quantify violin geometry. This project aims to develop a set of algorithmic tools in order to compute an adequate mathematical representation that is able to describe the 3D shape of an object (violin) acquired by photogrammetry. Based on that representation, machine learning techniques will be applied in order to perform clustering and classification with the aim of quantifying their geometric characteristics, their possible anomalies and, if applicable, their original morphology. From an engineering point of view, the development of 3D analyses of objects as complex as violins will allow the growth of methods and knowledge potentially applicable to many fields outside of organology.
Title: Biomechanics of tactile feedback during dexterous manipulation of objects
Speaker: Donatien Doumont(UCLouvain/INMA)
Abstract: The biomechanics of the skin and underlying tissues play a fundamental role in the sense of touch. Indeed, when the finger comes in contact with an object, the deformations induced within these living tissues are translated by the mechanoreceptors (or tactile sensors) into neural signals, which are then interpreted by the brain to generate the appropriate motor response. This project is aimed at better understanding the feedback provided by the tactile afferents during interactions with objects. To that end, we will combine passive stimulation tasks mimicking the mechanical interactions during object manipulation and active manipulation tasks while monitoring the activity of tactile afferents (only in the passive case) and several key biomechanical parameters such as external forces, surface skin deformation, and skin moisture. This combination of techniques will enable us to gain insight into the essential aspects of tactile feedback that enable stable and dexterous control of fingertip forces during manipulation. Such insight will have direct implications for the development of sensorized bionic hands.
Title: Simulation of the stride-to-stride variability of patients with Parkinson’s disease
Speaker: Clémence Vandamme(UCLouvain/INMA)
Abstract:Gait is a complex mechanism involving several neural structures such as the motor cortex, the cerebellum and basal ganglia. It also requires the integration of sensorial feedback from visual, vestibular and peripheral receptors. The complex coordination of these structures is not yet fully understood. In particular, recent studies revealed that gait fluctuations exhibit Long Range Autocorrelation (LRA) such that fluctuations at any given moment are statistically related to those that occur over many different time scales. Up to now, the physiological origin of LRA remains unclear. Moreover, this specific correlation structure is altered in elderly and in people suffering from neurodegenerative diseases, such as Parkinson’s disease. Consequently, it is thought that an optimal level of LRA is a marker of stable gait and high adaptive capabilities. Clinically, LRA indicators could complement the diagnostic of Parkinson's disease and constitutes a precautious marker of risk of fall. This project aims to provide a comprehensive model that includes a continuous-time control of the biomechanics plant for both a healthy population and Parkinson’s patients, as it is an important challenge for fundamental understanding as well as for clinical applications.
Title: Influence of vestibular, visual and somesthesic inputs on dexterous manipulation
Speaker: Simon Vandergooten(UCLouvain/INMA)
Abstract:Since birth we evolve in a steady gravitational environment in which our brain has learned to manipulate objects. For example, we are able to unconsciously optimize the force with which we squeeze objects when handling them. In order to better characterize the influence of gravity on upper limb movement kinematics as well as on the dynamics of prehension (i.e. the act of grasping), several parabolic flight campaigns as well as ground-based experiments were performed. These experiments ultimately led to GRIP, an experiment carried out in the International Space Station that focuses on long-term motor adaptation to microgravity, in particular during dexterous object manipulation. On the ground, our current setup comprising a rotating chair, four motion-tracking cameras and the grip-lift manipulandum allows us to perform control experiments that are necessary for the interpretation of and analysis of GRIP data. With this project, we aim to study how visual, vestibular and somatosensory feedbacks interact and affect the planning of movement trajectory, the perception of the verticality and the anticipative force mechanisms underlying finger-arm coordination.
|26/10/2021 (14h) [Location: Euler building (room a.002)] |
Title: Employing neural network based control models to understand how brain generates movements
Speaker: Hari Teja Kalidindi (ICTEAM/IoNS)
Abstract:Biological agents display impressive abilities to move under dynamic environmental conditions. Distributed regions in the central nervous system coordinate to generate the motor commands suitable for a given behavioral goal. However, the neural computations that underlie even the simplest of the movements are still equivocal. Historically, one of the main sources of the dispute has been due to the emphasis on describing the motor-encoding in isolated brain regions, by ignoring the effect of key control elements such as sensory-feedback, prediction and the physics of the body under control. In this talk, we demonstrate the utility of mechanistic neural networks as normative models to study neural computations underlying movement control. Particularly, we emulate recent observations from the primate motor cortex - the brain region that is implicated for generating movements.
Title: Sensing of low-rank plus sparse matrices
Speaker: Simon Vary (ICTEAM)
Abstract:Expressing a matrix as the sum of a low-rank matrix plus a sparse matrix is a flexible model capturing global and local features in data. This model is the foundation of robust principle component analysis, and popularized by dynamic-foreground/static-background separation amongst other applications. In this talk we develop guarantees showing that rank-r plus sparsity-s matrices can be recovered by computationally tractable methods from p=O(r(m+n-r)+s)log(mn/s) linear measurements. We establish that the restricted isometry constants for the low-rank plus sparse matrix set remain bounded independent of the problem size provided p/mn, s/p, and r(m+n-r)/p remain fixed. The developed theory and algorithms also apply to the fully observed case of Robust PCA.
|19/10/2021 (14h) [Location: Euler building (room a.002)]|
David Wozabal (Technical University of Munich)
Multi-Stage Stochastic Programming for AC Optimal Power Flow Problems
Abstract:We propose the first computationally tractable framework to solve multi-stage stochastic optimal power flow (OPF) problems in alternating current (AC) power systems. To this end, we use recent results on convex semi-definite programming (SDP) relaxations of OPF problems in order to adapt the SDDP algorithm for problems with a Markovian structure employing scenario lattices to discretize the underlying randomness. We show that the usual SDDP lower bound remains valid and that the algorithm converges to a globally optimal solution of the stochastic AC-OPF problem as long as the SDP relaxations are tight. In the last part of the paper, we set up an extensive case study demonstrating the practical viability of our approach. In particular, we use the IEEE RTS-GMLC network to set up a storage sitting, sizing, and operation problem under uncertainty about demand and renewable generation. We show that the convex SDP relaxation of the stochastic problem is tight except in very rare cases. Furthermore, we demonstrate that after a reasonable number of iterations the algorithm finds a policy with a relatively small SDDP optimality gap that yields a significant added value over rolling deterministic planning.
|12/10/2021 (14h) [Location: Euler building (room a.002)]|
Vivian De Smedt (PSI Metals Belgium)
Optimization in Steel Industry
Abstract:What are some of the optimization problems of the steel industry? We will discuss the differences between two approaches of the optimization problems: black box vs. expert systems. The size of the solution space being very large some heuristic approaches are needed. We will discuss how to select the heuristics and what are the strong and weak points of the different approaches. The optimization project have to be integrated into workflow that less automatic and have its strong points. We will discuss how to take into account this aspect of the problematic in the conception of the solution
|05/10/2021 (14h) [Location: Euler building (room a.002)]|
Guillaume Drion (Université de Liège)
Neuromorphic control principles
Abstract:Owing to the recent advances in control engineering and machine learning, and combined with the remarkable improvement of sensors, actuators and computer power, modern high-performance computing systems far surpass human performance in a plethora of complex tasks. However, they consume megawatts of power, are optimized for specific tasks, and are hardly portable. In sharp contrast, biological brains are energy-efficient, are polyvalent, and show impressive adaptive capabilities in uncertain environments. These main differences between brains and current computing architectures are a crucial bottleneck for the expansion of automation in modern society. This situation has led many universities and world-leading companies, which include IBM, Intel, IMEC, or Thales, to investigate novel, brain-inspired computing systems and artificial intelligence technologies, an approach called neuromorphic computing. In this talk, we will approach the design of neuromorphic computing systems from a control engineering approach. We will first study how excitability, a key property of neuronal signaling, can be analyzed and designed following a nonlinear, multiscale feedback approach. Secondly, we will see how such neuromorphic systems can be controlled through loop shaping via a physiological mechanism called neuromodulation. Finally, we will exploit these brain-inspired mechanisms to improve artificial neural network adaptivity and long-term memory on the one hand, and design robust and controllable neuromorphic electronic systems on the other hand
|28/9/2021 (14h) [Location: Euler building (room a.002) ]|
Josh Taylor (University of Toronto)
Convex Optimization of Bioprocesses
Abstract:In this talk, we begin by optimizing the gradostat, in which several chemostats are interconnected by mass flow and diffusion. The gradostat is of interest both as a classical nonlinear system and because the basic network structure and nonlinearities appear in a wide variety of bioprocesses, including wastewater treatment. We formulate a convex relaxation of the gradostat. The relaxation is exact under several conditions, for instance, if the gradostat is outflow connected and its flow matrix is irreducible. When the microbial growth in the bioreactors is described by the Monod or Contois functions, the relaxation is a second-order cone program, which can be solved at scales of over 10^5 variables in minutes with industrial software. We also discuss how to extend the work to a general class of bioprocesses, and present an example based on wastewater treatment.
|21/9/2021 (14h) [Location: Euler building (room a.002) ]|
Balázs Gerencser ( Eötvös Loránd University, Hungary)
From theoretical to computable convergence rate of push-sum for consensus
Abstract:We know that reaching average consensus using only local communication along a network is a fundamental building block in the area of distributed computing, leading to applications such as sensor fusion and distributed optimization. We currently analyze schemes based on one such protocol, push-sum, and our target is to understand their convergence speed. We prove a bound on the almost sure convergence rate that is also computable for a class of push-sum algorithms. This extends the works of Iutzeler, Ciblat and Hachem (2013) on similar bounds but in a more restrictive setup and conclusion, and complements the results of Gerencsér and Gerencsér (2019) identifying the exact convergence rate but providing no computable access or approximation