CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python

DOI10.1088/1361-6420/AD22E8arXiv2305.16951WikidataQ129283736 ScholiaQ129283736MaRDI QIDQ6149902FDOQ6149902

Authors: Nicolai A. B. Riis, Babak Maboudi Afkham, Felipe Uribe, Per Christian Hansen, Jakob S. Jørgensen 

Publication date: 5 March 2024

Published in: Inverse Problems (Search for Journal in Brave)

Abstract: Inverse problems, particularly those governed by Partial Differential Equations (PDEs), are prevalent in various scientific and engineering applications, and uncertainty quantification (UQ) of solutions to these problems is essential for informed decision-making. This second part of a two-paper series builds upon the foundation set by the first part, which introduced CUQIpy, a Python software package for computational UQ in inverse problems using a Bayesian framework. In this paper, we extend CUQIpy's capabilities to solve PDE-based Bayesian inverse problems through a general framework that allows the integration of PDEs in CUQIpy, whether expressed natively or using third-party libraries such as FEniCS. CUQIpy offers concise syntax that closely matches mathematical expressions, streamlining the modeling process and enhancing the user experience. The versatility and applicability of CUQIpy to PDE-based Bayesian inverse problems are demonstrated on examples covering parabolic, elliptic and hyperbolic PDEs. This includes problems involving the heat and Poisson equations and application case studies in electrical impedance tomography (EIT) and photo-acoustic tomography (PAT), showcasing the software's efficiency, consistency, and intuitive interface. This comprehensive approach to UQ in PDE-based inverse problems provides accessibility for non-experts and advanced features for experts.

Full work available at URL: https://arxiv.org/abs/2305.16951

zbMATH Keywords

softwareBayesian statisticsuncertainty quantificationpartial differential equationsprobabilistic programmingcomputational imaging

Mathematics Subject Classification ID

Statistics (62-XX) Theory of computing (68Qxx)

Cites Work

• \texttt{Pyclaw}: accessible, extensible, scalable tools for wave propagation problems
• Automated solution of differential equations by the finite element method. The FEniCS book
• Firedrake, automating the finite element method by composing abstractions
• hIPPYlib
• The no-U-turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo
• Statistical and computational inverse problems.
• Title not available (Why is that?)
• Existence and Uniqueness for Electrode Models for Electric Current Computed Tomography
• A first course in finite elements. With CD-ROM.
• Discrete inverse problems. Insight and algorithms.
• Numerical approximation of partial differential equations
• MCMC methods for functions: modifying old algorithms to make them faster
• Approximation of Bayesian Inverse Problems for PDEs
• A computational framework for infinite-dimensional Bayesian inverse problems. I: The linearized case, with application to global seismic inversion
• The Bayesian formulation of EIT: analysis and algorithms
• A Bayesian level set method for geometric inverse problems
• A computational framework for infinite-dimensional Bayesian inverse problems. II: stochastic Newton MCMC with application to ice sheet flow inverse problems
• Uncertainty Quantification of Inclusion Boundaries in the Context of X-Ray Tomography
• Iterative updating of model error for Bayesian inversion
• An introduction to data analysis and uncertainty quantification for inverse problems
• Detecting stochastic inclusions in electrical impedance tomography
• PDE-constrained optimization in medical image analysis
• Learning physics-based models from data: perspectives from inverse problems and model reduction
• BEYOND BLACK-BOXES IN BAYESIAN INVERSE PROBLEMS AND MODEL VALIDATION: APPLICATIONS IN SOLID MECHANICS OF ELASTOGRAPHY
• CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python

Cited In (2)

• CUQIpy: I. Computational uncertainty quantification for inverse problems in Python
• CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python

This page was built for publication: CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6149902)