Optimal integration of highly oscillating functions in general form

The development of information technology contributes to the improvement of mathematical models of phenomena and processes in many scientific areas of the technical direction. In particular, modern methods of digital signal and image processing use algorithms with new information operators. Cubature formulas are constructed for the approximate calculation of integrals of highly oscillating functions of many variables for various types of data. The paper deals with the estimation of the error in the numerical integration of highly oscillating functions of a general form on the class of differentiable functions of three variables in the case when information about the functions is given to their traces on the corresponding planes. The results obtained make it possible to research the quality of cubature formulas for the approximate calculation of triple integrals of highly oscillating functions of a general form.

Optimal Numerical Integration

Introduction. In many applied problems, such as statistical data processing, digital filtering, computed tomography, pattern recognition, and many others, there is a need for numerical integration, moreover, with a given (often quite high) accuracy. Classical quadrature formulas cannot always provide the required accuracy, since, as a rule, they do not take into account the oscillation of the integrand. In this regard, the development of methods for constructing optimal in accuracy (and close to them) quadrature formulas for the integration of rapidly oscillating functions is rather important and topical problem of computational mathematics. The purpose of the article is to use the example of constructing optimal in accuracy (and close to them) quadrature formulas for calculating integrals for integrands of various degrees of smoothness and for oscillating factors of different types and constructing a priori estimates of their total error, as well as applying to them of the theory of testing the quality of algorithms-programs to create a theory of optimal numerical integration. Results. The optimal in accuracy (and close to them) quadrature formulas for calculating the Fourier transform, wavelet transforms, and Bessel transform were constructed both in the classical formulation of the problem and for interpolation classes of functions corresponding to the case when the information operator about the integrand is given by a fixed table of its values. The paper considers a passive pure minimax strategy for solving the problem. Within the framework of this strategy, we used the method of “caps” by N. S. Bakhvalov and the method of boundary functions developed at the V.M. Glushkov Institute of Cybernetics of the NAS of Ukraine. Great attention is paid to the quality of the error estimates and the methods to obtain them. The article describes some aspects of the theory of algorithms-programs testing and presents the results of testing the constructed quadrature formulas for calculating integrals of rapidly oscillating functions and estimates of their characteristics. The problem of determining the ranges of admissible values of control parameters of programs for calculating integrals with the required accuracy, as well as their best values for integration with the minimum possible error, is considered for programs calculating a priori estimates of characteristics. Conclusions. The results obtained make it possible to create a theory of optimal integration, which makes it possible to reasonably choose and efficiently use computational resources to find the value of the integral with a given accuracy or with the minimum possible error. Keywords: quadrature formula, optimal algorithm, interpolation class, rapidly oscillating function, quality testing.

Three-nucleon bound state calculations using the three dimensional formalism

The traditional method of carrying out few-nucleon calculations is based on the angular momentum decomposition of operators relevant to the calculation. Expressing operators using a finite-sized partial wave basis enables the calculations to be carried out using a small amount of numerical work. Unfortunately, certain calculations that involve higher energies or long range potentials, require including a large number of partial waves in order to get converged results. This is problematic because such an approach requires a numerical implementation of heavily oscillating functions. Modern computers made it possible to carry out few-nucleon calculations without using angular momentum decomposition and instead to work directly with the three dimensional degrees of freedom of the nucleons. In this paper we briefly describe the, so called 3D approach and present preliminary results related to the 3He bound state obtained within this formalism.