CLERC, M. - KERIVEN, R. - FAUGERAS, O. - KYBIC, J. - PAPADOPOULO, Th.: The Fast Multipole Method for the direct E/MEG problem. In Proceedings of the First 2002 IEEE International Symposium on Biomedical Imaging: Macro to Nano (ISBI'02). New York: IEEE Press, 2002, p. 1023-1026. ISBN 0-7803-7585-8. DOI: 10.1109/ISBI.2002.1029438
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.