Abstract
In this paper, we consider the interval bounded generalized trust region subproblem (GTRS) which is the problem of minimizing a general quadratic function subject to an upper and lower bounded general quadratic constraint. Under the assumption that two matrices from the objective and the constraint functions can be simultaneously diagonalizable via congruence, a diagonalization-based algorithm is introduced to solve it by showing that GTRS is indeed equivalent to a linearly constrained convex univariate problem. Some numerical experiments are given to show the effectiveness of the proposed method and to compare it with the extended Rendl–Wolkowicz algorithm due to Pong and Wolkowicz.
Similar content being viewed by others
Notes
GTRS_ codes are available from http://www.math.uwaterloo.ca/~ptingkei/GTRS_codes/.
eigifp is a MATLAB program for computing a few algebraically smallest or largest eigenvalues and their corresponding eigenvectors of the generalized eigenvalue problem, available from http://www.ms.uky.edu/~qye/software.html.
References
Ben-Tal A, Teboulle M (1996) Hidden convexity in some nonconvex quadratically constrained quadratic programming. Math Program 72(1):51–63
Burer S, Yang B (2013) The trust region subproblem with non-intersecting linear constraints. Math Program 40:1–12
Beck A, Ben-Tal A, Teboulle M (2006) Finding a global optimal solution for a quadratically constrained fractional quadratic problem with applications to the regularized total least squares. SIAM J Matrix Anal Appl 28(2):425–445
Conn AR, Gould NI, Toint PL (2000) Trust region methods. SIAM, Philadelphia
Celis MR, Dennis JE, Tapia AR (1985) A trust region strategy for nonlinear equality constrained optimization. Numer Optim 1984:71–82
Fletcher R (1987) Practical methods of optimization. Wiley, New York
Fortin C, Wolkowicz H (2004) The trust region subproblem and semidefinite programming. Optim Methods Softw 19(1):41–67
Fu M, Luo ZQ, Ye Y (1998) Approximation algorithms for quadratic programming. J Combin Optim 2(1):29–50
Golub GH, Ye Q (2002) An inverse free preconditioned Krylov subspace method for symmetric generalized eigenvalue problems. SIAM J Sci Comput 24(1):312–334
Gould NI, Lucidi S, Roma M, Toint PL (1999) Solving the trust-region subproblem using the Lanczos method. SIAM J Optim 9(2):504–525
Gould NI, Robinson DP, Thorne HS (2010) On solving trust-region and other regularised subproblems in optimization. Math Program Comput 2(1):21–57
Golub GH, Van Loan CF (1980) An analysis of the total least squares problem. SIAM J Numer Anal 17(6):883–893
Hsia Y, Lin GX, Sheu RL (2014) A revisit to quadratic programming with one inequality quadratic constraint via matrix pencil. Pac J Optim 10(3):461–481
Hansen PC (1994) Regularization tools, a Matlab package for analysis of discrete regularization problems. Numer Algo 6:135
Jin Q, Fang SC, Xing W (2010) On the global optimality of generalized trust region subproblems. Optimization 59(8):1139–1151
Lancaster P, Rodman L (2005) Canonical forms for hermitian matrix pairs under strict equivalence and congruence. SIAM Rev 47(3):407–443
Moré JJ, Sorensen DC (1983) Computing a trust region step. SIAM J Sci Stat Comput 4(3):553–572
Moré JJ (1993) Generalizations of the trust region problem. Optim Methods Softw 2(3–4):189–209
Pong TK, Wolkowicz H (2014) The generalized trust region subproblem. Comput Optim Appl 58(2):273–322
Rendl F, Wolkowicz H (1997) A semidefinite framework for trust region subproblems with applications to large scale minimization. Math Program 77(1):273–299
Stern RJ, Wolkowicz H (1995) Indefinite trust region subproblems and nonsymmetric eigenvalue perturbations. SIAM J Optim 5(2):286–313
Sturm JF, Zhang S (2003) On cones of nonnegative quadratic functions. Math Oper Res 28(2):246–267
Salahi M, Fallahi S, Terlaky T (2016) Minimizing an indefinite quadratic function subject to a single indefinite quadratic constraint, Technical Report, University of Guilan
Sima DM, Van Huffel S, Golub GH (2004) Regularized total least squares based on quadratic eigenvalue problem solvers. BIT Numer Math 44(4):793–812
Van Huffel S, Vandewalle J (1991) The total least squares problem: computational aspects and analysis. Front Appl Math 9:1–87 (SIAM, Philadelphia)
Wang S, Xia Y (2015) Strong duality for generalized trust region subproblem: S-lemma with interval bounds. Optim Lett 9:1063–1073
Ye Y, Zhang S (2003) New results on quadratic minimization. SIAM J Optim 14(1):245–267
Acknowledgments
The authors would like to thank the reviewer for his/her useful comments and Prof. Wolkowicz for his suggestions on the early version of this work. The first author also would like to thank University of Guilan for supporting his sabbatical leave 2015–2016, hosted by Prof. H. Wolkowicz at the Department of Combinatorics and Optimization, University of Waterloo.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Natasa Krejic.
Rights and permissions
About this article
Cite this article
Salahi, M., Taati, A. An efficient algorithm for solving the generalized trust region subproblem. Comp. Appl. Math. 37, 395–413 (2018). https://doi.org/10.1007/s40314-016-0349-1
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40314-016-0349-1