Microwave Tomographic Imaging of Cerebrovascular Accidents by Using High-Performance Computing

The motivation of this work is the detection of cerebrovascular accidents by microwave tomographic imaging. This requires the solution of an inverse problem relying on a minimization algorithm (for example, gradient-based), where successive iterations consist in repeated solutions of a direct problem. The reconstruction algorithm is extremely computationally intensive and makes use of efficient parallel algorithms and high-performance computing. The feasibility of this type of imaging is conditioned on one hand by an accurate reconstruction of the material properties of the propagation medium and on the other hand by a considerable reduction in simulation time. Fulfilling these two requirements will enable a very rapid and accurate diagnosis. From the mathematical and numerical point of view, this means solving Maxwell's equations in time-harmonic regime by appropriate domain decomposition methods, which are naturally adapted to parallel architectures.

[1]  Frédéric Nataf,et al.  High performance domain decomposition methods on massively parallel architectures with freefem++ , 2012, J. Num. Math..

[2]  Xiao-Chuan Cai,et al.  A Restricted Additive Schwarz Preconditioner for General Sparse Linear Systems , 1999, SIAM J. Sci. Comput..

[3]  Martin J. Gander,et al.  Optimized Multiplicative, Additive, and Restricted Additive Schwarz Preconditioning , 2007, SIAM J. Sci. Comput..

[4]  B.D. Van Veen,et al.  An overview of ultra-wideband microwave imaging via space-time beamforming for early-stage breast-cancer detection , 2005, IEEE Antennas and Propagation Magazine.

[5]  A. Abbosh,et al.  Novel Preprocessing Techniques for Accurate Microwave Imaging of Human Brain , 2013, IEEE Antennas and Wireless Propagation Letters.

[6]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[7]  Michel Barlaud,et al.  Gradient and Newton-Kantorovich Methods for Microwave Tomography , 1997 .

[8]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[9]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2004, Future Gener. Comput. Syst..

[10]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[11]  J.C. Lin,et al.  Microwave imaging of cerebral edema , 1982, Proceedings of the IEEE.

[12]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2002, Future Gener. Comput. Syst..

[13]  Bernhard Seiser,et al.  Electromagnetic tomography for brain imaging: From virtual to human brain , 2014, 2014 IEEE Conference on Antenna Measurements & Applications (CAMA).

[14]  Frédéric Hecht,et al.  Scalable domain decomposition preconditioners for heterogeneous elliptic problems , 2014, Sci. Program..

[15]  Eric de Sturler,et al.  Recycling Krylov Subspaces for Sequences of Linear Systems , 2006, SIAM J. Sci. Comput..

[16]  Aurora Torrente,et al.  Brain Stroke Detection by Microwaves Using Prior Information from Clinical Databases , 2013 .

[17]  D. Corfield,et al.  Microwave Tomography for Brain Imaging: Feasibility Assessment for Stroke Detection , 2008 .

[18]  Amin M. Abbosh,et al.  Fast Frequency-Based Multistatic Microwave Imaging Algorithm With Application to Brain Injury Detection , 2016, IEEE Transactions on Microwave Theory and Techniques.

[19]  Frédéric Nataf,et al.  Scalable domain decomposition preconditioners for heterogeneous elliptic problems , 2013, 2013 SC - International Conference for High Performance Computing, Networking, Storage and Analysis (SC).

[20]  Ralf Hiptmair,et al.  Multilevel solution of the time‐harmonic Maxwell's equations based on edge elements , 1999 .

[21]  Jean Roman,et al.  SCOTCH: A Software Package for Static Mapping by Dual Recursive Bipartitioning of Process and Architecture Graphs , 1996, HPCN Europe.

[22]  Victorita Dolean,et al.  An introduction to domain decomposition methods - algorithms, theory, and parallel implementation , 2015 .

[23]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..