{"id":"arxiv/2606.28184","url":"https://arxiv.org/abs/2606.28184","tags":["math.NA","physics.comp-ph"],"score":0,"title":"A fast sum-of-Gaussians algorithm for the high-dimensional fractional Fokker-Planck equation","domain":"math.NA","source":"arxiv","authors":["Shidong Jiang","Dong Wang","Qi Zhou"],"reviews":0,"subject":"math.NA","summary":"We present a fast, high-order algorithm for the free-space fractional Fokker-Planck equation (FFPE) in arbitrary spatial dimension. Its fundamental solution, corresponding to a Dirac-delta initial condition, is obtained from the explicit Fourier representation by applying a sum-of-Gaussians (SOG) approximation to the nonseparable stretched exponential, using its complete monotonicity as the Laplace transform of a one-sided $α$-stable density. Each Gaussian term is an ordinary heat kernel and therefore factorizes across spatial coordinates. On a tensor-product grid, the separated form can be…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-06-26","reactions":{}}