Kernel Methods for Bayesian Elliptic Inverse Problems on Manifolds

Abstract

This paper investigates the formulation and implementation of Bayesian inverse problems to learn input parameters of partial differential equations (PDEs) defined on manifolds. Specifically, we study the inverse problem of determining the diffusion coefficient of a second-order elliptic PDE on a closed manifold from noisy measurements of the solution. Inspired by manifold learning techniques, we approximate the elliptic differential operator with a kernel-based integral operator that can be discretized via Monte-Carlo without reference to the Riemannian metric. The resulting computational method is mesh-free and easy to implement, and can be applied without full knowledge of the underlying manifold, provided that a point cloud of manifold samples is available. We adopt a Bayesian perspective to the inverse problem, and establish an upper-bound on the total variation distance between the true posterior and an approximate posterior defined with the kernel forward map. Supporting numerical results show the effectiveness of the proposed methodology.

References in corpus (5)

• Galerkin Approximation of Dynamical Quantities using Trajectory Data
• Convergence of the Point Integral method for Poisson equation on point cloud
• Greedy inference with structure-exploiting lazy maps
• Approximating solutions of linear elliptic PDE's on a smooth manifold using local kernels
• On the Consistency of Graph-based Bayesian Learning and the Scalability of Sampling Algorithms