On the Solution of the Eigenvalue Complementarity Problem by a Line Search Filter-SQP Algorithm

In this paper, the Eigenvalue Complementarity Problem (EiCP) with real symmetric matrices is addressed, which appears in the study of contact problem in mechanics. We discuss a quadratic programming formulation to the problem. The resulting problems are nonlinear programs that can be solved by a line search filter-SQP algorithm.


Introduction
The Eigenvalue Complementarity Problem (EiCP) appears in the study of static equilibrium states of finite dimensional mechanical systems with unilateral frictional contact [1] [2] and takes the following form: where B is a positive definite matrix.Note that any solution with 0 ω = correspond to a solution of the (GEiCP).
The Generalized Eigenvalue Complementarity Problem (GEiCP) is ( ) ( ) T 0 GEiCP : Find 0, 0,such that 0 0 0 where B is a positive definite matrix, is given, and . The (EiCP) is clearly the particular case of the (GEiCP) J with For any solution ( ) , x λ ω , the value λ is called a (general) complementary eigenvalue of ( ) , A B , and x is a corresponding (general) complementary eigenvector [3].
In this paper, the symmetric Eigenvalue Complementarity Problem is studied.
Some properties of EiCP are derived, including a necessary and sufficient condition for solvability.When this condition is verified, we can easily obtain an initial point.For the design of an initialization algorithm, we discuss a quadratic programming formulation and a line search filter-SQP algorithm is introduced to solve the (EiCP).

