Bounded-Variable Least-Squares: an Algorithm and Applications

The Fortran subroutine BVLS (bounded variable least-squares) solves linear least-squares problems with upper and lower bounds on the variables, using an active set strategy. The unconstrained least-squares problems for each candidate set of free variables are solved using the QR decomposition. BVLS has a “warm-start” feature permitting some of the variables to be initialized at their upper or lower bounds, which speeds the solution of a sequence of related problems. Such sequences of problems arise, for example, when BVLS is used to find bounds on linear functionals of a model constrained to satisfy, in an approximate lp-norm sense, a set of linear equality constraints in addition to upper and lower bounds. We show how to use BVLS to solve that problem when p = 1, 2, or ∞, and to solve minimum l1 and l∞ fitting problems. FORTRAN 77 code implementing BVLS is available from the statlib gopher at Carnegie Mellon University.

[1]  W. R. Burrus,et al.  MATHEMATICAL PROGRAMMING AND THE NUMERICAL SOLUTION OF LINEAR EQUATIONS. , 1976 .

[2]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[3]  Kathryn A. Whaler,et al.  Numerical methods for establishing solutions to the inverse problem of electromagnetic induction , 1981 .

[4]  P. Stark,et al.  Strict Bounds on Seismic Velocity in the Spherical Earth , 1986 .

[5]  Smooth profiles from τ (p) and X(p) data , 1987 .

[6]  Philip B. Stark,et al.  Velocity bounds from statistical estimates of τ (p) and X (p) , 1987 .

[7]  Philip B. Stark,et al.  Correction to “Velocity Bounds From Statistical Estimates of τ(p) and X(p)” by Philip B. Stark and Robert L. Parker , 1988 .

[8]  Parker,et al.  Test of Newton's inverse-square law in the Greenland ice cap. , 1989, Physical review letters.

[9]  R. L. Parker,et al.  An analysis of geophysical experiments to test Newton's law of gravity , 1989, Nature.

[10]  John A. Hildebrand,et al.  A seafloor and sea surface gravity survey of Axial Volcano , 1990 .

[11]  Parker,et al.  Submarine measurement of the Newtonian gravitational constant. , 1991, Physical review letters.

[12]  Robert L. Parker,et al.  A theory of ideal bodies for seamount magnetism , 1991 .

[13]  P. Stark.  Inference in infinite-dimensional inverse problems: Discretization and duality , 1992 .

[14]  K. Hudnut,et al.  Extremal bounds on earthquake movement from geodetic data - Application to the Landers earthquake , 1994 .

[15]  P. Stark,et al.  Uncertainties for two-dimensional models of solar rotation from helioseismic eigenfrequency splitting , 1995 .

[16]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[17]  William H. Press,et al.  Numerical recipes in C , 2002 .