Abstract
Optimality (or KKT) systems arise as primal-dual stationarity conditions for constrained optimization problems. Under suitable constraint qualifications, local minimizers satisfy KKT equations but, unfortunately, many other stationary points (including, perhaps, maximizers) may solve these nonlinear systems too. For this reason, nonlinear-programming solvers make strong use of the minimization structure and the naive use of nonlinear-system solvers in optimization may lead to spurious solutions. Nevertheless, in the basin of attraction of a minimizer, nonlinear-system solvers may be quite efficient. In this paper quasi-Newton methods for solving nonlinear systems are used as accelerators of nonlinear-programming (augmented Lagrangian) algorithms, with equality constraints. A periodically-restarted memoryless symmetric rank-one (SR1) correction method is introduced for that purpose. Convergence results are given and numerical experiments that confirm that the acceleration is effective are presented.
Similar content being viewed by others
References
Andreani, R., Birgin, E.G., Martínez, J.M., Schuverdt, M.L.: Augmented Lagrangian methods with general lower-level constraints. SIAM J. Optim. (2007, to appear)
Andreani, R., Birgin, E.G., Martínez, J.M., Schuverdt, M.L.: On Augmented Lagrangian methods with general lower-level constraints. Technical Report MCDO-050303, Department of Applied Mathematics, UNICAMP, Brazil (2005), available in Optimization On Line and in www.ime.usp.br/~egbirgin
Andreani, R., Martínez, J.M., Schuverdt, M.L.: On the relation between the constant positive linear dependence condition and quasinormality constraint qualification. J. Optim. Theory Appl. 125, 473–485 (2005)
Bertsekas, D.P.: Nonlinear Programming, 2nd edn. Athena Scientific, Belmont (1999)
Birgin, E.G., Castillo, R., Martínez, J.M.: Numerical comparison of augmented Lagrangian algorithms for nonconvex problems. Comput. Optim. Appl. 31, 31–56 (2005)
Birgin, E.G., Krejić, N., Martínez, J.M.: Globally convergent inexact quasi-Newton methods for solving nonlinear systems. Numer. Algorithms 32, 249–260 (2003)
Birgin, E.G., Martínez, J.M.: Large-scale active-set box-constrained optimization method with spectral projected gradients. Comput. Optim. Appl. 23, 101–125 (2002)
Birgin, E.G., Martínez, J.M., Raydan, M.: Nonmonotone spectral projected gradient methods on convex sets. SIAM J. Optim. 10, 1196–1211 (2000)
Birgin, E.G., Martínez, J.M., Raydan, M.: Algorithm 813: SPG-Software for convex-constrained optimization. ACM Trans. Math. Softw. 27, 340–349 (2001)
Birgin, E.G., Martínez, J.M., Raydan, M.: Inexact spectral projected gradient methods on convex sets. IMA J. Numer. Anal. 23, 539–559 (2003)
Conn, A.R., Gould, N.I.M., Toint, P.L.: A globally convergent augmented Lagrangian algorithm for optimization with general constraints and simple bounds. SIAM J. Numer. Anal. 28, 545–572 (1991)
Conn, A.R., Gould, N.I.M., Toint, P.L.: Convergence of quasi-Newton matrices generated by the symmetric rank one update. Math. Program. A 50, 177–195 (1991)
Conn, A.R., Gouldm, N.I.M., Toint, P.L.: LANCELOT: A Fortran Package for Large Scale Nonlinear Optimization. Springer, Berlin (1992)
Curtis, A.R., Powell, J.D., Reid, J.K.: On the estimation of sparse Jacobian matrices. J. Inst. Math. Appl. 13, 117–119 (1974)
Dennis, J.E. Jr., Schnabel, R.B.: Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Classics in Applied Mathematics, vol. 16. SIAM, Philadelphia (1996)
Deuflhard, P.: Newton Techniques for Highly Nonlinear Problems—Theory and Algorithms. Academic Press, San Diego (1996)
Deuflhard, P., Freund, R., Walter, A.: Fast secant methods for the iterative solution of large nonsymmetric linear systems. Impact Comput. Sci. Eng. 2, 244–276 (1990)
Diniz-Ehrhardt, M.A., Gomes-Ruggiero, M.A., Lopes, V.L.R., Martínez, J.M.: Discrete Newton’s method with local variations for solving large-scale nonlinear systems. Optimization 52, 417–440 (2003)
Dostál, Z.: Inexact semimonotonic augmented Lagrangians with optimal feasibility convergence for convex bound and equality constrained quadratic programming. SIAM J. Numer. Anal. 43, 96–115 (2005)
Dostál, Z., Friedlander, A., Santos, S.A.: Augmented Lagrangian with adaptive precision control for quadratic programming with simple bounds and equality constraints. SIAM J. Optim. 13, 1120–1140 (2003)
Dostál, Z., Gomes, F.A.M., Santos, S.A.: Duality based domain decomposition with natural coarse space for variational inequalities. J. Comput. Appl. Math. 126, 397–415 (2000)
Fletcher, R.: Practical Methods of Optimization. Academic Press, London (1987)
Friedlander, A., Gomes-Ruggiero, M.A., Kozakevich, D.N., Martínez, J.M., Santos, S.A.: Solving nonlinear systems of equations by means of quasi-Newton methods with a nonmonotone strategy. In: Optimization Methods and Software, vol. 8, pp. 25–51 (1997)
Gertz, E.M., Gill, P.E.: A primal-dual trust region algorithm for nonlinear optimization. Math. Program. 100, 49–94 (2004)
Gill, P.E., Murray, W., Saunders, M.A.: SNOPT: An SQP algorithm for large-scale constrained optimization. SIAM Rev. 47, 99–131 (2005)
Goldfarb, D., Toint, P.L.: Optimal estimation of Jacobian and Hessian matrices that arise in finite difference calculations. In: Mathematics of Computation, vol. 43, pp. 69–88 (1984)
Gomes-Ruggiero, M.A., Martínez, J.M., Moretti, A.C.: Comparing algorithms for solving sparse nonlinear systems of equations. SIAM J. Sci. Stat. Comput. 13, 459–483 (1992)
Gould, N.I.M., Orban, D., Toint, P.L.: GALAHAD: a library of thread-safe Fortran 90 packages for large-scale nonlinear optimization. ACM Trans. Math. Softw. 29, 353–372 (2003)
Hager, W.W., Zhang, H.: A new active set algorithm for box constrained optimization. Department of Mathematics, University of Florida, 4 July 2005
Hager, W.W., Zhang, H.: Algorithm 851: CG-DESCENT, a conjugate gradient method with guaranteed descent. ACM Trans. Math. Softw. 32, 113–137 (2006)
Hestenes, M.R.: Multiplier and gradient methods. J. Optim. Theory Appl. 4, 303–320 (1969)
Hock, W., Schittkowski, K.: Test Examples for Nonlinear Programming Codes. Lecture Notes in Economic and Mathematical Systems, vol. 187. Springer, New York (1981)
Kelley, C.T.: Iterative Methods for Linear and Nonlinear Equations. SIAM, Philadelphia (1995)
Kelley, C.T., Sachs, E.W.: Local convergence of the Symmetric Rank-One iteration. Comput. Optim. Appl. 9, 43–63 (1998)
Khalfan, H., Byrd, R.H., Schnabel, R.B.: A theoretical and experimental study of the symmetric rank one update. SIAM J. Optim. 3, 1–24 (1993)
Kozakevich, D.N., Martínez, J.M., Santos, S.A.: Solving nonlinear systems of equations with simple constraints. Comput. Appl. Math. 16, 215–235 (1997)
Krejić, N., Martínez, J.M., Mello, M.P., Pilotta, E.A.: Validation of an augmented Lagrangian algorithm with a Gauss-Newton Hessian approximation using a set of hard-spheres problems. Comput. Optim. Appl. 16, 247–263 (2000)
Li, D.-H., Fukushima, M.: A derivative-free line search and global convergence of Broyden-like method for nonlinear equations. Optim. Methods Softw. 13, 181–201 (2000)
Lukšan, L., Vlček, J.: Sparse and partially separable test problems for unconstrained and equality constrained optimization. Technical Report 767, Institute of Computer Science, Academy of Sciences of the Czech Republic (1999)
Martínez, J.M.: Practical quasi-Newton methods for solving nonlinear systems. J. Comput. Appl. Math. 124, 97–122 (2000)
Martínez, J.M., Pilotta, E.A.: Inexact restoration methods for nonlinear programming: advances and perspectives. In: Optimization and Control with Applications, pp. 271–292. Springer, New York (2005)
Martínez, J.M., Qi, L.: Inexact Newton methods for solving nonsmooth equations. J. Comput. Appl. Math. 60, 127–145 (1995)
Nocedal, J., Wright, S.J.: Numerical Optimization. Springer, New York (1999)
Powell, M.J.D.: A method for nonlinear constraints in minimization problems. In: Fletcher, R. (ed.) Optimization, pp. 283–298. Academic Press, New York (1969)
Qi, L., Sun, J.: A nonsmooth version of Newton method. Math. Program. 58, 353–367 (1993)
Qi, L., Wei, Z.: On the constant positive linear dependence condition and its application to SQP methods. SIAM J. Optim. 10, 963–981 (2000)
Rockafellar, R.T.: Augmented Lagrange multiplier functions and duality in nonconvex programming. SIAM J. Control 12, 268–285 (1974)
Shamanski, V.E.: A modification of Newton’s method. Ukrain Mat. Z. 19, 133–138 (1967)
Author information
Authors and Affiliations
Corresponding author
Additional information
This work was supported by FAPESP, CNPq, PRONEX-Optimization (CNPq / FAPERJ), FAEPEX, UNICAMP.
Rights and permissions
About this article
Cite this article
Ferreira-Mendonça, L., Lopes, V.L.R. & Martínez, J.M. Quasi-Newton acceleration for equality-constrained minimization. Comput Optim Appl 40, 373–388 (2008). https://doi.org/10.1007/s10589-007-9090-4
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10589-007-9090-4