Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

A fast parallel solver for the forward problem in electrical impedance tomography

Jehl, Markus ; Dedner, Andreas ; Betcke, Timo ; Aristovich, Kirill ; Klöfkorn, Robert LU orcid and Holder, David (2015) In IEEE Transactions on Biomedical Engineering 62(1). p.37-126
Abstract

Electrical impedance tomography (EIT) is a noninvasive imaging modality, where imperceptible currents are applied to the skin and the resulting surface voltages are measured. It has the potential to distinguish between ischaemic and haemorrhagic stroke with a portable and inexpensive device. The image reconstruction relies on an accurate forward model of the experimental setup. Because of the relatively small signal in stroke EIT, the finite-element modeling requires meshes of more than 10 million elements. To study the requirements in the forward modeling in EIT and also to reduce the time for experimental image acquisition, it is necessary to reduce the run time of the forward computation. We show the implementation of a parallel... (More)

Electrical impedance tomography (EIT) is a noninvasive imaging modality, where imperceptible currents are applied to the skin and the resulting surface voltages are measured. It has the potential to distinguish between ischaemic and haemorrhagic stroke with a portable and inexpensive device. The image reconstruction relies on an accurate forward model of the experimental setup. Because of the relatively small signal in stroke EIT, the finite-element modeling requires meshes of more than 10 million elements. To study the requirements in the forward modeling in EIT and also to reduce the time for experimental image acquisition, it is necessary to reduce the run time of the forward computation. We show the implementation of a parallel forward solver for EIT using the Dune-Fem C++ library and demonstrate its performance on many CPU's of a computer cluster. For a typical EIT application a direct solver was significantly slower and not an alternative to iterative solvers with multigrid preconditioning. With this new solver, we can compute the forward solutions and the Jacobian matrix of a typical EIT application with 30 electrodes on a 15-million element mesh in less than 15 min. This makes it a valuable tool for simulation studies and EIT applications with high precision requirements. It is freely available for download.

(Less)
Please use this url to cite or link to this publication:
author
; ; ; ; and
publishing date
type
Contribution to journal
publication status
published
subject
keywords
Algorithms, Diagnosis, Computer-Assisted/methods, Electric Impedance, Humans, Models, Biological, Numerical Analysis, Computer-Assisted, Plethysmography, Impedance/methods, Software
in
IEEE Transactions on Biomedical Engineering
volume
62
issue
1
pages
12 pages
publisher
IEEE - Institute of Electrical and Electronics Engineers Inc.
external identifiers
  • pmid:25069109
  • scopus:84919934580
ISSN
1558-2531
DOI
10.1109/TBME.2014.2342280
language
English
LU publication?
no
id
35a7ea3e-eb92-442c-9b89-ca6a33c04626
date added to LUP
2021-02-10 14:05:36
date last changed
2024-05-30 07:32:41
@article{35a7ea3e-eb92-442c-9b89-ca6a33c04626,
  abstract     = {{<p>Electrical impedance tomography (EIT) is a noninvasive imaging modality, where imperceptible currents are applied to the skin and the resulting surface voltages are measured. It has the potential to distinguish between ischaemic and haemorrhagic stroke with a portable and inexpensive device. The image reconstruction relies on an accurate forward model of the experimental setup. Because of the relatively small signal in stroke EIT, the finite-element modeling requires meshes of more than 10 million elements. To study the requirements in the forward modeling in EIT and also to reduce the time for experimental image acquisition, it is necessary to reduce the run time of the forward computation. We show the implementation of a parallel forward solver for EIT using the Dune-Fem C++ library and demonstrate its performance on many CPU's of a computer cluster. For a typical EIT application a direct solver was significantly slower and not an alternative to iterative solvers with multigrid preconditioning. With this new solver, we can compute the forward solutions and the Jacobian matrix of a typical EIT application with 30 electrodes on a 15-million element mesh in less than 15 min. This makes it a valuable tool for simulation studies and EIT applications with high precision requirements. It is freely available for download.</p>}},
  author       = {{Jehl, Markus and Dedner, Andreas and Betcke, Timo and Aristovich, Kirill and Klöfkorn, Robert and Holder, David}},
  issn         = {{1558-2531}},
  keywords     = {{Algorithms; Diagnosis, Computer-Assisted/methods; Electric Impedance; Humans; Models, Biological; Numerical Analysis, Computer-Assisted; Plethysmography, Impedance/methods; Software}},
  language     = {{eng}},
  number       = {{1}},
  pages        = {{37--126}},
  publisher    = {{IEEE - Institute of Electrical and Electronics Engineers Inc.}},
  series       = {{IEEE Transactions on Biomedical Engineering}},
  title        = {{A fast parallel solver for the forward problem in electrical impedance tomography}},
  url          = {{http://dx.doi.org/10.1109/TBME.2014.2342280}},
  doi          = {{10.1109/TBME.2014.2342280}},
  volume       = {{62}},
  year         = {{2015}},
}