Abstract
In this paper, we introduce a Sequential Partial Linearization (SPL) algorithm for finding a solution of the symmetric Eigenvalue Complementarity Problem (EiCP). The algorithm can also be used for the computation of a stationary point of a standard fractional quadratic program. A first version of the SPL algorithm employs a line search technique and possesses global convergence to a solution of the EiCP under a simple condition related to the minimum eigenvalue of one of the matrices of the problem. Furthermore, it is shown that this condition is verified for a simpler version of the SPL algorithm that does not require a line search technique. The main computational effort of the SPL algorithm is the solution of a strictly convex standard quadratic problem, which is efficiently solved by a finitely convergent block principal pivoting algorithm. Numerical results of the solution of test problems from different sources indicate that the SPL algorithm is in general efficient for the solution of the symmetric EiCP in terms of the number of iterations, accuracy of the solution and total computational effort.
Similar content being viewed by others
References
Aragón Artacho, F.J., Fleming, R.M.T., Vuong, P.T.: Accelerating the DC algorithm for smooth functions. Math. Program. 169, 95–118 (2018)
Adly, S., Rammal, H.: A new method for solving second-order cone eigenvalue complementarity problem. J. Optim. Theory Appl. 165, 563–585 (2015)
Brás, C., Fukushima, M., Júdice, J.J., Rosa, S.: Variational inequality formulation for the asymmetric eigenvalue complementarity problem and its solution by means of a gap function. Pacific J. Optim. 8, 197–215 (2012)
Brás, C., Iusem, A.N., Júdice, J.J.: On the quadratic eigenvalue complementarity problem. J. Global Optim. 66, 153–171 (2016)
Brás, C., Fischer, A., Júdice, J.J., Schönefeld, K., Seifert, S.: A block active set algorithm with spectral choice line search for the symmetric eigenvalue complementarity problem. Appl. Math. Comput. 294, 36–48 (2017)
Chen, Z.M., Qi, L.Q.: A semismooth Newton method for tensor eigenvalue complementarity problem. Comput. Optim. Appl. 165, 583–585 (2015)
DIMACS: Second DIMACS Challenge. Test instances available at. http://dimacs.rutgers.edu/challenges
Fernandes, R., Júdice, J.J., Trevisan, V.: Complementary eigenvalues of graphs. Linear Algebra Appl. 527, 216–231 (2017)
Fukushima, M., Mine, H.: A generalized proximal point algorithm for certain nonconvex minimization problems. Int. J. Syst. Sci. 12, 989–1000 (1981)
Fun, J., Nie, J., Zhou, A.: Tensor eigenvalue complementarity problems. Math. Program. Ser. A 170, 507–539 (2018)
Holubova, G., Nocedal, P.: A note on a relation between the Fucik spectrum and Pareto eigenvalues. J. Math. Anal. Appl. 427, 618–628 (2015)
Iusem, A.N., Júdice, J.J., Sessa, V., Sarabando, P.: Splitting methods for the Eigenvalue Complementarity Problem. Optim. Methods Software 34, 1184–1212 (2019)
Iusem, A.N., Júdice, J.J., Sessa, V., Sherali, H.: The second-order cone quadratic eigenvalue complementarity problem. Pacific J. Optim. 13, 475–500 (2017)
Júdice, J., Fukushima, M., Iusem, A.N., Martinez, M., Sessa, V.: An alternating direction method of multipliers for the eigenvalue complementarity problem. Optim. Methods Software. 1, 1 (2020). https://doi.org/10.1080/10556788.2020.1734804
Júdice, J.J., Pires, F.: A block principal pivoting algorithm for large-scale strictly monotone linear complementarity problems. Comput. Oper. Res. 21, 587–596 (1994)
Júdice, J.J., Sherali, H.D., Ribeiro, I., Rosa, S.: On the asymmetric eigenvalue complementarity problem. Optim. Methods Softw. 24, 549–586 (2009)
Júdice, J.J., Raydan, M., Rosa, S., Santos, S.: On the solution of the symmetric complementarity problem by the spectral projected gradient method. Numer. Algorithms 44, 391–407 (2008)
Le Thi, H.A., Pham Dinh, T.: The DC (difference of convex functions) programming and DCA revisited with DC models of real world nonconvex optimization problems. Annal. Oper. Res. 133, 23–46 (2005)
Le Thi, H.A., Moeini, M., Pham, D.T., Júdice, J.J.: A DC programming approach for solving the symmetric eigenvalue complementarity problem. Comput. Optim. Appl. 51, 1097–1117 (2012)
Mine, H., Fukushima, M.: A minimization method for the sum of a convex function and a continuously differentiable function. J. Optim. Theory Appl. 33, 9–23 (1981)
Murty, K.: Linear Complementary. Linear and Nonlinear Programming. Heldermann Verlag, Berlin (1988)
Nocedal, J., Wright, S.: Numerical Optimization. Springer, New York (2006)
de Oliveira, W., Tcheou, M.: An inertial algorithm for DC programming. Set-Valued Variation. Anal. 27, 895–919 (2019)
Pinto da Costa, A., Martins, J., Figueiredo, I., Júdice, J.J.: The directional instability problem in systems with frictional contact. Comput. Methods Appl. Mech. Eng. 193, 357–384 (2004)
Pinto da Costa, A., Seeger, A.: Cone constrained eigenvalue problems, theory and algorithms. Comput. Optim. Appl. 45, 25–57 (2010)
Pinto da Costa, A., Seeger, A., Simões, F.M.F.: Complementarity eigenvalue problems for nolinear matrix pencils. Appl. Math. Comput. 312, 134–148 (2017)
Queiroz, M., Júdice, J., Humes, C.: The symmetric eigenvalue complementarity problem. Math. Comput. 73, 1849–1863 (2003)
Seeger, A.: Eigenvalue analysis of equilibrium processes defined by linear complementarity conditions. Linear Algebra Appl. 294, 1–14 (1999)
Seeger, A.: Complementarity eigenvalue analysis of connected graphs. Linear Algebra Appl. 543, 205–225 (2018)
Seeger, A., Torki, M.: On eigenvalues induced by a cone constraint. Linear Algebra Appl. 372, 181–206 (2003)
Seeger, A., Torki, M.: Local minima of quadratic forms on convex cones. J. Global Optim. 44, 1–28 (2009)
Seeger, A., Vicente-Perez, J.: On cardinality of Pareto spectra. Electron. Linear Algebra J. 22, 758–766 (2011)
Seeger, A.: Quadratic eigenvalue problems under conic constraints. SIAM J. Matrix Anal. Appl. 32, 700–721 (2011)
Tseng, P.: Approximation accuracy, gradient methods, and error bound for structured convex optimization. Math. Program. 125, 263–295 (2010)
Toh, K.C., Yun, S.: An accelerated proximal gradient algorithm for nuclear norm regularized linear least squares problems. Pacific J. Optim. 6, 615–640 (2010)
Wächter, A., Biegler, L.T.: On the implementation of a primal-dual interior point filter line search algorithm for large-scale nonlinear programming. Math. Program. 106, 25–57 (2006)
Yu, G., Song, Y., Xu, Y.: Spectral projected gradient methods for generalized tensor eigenvalue complementarity problems. Numer. Algorithms 80, 1181–1201 (2019)
Zhang, L., Shen, C., Yang, M., Júdice, J.J.: A Lanczos method for large-scale extreme Lorentz eigenvalue problems. SIAM J. Matrix Anal. Appl. 39, 611–631 (2018)
Zhou, Y.H., Gowda, M.S.: On the finiteness of the cone spectrum of certain linear transformations on Euclidean Jordan algebras. Linear Algebra Appl. 431, 772–782 (2009)
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The research of J. Júdice is funded by FCT/MCTES through national funds and when applicable co-funded EU funds under the project UIDB/EEA/50008/2020.
W. de Oliveira acknowledges financial support from the Gaspard-Monge program for Optimization and Operations Research (PGMO) Project “Models for planning energy investment under uncertainty”.
Rights and permissions
About this article
Cite this article
Fukushima, M., Júdice, J., de Oliveira, W. et al. A sequential partial linearization algorithm for the symmetric eigenvalue complementarity problem. Comput Optim Appl 77, 711–728 (2020). https://doi.org/10.1007/s10589-020-00226-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10589-020-00226-7