2024/11/18 by Laura Bazahica, Vesa Kaarnioja, Bazahica, Laura +3 · 1 citation
Decision Sciences · #FOS: Mathematics #Numerical Analysis (math.NA) #Probabilistic and Robust Engineering Design
paper · pdf · doi:10.48550/arxiv.2411.11538
openalex publication_date 2024/11/18 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
The theoretical development of quasi-Monte Carlo (QMC) methods for uncertainty quantification of partial differential equations (PDEs) is typically centered around simplified model problems such as elliptic PDEs subject to homogeneous zero Dirichlet boundary conditions. In this paper, we present a theoretical treatment of the application of randomly shifted rank-1 lattice rules to electrical impedance tomography (EIT). EIT is an imaging modality, where the goal is to reconstruct the interior conductivity of an object based on electrode measurements of current and voltage taken at the boundary of the object. This is an inverse problem, which we tackle using the Bayesian statistical inversion paradigm. As the reconstruction, we consider QMC integration to approximate the unknown conductivity given current and voltage measurements. We prove under moderate assumptions placed on the parameterization of the unknown conductivity that the QMC approximation of the reconstructed estimate has a dimension-independent, faster-than-Monte Carlo cubature convergence rate. Finally, we present numerical results for examples computed using simulated measurement data.