The fast multipole method for the direct E/MEG problem
نویسندگان
چکیده
Reconstructing neuronal activity from MEG and EEG measurements requires the accurate calculation of the electromagnetic field inside the head. The boundary element formulation of this problem leads to a dense linear system which is too large to be solved directly. We propose to accelerate the computations via the fast multipole method. This method approximates the electromagnetic interaction between surface elements by performing multipole expansions at a coarse resolution. It significantly reduces the computational complexity of the matrix-vector products needed for the iterative solution of the linear system, and avoids the storage of its matrix. We describe the single-level fast multipole method and present several experiments demonstrating its accuracy and performance.
منابع مشابه
Fast multipole acceleration of the MEG/EEG boundary element method.
The accurate solution of the forward electrostatic problem is an essential first step before solving the inverse problem of magneto- and electroencephalography (MEG/EEG). The symmetric Galerkin boundary element method is accurate but cannot be used for very large problems because of its computational complexity and memory requirements. We describe a fast multipole-based acceleration for the sym...
متن کاملA Comparative Study of Multipole and Empirical Relations Methods for Effective Index and Dispersion Calculations of Silica-Based Photonic Crystal Fibers
In this paper, we present a solid-core Silica-based photonic crystal fiber (PCF) composed of hexagonal lattice of air-holes and calculate the effective index and chromatic dispersion of PCF for different physical parameters using the empirical relations method (ERM). These results are compared with the data obtained from the conventional multipole method (MPM). Our simulation results reveal tha...
متن کاملA New Guideline for the Allocation of Multipoles in the Multiple Multipole Method for Two Dimensional Scattering from Dielectrics
A new guideline for proper allocation of multipoles in the multiple multipole method (MMP) is proposed. In an ‘a posteriori’ approach, subspace fitting (SSF) is used to find the best location of multipole expansions for the two dimensional dielectric scattering problem. It is shown that the best location of multipole expansions (regarding their global approximating power) coincides with the med...
متن کاملFast Multipole Method for Accelerating the Evaluation of Splines
We consider the problem of interpolating scattered data using spline methods and present a general framework of using the multipole method to accelerate the evaluation of splines. The method depends on a tree-data structure and two hierarchical approximations: an upward multipole expansion approximation and a downward local Taylor series approximation. We also illustrate the performance of the ...
متن کاملAnalytical integration of the moments in the diagonal form fast multipole boundary element method for 3-D acoustic wave problems
A diagonal form fast multipole boundary element method (BEM) is presented in this paper for solving 3-D acoustic wave problems based on the Burton–Miller boundary integral equation (BIE) formulation. Analytical expressions of the moments in the diagonal fast multipole BEM are derived for constant elements, which are shown to be more accurate, stable and efficient than those using direct numeric...
متن کامل