Physics‐constrained non‐Gaussian probabilistic learning on manifolds

Christian Soize (1) , Roger Ghanem (2)

1 MSME - Laboratoire de Modélisation et Simulation Multi Echelle
2 USC - University of Southern California

An extension of the probabilistic learning on manifolds (PLoM), recently introduced by the authors, has been presented: In addition to the initial data set given for performing the probabilistic learning, constraints are given, which correspond to statistics of experiments or of physical models. We consider a non-Gaussian random vector whose unknown probability distribution has to satisfy constraints. The method consists in constructing a generator using the PLoM and the classical Kullback-Leibler minimum cross-entropy principle. The resulting optimization problem is reformulated using Lagrange multipliers associated with the constraints. The optimal solution of the Lagrange multipliers is computed using an efficient iterative algorithm. At each iteration, the Markov chainMonte Carlo algorithm developed for the PLoM is used, consisting in solving an Itô stochastic differential equation that is projected on a diffusion-maps basis. The method and the algorithm are efficient and allow the construction of probabilistic models for high-dimensional problems from small initial data sets and for which an arbitrary number of constraints are specified. The first application is sufficiently simple in order to be easily reproduced. The second one is relative to a stochastic elliptic boundary value problem in high dimension.

Mots clés: probabilistic learning; Kullback-Leibler; data driven; statistical constraints; uncertainty quantification; machine learning

Domaines: Statistiques [math.ST]; Machine Learning [stat.ML]; Probabilités [math.PR]; Sciences de l'ingénieur [physics]

Christian Soize, Roger Ghanem. Physics‐constrained non‐Gaussian probabilistic learning on manifolds. International Journal for Numerical Methods in Engineering, 2020, 121 (1), pp.110-145. ⟨10.1002/nme.6202⟩. ⟨hal-02341991⟩