5,085 research outputs found
Higher order quasi-Monte Carlo for Bayesian shape inversion
In this article, we consider a Bayesian approach towards data assimilation and uncertainty quantification in diffusion problems on random domains. We provide a rigorous analysis of parametric regularity of the posterior distribution given that the data exhibit only limited smoothness. Moreover, we present a dimension truncation analysis for the forward problem, which is formulated in terms of the domain mapping method. Having these novel results at hand, we shall consider as a practical example Electrical Impedance Tomography in the regime of constant conductivities. We are interested in computing moments, in particular expectation and variance, of the contour of an unknown inclusion, given perturbed surface measurements. By casting the forward problem into the framework of elliptic diffusion problems on random domains, we can directly apply the presented analysis. This straightforwardly yields parametric regularity results for the system response and for the posterior measure, facilitating the application of higher order quadrature methods for the approximation of moments of quantities of interest. As an example of such a quadrature method, we consider here recently developed higher order quasi-Monte Carlo methods. To solve the forward problem numerically, we employ a fast boundary integral solver. Numerical examples are provided to illustrate the presented approach and validate the theoretical findings
Recent advances in higher order quasi-Monte Carlo methods
In this article we review some of recent results on higher order quasi-Monte
Carlo (HoQMC) methods. After a seminal work by Dick (2007, 2008) who originally
introduced the concept of HoQMC, there have been significant theoretical
progresses on HoQMC in terms of discrepancy as well as multivariate numerical
integration. Moreover, several successful and promising applications of HoQMC
to partial differential equations with random coefficients and Bayesian
estimation/inversion problems have been reported recently. In this article we
start with standard quasi-Monte Carlo methods based on digital nets and
sequences in the sense of Niederreiter, and then move onto their higher order
version due to Dick. The Walsh analysis of smooth functions plays a crucial
role in developing the theory of HoQMC, and the aim of this article is to
provide a unified picture on how the Walsh analysis enables recent developments
of HoQMC both for discrepancy and numerical integration
Bayesian Coronal Seismology
In contrast to the situation in a laboratory, the study of the solar
atmosphere has to be pursued without direct access to the physical conditions
of interest. Information is therefore incomplete and uncertain and inference
methods need to be employed to diagnose the physical conditions and processes.
One of such methods, solar atmospheric seismology, makes use of observed and
theoretically predicted properties of waves to infer plasma and magnetic field
properties. A recent development in solar atmospheric seismology consists in
the use of inversion and model comparison methods based on Bayesian analysis.
In this paper, the philosophy and methodology of Bayesian analysis are first
explained. Then, we provide an account of what has been achieved so far from
the application of these techniques to solar atmospheric seismology and a
prospect of possible future extensions.Comment: 19 pages, accepted in Advances in Space Researc
- …