A Bayesian multifidelity method using a fast parallel GPU eikonal solver enabled patient-specific uncertainty quantification of cardiac activation sequences in less than 5 minutes.
A novel Bayesian multifidelity computational approach significantly reduces the time required for patient-specific uncertainty quantification in cardiac electrophysiology to under 5 minutes, meeting clinical time constraints.
We present a fast, patient-specific methodology for uncertainty quantification in electrophysiology, aimed at meeting the time constraints of clinical practitioners. We focus on computing the statistics of the activation map, given the uncertainties associated with the conductivity tensor modeling the fiber orientation in the heart. We use a fast parallel solution method implemented on a graphics processing unit for the eikonal approximation, in order to compute the activation map and to sample the random fiber field with correlation on the basis of geodesic distances. While this enables to perform uncertainty quantification studies with a manageable computational effort, the required time frame still exceeds clinically suitable time expectations. In order to reduce it further by 2 orders of magnitude, we rely on Bayesian multifidelity methods. In particular, we propose a low-fidelity model that is patient-specific and free from the additional training cost associated with reduced models. This is achieved by a sound physics-based simplification of the full eikonal model. The low-fidelity output is then corrected by the standard multifidelity framework. In practice, the complete procedure only requires approximately 100 new runs of our eikonal graphics processing unit solver for producing the sought estimates and their associated credible intervals, enabling a full online analysis in less than 5 minutes.
Quaglino et al. (Mon,) conducted a other in cardiac electrophysiology. Bayesian multifidelity methods with GPU eikonal solver was evaluated on computational time for uncertainty quantification. A Bayesian multifidelity method using a fast parallel GPU eikonal solver enabled patient-specific uncertainty quantification of cardiac activation sequences in less than 5 minutes.