Skip to content Skip to navigation
University of Warwick
  • Study
  • |
  • Research
  • |
  • Business
  • |
  • Alumni
  • |
  • News
  • |
  • About

University of Warwick
Publications service & WRAP

Highlight your research

  • WRAP
    • Home
    • Search WRAP
    • Browse by Warwick Author
    • Browse WRAP by Year
    • Browse WRAP by Subject
    • Browse WRAP by Department
    • Browse WRAP by Funder
    • Browse Theses by Department
  • Publications Service
    • Home
    • Search Publications Service
    • Browse by Warwick Author
    • Browse Publications service by Year
    • Browse Publications service by Subject
    • Browse Publications service by Department
    • Browse Publications service by Funder
  • Help & Advice
University of Warwick

The Library

  • Login
  • Admin

Geometric MCMC for infinite-dimensional inverse problems

Tools
- Tools
+ Tools

Beskos, Alexandros, Girolami, Mark, Lan, Shiwei, Farrell, Patrick E. and Stuart, A. M. (2017) Geometric MCMC for infinite-dimensional inverse problems. Journal of Computational Physics, 335 . pp. 327-351. doi:10.1016/j.jcp.2016.12.041

[img]
Preview
PDF
WRAP-geometric-MCMC-infinite-dimensional-inverse-Stuart-2017.pdf - Accepted Version - Requires a PDF viewer.
Available under License Creative Commons Attribution Non-commercial No Derivatives 4.0.

Download (4Mb) | Preview
Official URL: http://dx.doi.org/10.1016/j.jcp.2016.12.041

Request Changes to record.

Abstract

Bayesian inverse problems often involve sampling posterior distributions on infinite-dimensional function spaces. Traditional Markov chain Monte Carlo (MCMC) algorithms are characterized by deteriorating mixing times upon mesh-refinement, when the finite-dimensional approximations become more accurate. Such methods are typically forced to reduce step-sizes as the discretization gets finer, and thus are expensive as a function of dimension. Recently, a new class of MCMC methods with mesh-independent convergence times has emerged. However, few of them take into account the geometry of the posterior informed by the data. At the same time, recently developed geometric MCMC algorithms have been found to be powerful in exploring complicated distributions that deviate significantly from elliptic Gaussian laws, but are in general computationally intractable for models defined in infinite dimensions. In this work, we combine geometric methods on a finite-dimensional subspace with mesh-independent infinite-dimensional approaches. Our objective is to speed up MCMC mixing times, without significantly increasing the computational cost per step (for instance, in comparison with the vanilla preconditioned Crank–Nicolson (pCN) method). This is achieved by using ideas from geometric MCMC to probe the complex structure of an intrinsic finite-dimensional subspace where most data information concentrates, while retaining robust mixing times as the dimension grows by using pCN-like methods in the complementary subspace. The resulting algorithms are demonstrated in the context of three challenging inverse problems arising in subsurface flow, heat conduction and incompressible flow control. The algorithms exhibit up to two orders of magnitude improvement in sampling efficiency when compared with the pCN method.

Item Type: Journal Article
Subjects: Q Science > QA Mathematics
Divisions: Faculty of Science > Mathematics
Library of Congress Subject Headings (LCSH): Algorithms, Markov processes, Monte Carlo method, Hilbert space
Journal or Publication Title: Journal of Computational Physics
Publisher: Academic Press Inc. Elsevier Science
ISSN: 0021-9991
Official Date: 15 April 2017
Dates:
DateEvent
15 April 2017Published
13 December 2016Accepted
Volume: 335
Page Range: pp. 327-351
DOI: 10.1016/j.jcp.2016.12.041
Status: Peer Reviewed
Publication Status: Published
Access rights to Published version: Restricted or Subscription Access
RIOXX Funder/Project Grant:
Project/Grant IDRIOXX Funder NameFunder ID
UNSPECIFIEDLeverhulme Trusthttp://dx.doi.org/10.13039/501100000275
EP/K034154/1Engineering and Physical Sciences Research Councilhttp://dx.doi.org/10.13039/501100000266
EP/J016934/2Engineering and Physical Sciences Research Councilhttp://dx.doi.org/10.13039/501100000266
EP/K030930/1Engineering and Physical Sciences Research Councilhttp://dx.doi.org/10.13039/501100000266
EP/M019721/1 Engineering and Physical Sciences Research Councilhttp://dx.doi.org/10.13039/501100000266
UNSPECIFIEDNorges Forskningsrådhttp://dx.doi.org/10.13039/501100005416
UNSPECIFIEDOffice of Naval Researchhttp://dx.doi.org/10.13039/100000006

Request changes or add full text files to a record

Repository staff actions (login required)

View Item View Item

Downloads

Downloads per month over past year

View more statistics

twitter

Email us: wrap@warwick.ac.uk
Contact Details
About Us