A new interior-point approach for large two-stage stochastic problems

Two-stage stochastic models give rise to very large optimization problems. Several approaches havebeen devised for efficiently solving them, including interior-point methods (IPMs). However, usingIPMs, the linking columns associated to first-stage decisions cause excessive fill-in for the solutionof...

ver descrição completa

Detalhes bibliográficos
Autores: Castro Pérez, Jordi|||0000-0003-3573-4568, Lama Zubirán, Paula de la|||0000-0001-5735-9581
Formato: informe técnico
Fecha de publicación:2020
País:España
Recursos: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/184307
Acesso em linha:https://hdl.handle.net/2117/184307
Access Level:acceso abierto
Palavra-chave:Interior-point methods
Stochastic optimization
Structured problems
Large-scale optimization
Classificació AMS::90 Operations research, mathematical programming::90C Mathematical programming
Àrees temàtiques de la UPC::Matemàtiques i estadística::Investigació operativa::Optimització
Descrição
Resumo:Two-stage stochastic models give rise to very large optimization problems. Several approaches havebeen devised for efficiently solving them, including interior-point methods (IPMs). However, usingIPMs, the linking columns associated to first-stage decisions cause excessive fill-in for the solutionof the normal equations. This downside is usually alleviated if variable splitting is applied to first-stage variables. This work presents a specialized IPM that applies variable splitting and exploits thestructure of the deterministic equivalent of the stochastic problem. The specialized IPM combinesCholesky factorizations and preconditioned conjugate gradients for solving the normal equations.This specialized IPM outperforms other approaches when the number of first-stage variables is largeenough. This paper provides computational results for two stochastic problems: (1) a supply chainsystem and (2) capacity expansion in an electric system. Both linear and convex quadratic formu-lations were used, obtaining instances of up to 38 million variables and six million constraints. Thecomputational results show that our procedure is more efficient than alternative state-of-the-art IPMimplementations (e.g., CPLEX) and other specialized solvers for stochastic optimization