Gunther Schmidt

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

In the present paper we study high-order cubature formulas for the computation of advection-diffusion potentials over boxes. By using the basis functions introduced in the theory of approximate approximations, the cubature of a potential is reduced to the quadrature of one dimensional integrals. For densities with separated approximation, we derive a tensor product representation of the integral operator which admits efficient cubature procedures in very high dimensions. Numerical tests show that these formulas are accurate and provide approximation of order $O(h^6)$ up to dimension $10^8$.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

We propose a fast method for high order approximation of potentials of the Helmholtz type operator Delta+kappa^2 over hyper-rectangles in R^n. By using the basis functions introduced in the theory of approximate approximations, the cubature of a potential is reduced to the quadrature of one-dimensional integrals with separable integrands. Then a separated representation of the density, combined with a suitable quadrature rule, leads to a tensor product representation of the integral operator. Numerical tests show that these formulas are accurate and provide approximations of order 6 up to dimension 100 and kappa^2=100.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

In this paper we derive approximate quasi-interpolants when the values of a function $u$ and of some of its derivatives are prescribed at the points of a uniform grid. As a byproduct of these formulas we obtain very simple approximants which provide high order approximations for solutions to elliptic differential equations with constant coefficients.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

The paper is devoted to the efficient computation of high-order cubature formulas for volume potentials obtained within the framework of approximate approximations. We combine this approach with modern methods of structured tensor product approximations. Instead of performing high-dimensional discrete convolutions the cubature of the potentials can be reduced to a certain number of one-dimensional convolutions leading to a considerable reduction of computing resources. We propose one-dimensional integral representions of high-order cubature formulas for n-dimensional harmonic and Yukawa potentials, which allow low rank tensor product approximations.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

In this paper we propose fast solution methods for the Cauchy problem for the multidimensional Schrödinger equation. Our approach is based on the approximation of the data by the basis functions introduced in the theory of approximate approximations. We obtain high-order approximations also in higher dimensions up to a small saturation error, which is negligible in computations, and we prove error estimates in mixed Lebesgue spaces for the inhomogeneous equation. The proposed method is very efficient in high dimensions if the densities allow separated representations. We illustrate the efficiency of the procedure on different examples, up to approximation order 6 and space dimension 200.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

We derive new formulas for the high dimensional biharmonic potential acting on Gaussians or Gaussians times special polynomials. These formulas can be used to construct accurate cubature formulas of an arbitrary high order which are fast and effective also in very high dimensions. Numerical tests show that the formulas are accurate and provide the predicted approximation rate (O(h^8)) up to the dimension 10^7.

Flavia Lanzara, Vladimir Maz'ya, Gunther Schmidt

A fast method of an arbitrary high order for approximating volume potentials is proposed, which is effective also in high dimensional cases. Basis functions introduced in the theory of approximate approximations are used. Results of numerical experiments, which show approximation order O(h^8) for the Newton potential in high dimensions, for example, for n= 200 000, are provided. The computation time scales linearly in the space dimension. New one-dimensional integral representations with separable integrands of the potentials of advection-diffusion and heat equations are obtained.