This is a Julia interface to the Fast Multipole Method (FMM) library FMMLIB2D by Leslie Greengard and Zydrunas Gimbutas. Currently v1.2.3 is used.
Documentation for the library can be found in the FMMLIB2D User's Guide.
This package currently provides interfaces to the FMM's for Laplace (real and complex), Helmholtz, and complex sums:
rfmm2dpartself, rfmm2dparttarg, lfmm2dpartself, lfmm2dparttarg, hfmm2dparttarg, zfmm2dparttarg
The most convenient way of calling them is through the Julia interfaces with keyword arguments, e.g.
x = rand(2, 10)
y = rand(2, 20)
q = rand(10) + 1im*rand(10)
U = lfmm2d(source=x, charge=q, target=y, ifgradtarg=true, tol=1e-9)U = rfmm2d(source::Array{Float64} = ...,
target::Array{Float64} = ...,
charge::Array{Float64} = ...,
dipstr::Array{Float64} = ...,
dipvec::Array{Float64} = ...,
tol::Float64 = ...,
ifpot::Bool = ...,
ifgrad::Bool = ...,
ifhess::Bool = ...,
ifpottarg::Bool = ...,
ifgradtarg::Bool = ...,
ifhesstarg::Bool = ...,
)Output format:
U.pot (Nsrc)
U.grad (2,Nsrc)
U.hess (3,Nsrc)
U.pottarg (Ntrg)
U.gradtarg (2,Ntrg)
U.hesstarg (3,Ntrg)U = lfmm2d(source::Array{Float64} = ...,
target::Array{Float64} = ...,
charge::Array{ComplexF64} = ...,
dipstr::Array{ComplexF64} = ...,
dipvec::Array{Float64} = ...,
tol::Float64 = ...,
ifpot::Bool = ...,
ifgrad::Bool = ...,
ifhess::Bool = ...,
ifpottarg::Bool = ...,
ifgradtarg::Bool = ...,
ifhesstarg::Bool = ...,
)Output format:
U.pot (Nsrc)
U.grad (2,Nsrc)
U.hess (3,Nsrc)
U.pottarg (Ntrg)
U.gradtarg (2,Ntrg)
U.hesstarg (3,Ntrg)U = hfmm2d(source::Array{Float64} = ...,
target::Array{Float64} = ...,
charge::Array{ComplexF64} = ...,
dipstr::Array{ComplexF64} = ...,
dipvec::Array{Float64} = ...,
tol::Float64 = ...,
zk::ComplexF64 = ...,
ifpot::Bool = ...,
ifgrad::Bool = ...,
ifhess::Bool = ...,
ifpottarg::Bool = ...,
ifgradtarg::Bool = ...,
ifhesstarg::Bool = ...,
)Output format:
U.pot (Nsrc)
U.grad (2,Nsrc)
U.hess (3,Nsrc)
U.pottarg (Ntrg)
U.gradtarg (2,Ntrg)
U.hesstarg (3,Ntrg)U = zfmm2d(source::Array{ComplexF64} = ...,
target::Array{ComplexF64} = ...,
dipstr::Array{ComplexF64} = ...,
tol::Float64 = 1e-15,
ifpot::Bool = true,
ifgrad::Bool = false,
ifhess::Bool = false,
ifpottarg::Bool = true,
ifgradtarg::Bool = false,
ifhesstarg::Bool = false
)Output format:
U.pot (Nsrc)
U.grad (Nsrc)
U.hess (Nsrc)
U.pottarg (Ntrg)
U.gradtarg (Ntrg)
U.hesstarg (Ntrg)