Search

High-dimensional Lagrange interpolation plays a pivotal role in finite element methods, where ensuring the unisolvence and symmetry of its interpolation space and nodes set is crucial. In this paper, we leverage group action and group representation theories to precisely delineate the conditions for unisolvence. We establish a necessary condition for unisolvence: the symmetry of the interpolation nodes set is determined by the given interpolation space. Our findings not only contribute to a deeper theoretical understanding but also promise practical benefits by reducing the computational overhead associated with identifying appropriate interpolation nodes., Comment: 18 pages

In the present paper, we consider one-hidden layer ANNs with a feedforward architecture, also referred to as shallow or two-layer networks, so that the structure is determined by the number and types of neurons. The determination of the parameters that define the function, called training, is done via the resolution of the approximation problem, so by imposing the interpolation through a set of specific nodes. We present the case where the parameters are trained using a procedure that is referred to as Extreme Learning Machine (ELM) that leads to a linear interpolation problem. In such hypotheses, the existence of an ANN interpolating function is guaranteed. The focus is then on the accuracy of the interpolation outside of the given sampling interpolation nodes when they are the equispaced, the Chebychev, and the randomly selected ones. The study is motivated by the well-known bell-shaped Runge example, which makes it clear that the construction of a global interpolating polynomial is accurate only if trained on suitably chosen nodes, ad example the Chebychev ones. In order to evaluate the behavior when growing the number of interpolation nodes, we raise the number of neurons in our network and compare it with the interpolating polynomial. We test using Runge's function and other well-known examples with different regularities. As expected, the accuracy of the approximation with a global polynomial increases only if the Chebychev nodes are considered. Instead, the error for the ANN interpolating function always decays and in most cases we observe that the convergence follows what is observed in the polynomial case on Chebychev nodes, despite the set of nodes used for training.

The Sinc approximation applied to double-exponentially decaying functions is referred to as the DE-Sinc approximation. Because of its high efficiency, this method has been used in various applications. In the Sinc approximation, the mesh size and truncation numbers should be optimally selected to achieve its best performance. However, the standard selection formula has only been "near-optimally" selected because the optimal formula of the mesh size cannot be expressed in terms of elementary functions of truncation numbers. In this study, we propose two improved selection formulas. The first one is based on the concept by an earlier research that resulted in a better selection formula for the double-exponential formula. The formula performs slightly better than the standard one, but is still not optimal. As a second selection formula, we introduce a new parameter to propose truly optimal selection formula. We provide explicit error bounds for both selection formulas. Numerical comparisons show that the first formula gives a better error bound than the standard formula, and the second formula gives a much better error bound than the standard and first formulas., Comment: Keywords: Sinc approximation, double-exponential transformation, error bound, mesh size, truncation number

Many applications of geometry modeling and computer graphics necessite accurate curvature estimations of curves on the plane or on manifolds. In this paper, we define the notion of the discrete geodesic curvature of a geodesic polygon on a smooth surface. We show that, when a geodesic polygon P is closely inscribed on a $C^2$-regular curve, the discrete geodesic curvature of P estimates the geodesic curvature of C. This result allows us to evaluate the geodesic curvature of discrete curves on surfaces. In particular, we apply such result to planar and spherical 4-point angle-based subdivision schemes. We show that such schemes cannot generate in general $G^2$-continuous curves. We also give a novel example of $G^2$-continuous subdivision scheme on the unit sphere using only points and discrete geodesic curvature called curvature-based 6-point spherical scheme.

An important component of a number of computational modeling algorithms is an interpolation method that preserves the positivity of the function being interpolated. This report describes the numerical testing of a new positivity-preserving algorithm that is designed to be used when interpolating from a solution defined on one grid to different spatial grid. The motivating application for this work was a numerical weather prediction (NWP) code that uses a spectral element mesh discretization for its dynamics core and a cartesian tensor product mesh for the evaluation of its physics routines. This coupling of spectral element mesh, which uses nonuniformly spaced quadrature/collocation points, and uniformly-spaced cartesian mesh combined with the desire to maintain positivity when moving between these meshes necessitates our work. This new approach is evaluated against several typical algorithms in use on a range of test problems in one or more space dimensions. The results obtained show that the new method is competitive in terms of observed accuracy while at the same time preserving the underlying positivity of the functions being interpolated., Comment: 57 pages, 15 figures

A collection of algorithms in object-oriented MATLAB is described for numerically computing with smooth functions defined on the unit ball in the Chebfun software. Functions are numerically and adaptively resolved to essentially machine precision by using a three-dimensional analogue of the double Fourier sphere method to form "ballfun" objects. Operations such as function evaluation, differentiation, integration, fast rotation by an Euler angle, and a Helmholtz solver are designed. Our algorithms are particularly efficient for vector calculus operations, and we describe how to compute the poloidal-toroidal and Helmholtz--Hodge decomposition of a vector field defined on the ball., Comment: 23 pages, 9 figures