The Symmetric Eigenvalue Complementarity Problem
In this section, the symmetric EiCP is considered, where the matrices A and B are both symmetric.In this case, the EiCP is closely related to the classical Eigenvalue Problem , and some properties are derived.
Proposition 2.1.Proposition The solvability of the (GEiCP) J is in general a NP-complete problem.
It follows that solving the (GEiCP) is in general a NP-hard problem.Despite this fact, for some classes of matrices the solvability of the corresponding (EiCP) can be answered easily.
Since the set of complementary eigenvectors of a given complementary eigenvalue is a cone, there is no loss of generality in restricting the problem to finding solutions satisfying  : 0; , : 0, 0 and 0; A is a S-matrix ( Then a point x X ∈ such that T 0 x Ax > can be easily obtained and the corresponding (EiCP) is solvable.
An equivalent way to formulate the (EiCP) is through the quadratic formulation ( ) Theorem 2.1.If A is strictly copositive and 0 x ≥ is a stationary point of (QB), then the pair x , is a solution of (EiCP).

A Line Search Filter-SQP Algorithm
In this section a line search filter-SQP algorithm is introduced to solve the formulation of the previous section.As we all know, sequential Quadratic Programming (SQP) algorithm is one of the most efficient methods for the numerical solution of constrained nonlinear optimization problems [6].In the previous methods, penalty function is always used as a merit function to test the acceptability of the iterate points.However, there are several difficulties associated with the use of penalty function and in particular with the choice of the penalty parameter.Too low a choice may result in the loss of an optimal solution; on the other hand, too large a choice damps out the effect of the objective function.Fletcher and Leyffer proposed a filter method as an alternative to the traditional merit function for solving the nonlinear optimization [7].The significant advantage of the filter method is that it does not need to estimate the penalty parameter which could be difficult to obtain.Recently, the filter method was developed in [8] [9].
Considering a suitable continuously differentiable merit function ( ) f x , it is possible to reduce the (NLP) to the following nonlinear program ( ) ( ) ( ) where the objective function  and the constraints In SQP methods, at each iteration the search direction is generally obtained by solving the subproblem as follows: However, the previous QP subproblem has a serious shortcoming that constraints in (QP 1 ) may be inconsistent.To overcome this disadvantage, much attention has been paid [10] [11] [12] [13].In addition, Liu and Li [12] and Liu and Zeng [13] proposed SQP algorithms with cautious update criteria, which can be considered as modifications of the SQP algorithm given in [11].In our algorithm, stimulated by [10] [11] [12] [13], the quadratic subproblem (QP 1 ) is replaced by the following problem: , 0 where k b is a positive parameter.Clearly, this subproblem is always consistent and is a convex programming if k B is positive semidefinite.The KKT condition for the subproblem ( ) ( ) In this paper, the Lagrangian function value instead of the objective function value is used in the filter together with an appropriate infeasibility measure.
( ) ( ) ( ) and we define the constraints violation function ( ) where λ is the Lagrangian multiplier corresponding to NLP satisfied at x and In the line search filter technique, for fixed constants Similar to the traditional strategy of the filter method, to avoid the convergence to a feasible point but not an optimal solution, we consider the following L-type [14] [15] switching condition: where 0, 0, 1 then the Armijo-type reduction condition is employed as follows: ( where 1 0, 2 For the sake of a simplified notation, the filter is defined in this paper not as a list but as a set , h L pairs that are prohibited in iteration k.We say that a trial point ( ) x α is acceptable to the filter if its ( ) , h L pair does not in k F .Given a 0 F , throughout the algorithm ,the filter is then augmented in some iterations after the new iterate point x + has been accepted.We use the updating formula , In the situation where the trial step size is too small to guarantee sufficient reduction as defined by (10), the method switches to a feasibility restoration phrase, whose purpose is to find a new iterate point that satisfies (10) and is also acceptable to the current filter by trying to decrease the constraint violation.In order to detect the situation where no admissible step size can be found and the restoration phase has to be invoked, we approximate a minimum desired step size using linear models of involved functions.For this, we define that δ ω ω are form (11) and , h L γ γ are form (10).The proposed algorithm switches to the feasibility restoration phase when ( ) , and 0, 1 Step 2. Solve the subproblem ( ) to obtain the solution ( ) and let ( ) Step 3. If 0 and , F and go to Step 5.5; otherwise, go to Step 5.4.
Step 5. Step 6. Augment filter if it is necessary.If k is not an L-type iteration, augment the filter using (13); otherwise leave the filter unchanged; that is, set x + from a feasible restoration phrase.Set

Global Convergence of Algorithm
In this section, we will show the proposed algorithm is well defined and globally convergent under some mild donditions.Let x be a feasible point of NLP (1) and λ be the multiplies corresponding to the nonlinear constraints satisfied at x .Denote the ith component of λ by i λ .Then throughout this paper, we always assume that them following assumptions hold: Assumptions 1) The function ( ) , L x λ and ( ) g x are twice continuously differentiable.
2) The sequence { } k x remains in a compact and convex subset n X ∈  .
3) The strict Mangasarian-Fromowitz constraint qualification (SMFCQ) condition holds when x is a feasible point of the NLP.Then for x , the . And there exist two constants 0 q p ≥ > such that the matrices sequence { } , shows that the search direction, the optimal solution of the subproblem is a sufficient descent direction for Lagrangian function at points that are nonoptimal and sufficiently close to feasible point.And the next lemma shows that no ( ) , h L -pair corresponding to a feasible point is ever including in the filter.  .Then there exists a constant 0 α > so that for all i k and α α < , we have The Lemma 4.5 shows that there exists a step length bounded away from zero so that the Armijo condition for the Lagrangian function is satisfied under some conditions.The next lemma shows that there is no cycle between Step 3 and Step 7 in above algorithm.The proof of the Lemma 4.6 is similar to Lemma 4 in [16].
Lemma 4.6.Suppose that Assumptions G hold.Then the inner loop terminates in a finite number of iterations.
Lemma 4.7.Suppose that the filter is augmented only a finite number of times; that is, There the proof is stated for slightly different circumstances, but it is easy to verify that it is still valid in our context.And we now show the global convergence result under some mild conditions.The proof of the next two theorems is similar to the Theorem 10 and Theorem 11 in [17].Theorem 4.2.Suppose that all stated assumptions hold.Furthermore, assume that { } k x is an infinite sequence generated by Algorithm and < ∞ L .Then every limit point is a KKT point.
Theorem 4.3.Suppose that all stated assumptions hold.Furthermore, assume that { } k x is an infinite sequence generated by Algorithm and = ∞ L .Then there exists at least one accumulation point which is a KKT point.

