Séminaire de mathématiques appliquées (archives)

Nom de l'orateur
Jerome Stenger
Etablissement de l'orateur
EDF et IMT
Lieu de l'exposé
Salle des séminaires
Date et heure de l'exposé

La quantification d'incertitudes a pour but d'évaluer l'impact d'un manque de connaissance des paramètres d'entrées (considérés aléatoires) sur les résultats d'une expérience numérique. Dans ce travail, nous prenons en compte un second niveau d'incertitude qui affecte le choix du modèle probabiliste des paramètres d'entrées. Nous évaluons les bornes d'une quantité d'intérêt sur l'ensemble des mesures de probabilités uniquement définies par leur bornes et certains de leurs moments. Du fait du grand nombre de contraintes, l'optimisation numérique est complexe. Nous montrons que le problème d'optimisation peut se paramétriser sur les points extrémaux de cet espace de mesures de probabilité contraintes.

Nom de l'orateur
Lionel Riou-Durand
Etablissement de l'orateur
University of Warwick
Lieu de l'exposé
Salle des séminaires
Date et heure de l'exposé

Sampling approximations for high dimensional statistical models often rely on so-called gradient-based MCMC algorithms. It is now well established that these samplers scale better with the dimension than other state of the art MCMC samplers, but are also more sensitive to tuning [5]. Among these, Hamiltonian Monte Carlo is a widely used sampling method shown to achieve gold standard d^{​​​​​1/4}​​​​​ scaling with respect to the dimension [1]. However it is also known that its efficiency is quite sensible to the choice of integration time, see e.g. [4], [2]. This problem is related to periodicity in the autocorrelations induced by the deterministic trajectories of Hamiltonian dynamics.

Nom de l'orateur
Raed Blel
Etablissement de l'orateur
ENPC
Lieu de l'exposé
Salle des séminaires
Date et heure de l'exposé

The main focus of this article is to provide a mathematical study of the algorithm proposed in [6] where the authors proposed a variance reduction technique for the computation of parameter-dependent expectations using a reduced basis paradigm. We study the effect of Monte-Carlo sampling on the the- oretical properties of greedy algorithms. In particular, using concentration inequalities for the empirical measure in Wasserstein distance proved in [14], we provide sufficient conditions on the number of samples used for the computation of empirical variances at each iteration of the greedy procedure to guarantee that the resulting method algorithm is a weak greedy algorithm with high probability.

Nom de l'orateur
Camilla Fiorini
Etablissement de l'orateur
M2N, Conservatoire National des Arts et Métiers
Lieu de l'exposé
Salle des séminaires
Date et heure de l'exposé

In this work we consider the surface quasi-geostrophic (SQG) system under location uncertainty (LU) and propose a Milstein-type scheme for these equations, which is then used in a multi-step method. The LU framework, is based on the decomposition of the Lagrangian velocity into two components: a large-scale smooth component and a small-scale stochastic one. This decomposition leads to a stochastic transport operator, and one can, in turn, derive the stochastic LU version of every classical fluid-dynamics system.

SQG in particular consists of one partial differential equation, which models the stochastic transport of the buoyancy, and an operator which relies the velocity and the buoyancy.

Nom de l'orateur
Michael Fanuel
Etablissement de l'orateur
Université de Lille
Lieu de l'exposé
Date et heure de l'exposé

Determinantal Point Processes (DPPs) elegantly model repulsive point patterns. A natural problem is the estimation of a DPP given a few samples. Parametric and nonparametric inference methods have been studied in the finite case, i.e. when the point patterns are sampled in a finite ground set. In the continuous case, several parametric methods have been proposed but nonparametric methods have received little attention. In this talk, we discuss a nonparametric approach for continuous DPP estimation leveraging recent advances in kernel methods. We show that a restricted version of this maximum likelihood (MLE) problem falls within the scope of a recent representer theorem for nonnegative functions in a Reproducing Kernel Hilbert Space.

Nom de l'orateur
Alexandre Poulain
Etablissement de l'orateur
Simula Research Laboratory, Oslo
Lieu de l'exposé
Zoom Planet
Date et heure de l'exposé

The Cahn-Hilliard equation, arising from physics, describes the phase separation occurring in a material during a sudden cooling process and is the subject of many pieces of research [2]. An interesting application of this equation is its capacity to model cell populations undergoing attraction and repulsion effects. For this application, we consider a variant of the Cahn-Hilliard equation with a single-well potential and a degenerate mobility. This particular form introduces numerous di culties especially for numerical simulations. We propose a relaxation of the equation to tackle these issues and analyze the resulting system. Interestingly, this relaxed version of the degenerate Cahn-Hilliard equation bears some similarity with a nonlinear Keller-Segel model.

Nom de l'orateur
Titouan Vayer
Etablissement de l'orateur
ENS Lyon
Lieu de l'exposé
Salle de séminaire
Date et heure de l'exposé

Abstract: Nowadays large-scale machine learning faces a number of fundamental computational challenges, triggered by the high dimensionality of modern data and the increasing availability of very large training collections. These data can also be of a very complex nature, such as such as those described by the graphs that are integral to many application areas. In this talk I will present some solutions to these problems. I will introduce the Compressive Statistical Learning (CSL) theory, a general framework for resource-efficient large scale learning in which the training data is summarized in a small single vector (called sketch) that captures the information relevant to the learning task.

Nom de l'orateur
Josselin Massot
Etablissement de l'orateur
IRMAR - Université de Rennes 1
Lieu de l'exposé
Salle Eole
Date et heure de l'exposé

Dans cet exposé, nous étudierons un plasma électronique où les particules peuvent être distribuées en deux populations distinctes, froides et chaudes, menant au modèle de Vlasov-Maxwell hybride fluide/cinétique linéarisé, restreint ici à 1 dimension en espace et 3 en vitesse. Notre objectif sera de proposer deux méthodes numériques pour résoudre ce modèle. La première est basée sur la structure hamiltonienne du système, et la seconde utilise un intégrateur exponentielle (ou méthode de Lawson), permettant facilement de monter en ordre en retirant une contrainte de stabilité provenant de la partie linéaire du problème. Nous étudierons ensuite la possibilité d'approximer l'exponentielle de la partie linéaire lorsqu'il n'est pas possible de déterminer celle-ci formellement.

Nom de l'orateur
Philipp Trunschke
Etablissement de l'orateur
Technische Universität Berlin
Lieu de l'exposé
Zoom
Date et heure de l'exposé

We consider best approximation problems in a nonlinear subset of a Banach space of functions. The norm is assumed to be a generalization of the L2-norm for which only a weighted Monte Carlo estimate can be computed. We establish error bounds for the empirical best approximation error in this general setting and use these bounds to derive a new, sample efficient algorithm for the model set of low-rank tensors. The viability of this algorithm is demonstrated by recovering quantities of interest for a classical random partial differential equation.