Exploiting multilevel Toeplitz structures in high dimensional nonlocal diffusion
Abstract
We present a finite element implementation for the steady-state nonlocal Dirichlet problem with homogeneous volume constraints. Here, the nonlocal diffusion operator is defined as integral operator characterized by a certain kernel function. We assume that the domain is an arbitrary -dimensional hyperrectangle and the kernel is translation invariant. Under these assumptions, we carefully analyze the structure of the stiffness matrix resulting from a continuous Galerkin method with multilinear elements and exploit this structure in order to cope with the curse of dimensionality associated to nonlocal problems. For the purpose of illustration we choose a particular kernel, which is related to space-fractional diffusion and present numerical results in 1d, 2d and for the first time also in 3d.
Cite
@article{arxiv.1708.02526,
title = {Exploiting multilevel Toeplitz structures in high dimensional nonlocal diffusion},
author = {Christian Vollmann and Volker Schulz},
journal= {arXiv preprint arXiv:1708.02526},
year = {2019}
}