Quadratic regularizations in an interior-point method for primal block-angular problems

One of the most efficient interior-point methods for some classes of primal block-angular problems solves the normal equations by a combination of Cholesky factorizations and preconditioned conjugate gradient for, respectively, the block and linking constraints. Its efficiency depends on the spectra...

Descripción completa

Detalles Bibliográficos
Autores: Castro Pérez, Jordi|||0000-0003-3573-4568, Cuesta Andrea, Jordi
Tipo de recurso: artículo
Fecha de publicación:2011
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/14412
Acceso en línea:https://hdl.handle.net/2117/14412
https://dx.doi.org/10.1007/s10107-010-0341-2
Access Level:acceso abierto
Palabra clave:Programming (Mathematics)
Programació (Matemàtica)
Classificació AMS::90 Operations research, mathematical programming::90C Mathematical programming
Àrees temàtiques de la UPC::Matemàtiques i estadística::Investigació operativa::Optimització
Descripción
Sumario:One of the most efficient interior-point methods for some classes of primal block-angular problems solves the normal equations by a combination of Cholesky factorizations and preconditioned conjugate gradient for, respectively, the block and linking constraints. Its efficiency depends on the spectral radius—in [0,1)— of a certain matrix in the definition of the preconditioner. Spectral radius close to 1 degrade the performance of the approach. The purpose of this work is twofold. First, to show that a separable quadratic regularization term in the objective reduces the spectral radius, significantly improving the overall performance in some classes of instances. Second, to consider a regularization term which decreases with the barrier function, thus with no need for an extra parameter. Computational experience with some primal block-angular problems confirms the efficiency of the regularized approach. In particular, for some difficult problems, the solution time is reduced by a factor of two to ten by the regularization term, outperforming state-of-the-art commercial solvers.