symp

On Hessian- and Jacobian-free SQP methods - a total quasi-Newton scheme with compact storage

symp

In this paper we describe several modifications to reduce the memory requirement of the total quasi-Newton method proposed by Andreas Griewank et al. The idea is based on application of the compact representation formulae for the wellknown BFGS and SR1 update for unconstrained optimization. It is shown how these definitions can be extended to a total quasi-Newton approach for the constrained case. A brief introduction to the limited-memory approach is described in the present paper using an updated null-space factorization for the KKT system as well as an efficient numerical implementation of the null-space method in which the null-space representation is not stored directly. It can be proven that the number of operations per iteration is bounded by a bilinear order O(n⋅max(m,l)) instead of a cubic order O(m⋅n2) for standard SQP methods. Here ndenotes the number of variables, m the maximal number of active constraints, and lthe user-selected number of stored update vectors.

Next Post