Implementing Generating Set Search Methods for Linearly Constrained Minimization

We discuss an implementation of a derivative-free generating set search method for linearly constrained minimization with no assumption of nondegeneracy placed on the constraints. The convergence guarantees for generating set search methods require that the set of search directions possesses certain geometrical properties that allow it to approximate the feasible region near the current iterate. In the hard case, the calculation of the search directions corresponds to finding the extreme rays of a cone with a degenerate vertex at the origin, a difficult problem. We discuss here how state-of-the-art computational geometry methods make it tractable to solve this problem in connection with generating set search. We also discuss a number of other practical issues of implementation, such as the careful treatment of equality constraints and the desirability of augmenting the set of search directions beyond the theoretically minimal set. We illustrate the behavior of the implementation on several problems from the CUTEr test suite. We have found it to be successful on problems with several hundred variables and linear constraints.

[1]  Robert Michael Lewis,et al.  Pattern Search Algorithms for Bound Constrained Minimization , 1999, SIAM J. Optim..

[2]  Komei Fukuda,et al.  Double Description Method Revisited , 1995, Combinatorics and Computer Science.

[3]  Marco Sciandrone,et al.  A Derivative-Free Algorithm for Bound Constrained Optimization , 2002, Comput. Optim. Appl..

[4]  D. Avis lrs : A Revised Implementation of the Rev rse Search Vertex Enumeration Algorithm , 1998 .

[5]  Jorge J. Moré,et al.  Testing Unconstrained Optimization Software , 1981, TOMS.

[6]  Nicholas I. M. Gould,et al.  Trust Region Methods , 2000, MOS-SIAM Series on Optimization.

[7]  D. Avis A Revised Implementation of the Reverse Search Vertex Enumeration Algorithm , 2000 .

[8]  Tamara G. Kolda,et al.  Optimization by Direct Search: New Perspectives on Some Classical and Modern Methods , 2003, SIAM Rev..

[9]  H. Raiffa,et al.  3. The Double Description Method , 1953 .

[10]  W. Yu POSITIVE BASIS AND A CLASS OF DIRECT SEARCH TECHNIQUES , 1979 .

[11]  Paul Tseng,et al.  Objective-derivative-free methods for constrained optimization , 2002, Math. Program..

[12]  I M GouldNicholas,et al.  CUTEr and SifDec , 2003 .

[13]  Robert Michael Lewis,et al.  Pattern Search Methods for Linearly Constrained Minimization , 1999, SIAM J. Optim..

[14]  S. Nash,et al.  Linear and Nonlinear Programming , 1987 .

[15]  O. H. Brownlee,et al.  ACTIVITY ANALYSIS OF PRODUCTION AND ALLOCATION , 1952 .

[16]  Tamara G. Kolda,et al.  Stationarity Results for Generating Set Search for Linearly Constrained Optimization , 2006, SIAM J. Optim..

[17]  V. Torczon,et al.  RANK ORDERING AND POSITIVE BASES IN PATTERN SEARCH ALGORITHMS , 1996 .

[18]  Nicholas I. M. Gould,et al.  Convergence Properties of an Augmented Lagrangian Algorithm for Optimization with a Combination of General Equality and Linear Constraints , 1996, SIAM J. Optim..

[19]  David Avis,et al.  A pivoting algorithm for convex hulls and vertex enumeration of arrangements and polyhedra , 1991, SCG '91.

[20]  Philip E. Gill,et al.  Practical optimization , 1981 .

[21]  E. H. Zarantonello Projections on Convex Sets in Hilbert Space and Spectral Theory: Part I. Projections on Convex Sets: Part II. Spectral Theory , 1971 .

[22]  Chandler Davis THEORY OF POSITIVE LINEAR DEPENDENCE. , 1954 .

[23]  David Avis,et al.  A pivoting algorithm for convex hulls and vertex enumeration of arrangements and polyhedra , 1992, Discret. Comput. Geom..

[24]  C. G. Broyden A Class of Methods for Solving Nonlinear Simultaneous Equations , 1965 .

[25]  Nicholas I. M. Gould,et al.  CUTE: constrained and unconstrained testing environment , 1995, TOMS.

[26]  David Avis,et al.  How good are convex hull algorithms? , 1995, SCG '95.

[27]  Nicholas I. M. Gould,et al.  CUTEr and SifDec: A constrained and unconstrained testing environment, revisited , 2003, TOMS.