N. Dyn, O. Kounchev, D. Levin et al.
We investigate non-stationary orthogonal wavelets based on a non-stationary interpolatory subdivision scheme reproducing a given set of exponentials. The construction is analogous to the construction of Daubechies wavelets using the subdivision scheme of Deslauriers-Dubuc. The main result is the smoothness of these Daubechies type wavelets.
O. Kounchev, H. Render
Unlike the classical polynomial case there has not been invented up to very recently a tool similar to the Bernstein-Bezier representation which would allow us to control the behavior of the exponential polynomials. The exponential analog to the classical Bernstein polynomials has been introduced in a recent authors' paper which appeared in Constructive Approximations, and this analog retains all basic properties of the classical Bernstein polynomials. The main purpose of the present paper is to contribute in this direction, by proving some important properties of the "Bernstein exponential operator" which has been introduced. We also fix our attention upon some special type of exponential polynomials which are particularly important for the further development of theory of representation of Multivariate data.
O. Kounchev, H. Render, G. Simeonov et al.
Interpolation and smoothing using cubic and generalized splines are fundamental tools in data analysis and statistical modeling. Recently, fast computational algorithms were developed for natural $L$-splines of order four, which arise as piecewise solutions to the differential operator $L_ξ^2 = (\frac{d^2}{dt^2} - ξ^2)^2$. In this paper, we extend this mathematical framework to the important case of clamped (or complete) boundary conditions, where the first derivatives at the interval endpoints are prescribed. We explicitly construct the governing linear system for the interpolation problem and mathematically prove that the resulting tridiagonal matrix is strictly row diagonally dominant, thereby guaranteeing its invertibility and the numerical stability of the fast algorithm. The proposed method is implemented in MATLAB. Furthermore, the developed clamped $L$-splines provide a foundation for constructing multivariate clamped polysplines, which serve as a promising alternative to Physics-Informed Neural Networks (PINNs) for solving partial differential equations in Mathematical Physics.
O. Kounchev, H. Render
We introduce a new type of cubature formula for the evaluation of an integral over the disk with respect to a weight function. The method is based on an analysis of the Fourier series of the weight function and a reduction of the bivariate integral into an infinite sum of univariate integrals. Several experimental results show that the accuracy of the method is superior to standard cubature formula on the disk. Error estimates provide the theoretical basis for the good performance of the new algorithm.
O. Kounchev, H. Render
This paper is a second part of our study of the Discrete Polyharmonic Cubature Formulas on the disc. It completes our study and provides a satisfactory cubature formula in terms of precision and number of evaluation points (coefficient of efficiency of the formula).