Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
36/2008

Fast and Accurate Tensor Approximation of Multivariate Convolution with Linear Scaling in Dimension

Boris N. Khoromskij

Abstract

In the present paper we present the tensor-product approximation of multi-dimensional convolution transform discretized via collocation-projection scheme on the uniform or composite refined grids. Examples of convolving kernels are given by the classical Newton, Slater (exponential) and Yukawa potentials, ${1}/{\|x\|} $, $e^{-\lambda \|x\|} $ and ${e^{-\lambda\|x\| }}/{\|x\|} $ with $x \in \mathbb{R}^{d} $.

For piecewise constant elements on the uniform grid of size $n^d$, we prove the quadratic convergence $O(h^2) $ in the mesh parameter $h=1/n$, and then justify the Richardson extrapolation method on a sequence of grids that improves the order of approximation up to $O(h^3) $.

The fast algorithm of complexity $O(d R_1 R_2 n \log n)$ is described for tensor-product convolution on the uniform/composite grids of size $n^d$, where $R_1, R_2$ are tensor ranks of convolving functions. We also present the tensor-product convolution scheme in the two-level Tucker-canonical format and discuss the consequent rank reduction strategy.

Finally, we give numerical illustrations confirming:
(a) the approximation theory for convolution schemes of order $O(h^2) $ and $O(h^3) $;
(b) linear-logarithmic scaling of 1D discrete convolution on composite grids;
(c) linear-logarithmic scaling in $n$ of our tensor-product convolution method on $n \times n \times n$ grid in the range $n\leq 16384 $.

Received:
Apr 22, 2008
Published:
May 6, 2008
MSC Codes:
65F30, 65F50, 65F35
Keywords:
multi-dimensional convolution, Tucker decomposition, composite grids, Richardson extrapolation

Related publications

inJournal
2010 Repository Open Access
Boris N. Khoromskij

Fast and accurate tensor approximation of a multivariate convolution with linear scaling in dimension

In: Journal of computational and applied mathematics, 234 (2010) 11, Special Issue, pp. 3122-3139