A new source of structured singular value decomposition problems.
Necessity of computing large sparse Hessian matrices gave birth to many methods for their effective approximation by differences of gradients. We adopt the so-called direct methods for this problem that we faced when developing programs for nonlinear optimization. A new approach used in the frame of symmetric sequential coloring is described. Numerical results illustrate the differences between this method and the popular Powell-Toint method.
This paper is devoted to Hermite interpolation with Chebyshev-Lobatto and Chebyshev-Radau nodal points. The aim of this piece of work is to establish the rate of convergence for some types of smooth functions. Although the rate of convergence is similar to that of Lagrange interpolation, taking into account the asymptotic constants that we obtain, the use of this method is justified and it is very suitable when we dispose of the appropriate information.
A numerical method of fitting a multiparameter function, non-linear in the parameters which are to be estimated, to the experimental data in the norm (i.e., by minimizing the sum of absolute values of errors of the experimental data) has been developed. This method starts with the least squares solution for the function and then minimizes the expression , where is the error of the -th experimental datum, starting with an comparable with the root-mean-square error of the least squares solution...
In current textbooks the use of Chebyshev nodes with Newton interpolation is advocated as the most efficient numerical interpolation method in terms of approximation accuracy and computational effort. However, we show numerically that the approximation quality obtained by Newton interpolation with Fast Leja (FL) points is competitive to the use of Chebyshev nodes, even for extremely high degree interpolation. This is an experimental account of the analytic result that the limit distribution of FL...
There are two grounds the spline theory stems from - the algebraic one (where splines are understood as piecewise smooth functions satisfying some continuity conditions) and the variational one (where splines are obtained via minimization of some quadratic functionals with constraints). We use the general variational approach called smooth interpolation introduced by Talmi and Gilat and show that it covers not only the cubic spline and its 2D and 3D analogues but also the well known tension spline...
We present the error analysis of Lagrange interpolation on triangles. A new a priori error estimate is derived in which the bound is expressed in terms of the diameter and circumradius of a triangle. No geometric conditions on triangles are imposed in order to get this type of error estimates. To derive the new error estimate, we make use of the two key observations. The first is that squeezing a right isosceles triangle perpendicularly does not reduce the approximation property of Lagrange interpolation....
In signal and image processing as well as in numerical solution of differential equations, wavelets with short support and with vanishing moments are important because they have good approximation properties and enable fast algorithms. A B-spline of order is a spline function that has minimal support among all compactly supported refinable functions with respect to a given smoothness. And recently, B. Han and Z. Shen constructed Riesz wavelet bases of with vanishing moments based on B-spline...