A partially reflecting random walk on spheres algorithm for electrical impedance tomography

Sylvain Maire 1 Martin Simon 2
1 TOSCA - TO Simulate and CAlibrate stochastic models
CRISAM - Inria Sophia Antipolis - Méditerranée , IECL - Institut Élie Cartan de Lorraine : UMR7502
Abstract : In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias and the variance of the new estimator both theoretically and experimentally. In a second step, the variance is considerably reduced via a novel control variate conditional sampling technique.
Type de document :
Article dans une revue
Journal of Computational Physics, Elsevier, 2015, 303, 〈10.1016/j.jcp.2015.10.005〉
Liste complète des métadonnées

Littérature citée [44 références]  Voir  Masquer  Télécharger

https://hal.inria.fr/hal-01253538
Contributeur : Sylvain Maire <>
Soumis le : mercredi 13 janvier 2016 - 20:20:11
Dernière modification le : samedi 27 janvier 2018 - 01:30:41
Document(s) archivé(s) le : jeudi 10 novembre 2016 - 22:41:09

Fichier

PRWOS.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Collections

Citation

Sylvain Maire, Martin Simon. A partially reflecting random walk on spheres algorithm for electrical impedance tomography. Journal of Computational Physics, Elsevier, 2015, 303, 〈10.1016/j.jcp.2015.10.005〉. 〈hal-01253538〉

Partager

Métriques

Consultations de la notice

215

Téléchargements de fichiers

78