We are concerned with the solution of the bound constrained minimization problem {minf(x), la parts per thousand currency signxa parts per thousand currency signu}. For the solution of this problem we propose an active set method that combines ideas from projected and nonmonotone Newton-type methods. It is based on an iteration of the form x (k+1)=[x (k) +alpha (k) d (k) ](a (TM)-), where alpha (k) is the steplength, d (k) is the search direction and [a <...](a (TM)-) is the projection operator on the set [l,u]. At each iteration a new formula to estimate the active set is first employed. Then the components of d (k) corresponding to the free variables are determined by a truncated Newton method, and the components of d (k) corresponding to the active variables are computed by a Barzilai-Borwein gradient method. The steplength alpha (k) is computed by an adaptation of the nonmonotone stabilization technique proposed in Grippo et al. (Numer. Math. 59:779-805, 1991). The method is a feasible one, since it maintains feasibility of the iterates x (k) , and is well suited for large-scale problems, since it uses matrix-vector products only in the truncated Newton method for computing . We prove the convergence of the method, with superlinear rate under usual additional assumptions. An extensive numerical experimentation performed on an algorithmic implementation shows that the algorithm compares favorably with other widely used codes for bound constrained problems.

An active set feasible method for large-scale minimization problems with bound constraints / DE SANTIS, MARIANNA; DI PILLO, Gianni; LUCIDI, Stefano. - In: COMPUTATIONAL OPTIMIZATION AND APPLICATIONS. - ISSN 0926-6003. - STAMPA. - 53:(2012), pp. 395-423. [10.1007/s10589-012-9506-7]

### An active set feasible method for large-scale minimization problems with bound constraints

#### Abstract

We are concerned with the solution of the bound constrained minimization problem {minf(x), la parts per thousand currency signxa parts per thousand currency signu}. For the solution of this problem we propose an active set method that combines ideas from projected and nonmonotone Newton-type methods. It is based on an iteration of the form x (k+1)=[x (k) +alpha (k) d (k) ](a (TM)-), where alpha (k) is the steplength, d (k) is the search direction and [a <...](a (TM)-) is the projection operator on the set [l,u]. At each iteration a new formula to estimate the active set is first employed. Then the components of d (k) corresponding to the free variables are determined by a truncated Newton method, and the components of d (k) corresponding to the active variables are computed by a Barzilai-Borwein gradient method. The steplength alpha (k) is computed by an adaptation of the nonmonotone stabilization technique proposed in Grippo et al. (Numer. Math. 59:779-805, 1991). The method is a feasible one, since it maintains feasibility of the iterates x (k) , and is well suited for large-scale problems, since it uses matrix-vector products only in the truncated Newton method for computing . We prove the convergence of the method, with superlinear rate under usual additional assumptions. An extensive numerical experimentation performed on an algorithmic implementation shows that the algorithm compares favorably with other widely used codes for bound constrained problems.
##### Scheda breve Scheda completa Scheda completa (DC)
2012
53
395
423
DE SANTIS, MARIANNA; DI PILLO, Gianni; LUCIDI, Stefano
File in questo prodotto:
File
VE_2012_11573-485376.pdf

Accesso chiuso

Licenza: Tutti i diritti riservati
Dimensione 914.45 kB
COAP_BOX-NT-2012.pdf

Accesso chiuso

Licenza: Tutti i diritti riservati
Dimensione 914.45 kB
Utilizza questo identificatore per citare o creare un link a questa risorsa: `https://hdl.handle.net/2158/1350102`