COOKIES: By using this website you agree that we can place Google Analytics Cookies on your device for performance monitoring. |
University of Cambridge > Talks.cam > Isaac Newton Institute Seminar Series > Fast sampling from parameterised Gaussian random fields
Fast sampling from parameterised Gaussian random fieldsAdd to your list(s) Download to your calendar using vCal
If you have a question about this talk, please contact INI IT. This talk has been canceled/deleted Gaussian random fields are popular models for spatially varying uncertainties, arising for instance in geotechnical engineering, hydrology or image processing. A Gaussian random field is fully characterised by its mean function and covariance operator. In more complex models these can also be partially unknown. Then we need to handle a family of Gaussian random fields indexed with hyperparameters. Sampling for a fixed configuration of hyperparameters is already very expensive, as it requires a Cholesky or spectral decomposition of the discretised covariance operator, which is in general a large dense matrix. Sampling from multiple configurations increases the total computational cost severely. In this report we construct a reduced basis surrogate for parameterised Karhunen-Loève expansions – upon which our sampling procedure relies. The reduced basis is built using snapshots of Karhunen-Loève eigenvectors. In particular, we consider Matern type covariance operators with unknown correlation length and standard deviation. We suggest a linearisation of the covariance function and describe the associated online-offline decomposition. In numerical experiments, we investigate the approximation error of the reduced eigenpairs. As an application, we consider forward uncertainty propagation and Bayesian inversion in an elliptic PDE , where the log of the diffusion coefficient is a parameterised Gaussian random field. We discretise PDE operators and covariance operators with finite elements. In the Bayesian inverse problem we employ a Markov Chain Monte Carlo method on the reduced space to generate samples from the posterior measure. All numerical experiments are done in 2D space with non-separable covariance operators on finite element grids with tens of thousands of degrees of freedom. This talk is part of the Isaac Newton Institute Seminar Series series. This talk is included in these lists:This talk is not included in any other list Note that ex-directory lists are not shown. |
Other listsIsotope Coffee: Geochemistry and Petrology Seminars Department of Earth Sciences Physics of Living Matter lectures https://data.mendeley.com/datasets?... Traduire cette page N Boudemagh. N Boudemagh. Contribution: PhD, network ASSET MANAGEMENT. 07 Nov 2016 in: Smart Transportation. aPPLIED MATHEMA. Viewed.Other talksParkinson's Rehabilitation using interactive Dance Technology Assessing the Impact of Open IP in Emerging Technologies CANCELLED First year PhD student fieldwork seminar South American Opuntioids Hide and seek: medieval creatures on the manuscript page MOVED TO 28 JUNE 2018 It takes two to tango:platelet collagen receptor GPVI-dimer in thrombosis and clinical implications |