MiS Preprint Repository

Delve into the future of research at MiS with our preprint repository. Our scientists are making groundbreaking discoveries and sharing their latest findings before they are published. Explore repository to stay up-to-date on the newest developments and breakthroughs.

MiS Preprint

Convolution of $hp$-functions on locally refined grids - extended version

Wolfgang Hackbusch


Usually, the fast evaluation of a convolution integral $\int_{\mathbb{R}}f(y)g(x-y)\mathrm{d}y$ requires that the functions $f,g$ have a simple structure based on an equidistant grid in order to apply the fast Fourier transform. Here we discuss the efficient performance of the convolution of $hp$-functions in certain locally refined grids. More precisely, the convolution result is projected into some given $hp$-space (Galerkin approximation). The overall cost is $\mathcal{O}(p^{2}N\log N),$ where $N$ is the sum of the dimensions of the subspaces containing $f$, $g$ and the resulting function, while $p$ is the maximal polynomial degree.

Apr 11, 2007
Apr 11, 2007
MSC Codes:
44A35, 42A55
convolution integral, hp-finite elements, non-uniform grids, local refinement

Related publications

2009 Repository Open Access
Wolfgang Hackbusch

Convolution of hp-functions on locally refined grids

In: IMA journal of numerical analysis, 29 (2009) 4, pp. 960-985