A Caltech Library Service

Parallel inverse-problem solver for time-domain optical tomography with perfect parallel scaling

Gaggioli, E. L. and Bruno, Oscar P. (2022) Parallel inverse-problem solver for time-domain optical tomography with perfect parallel scaling. Journal of Quantitative Spectroscopy and Radiative Transfer, 290 . Art. No. 108300. ISSN 0022-4073. doi:10.1016/j.jqsrt.2022.108300.

[img] PDF - Submitted Version
See Usage Policy.


Use this Persistent URL to link to this item:


This paper presents an efficient parallel radiative transfer-based inverse-problem solver for time-domain optical tomography. The radiative transfer equation provides a physically accurate model for the transport of photons in biological tissue, but the high computational cost associated with its solution has hindered its use in time-domain optical-tomography and other areas. In this paper this problem is tackled by means of a number of computational and modeling innovations, including (1) A spatial parallel-decomposition strategy with perfect parallel scaling for the forward and inverse problems of optical tomography on parallel computer systems; and, (2) A Multiple Staggered Source method (MSS) that solves the inverse transport problem at a computational cost that is independent of the number of sources employed, and which significantly accelerates the reconstruction of the optical parameters: a six-fold MSS acceleration factor is demonstrated in this paper. Finally, this contribution presents (3) An intuitive derivation of the adjoint-based formulation for evaluation of functional gradients, including the highly-relevant general Fresnel boundary conditions—thus, in particular, generalizing results previously available for vacuum boundary conditions. Solutions of large and realistic 2D inverse problems are presented in this paper, which were produced on a 256-core computer system. The combined parallel/MSS acceleration approach reduced the required computing times by several orders of magnitude, from months to a few hours.

Item Type:Article
Related URLs:
URLURL TypeDescription Paper
Gaggioli, E. L.0000-0001-6939-6347
Bruno, Oscar P.0000-0001-8369-3014
Additional Information:© 2022 Elsevier. Received 17 February 2022, Revised 30 May 2022, Accepted 21 June 2022, Available online 24 June 2022, Version of Record 30 June 2022. This work was supported by NSF, DARPA and AFOSR through contracts DMS-2109831 and HR00111720035 and FA9550-21-1-0373, and by the NSSEFF Vannevar Bush Fellowship under contract number N00014-16-1-2808. ELG acknowledge financial support from CONICET. The codes utilized in this paper are available from the first author upon reasonable request. CRediT authorship contribution statement. E.L. Gaggioli: Conceptualization, Methodology, Software, Validation, Investigation, Visualization, Writing – original draft. Oscar P. Bruno: Conceptualization, Methodology, Validation, Investigation, Resources, Writing – original draft, Supervision, Funding acquisition. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Funding AgencyGrant Number
Defense Advanced Research Projects Agency (DARPA)HR00111720035
Air Force Office of Scientific Research (AFOSR)FA9550-21-1-0373
Vannever Bush Faculty FellowshipUNSPECIFIED
National Security Science and Engineering Faculty FellowshipN00014-16-1-2808
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)UNSPECIFIED
Record Number:CaltechAUTHORS:20220705-346498000
Persistent URL:
Official Citation:E.L. Gaggioli, Oscar P. Bruno, Parallel inverse-problem solver for time-domain optical tomography with perfect parallel scaling, Journal of Quantitative Spectroscopy and Radiative Transfer, Volume 290, 2022, 108300, ISSN 0022-4073,
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:115323
Deposited By: George Porter
Deposited On:07 Jul 2022 23:14
Last Modified:07 Jul 2022 23:14

Repository Staff Only: item control page