Fast computation of elastic and hydrodynamic potentials using approximate approximations

We propose fast cubature formulas for the elastic and hydrodynamic potentials based on the approximate approximation of the densities with Gaussian and related functions. For densities with separated representation, we derive a tensor product representation of the integral operator which admits efficient cubature procedures. We obtain high order approximations up to a small saturation error, which is negligible in computations. Results of numerical experiments which show approximation order O(h2M)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathscr {O}(h^{2M})$$\end{document}, M=1,2,3,4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M=1,2,3,4$$\end{document}, are provided.


Introduction
In this paper we consider volume potentials which arise in the solution of the threedimensional problems in elasticity and hydrodynamics.
The concept of approximate approximations and first related results were published by Maz'ya [14,15]. Various aspects of a general theory of these approximations were systematically investigated (cf. [17] and the references therein). This approximation procedure employs linear combination of smooth and rapidly decaying functions to approximate a given function with high order within a prescribed accuracy. Although in general the method does not converge in a rigorous sense, even simple quasiinterpolants approximate functions very accurately. Additionally, this approximation process has considerable advantages due to the great flexibility in the choice of approximating functions. This flexibility makes it possible to get simple formulas for the approximation of various integral and pseudo-differential operators. Example of those cubature formulas for important integral operators of mathematical physics on R 3 were studied in [16]. By combining semi-analytic cubature formulas for volume potentials based on approximate approximations with the strategy of separated representations, we derive a method for approximating volume potentials, which is accurate and fast also in the multidimensional case and provides approximation formulas of high order. This approach was introduced for the first time in [9] to obtain fast cubature of high-dimensional harmonic potentials. Fast cubature formulas of advection-diffusion potentials in the frame of approximate approximations have been obtained in [10]. In [11] the procedure has been extended to parabolic problems of second order and in [12] to the biharmonic operator.
We describe the idea for the elastic potential (1.4). The hydrodynamic potential will be considered in Sect. 4. Our approximation formulas are based on replacing the density f of the integral operator (1.4) by quasi-interpolants of the form with fixed positive parameters h and D and with some generating function η sufficiently smooth and of rapid decay chosen such that the integrals can be effectively determined. Then the sum provides an approximation formula for Γ (k, ) ( f ), k, = 1, 2, 3.
If the generating function η satisfies the moment conditions of order N then at any point x . A proper choice of the parameter D allows to make the saturation error ε 0 (η, D) as small as necessary, e.g., less than the precision machine. The cubature error (1.12) Therefore it is very efficient to precompute and store the values Γ (k, ) (η 2M )((j − m)/ √ D), which can be used for any gridsize h. However, due to the operation number proportional to h −3 , the convolution (1.12) is not practical. We propose a method which reduces the computational effort. As basis functions, we choose (1.13) where H k are the Hermite polynomials The functions (1.13) satisfy the moment conditions (1.10) and consequently give rise to approximation formulas of order 2M, plus the saturation error O(h 2 e −Dπ 2 ). Moreover, the action of three-dimensional elastic and hydrodynamic potentials onto these basis functions allows one-dimensional integral representations. For the special case of Gaussian functions, which may be even orthotropic or anisotropic, this was obtained in [17, pp. 131-136]. For example, In this paper we obtain new one-dimensional integral representations for the threedimensional integrals Γ (k, ) and Ψ (k, ) applied to the tensor product functions (1.13). The integrand is separated, i.e., it is a product of functions depending only on one of the variables. We show how a suitable quadrature of the one-dimensional integrals combined with the approximation of a separated representation of the density leads to a separated approximation of the integral operator. Since the separated representation of the density can be approximated with high order and only one-dimensional operations are used, the resulting method is very efficient and can provide high order approximations.
The outline of the paper is as follows. In Sect. 2 we derive one-dimensional integral representations of the three-dimensional elastic potential applied to the functions (1.13). In Sect. 3, for densities with separated representation, we derive a tensor product representation of the integral operator which admits efficient one-dimensional operations. We provide numerical tests, showing that these formulas are efficient and provide approximation of order O(h 2M ), M = 1, 2, 3, 4. In Sect. 4 we apply the procedure to the hydrodynamic potentials (1.6) and (1.8), and provide approximation up to order O(h 8 ).

Elastic potential
where L denotes the harmonic potential and It remains to determine a one-dimensional integral representation for I k (e −|·| 2 ).
Thus, from the formulas (2.1), (2.4) and (2.5) we obtain Let M > 1. We are looking for a one-dimensional integral representation with separated integrand of the potential Γ (k, ) with the basis functions (1.13). Keeping in mind (2.1), we write A one-dimensional integral representation for the harmonic potential applied to η 2M in (1.13) is known ([9, p. 894]) is a polynomial in x whose coefficients depend on t, defined by

9)
and, for k = , and integrating by parts we get Therefore, by using (2.5), we derive that Then it is easy to verify that Q Hence, (2.14) Let Q (1) M (x, t) and Q (2) M (x, t) be polynomials in x such that
Let M ≥ 1. We replace f = ( f 1 , f 2 , f 3 ) in (1.6) with the quasi-interpolant (1.9) and the generating functions (1.13). Then an approximations of Ψ (k, ) where rewriting (1.7) as and combining it with (2.2) and (2.3), we have In view of (2.7), (2.9) and (2.10), we can write Ψ (k, ) (η 2M ) in the form An interesting feature of the cubature formulas based on the approximate approximation of the density is that the gradient of L f is approximated by the gradient of (4.4) Moreover, by the smoothing properties of the harmonic potential and its gradient, the corresponding saturation errors converge with the rate h 2 and h, respectively, as h tends to 0 (cf. [17,).
In Table 7 we report on the absolute errors and approximation rates for the value We conclude this section by showing accuracy and convergence rate of formula (4.5) for the approximation of P in (4.3) (see Table 8). We assume f 1 = e −|x| 2 (3 − 2|x| 2 ), f 2 = f 3 = 0, which provides, keeping in mind that L ( f 1 ) = e −|x| 2 /2, the exact solution P(x) = x 1 e −|x| 2 .
For the calculations we choose the parameter ν = 2, the quadrature rule with α = 5 and β = 6 in (3.1), τ = 0.003 and 250 summands in the quadrature sum. We choose D = 4 to have the saturation error comparable with the double precision rounding errors. The numerical results confirm the h 2M convergence of the cubature formula when M = 1, 2, 3, 4.