In this paper, we propose an interpolation formula for periodic functions. This formula can be regarded as an analog of the Sinc approximation, which is an interpolation formula for functions defined on the entire infinite interval. Theoretical error analysis and numerical examples show that the proposed formula converges exponentially for analytic periodic functions., Comment: 9 pages, 1 figure I have decided to withdraw the paper because the contents of the paper were already proposed by other researchers

We propose a novel approach to allocating resources for expensive simulations of high fidelity models when used in a multifidelity framework. Allocation decisions that distribute computational resources across several simulation models become extremely important in situations where only a small number of expensive high fidelity simulations can be run. We identify this allocation decision as a problem in optimal subset selection, and subsequently regularize this problem so that solutions can be computed. Our regularized formulation yields a type of group lasso problem that has been studied in the literature to accomplish subset selection. Our numerical results compare performance of algorithms that solve the group lasso problem for algorithmic allocation against a variety of other strategies, including those based on classical linear algebraic pivoting routines and those derived from more modern machine learning-based methods. We demonstrate on well known synthetic problems and more difficult real-world simulations that this group lasso solution to the relaxed optimal subset selection problem performs better than the alternatives., Comment: 27 pages, 10 figures

We study the problem of reconstructing a function on a manifold satisfying some mild conditions, given data on the values and some derivatives of the function at arbitrary points on the manifold. While the problem of finding a polynomial of two variables with total degree $\le n$ given the values of the polynomial and some of its derivatives at exactly the same number of points as the dimension of the polynomial space is sometimes impossible, we show that such a problem always has a solution in a very general situation if the degree of the polynomials is sufficiently large. We give estimates on how large the degree should be, and give explicit constructions for such a polynomial even in a far more general case. As the number of sampling points at which the data is available increases, our polynomials converge to the target function on the set where the sampling points are dense. Numerical examples in single and double precision show that this method is stable and of high-order.

In this paper, I generalize a previous one about hilbertian kernels and approximation theory

I want to prove that all classical techniques of interpolation and approximation as Lagrange, Taylor, Hermite interpolations Beziers interpolants, Quasi interpolants, Box splines and others (radial splines, simplicial splines) are derived from a \textbf{unique} simple hilbertian scheme. For sake of simplicity, we shall consider only elementary examples which could be easily generalized.

Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent, however, for the data sets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large data sets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs., Comment: Readers might also like a follow up work to this paper, recently published in Engineering Analysis and Boundary Elements. arXiv:1606.03258, Computational Geoscience, 2018

We derive analyticity criteria for explicit error bounds and an exponential rate of convergence of the magic point empirical interpolation method introduced by Barrault et al. (2004). Furthermore, we investigate its application to parametric integration. We find that the method is well-suited to Fourier transforms and has a wide range of applications in such diverse fields as probability and statistics, signal and image processing, physics, chemistry and mathematical finance. To illustrate the method, we apply it to the evaluation of probability densities by parametric Fourier inversion. Our numerical experiments display convergence of exponential order, even in cases where the theoretical results do not apply.

A pattern of interpolation nodes on the disk is studied, for which the interpolation problem is theoretically unisolvent, and which renders a minimal numerical condition for the collocation matrix when the standard basis of Zernike polynomials is used. It is shown that these nodes have an excellent performance also from several alternative points of view, providing a numerically stable surface reconstruction, starting from both the elevation and the slope data. Sampling at these nodes allows for a more precise recovery of the coefficients in the Zernike expansion of a wavefront or of an optical surface., Comment: 21 pages, 7 figures. Submitted to Appl. Math. Comput

In this work we develop a dynamically adaptive sparse grids (SG) method for quasi-optimal interpolation of multidimensional analytic functions defined over a product of one dimensional bounded domains. The goal of such approach is to construct an interpolant in space that corresponds to the "best $M$-terms" based on sharp a priori estimate of polynomial coefficients. In the past, SG methods have been successful in achieving this, with a traditional construction that relies on the solution to a Knapsack problem: only the most profitable hierarchical surpluses are added to the SG. However, this approach requires additional sharp estimates related to the size of the analytic region and the norm of the interpolation operator, i.e., the Lebesgue constant. Instead, we present an iterative SG procedure that adaptively refines an estimate of the region and accounts for the effects of the Lebesgue constant. Our approach does not require any a priori knowledge of the analyticity or operator norm, is easily generalized to both affine and non-affine analytic functions, and can be applied to sparse grids build from one dimensional rules with arbitrary growth of the number of nodes. In several numerical examples, we utilize our dynamically adaptive SG to interpolate quantities of interest related to the solutions of parametrized elliptic and hyperbolic PDEs, and compare the performance of our quasi-optimal interpolant to several alternative SG schemes.

