L2CXFT
A Fortran subroutine applies the method of {\itI. C. Demetriou} and {\itM. J. D. Powell} [IMA J. Numer. Anal. 11, No. 3, 433--448 (1991; Zbl 726.65009)] to restore convexity in $n$ measurements of a convex function contaminated by random errors. The method minimizes the sum of the squares of the errors, subject to nonnegativity of second divided differences, in two phases. \parFirst, an approximation close to the optimum is derived in $O(n)$ operations. Then, this approximation is used as the starting point of a dual-feasible quadratic programming algorithm that completes the calculation of the optimum. The constraints allow $B$-splines to be used, which reduce the problem to an equivalent one with fewer variables where the knots of the splines are determined automatically from the data points due to the constraint equations. The subroutine benefits from this reduction, since common submatrices that occur during the calculation are updated suitably. Iterative refinement improves the accuracy of some calculations when round-off errors accumulate. \parThe subroutine has been applied to a variety of data having substantial differences and has performed fast and stably even for small data spacing, large $n$, and single-precision arithmetic. Driver programs and examples with output are provided to demonstrate the use of the subroutine.
(Source: http://dl.acm.org/)
This software is also peer reviewed by journal TOMS.
This software is also peer reviewed by journal TOMS.
Keywords for this software
References in zbMATH (referenced in 8 articles , 1 standard article )
Showing results 1 to 8 of 8.
Sorted by year (- Vassiliou, E.E.; Demetriou, I.C.: A linearly distributed lag estimator with $r$-convex coefficients (2010)
- Georgiadou, Sofia A.; Demetriou, Ioannis C.: A computational method for the Karush-Kuhn-Tucker test of convexity of univariate observations and certain economic applications (2008)
- Demetriou, I.C.: L2CXCV: a Fortran 77 package for least squares convex/concave data smoothing (2006)
- Demetriou, Ioannis C.; Lipitakis, Elias A.; Vassiliou, Evangelos E.: Efficient use of Toeplitz matrices for least squares data fitting by nonnegative differences (2005)
- Demetriou, I.C.: Signs of divided differences yield least squares data fitting with constrained monotonicity or convexity (2002)
- Demetriou, I.C.; Lipitakis, E.A.: Certain positive definite submatrices that arise from binomial coefficient matrices (2001)
- Demetriou, I.C.: CXFTV2: A Fortran subroutine for the discrete least squares convex approximation (1997)
- Demetriou, I.C.: Algorithm 742: L2CXFT: A Fortran subroutine for least-squares data fitting with nonnegative second divided differences (1995)