Large-Scale Optimization with Linear Equality Constraints Using Reduced Compact Representation
Abstract
For optimization problems with linear equality constraints, we prove that the (1,1) block of the inverse KKT matrix remains unchanged when projected onto the nullspace of the constraint matrix. We develop reduced compact representations of the limited-memory inverse BFGS Hessian to compute search directions efficiently when the constraint Jacobian is sparse. Orthogonal projections are implemented by a sparse QR factorization or a preconditioned LSQR iteration. In numerical experiments two proposed trust-region algorithms improve in computation times, often significantly, compared to previous implementations of related algorithms and compared to IPOPT.
- Publication:
-
SIAM Journal on Scientific Computing
- Pub Date:
- February 2022
- DOI:
- 10.1137/21M1393819
- arXiv:
- arXiv:2101.11048
- Bibcode:
- 2022SJSC...44A.103B
- Keywords:
-
- Mathematics - Optimization and Control;
- Mathematics - Numerical Analysis;
- Statistics - Computation;
- 68Q25;
- 68R10;
- 68U05
- E-Print:
- SIAM Journl. Sci. Comput. 44(1), 2022