Positive polynomial operator that approximates Urison operator, when integration domain is a "regular triangle" is investigated. We obtain Bernstein Polynomials as a particular case.

Radial basis functions provide highly useful and flexible interpolants to multivariate functions. Further, they are beginning to be used in the numerical solution of partial differential equations. Unfortunately, their construction requires the solution of a dense linear system. Therefore much attention has been given to iterative methods. In this paper, we present a highly efficient preconditioner for the conjugate gradient solution of the interpolation equations generated by gridded data. Thus our method applies to the corresponding Toeplitz matrices. The number of iterations required to achieve a given tolerance is independent of the number of variables.

In this paper, we consider methods to compute the coefficients of interpolants relative to a basis of polynomials satisfying a three-term recurrence relation. Two new algorithms are presented: the first constructs the coefficients of the interpolation incrementally and can be used to update the coefficients whenever a nodes is added to or removed from the interpolation. The second algorithm, which constructs the interpolation coefficients by decomposing the Vandermonde-like matrix iteratively, can not be used to update or downdate an interpolation, yet is more numerically stable than the first algorithm and is more efficient when the coefficients of multiple interpolations are to be computed over the same set of nodes., Comment: 18 pages, submitted to the Journal of Scientific Computing.

For the last almost three decades, since the famous Buchberger-M\"oller(BM) algorithm emerged, there has been wide interest in vanishing ideals of points and associated interpolation polynomials. Our paradigm is based on the theory of bivariate polynomial interpolation on cartesian point sets that gives us related degree reducing interpolation monomial and Newton bases directly. Since the bases are involved in the computation process as well as contained in the final output of BM algorithm, our paradigm obviously simplifies the computation and accelerates the BM process. The experiments show that the paradigm is best suited for the computation over finite prime fields that have many applications., Comment: 24 pages, 7 figures, submitted to JCAM

The development of high-degree interpolation polynomials which use the values of the function and its subsequent derivatives is reformulated. Also, we present a variant of new formula in barycentric form., Comment: This paper has been withdrawn by the author due to a crucial sign error in main formula

This paper is devoted to the convergence and stability analysis of a class of nonlinear subdivision schemes and associated multi-resolution transforms. These schemes are defined as a perturbation of a linear subdivision scheme. Assuming a contractivity property, stability and convergence are derived. These results are then applied to various schemes such as uncentered interpolatory linear scheme, WENO scheme [13], Power-P scheme [16] and a non linear scheme using local spherical coordinates [18]., Comment: 25 pages, 4 figures

Through introducing a new iterative formula for divided differnce using Neville's and Aitken's algorithms,we study new iterative methods for interpolation,numerical differentiation and numerical integration formulas with arbitrary order of accuracy for evanly or unevanly spaced data. Basic computer algorithms for new methods are given, Comment: 11 pages, comments are invited

We consider the problem of approximating a multiple-input multiple-output (MIMO) $p\times m$ rational transfer function $H(s)$ of high degree by another $p\times m$ rational transfer function $\hat H(s)$ of much smaller degree, so that the ${\cal H}_2$ norm of the approximation error is minimized. We characterize the stationary points of the ${\cal H}_2$ norm of the approximation error by tangential interpolation conditions and also extend these results to the discrete-time case. We analyze whether it is reasonable to assume that lower-order models can always be approximated arbitrarily closely by imposing only first-order interpolation conditions. Finally, we analyze the ${\cal H}_2$ norm of the approximation error for a simple case in order to illustrate the complexity of the minimization problem., Comment: Paper ID sheet: http://www.inma.ucl.ac.be/~absil/Publi/H2-modred-mimo.htm

We consider regular polynomial interpolation algorithms on recursively defined sets of interpolation points which approximate global solutions of arbitrary well-posed systems of linear partial differential equations. Convergence of the 'limit' of the recursively constructed family of polynomials to the solution and error estimates are obtained from a priori estimates for some standard classes of linear partial differential equations, i.e. elliptic and hyperbolic equations. Another variation of the algorithm allows to construct polynomial interpolations which preserve systems of linear partial differential equations at the interpolation points. We show how this can be applied in order to compute higher order terms of WKB-approximations of fundamental solutions of a large class of linear parabolic equations. The error estimates are sensitive to the regularity of the solution. Our method is compatible with recent developments for solution of higher dimensional partial differential equations, i.e. (adaptive) sparse grids, and weighted Monte-Carlo, and has obvious applications to mathematical finance and physics., Comment: 28 pages