Numerical Experiments
In this section, some computational experience is presented to illustrate the efficiency of the algorithms described in this paper for the solution of symmetric , where i e is the vector i of the canonical basis.
The parameters used in the algorithm are follows:  Based on the results, we claim that this algorithm is a efficient procedure to solve Symmetric Eigenvalue Complementarity Problems.

Conclusion
In this paper, the Eigenvalue Complementarity Problem with real symmetric matrices have been considered, and we reformulate the EiCP to a NLP and use a line search filter-SQP algorithm to solve them.The numerical experiments show that this method is a promising method for solving the EiCP.Actually, we can solve the problem as a quadratically constrained quadratic programming problem not a NLP.And the design of efficient procedure for solving the quadratically constrained quadratic programming subproblems is also one of our current research areas.

2 1x
= , which replaces the constraint 0 x ≠ .In case of the (EiCP) the linear constraint B are both symmetric, the (EiCP) is closely related to the classical Eigenvalue Problem[4] [5].The complementary condition T algorithm for solving the inequality constrained optimization problem (NLP) can be stated as follows: Algorithm Step 1. Initialization: choose an initial point 0 n x ∈  , an initial filter

2 .
go to the feasible restoration phrase in Step 9. Compute a new trial point.If the trial step size to the filter: if the trial point the trial step size and go to Step 5.4.Step 5.3.Check sufficient decrease with respect to current iterate point: Step 5.3.1.Case I: condition (11) holds.If the Armijo condition (12) holds, accept the trial step (that is, an L-type iteration), set ( )

2 with k replaced by 1 k
+ .Step 9. Obtain a new point 1 k

+
, and go to Step 2. Remark 1.The mechanisms of the filter could ensure that Q. Yu et al.

Remark 2 .
The feasibility restoration phrase in Step 9 could be any iterative algorithm with the goal of finding a less infeasible point.Remark 3. Step 3 -Step 7 which are call the inner loop are backtracking line search whose goal is to find the k α

2 . 3 .
to zero, and all limit points of { } k x are feasible.Lemma 4.1.Suppose that Assumptions hold and that the filter is augmented only finite number of times.Then Suppose that Assumptions hold and let { } k x be the sequence of iterates generated by Algorithm so that the filter is augmented in every iteration kproof of the Lemma 4.1 and Lemma 4.2 are similar to Lemma 5 in [14] by replacing ( ) f x by ( ) , L x λ , we omit it.After we have the previous two lemmas, we can get the following theorem.That is to say that the limit point of { } k x is feasible.The proof of the next theorem is similar to the Theorem 1 Suppose that Assumptions hold and let x be a feasible point of NLP at which SMFCQ condition holds, but not a KKT point.If { } k x is a subsequence of iterate points for which ( ) 0 et al.

Lemma 4 . 5 .
Suppose that Assumptions hold.Let { } b .The test results are given in Table, the notations mean follows: No: the number of the problems, m: the dimension of the vector; T: the total CPU time in seconds for solving the problem; OPT: the final value of the objective function; λ : the eigenvalue of the EiCP.
Proposition 2.2.shows a necessary and sufficient condition for solvability.For the practical solution of the symmetric (EiCP), finding an initial point is a NP-hard problem.But for some class of matrices which an initial solution for the (OEiCP) can be easily obtained.
Analogously to the classical case, equilibrium points of the Rayleigh quotient in the nonnegative orthant with ( ) 0 x λ > are solutions of the (EiCP).This is the main result concerning he practical solution of the symmetric (EiCP).Q. Yu et al.Journal of Applied Mathematics and Physics Proposition 2.2.(EiCP) is solvable if and only if there exists some As described in Proposition 2.3, the initial solution 0x can be chosen by one several processes.In particular if A has at least one diagonal element 0 according to the procedure described in[18].The scaling is important because the matrices that we are using are badly conditioned, and without this tool some of the problems cannot be solved.Q. Yu et al.DOI: 10.4236/jamp.2017.5101681994 Journal of Applied Mathematics and Physics