Dear Kaustians, We are excited to announce the upcoming Stochastic Numerics and Statistical Learning: Theory and Applications Workshop 2024, taking place at KAUST, Building 9, from May 19-30, 2024. Following the highly successful 2022 and 2023 editions, this year's workshop promises to be another engaging and insightful event for researchers, faculty members, and students interested in stochastic algorithms, statistical learning, optimization, and approximation. The 2024 workshop aims to build on the achievements of last year's event, which featured 28 talks, two mini-courses, and two poster sessions, attracting over 150 participants from various universities and research institutes. In previous two years attendees had the opportunity to learn from through insightful talks, interactive mini-courses, and vibrant poster sessions. This year, the workshop will once again showcase contributions that offer mathematical foundations for algorithmic analysis or highlight relevant applications. Confirmed speakers include renowned experts from institutions such as Ecole Polytechnique, EPFL, Université Pierre et Marie Curie - Paris VI, and Imperial College London, among others.

We will discuss the solution of eigenvalue problems associated with partial differential equations (PDE)s that can be written in the generalised form Ax = λMx, where the matrices A and/or M may depend on a scalar parameter. Parameter dependent matrices occur frequently when stabilised formulations are used for the numerical approximation of PDEs. With the help of classical numerical examples we will show that the presence of one (or both) parameters can produce unexpected results.

Hessian-dependent functionals play a pivotal role in a wide latitude of problems in mathematics. Arising in the context of differential geometry and probability theory, this class of problems find applications in the mechanics of deformable media (mostly in elasticity theory) and the modelling of slow viscous fluids. We study such functionals from three distinct perspectives.

In this work, we employ importance sampling (IS) techniques to track a small over-threshold probability of a running maximum associated with the solution of a stochastic differential equation (SDE) within the framework of ensemble Kalman filtering (EnKF).

In this work, we analyze the convergence rate of randomized quasi-Monte Carlo (RQMC) methods under Owen's boundary growth condition [Owen, 2006] via spectral analysis.

We present a theoretical analysis of the Weak Adversarial Networks (WAN) method, recently proposed in [1, 2], as a method for approximating the solution of partial differential equations in high dimensions and tested in the framework of inverse problems. In a very general abstract framework.

Analogies between financial markets and critical phenomena have long been observed empirically. So far, no convincing theory has emerged that can explain these empirical observations. Here, we take a step towards such a theory by modeling financial markets as a lattice gas.

Wave propagation and scattering problems are of huge importance in many applications in science and engineering - e.g., in seismic and medical imaging and more generally in acoustics and electromagnetics.

In this short course, we will introduce some elements in deriving the hp a posteriori error estimate for a high-order unfitted finite element method for elliptic interface problems. The key ingredient is an hp domain inverse estimate, which allows us to prove a sharp lower bound of the hp a posteriori error estimator.

Dear Kaustians, We are excited to announce the upcoming Stochastic Numerics and Statistical Learning: Theory and Applications Workshop 2023, taking place at KAUST, Building 9, from May 21 to June 1, 2023. Following the highly successful 2022 edition, this year's workshop promises to be another engaging and insightful event for researchers, faculty members, and students interested in stochastic algorithms, statistical learning, optimization, and approximation. The 2023 workshop aims to build on the achievements of last year's event, which featured 28 talks, two mini-courses, and two poster sessions, attracting over 150 participants from various universities and research institutes. In 2022, attendees had the opportunity to learn from through insightful talks, interactive mini-courses, and vibrant poster sessions. This year, the workshop will once again showcase contributions that offer mathematical foundations for algorithmic analysis or highlight relevant applications. Confirmed speakers include renowned experts from institutions such as Ecole Polytechnique, EPFL, Université Pierre et Marie Curie - Paris VI, and Imperial College London, among others.

This scientific meeting will concentrate on stochastic algorithms and their rigorous numerical analysis for various problems, including statistical learning, optimization, and approximation. Stochastic algorithms are valuable tools when addressing challenging computational problems.

This thesis studies novel and efficient computational sampling methods for applications in three types of stochastic inversion problems: seismic waveform inversion, filtering problems, and static parameter estimation.

In biochemically reactive systems with small copy numbers of one or more reactant molecules, stochastic effects dominate the dynamics. In the first part of this thesis, we design novel efficient simulation techniques, based on multilevel Monte Carlo methods and importance sampling, for a reliable and fast estimation of various statistical quantities for stochastic biological and chemical systems under the framework of Stochastic Reaction Networks (SRNs). In the second part of this thesis, we design novel numerical methods for pricing financial derivatives. Option pricing is usually challenging due to a combination of two complications: 1) The high dimensionality of the input space, and 2) The low regularity of the integrand on the input parameters. We address these challenges by using different techniques for smoothing the integrand to uncover the available regularity and improve quadrature methods' convergence behavior. We develop different ways of smoothing that depend on the characteristics of the problem at hand. Then, we approximate the resulting integrals using hierarchical quadrature methods combined with Brownian bridge construction and Richardson extrapolation.

Optimal experimental design for parameter estimation is a fast-growing area of research. Let us consider the experimental goal to be the inference of some attributes of a complex system using measurement data of some chosen system responses, and the optimal designs are those that maximize the value of measurement data. The value of data is quantified by the expected information gain utility, which measures the informativeness of an experiment. Often, a mathematical model is used that approximates the relationship between the system responses and the model parameters acting as proxies for the attributes of interest.

Hierarchical porous media that feature properties and processes at multiple scales arise in many engineering applications including the design of novel materials for energy storage devices. Microscopic (pore-scale) properties of the media impact their macroscopic (continuum- or Darcy-scale) counterparts and understanding the relationships between processes on these two scales is essential for informing engineering decision tasks.

I​n the field of uncertainty quantification, the effects of parameter uncertainties on scientific simulations may be studied by integrating or approximating a quantity of interest as a function over the parameter space. If this is done numerically, using regular grids with a fixed resolution, the required computational work increases exponentially with respect to the number of uncertain parameters -- a phenomenon known as the curse of dimensionality.

The probability that a sum of random variables (RVs) exceeds (respectively falls below) a given threshold, is often encountered in the performance analysis of wireless communication systems. Generally, a closed-form expression of the sum distribution does not exist and a naive Monte Carlo (MC) simulation is computationally expensive when dealing with rare events. An alternative approach is represented by the use of variance reduction techniques, known for their efficiency in requiring less computations for achieving the same accuracy requirement.

This work employs statistical and Bayesian techniques to analyze mathematical forward models with several sources of uncertainty. The forward models usually arise from phenomenological and physical phenomena and are expressed through regression-based models or partial differential equations (PDEs) associated with uncertain parameters and input data. One of the critical challenges in real-world applications is, for any proposed model, to estimate and quantify uncertainties of its unknown parameters using the available observations.

A feasible, nonparametric, probabilistic approach for modeling and quantifying model-form uncertainties associated with a High-Dimensional Computational Model (HDM) and/or a corresponding Hyperreduced Projection-based Reduced-Order Model (HPROM) [1,2] designed for the solution of computational mechanics problems is presented.

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Ms. Soumaya ElKantassi is a master student in Applied Mathematics and Computational Science and a member of the Stochastic Numerics Group at KAUST. She obtained her Bachelor's degree in Engineering from Tunisia Polytechnic School majoring in Economic and Scientific Management.

Most problems in engineering and natural sciences involve parametric equations in which the parameters are not known exactly due to measurement errors, lack of measurement data, or even intrinsic variability. In such problems, one objective is to compute point or aggregate values, called "quantities of interest".

Abstract:

Time inconsistent performance crit