A fast and accurate algorithm for solving a Bernstein-Vandermonde linear system is presented. The algorithm is derived by using results related to the bidiagonal decomposition of the inverse of a totally positive matrix by means of Neville elimination. The use of explicit expressions for the determinants involved in the process serves to make the algorithm both fast and accurate., Comment: 13 pages. We have extended the numerical experiments

The Euclidean Algorithm is the often forgotten key to rational approximation techniques, including Taylor, Lagrange, Hermite, osculating, cubic spline, Chebyshev, Pade and other interpolation schemes. A unified view of these various interpolation techniques is eloquently expressed in terms of the concept of the spectral basis of a factor ring of polynomials. When these methods are applied to the minimal polynomial of a matrix, they give a family of rational forms of functions of that matrix., Comment: 8 pages, 2 figures

Formulas for stable differentiation of piecewise-smooth functions are given. The data are noisy values of these functions. The locations of discontinuity points and the sizes of the jumps across these points are not assumed known, but found stably from the noisy data.

The scheme of divided differences is widely used in many approximation and interpolation problems. Computing the Newton coefficients of the interpolating polynomial is the first step of the Bj\"{o}rck and Pereyra algorithm for solving Vandermonde systems of equations (Cf. \cite{bjorck: 70}). Very often this algorithm produces very accurate solution. The problem of determining the Newton coefficients is intimately related with the problem of evaluation the Lagrange interpolating polynomial, which can be realized by many algorithms. For these reasons we use the uniform approach and analyze also Aitken's algorithm of the evaluation of an interpolating polynomial. We propose new algorithms that are always numerically stable with respect to perturbation in the function values and more accurate than the Aitken's algorithm and the scheme of divided differences, even for complex data., Comment: 19 pages, 4 figures

A new global basis of B-splines is defined in the space of generalized quadratic splines (GQS) generated by Merrien subdivision algorithm. Then, refinement equations for these B-splines and the associated corner-cutting algorithm are given. Afterwards, several applications are presented. First a global construction of monotonic and/or convex generalized splines interpolating monotonic and/or convex data. Second, convergence of sequences of control polygons to the graph of a GQS. Finally, a Lagrange interpolant and a quasi-interpolant which are exact on the space of affine polynomials and whose infinite norms are uniformly bounded independently of the partition., Comment: 2004-13

We prove the existence of Absolutely Minimizing Lipschitz Extensions by a method which differs from those used by G. Aronsson in general metrically convex compact metric spaces and R. Jensen in Euclidean spaces. Assuming Jensen's hypotheses, our method yields numerical schemes for computing, in euclidean $\mathbb R$, the solution of viscosity of equation $\Delta_\infty=0$ with Dirichlet's condition., Comment: 2004-02

We use simple sub-Riemannian techniques to prove that an arbitrary geometric p-rough path in the sense of Lyons (98) is the limit in sup-norm of a sequence of canonically lifted smooth paths, which are uniformly bounded in p-variation, clarifying the two different definitions of a geometric p-rough path, Lyons (98), Lyons/Qian (02). Our proofs are based on fine estimates in terms of control functions and are sufficiently general to include the case of Hoelder- and modulus-type regularity, Friz/Victoir (03). This allows us to extend a few classical results on Hoelder-spaces (Ciesielski, Musielak/Semadeni) and p-variation spaces (Wiener,Dudley) to the non-commutative setting necessary for the theory of rough paths.

The goal of this article is that of understanding how the oscillation and concentration effects developed by a sequence of functions in $\mathbb{R}^{d} $ are modified by the action of Sampling and Reconstruction operators on regular grids. Our analysis is performed in terms of Wigner and defect measures, which provide a quantitative description of the high frequency behavior of bounded sequences in $L^{2}(mathbb{R}^{d}) $. We actually present explicit formulas that make possible to compute such measures for sampled/reconstructed sequences. As a consequence, we are able to characterize sampling and reconstruction operators that preserve or filter the high-frequency behavior of specific classes of sequences. The proofs of our results rely on the construction and manipulation of Wigner measures associated to sequences of discrete functions.

This note proposes a method for pricing high-dimensional American options based on modern methods of multidimensional interpolation. The method allows using sparse grids and thus mitigates the curse of dimensionality. A framework of the pricing algorithm and the corresponding interpolation methods are discussed, and a theorem is demonstrated that suggests that the pricing method is less vulnerable to the curse of dimensionality. The method is illustrated by an application to rainbow options and compared to Least Squares Monte Carlo and other benchmarks., Comment: 12 pages, tables omitted