A Fast Multipole Method for the Bessel potential
The fast multipole method (FMM) proposed by Greengard and Rokhlin provides a method by which the O(N-squared) many-body problem can be reduced to O(N) complexity. In this project, a multipole method is developed to calculate the energy of a system of vortices in a high temperature superconductor, where the many-body interactions give rise to rich and complex physics. The method developed here is suitable for systems where the interactions are governed by a Bessel potential rather than the usual logarithmic potentials occurring in gravitational and electrostatic problems. We derive and apply vectorised forms of the Gegenbauer addition formulae in order to achieve the O(N) scaling associated with fast multipole methods.
Physical Systems and Engineering simulation: Superconductivity
Algorithms and computational methods: Multipole methods
Visualisation and data handling software: Gnuplot
Software Engineering Tools: SVN