As a sequel to our previous paper on extending the Fast Multipole Method (FMM) for charges inside a dielectric sphere [J. Comput. Phys. 223 (2007) 846-864], this paper further extends the FMM to the electrostatic calculation for charges inside a dielectric sphere immersed in an ionic solvent, a scenery with more relevance in biological applications. The key findings include two fourth-order mul...

The Fast Marching Method (FMM) is an efficient technique to numerically solve the eikonal equation. The parallelisation of the FMM is not easy because of its intrinsic sequential nature. In this paper we propose a novel approach to parallelise the FMM which is not based on a classic domain-decomposition procedure. Compared to other techniques in the field, our method is much simpler to implemen...

A new fast multipole formulation for the hypersingular BIE (HBIE) for 2D elasticity is presented in this paper based on a complex-variable representation of the kernels, similar to the formulation developed earlier for the conventional BIE (CBIE). A dual BIE formulation using a linear combination of the developed CBIE and HBIE is applied to analyze multi-domain problems with thin inclusions or ...

We develop a fast numerical method to approximate the solutions of a wide class of equations associated to the Shape From Shading problem. Our method, which is based on the control theory and the interfaces propagation, is an extension of the “Fast Marching Method” (FMM) [30,25]. In particular our method extends the FMM to some equations for which the solution is not systematically decreasing a...

A symmetric Galerkin boundary-element method is used for the solution of boundary-value problems with mixed boundary conditions of Dirichlet and Neumann type. As a model problem we consider the Laplace equation. When an iterative scheme is employed for solving the resulting linear system, the discrete boundary integral operators are realized by the fast multipole method. While the single-layer ...

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 be...