The modified absolute-value factorization norm for trust-region minimization

A trust-region method for unconstrained minimization, using a trustregion norm based upon a modified absolute-value factorization of the model Hessian, is proposed. It is shown that the resulting trust-region subproblem may be solved using a single factorization. In the convex case, the method reduces to a backtracking Newton linesearch procedure. The resulting software package is available as HSL_VF06 within the Harwell Subroutine Library. Numerical evidence shows that the approach is effective in the nonconvex case.

[1]  W. Burnside,et al.  Theory of equations , 1886 .

[2]  A note on the solution of quartic equations , 1960 .

[3]  J. Greenstadt On the relative efficiencies of gradient methods , 1967 .

[4]  J. Bunch,et al.  Direct Methods for Solving Symmetric Indefinite Systems of Linear Equations , 1971 .

[5]  M. D. Hebden,et al.  An algorithm for minimization using exact second derivatives , 1973 .

[6]  Philip E. Gill,et al.  Newton-type methods for unconstrained and linearly constrained optimization , 1974, Math. Program..

[7]  J. Bunch,et al.  Some stable methods for calculating inertia and solving symmetric linear systems , 1977 .

[8]  R. Fletcher Factorizing symmetric indefinite matrices , 1976 .

[9]  I. Duff,et al.  Direct Solution of Sets of Linear Equations whose Matrix is Sparse, Symmetric and Indefinite , 1979 .

[10]  Donald Goldfarb The Use of Negative Curvature in Minimization Algorithms , 1980 .

[11]  David M. author-Gay Computing Optimal Locally Constrained Steps , 1981 .

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

[13]  Iain S. Duff,et al.  MA27 -- A set of Fortran subroutines for solving sparse symmetric sets of linear equations , 1982 .

[14]  D. Sorensen Newton's method with a model trust region modification , 1982 .

[15]  Jorge J. Moré,et al.  Computing a Trust Region Step , 1983 .

[16]  S. Vajda Nonlinear Optimization 1981 , 1983 .

[17]  John K. Reid,et al.  The Multifrontal Solution of Indefinite Sparse Symmetric Linear , 1983, TOMS.

[18]  Elizabeth Eskow,et al.  A New Modified Cholesky Factorization , 1990, SIAM J. Sci. Comput..

[19]  P. Toint,et al.  Lancelot: A FORTRAN Package for Large-Scale Nonlinear Optimization (Release A) , 1992 .

[20]  Michael A. Saunders,et al.  Preconditioners for Indefinite Systems Arising in Optimization , 1992, SIAM J. Matrix Anal. Appl..

[21]  M. SIAMJ. STABILITY OF THE DIAGONAL PIVOTING METHOD WITH PARTIAL PIVOTING , 1995 .

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

[23]  John K. Reid,et al.  Exploiting zeros on the diagonal in the direct solution of indefinite sparse symmetric linear systems , 1996, TOMS.

[24]  Annick Sartenaer,et al.  Automatic Determination of an Initial Trust Region in Nonlinear Programming , 1997, SIAM J. Sci. Comput..

[25]  S. H. Cheng,et al.  A Modified Cholesky Algorithm Based on a Symmetric Indefinite Factorization , 1998, SIAM J. Matrix Anal. Appl..

[26]  John G. Lewis,et al.  Accurate Symmetric Indefinite Linear Equation Solvers , 1999, SIAM J. Matrix Anal. Appl..

[27]  Nicholas I. M. Gould,et al.  Solving the Trust-Region Subproblem using the Lanczos Method , 1999, SIAM J. Optim..

[28]  Chih-Jen Lin,et al.  Incomplete Cholesky Factorizations with Limited Memory , 1999, SIAM J. Sci. Comput..