Jan Kybic et al 2005 Phys. Med. Biol. 50 4695 doi:10.1088/0031-9155/50/19/018
Jan Kybic1,2, Maureen Clerc2, Olivier Faugeras2, Renaud Keriven2 and Théo Papadopoulo2
Show affiliationsThe 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 symmetric boundary element method (BEM). It creates a hierarchical structure of the elements and approximates far interactions using spherical harmonics expansions. The accelerated method is shown to be as accurate as the direct method, yet for large problems it is both faster and more economical in terms of memory consumption.
87.50.C- Static and low-frequency electric and magnetic fields effects
02.70.Dh Finite-element and Galerkin methods
02.60.Lj Ordinary and partial differential equations; boundary value problems
65N30 Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Issue 19 (7 October 2005)
Received 6 March 2005, in final form 17 June 2005
Published 21 September 2005
Jan Kybic et al 2005 Phys. Med. Biol. 50 4695
A Konin and R Raguotis 2000 J. Phys.: Condens. Matter 12 9163
James Atkinson et al 2008 J. Phys. A: Math. Theor. 41 142001
Miju Kang et al. 2009 ApJ 701 454
R. Indebetouw et al. 2005 ApJ 619 931
V I Korobov 2004 J. Phys. B: At. Mol. Opt. Phys. 37 2331
P W Forder 1999 J. Phys. B: At. Mol. Opt. Phys. 32 4581
Subarna Mitra et al 2006 J. Phys. D: Appl. Phys. 39 4228
Richard B Larson 2003 Rep. Prog. Phys. 66 1651
Xinglong Liu and Zhongwei Jiang 2009 Smart Mater. Struct. 18 125017