Wavelet compression of anisotropic integrodifferential operators on sparse tensor product spaces

Abstract

Abstract. For a class of anisotropic integrodifferential operators B arising as semigroup generators ofMarkov processes, we present a sparse tensor product wavelet compression scheme for the Galerkin finite element discretization of the corresponding integrodifferential equations Bu = f on [0,1] n with possibly large n. Under certain conditions on B, the scheme is of essentially optimal and dimension independent complexity O(h −1 | log h| 2(n−1) ) without corrupting the convergence or smoothness requirements of the original sparse tensor finite element scheme. If the conditions on B are not satisfied, the complexity can be bounded by O(h −(1+ε) ), where ε 1 tends to zero with increasing number of the wavelets' vanishing moments. Here h denotes the width of the corresponding finite element mesh. The operators under consideration are assumed to be of non-negative (anisotropic) order and admit a non-standard kernel κ(·, ·) that can be singular on all secondary diagonals. Practical examples of such operators from Mathematical Finance are given and some numerical results are presented.Mathematics Subject Classification. 47A20, 65F50, 65N12, 65Y20, 68Q25, 45K05, 65N30.

Space-time adaptive wavelet methods for parabolic evolution problems

“…These conditions are already implied by (8.7). For certain integrodifferential operators, s * -computability in case (B) has been investigated in [Rei08].…”

Numerical Analysis of Additive, Lévy and Feller Processes with Applications to Option Pricing

“…This can be achieved by defining an appropriate approximation matrix A J corresponding to a bilinear form a (u, v). The analysis of compression schemes for high-dimensional anisotropic Lévy type operators was done by [64,65,66]. Processes with state spaces in R 1 whose generators are Sobolev spaces of variable order were treated by [77].…”

Adaptive wavelet methods for solving operator equations: An overview

“…The key is that the sizes of the entries decay exponentially as function of the sum of the absolute differences in levels of the tensor product wavelets involved. Compressiblity of integrodifferential operators has been investigated in [Rei08].…”