The algebraic eigenvalue problem
The algebraic eigenvalue problem
A tool for the analysis of Quasi-Newton methods with application to unconstrained minimization
SIAM Journal on Numerical Analysis
On the limited memory BFGS method for large scale optimization
Mathematical Programming: Series A and B
Iterative solution methods
A semidefinite framework for trust region subproblems with applications to large scale minimization
Mathematical Programming: Series A and B
Computing an Eigenvector with Inverse Iteration
SIAM Review
LAPACK Users' guide (third ed.)
LAPACK Users' guide (third ed.)
Graph Partitioning and Continuous Quadratic Programming
SIAM Journal on Discrete Mathematics
Testing Unconstrained Optimization Software
ACM Transactions on Mathematical Software (TOMS)
Trust-region methods
A New Matrix-Free Algorithm for the Large-Scale Trust-Region Subproblem
SIAM Journal on Optimization
The Barzilai and Borwein Gradient Method for the Large Scale Unconstrained Minimization Problem
SIAM Journal on Optimization
Nonmonotone Spectral Projected Gradient Methods on Convex Sets
SIAM Journal on Optimization
Solving the Trust-Region Subproblem using the Lanczos Method
SIAM Journal on Optimization
Minimizing a Quadratic Over a Sphere
SIAM Journal on Optimization
Implementing a proximal algorithm for some nonlinear multicommodity flow problems
Networks - Special Issue on Multicommodity Flows and Network Design
Algorithm 873: LSTRS: MATLAB software for large-scale trust-region subproblems and regularization
ACM Transactions on Mathematical Software (TOMS)
A proximal subgradient projection algorithm for linearly constrained strictly convex problems
Optimization Methods & Software
Hi-index | 0.00 |
We present a new matrix-free method for the large-scale trust-region subproblem, assuming that the approximate Hessian is updated by the L-BFGS formula with m=1 or 2. We determine via simple formulas the eigenvalues of these matrices and, at each iteration, we construct a positive definite matrix whose inverse can be expressed analytically, without using factorization. Consequently, a direction of negative curvature can be computed immediately by applying the inverse power method. The computation of the trial step is obtained by performing a sequence of inner products and vector summations. Furthermore, it immediately follows that the strong convergence properties of trust region methods are preserved. Numerical results are also presented.