In the previous paper by the first and the third authors, we present six algorithms for determining
whether a given symmetric matrix is strictly copositive, copositive (but not strictly), or not copositive.
The algorithms for matrices of order are not guaranteed to produce an answer. It also shows that
for 1000 symmetric random matrices of order 8, 9, and 10 with unit diagonal and with positive entries all
being less than or equal to 1 and negative entries all being greater than or equal to , there are 8, 6, and 2 matrices remaing undetermined, respectively. In this paper we give two more algorithms for and our experiment shows that no such matrix of order 8 or 9 remains undetermined; and almost always no
such matrix of order 10 remains undetermined. We also do some discussion based on our experimental
results.
1. Introduction
Reference [1] gives six algorithms for determining whether a given symmetric matrix is strictly copositive, copositive (but not strictly), or not copositive. The algorithms for matrices of order 3, 4, 5, 6 or 7 are efficient. But for matrices of order , it cannot guarantee to produce an answer. Table of [1] shows that for 1000 symmetric random matrices of order with unit diagonal and with positive entries all being less than or equal to 1 and negative entries all being greater than or equal to , there are 8, 6, and 2 matrices remaining undetermined when , respectively. In this paper we continue our study as in [1] and give two algorithms for and our experiment shows that no such matrix of order 8 or 9 remains undetermined; and almost always no such matrix of order 10 remains undetermined. We also do some discussion based on our experimental results.
In this paper we use all the concepts and notations of [1, 2] without explanation. Our main theorems will give the necessary and sufficient conditions for symmetric matrices of order 8 or 9 to be (strictly) copositive.
Let be symmetric and be partitioned into
with , . As in [2], let
be the simplex of order ; and let
be the standard simplex of order whose vertices are all vertices of . It is proved in [2] that an symmetric matrix is copositive if and only if for all . Consider the polyhedron in is the given vector of dimension in (1.1)) which has some vertices being vertices of , and all the other vertices being in the hyperplane . It is known (see [2, Section and Lemma ]) that the polyhedron can be subdivided into simplices in such that , is a subsimplex of and if , and the vertices of are all vertices of . We mention this fact since that is subdivided into simplices .
Denote the vertices of by , then is a vertex of , or a common point of the line connecting two vertices of and the hyperplane and should be presented in the barycenter coordinates of . If is the th vertex of , then it is represented by the coordinate vector with a 1 in the th position and all 0's elsewhere; otherwise write to denote that it is the common point of line and the hyperplane . Each determines a matrix (see [2, Lemma ]), to simplify the notation we still write with or . For example, if share only one vertex with and the other vertices are , then
Lemma 1.1 (see [2]). Let be symmetric and partitioned as in (1.1) with , being copositive and is subdivided into simplices which determine matrices . Then is copositive if and only if , are all copositive (see [2, Lemma 3.1]); is strictly copositive if and only if , are all strictly copositive and and is strictly copositive (see [1]).
It is noticed from [2] that if the polyhedron contains vertices (coordinate vectors of the standard simplex not in the hyperplane , then contains exact vertices in the hyperplane , and that can be subdivided into simplices of dimension such that is a simplex of dimension for and is a simplex of dimension when .
Lemma 1.2 (see [1]). Let . If there are -triples of pairwise different vertices of satisfying the following two conditions: (i)each contains at least one coordinate vector vertex;(ii) has exactly vertices for , and has less than vertices when , then can be subdivided into simplices , where is the simplex whose vertices are the elements of .
These two lemmas are basic for proving Theorems , 2.6, 2.7, and 2.8 in [1]; they are also basic for proving Theorems 2.1 and 2.2 of this paper.
2. Main Theorems and Algorithms
The following two theorems give two algorithms for determining the copositivity of a given symmetric matrix of order 8 or 9. These two theorems can be proved by Lemma 1.1 and Lemma 1.2 following the same pattern as in [1].
Theorem 2.1. Let be symmetric and be partitioned as in (1.1) and , then at least one of the following cases must happen:
(a) If one principal submatrix of is not copositive, then is not copositive. Otherwise it holds that and is copositive.
(b) If then is copositive; if with and is strictly copositive, then is strictly copositive.
(c) If , then is copositive if and only if is copositive; is strictly copositive if and only if is strictly copositive and and is strictly copositive.
(d) If has exactly one negative entry: , then is copositive if and only if is copositive; is strictly copositive if and only if is strictly copositive, and and is strictly copositive, where
(e) If has exactly two negative entries:, and , then is copositive if and only if and are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(f) If has exactly three negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(g) If has exactly four negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(h) If has exactly five negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(i) If has exactly six negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
It is clear (see [1, Remark ]) that if is odd, then a copositive matrix must have a row with an even number of negative entries. In other words, if a symmetric matrix of odd order has row with an even number of negative entries, then some principal submatrices of it are not copositive. This fact will be used in Theorem 2.2.
Theorem 2.2. If is symmetric, then at least one of the following cases must happen:
(a) If one principal submatrix of is not copositive, then is not copositive.
Otherwise ( must have a row with an even number of negative entries and is copositive) find a row of which has exactly negative entries. If the th row does, then interchange the th row and column with the first row and column, and partition into (1.1) as in Theorem 2.1.
(b) If , then and is copositive; if with and is strictly copositive, then is strictly copositive.
(c) If , then , then is copositive if and only if is copositive; is strictly copositive if and only if is strictly copositive and and is strictly copositive.
(d) If , then has exactly two negative entries: , and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(e) If , then has exactly four negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where
(f) If , then has exactly six negative entries: and , then is copositive if and only if are all copositive; is strictly copositive if and only if are all strictly copositive and and is strictly copositive, where