Abstract:A variety of problems in computational physics and engineering require the convolution of the heat kernel (a Gaussian) with either discrete sources, densities supported on boundaries, or continuous volume distributions. We present a unified fast Gauss transform for this purpose in two dimensions, making use of an adaptive quad-tree discretization on a unit square which is assumed to contain all sources. Our implementation permits either free-space or periodic boundary conditions to be imposed, and is efficient for any choice of variance in the Gaussian.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: In the fields of computational physics and engineering, many problems require solving the Heat Equation, especially those involving the convolution of the Gaussian Kernel with discrete sources, densities on boundaries or continuous volume distributions. Although the traditional Fast Gauss Transform (FGT) can accelerate these calculations, it has limitations when dealing with different types of sources (such as discrete point sources, volume sources and densities on boundaries) in two - dimensional space, especially in adapting to different scales and boundary conditions.
To solve these problems, the author proposes an Adaptive Fast Gauss Transform algorithm (AFGT), which can efficiently handle different types of sources in two - dimensional space and is applicable to free - space or periodic boundary conditions. Specifically, the main objectives of this algorithm include:
1. **Improve computational efficiency**: By introducing the adaptive Quad - tree discretization method, the algorithm can automatically adjust the grid resolution according to the distribution of sources, thereby reducing unnecessary computational load.
2. **Handle multiple source types**: It can not only handle discrete point sources, but also volume sources and densities on boundaries, providing a unified solution.
3. **Adapt to different boundary conditions**: Whether it is free - space or periodic boundary conditions, this algorithm can efficiently handle them to ensure the accuracy of calculation results.
4. **Maintain insensitivity to Gaussian variance**: Regardless of how the variance of the Gaussian kernel changes, the algorithm can maintain efficient computational performance.
### Summary of Mathematical Formulas
- **Gaussian Kernel**:
\[
G(x, t)=\frac{e^{-\|x\|^{2} / 4t}}{(4\pi t)^{d / 2}}
\]
where \(d\) is the spatial dimension.
- **Initial heat potential and domain heat potential**:
\[
J[f](x, \Delta t)=\int_{\Omega} G(x - y, \Delta t) f(y) dy
\]
\[
V[F](x, \Delta t)=\int_{0}^{\Delta t} \int_{\Omega} G(x - y, \Delta t-\tau) F(y, \tau) dy d\tau
\]
- **Discrete Gaussian transform**:
\[
F(x_{i})=\sum_{j = 1}^{N} q_{j}\cdot e^{-\|x_{i}-y_{j}\|^{2} / \delta}
\]
where \(x_{i}\) is the target point, \(y_{j}\) is the source point, and \(q_{j}\) is the source intensity.
Through these improvements, this algorithm can significantly improve computational efficiency while maintaining high precision and is applicable to a wide range of physical and engineering problems.