Parallel implementation of the accelerated BEM approach for EMSI of the human brain

Akalin-Acar, Z.
Acar, C. E.
Gençer, Nevzat Güneri
Boundary element method (BEM) is one of the numerical methods which is commonly used to solve the forward problem (FP) of electro-magnetic source imaging with realistic head geometries. Application of BEM generates large systems of linear equations with dense matrices. Generation and solution of these matrix equations are time and memory consuming. This study presents a relatively cheap and effective solution for parallel implementation of the BEM to reduce the processing times to clinically acceptable values. This is achieved using a parallel cluster of personal computers on a local area network. We used eight workstations and implemented a parallel version of the accelerated BEM approach that distributes the computation and the BEM matrix efficiently to the processors. The performance of the solver is evaluated in terms of the CPU operations and memory usage for different number of processors. Once the transfer matrix is computed, for a 12,294 node mesh, a single FP solution takes 676 ms on a single processor and 72 ms on eight processors. It was observed that workstation clusters are cost effective tools for solving the complex BEM models in a clinically acceptable time.


Evaluation of multivariate adaptive non-parametric reduced-order model for solving the inverse electrocardiography problem: a simulation study
Onak, Onder Nazim; Serinağaoğlu Doğrusöz, Yeşim; Weber, Gerhard Wilhelm (Springer Science and Business Media LLC, 2019-05-01)
In the inverse electrocardiography (ECG) problem, the goal is to reconstruct the heart's electrical activity from multichannel body surface potentials and a mathematical model of the torso. Over the years, researchers have employed various approaches to solve this ill-posed problem including regularization, optimization, and statistical estimation. It is still a topic of interest especially for researchers and clinicians whose goal is to adopt this technique in clinical applications. Among the wide range of...
A Kalman filter-based approach to reduce the effects of geometric errors and the measurement noise in the inverse ECG problem
Aydin, Umit; Serinağaoğlu Doğrusöz, Yeşim (Springer Science and Business Media LLC, 2011-09-01)
In this article, we aimed to reduce the effects of geometric errors and measurement noise on the inverse problem of Electrocardiography (ECG) solutions. We used the Kalman filter to solve the inverse problem in terms of epicardial potential distributions. The geometric errors were introduced into the problem via wrong determination of the size and location of the heart in simulations. An error model, which is called the enhanced error model (EEM), was modified to be used in inverse problem of ECG to compens...
Non-destructive recognition of dielectric coated conducting objects by using WD type time-frequency transformation and PCA-based fusion
Sayan, Gönül (Wiley, 2013-07-01)
This article demonstrates the applications of a non-destructive electromagnetic target recognition method, called Wigner distribution-principal component analysis (WD-PCA) method, to dielectric coated conducting spheres. These spheres are chosen to be highly similar having the same overall size but slightly different permittivity and thickness values in coating layers. Four different classifiers are simulated by using the WD-PCA method for varying sizes of object libraries under different noise conditions. ...
Novel electromagnetic surface integral equations for highly accurate computations of dielectric bodies with arbitrarily low contrasts
Ergül, Özgür Salih (Elsevier BV, 2008-12-01)
We present a novel stabilization procedure for accurate surface formulations of electromagnetic scattering problems involving three-dimensional dielectric objects with arbitrarily low contrasts. Conventional surface integral equations provide inaccurate results for the scattered fields when the contrast of the object is low, i.e., when the electromagnetic material parameters of the scatterer and the host medium are close to each other. We propose a stabilization procedure involving the extraction of nonradi...
Implementation and comparison of reconstruction algorithms for magnetic resonance-electric impedance tomography (mr-eit)
Martin Lorca, Dario; Eyüboğlu, Behçet Murat; Department of Biomedical Engineering (2007)
In magnetic resonance electrical impedance tomography (MR-EIT), crosssectional images of a conductivity distribution are reconstructed. When current is injected to a conductor, it generates a magnetic field, which can be measured by a magnetic resonance imaging (MRI) scanner. MR-EIT reconstruction algorithms can be grouped into two: current density based reconstruction algorithms (Type-I) and magnetic flux density based reconstruction algorithms (Type-II). The aim of this study is to implement a series of r...
Citation Formats
Y. ATASEVEN, Z. Akalin-Acar, C. E. Acar, and N. G. Gençer, “Parallel implementation of the accelerated BEM approach for EMSI of the human brain,” MEDICAL & BIOLOGICAL ENGINEERING & COMPUTING, pp. 671–679, 2008, Accessed: 00, 2020. [Online]. Available: