to appear in Journal of Inequalities and Applications, 2016
An alternative approach for a distance inequality associated with the second-order cone and the circular cone
Xin-He Miao1
Department of Mathematics Tianjin University, China
Tianjin 300072, China
Yen-chi Roger Lin 2 Department of Mathematics National Taiwan Normal University
Taipei 11677, Taiwan
Jein-Shan Chen 3 Department of Mathematics National Taiwan Normal University
Taipei 11677, Taiwan.
April 28, 2016
(revised on November 10, 2016)
Abstract. It is known that the second-order cone and the circular cone have many anal- ogous properties. In particular, there exists an important distance inequality associated with second-order cone and circular cone. The inequality indicates that the distances of arbitrary points to the second-order cone and the circular cone are equivalent, which is crucial in analyzing the tangent cone and normal cone for the circular cone. In this paper, we provide an alternative approach to achieve the aforementioned inequality. Al- though the proof is a bit longer than the existing one, the new approach offers a way to clarify when the equality holds. Such clarification is helpful for further studying in the relationship between the second-order cone programming problems and the circular cone programming problems.
Keywords. Second-order cone, circular cone, projection, distance.
1E-mail: xinhemiao@tju.edu.cn.
2E-mail: yclin@math.ntnu.edu.tw
3Corresponding author. E-mail:jschen@math.ntnu.edu.tw.
1 Introduction
The circular cone [3, 14] is a pointed closed convex cone having hyperspherical sections orthogonal to its axis of revolution about which the cone is invariant to rotation. Let Lθ denote the circular cone in IRn, which is defined by
Lθ := x = (x1, x2) ∈ IR × IRn−1 | kxk cos θ ≤ x1
= x = (x1, x2) ∈ IR × IRn−1 | kx2k ≤ x1tan θ , (1) with k · k denoting the Euclidean norm and θ ∈ (0,π2). When θ = π4, the circular cone Lθ reduces to the well-known second-order cone (SOC) Kn [2, 6] (also called the Lorentz cone), i.e.,
Kn:=(x1, x2) ∈ IR × IRn−1 | kx2k ≤ x1 .
In particular, K1 is the set of nonnegative reals IR+. It is well known that the second- order cone Kn is a special kind of symmetric cones [4]. But when θ 6= π4, the circular cone Lθ is a non-symmetric cone [3, 13, 15].
In [14], Zhou and Chen showed that there is a special relationship between the SOC and the circular cone as follows:
x ∈ Lθ ⇐⇒ Ax ∈ Kn with A = tan θ 0 0 In−1
, (2)
where In−1is the (n − 1) × (n − 1) identity matrix. Based on the relationship (2) between the SOC and circular cone, Miao et al. [12] showed that circular cone complementar- ity problems can be transformed into the second-order cone complementarity problems.
Furthermore from the relationship (2), there hold
x ∈ int Lθ ⇐⇒ Ax ∈ int Kn and x ∈ bd Lθ ⇐⇒ Ax ∈ bd Kn.
Besides the relationship between second-order cone and circular cone, some topological structures play important roles in theoretical analysis for optimization problems. For example, the projection formula onto a cone facilitates designing algorithms for solving conic programming problems [1, 5, 7, 9, 16]; the distance formula from an element to a cone is an important factor in the approximation theory; the tangent cone and normal cone are crucial in analyzing the structure of the solution set for optimization problems [8, 10, 11]. From above illustrations, it raises an interesting question: What is the rela- tionship between second-order cone and circular cone regarding the projection formula, the distance formula, the tangent cone and normal cone, and so on? The issue of the tangent cone and normal cone has been studied in [14, Theorem 2.3]. In this paper, we focus on the other two issues.
More specifically, we provide an alternative approach to achieve an inequality which was obtained in [14, Theorem 2.2]. Although the proof is a bit longer than the existing
one, the new approach offers a way to clarify when the equality holds which is helpful for further studying in the relationship between the second-order cone programming prob- lems and the circular cone programming problems.
In order to study the relationship between second-order cone and circular cone, we need to recall some background materials. For any vector x = (x1, x2) ∈ IR × IRn−1, the spectral decomposition of x with respect to second-order cone is given by
x = λ1(x)u(1)x + λ2(x)u(2)x , (3) where λ1(x), λ2(x), u(1)x , and u(2)x are expressed as
λi(x) = x1+ (−1)ikx2k and u(i)x = 1 2
1
(−1)iw
i = 1, 2, (4) with w = kxx2
2k if x2 6= 0, or any vector in IRn−1 satisfying kwk = 1 if x2 = 0. In the setting of circular cone Lθ, Zhou and Chen [14] gave the following spectral decomposition of x ∈ IRn with respect to Lθ:
x = µ1(x)vx(1)+ µ2(x)vx(2), (5) where µ1(x), µ2(x), vx(1), and vx(2) are expressed as
µ1(x) = x1− kx2k cot θ,
µ2(x) = x1+ kx2k tan θ, and
v(1)x = 1+cot1 2θ
1
− cot θ · w
, v(2)x = 1+tan1 2θ
1
tan θ · w
,
(6)
with w = kxx2
2k if x2 6= 0, or any vector in IRn−1 satisfying kwk = 1 if x2 = 0. Moreover, λ1(x), λ2(x) and u(1)x , u(2)x are called the spectral values and the spectral vectors of x associated with Kn, whereas µ1(x), µ2(x) and vx(1), vx(2) are called the spectral values and the spectral vectors of x associated with Lθ, respectively.
To proceed, we denote x+ (resp. xθ+) the projection of x onto Kn (resp. Lθ); also we set a+ = max{a, 0} for any real number a ∈ IR. According to the spectral decompositions (3) and (5) of x, the expressions of x+ and xθ+ can be obtained explicitly, as stated in the following lemma.
Lemma 1.1 ([6, 14]). Let x = (x1, x2) ∈ IR×IRn−1have the spectral decompositions given as (3) and (5) with respect to SOC and circular cone, respectively. Then the following hold:
(a)
x+ = (x1− kx2k)+u(1)x + (x1+ kx2k)+u(2)x
=
x, if x ∈ Kn,
0, if x ∈ −(Kn)∗ = −Kn, u, otherwise,
where u =
" x
1+kx2k x1+kx22k
2 x2
kx2k
#
;
(b)
xθ+ = (x1− kx2k cot θ)+u(1)x + (x1+ kx2k tan θ)+u(2)x
=
x, if x ∈ Lθ,
0, if x ∈ −(Lθ)∗ = −Lπ
2−θ, v, otherwise,
where v =
" x1+kx2k tan θ 1+tan2θ
x1+kx2k tan θ
1+tan2θ tan θ
x2
kx2k
# .
Based on the expression of the projection xθ+onto Lθ in Lemma 1.1, it is easy to obtain that for any x = (x1, x2) ∈ IR × IRn−1, the explicit formula of projection of x ∈ IRn onto the dual cone L∗θ (denoted by (xθ)∗+) is as below:
(xθ)∗+ = (x1− kx2k tan θ)+u(1)x + (x1+ kx2k cot θ)+u(2)x
=
x, if x ∈ L∗θ = Lπ
2−θ, 0, if x ∈ −(L∗θ)∗ = −Lθ, ω, otherwise,
where ω =
" x1+kx2k cot θ
1+cot2θ
x
1+kx2k cot θ
1+cot2θ cot θ
x2
kx2k
# .
2 Main results
In this section, we give the main results of this paper.
Theorem 2.1. For any x ∈ IRn, let xθ+ and (xθ)∗+ be the projections of x onto the circular cone Lθ and its dual cone L∗θ, respectively. Let A be the matrix defined as in (2). Then the following hold:
(a) If Ax ∈ Kn, then (Ax)+= Axθ+.
(b) If Ax ∈ −Kn, then (Ax)+ = A(xθ)∗+ = 0.
(c) If Ax /∈ Kn ∪ (−Kn), then (Ax)+ = 1+tan2 2θA−1(xθ)∗+, where (xθ)∗+ retains its expression only in the case of x /∈ L∗θ∪ (−Lθ).
Theorem 2.2. Let x = (x1, x2) ∈ IR×IRn−1 have the spectral decompositions with respect to the SOC and the circular cone given as in (3) and (5), respectively. Then the following hold:
(a) dist(Ax, Kn) = q1
2(x1tan θ − kx2k)2−+ 12(x1tan θ + kx2k)2−. (b) dist(x, Lθ) =
q cot2θ
1+cot2θ(x1tan θ − kx2k)2−+ 1+tantan2θ2θ(x1cot θ + kx2k)2−, where (a)− = min{a, 0}.
Now, applying Theorem 2.2, we can obtain the relation on the distance formulas associated with the second-order cone and the circular cone. Note that when θ = π4, we know Lθ = Kn and Ax = x. Thus, it is obvious that dist(Ax, Kn) = dist(x, Lθ). In the following theorem, we only consider the case θ 6= π4.
Theorem 2.3. For any x = (x1, x2) ∈ IR × IRn−1, according to the expressions of the distance formulas dist(Ax, Kn) and dist(x, Lθ), the following hold.
(a) For θ ∈ (0,π4), we have
dist(Ax, Kn) ≤ dist(x, Lθ) ≤ cot θ · dist(Ax, Kn).
(b) For θ ∈ (π4,π2), we have
dist(x, Lθ) ≤ dist(Ax, Kn) ≤ tan θ · dist(x, Lθ).
3 Proofs of main results
3.1 Proof of Theorem 2.1
(a) If Ax ∈ Kn, by the relationship (2) between the SOC and the circular cone, we have x ∈ Lθ. Thus, it is easy to see that (Ax)+= Ax = Axθ+.
(b) If Ax ∈ −Kn, we know that −Ax ∈ Kn, which implies (Ax)+ = 0. Besides, combining with (2), we have −x ∈ Lθ, which leads to (xθ)∗+= 0. Hence, we have (Ax)+= A(xθ)∗+= 0.
(c) If Ax /∈ Kn∪ (−Kn), from Lemma 1.1(a), we have
(Ax)+ =
" x
1tan θ+kx2k x1tan θ+kx2 2k
2
x2
kx2k
#
= tan θ
2 (1 + cot2θ)
" x
1+kx2k cot θ 1+cot2θ x1+kx2k cot θ
1+cot2θ x2
kx2k
#
= 1 + tan2θ 2 · A−1
" x
1+kx2k cot θ 1+cot2θ x1+kx2k cot θ
1+cot2θ · cot θ · kxx2
2k
#
= 1 + tan2θ
2 · A−1(xθ)∗+. The proof is complete. 2
Remark 3.1. Here, we say a few more words for the part(c) in Theorem 2.1. Indeed, if Ax /∈ Kn∪ (−Kn), there are two cases for the element x ∈ IRn, i.e., x /∈ L∗θ ∪ (−Lθ) or x ∈ L∗θ. When x /∈ L∗θ∪(−Lθ), the relationship between (Ax)+and (xθ)∗+ is just as stated in Theorem 2.1(c), that is, (Ax)+= 1+tan2 2θA−1(xθ)∗+. However, when x ∈ L∗θ, we have (xθ)∗+= x. This implies that the relationship between (Ax)+ and (xθ)∗+ is not very clear.
Hence, the relation between (Ax)+ and (xθ)∗+ in Theorem 2.1(c) is a bit limited.
3.2 Proof of Theorem 2.2
(a) For any x = (x1, x2) ∈ IR × IRn−1, from the spectral decomposition (3) with respect to the SOC, we have Ax = (x1tan θ − kx2k)u(1)x + (x1tan θ + kx2k)u(2)x , where u(1)x and u(2)x
are given as in (4). It follows from Lemma 1.1(a) that (Ax)+ = (x1tan θ − kx2k)+u(1)x + (x1tan θ + kx2k)+u(2)x . Hence, we obtain the distance dist(Ax, Kn) as bellow:
dist(Ax, Kn) = kAx − (Ax)+k
=
(x1tan θ − kx2k)−u(1)x + (x1tan θ + kx2k)−u(2)x
= r1
2(x1tan θ − kx2k)2−+ 1
2(x1tan θ + kx2k)2−.
(b) For any x = (x1, x2) ∈ IR × IRn−1, from the spectral decomposition (5) with respect to circular cone and Lemma 1.1(b), with the same argument, it is easy to get that
dist(x, Lθ) =
x − xθ+
= s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2−+ tan2θ
1 + tan2θ(x1cot θ + kx2k)2−. 2
3.3 Proof of Theorem 2.3
(a) For θ ∈ (0,π4), we have 0 < tan θ < 1 < cot θ and Lθ ⊂ Kn ⊂ L∗θ. We discuss three cases according to x ∈ Lθ, x ∈ −L∗θ and x /∈ Lθ∪ (−L∗θ).
Case 1. If x ∈ Lθ, then Ax ∈ Kn which clearly yields dist(Ax, Kn) = dist(x, Lθ) = 0.
Case 2. If x ∈ −L∗θ, then x1cot θ ≤ −kx2k and dist(x, Lθ) = kxk =
q
x21+ kx2k2.
Under this case, there are two subcases for the element Ax. If x1cot θ ≤ x1tan θ ≤
−kx2k, i.e., Ax ∈ −Kn, it follows that dist(Ax, Kn) = kAxk =
q
x21tan2θ + kx2k2 ≤ q
x21+ kx2k2
= dist(x, Lθ) ≤ q
x21 + kx2k2cot2θ
= cot θ · q
x21tan2θ + kx2k2
= cot θ · dist(Ax, Kn),
where the first inequality holds since tan θ < 1 (it becomes equality only in the case of x = 0), and the second inequality holds since cot θ > 1 (it becomes equality only in the case of x2 = 0). On the other hand, if x1cot θ ≤ −kx2k < x1tan θ ≤ 0, we have
dist(Ax, Kn) = kAx − (Ax)+k
= r1
2(x1tan θ − kx2k)2−+ 1
2(x1tan θ + kx2k)2−
= r1
2(x1tan θ − kx2k)2 ≤ q
x21tan2θ + kx2k2 <
q
x21+ kx2k2
= dist(x, Lθ)
≤ q
x21 − x1kx2k cot θ
<
r1
2x21+1
2kx2k2cot2θ − x1kx2k cot θ
= cot θ · dist(Ax, Kn),
where the third inequality holds because kx2k ≤ −x1cot θ, and the fourth inequality holds since kx2k > −x1tan θ ≥ 0. Therefore, for the subcases of x ∈ −L∗θ, we can conclude that
dist(Ax, Kn) ≤ dist(x, Lθ) ≤ cot θ · dist(Ax, Kn), and dist(x, Lθ) = cot θ · dist(Ax, Kn) holds only in the case of x2 = 0.
Case 3. If x /∈ Lθ∪ (−L∗θ), then −kx2k tan θ < x1 < kx2k cot θ, which yields x1tan θ <
kx2k and x1cot θ > −kx2k. Thus, we have dist(x, Lθ) =
x − xθ+
= s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2−+ tan2θ
1 + tan2θ(x1cot θ + kx2k)2−
= s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2.
On the other hand, it follows from −kx2k tan θ < x1 < kx2k cot θ and θ ∈ (0,π4) that
−kx2k < −kx2k tan2θ < x1tan θ < kx2k.
This implies that
dist(Ax, Kn) = kAx − (Ax)+k
= r1
2(x1tan θ − kx2k)2−+ 1
2(x1tan θ + kx2k)2−
= r1
2(x1tan θ − kx2k)2. From this and θ ∈ (0,π4), we see that
dist(Ax, Kn) = r1
2(x1tan θ − kx2k)2 <
s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2
= dist(x, Lθ) =
r 2
1 + tan2θ r1
2(x1tan θ − kx2k)2
=
r 2
1 + tan2θdist(Ax, Kn).
Therefore, in these cases of x /∈ Lθ∪ (−L∗θ), we can conclude dist(Ax, Kn) < dist(x, Lθ) =
r 2
1 + tan2θdist(Ax, Kn), To sum up, from all the above and the fact that max{cot θ,q
2
1+tan2θ} = cot θ for θ ∈ (0,π4), we obtain that
dist(Ax, Kn) ≤ dist(x, Lθ) ≤ cot θ · dist(Ax, Kn).
(b) For θ ∈ (π4,π2), we have 0 < cot θ < 1 < tan θ and L∗θ ⊂ Kn ⊂ Lθ. Again we discuss the following three cases.
Case 1. If x ∈ Lθ, then Ax ∈ Kn which implies that dist(Ax, Kn) = dist(x, Lθ) = 0.
Case 2. If x ∈ −L∗θ, then x1cot θ ≤ −kx2k and dist(x, Lθ) = kxk =
q
x21+ kx2k2.
It follows from x1cot θ ≤ −kx2k and θ ∈ (π4,π2) that x1tan θ ≤ x1cot θ ≤ −kx2k, which leads to Ax ∈ −Kn. Hence, we have
dist(Ax, Kn) = kAxk = q
x21tan2θ + kx2k2.
With this, it is easy to verify that dist(Ax, Kn) ≥ dist(x, Lθ) for θ ∈ (π4,π2). Moreover, we note that
tan θ · dist(x, Lθ) = q
tan2θ(x21+ kx2k2) ≥ q
x21tan2θ + kx2k2 = dist(Ax, Kn).
Thus, it follows that
dist(x, Lθ) ≤ dist(Ax, Kn) ≤ tan θ · dist(x, Lθ), and dist(Ax, Kn) = tan θ · dist(x, Lθ) holds only in the case of x2 = 0.
Case 3. If x /∈ Lθ∪ (−L∗θ), then we have −kx2k tan θ < x1 < kx2k cot θ and dist(x, Lθ) =
x − xθ+
= s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2−+ tan2θ
1 + tan2θ(x1cot θ + kx2k)2−
= s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2.
Since −kx2k tan θ < x1 < kx2k cot θ, it follows immediately that −kx2k tan2θ < x1tan θ <
kx2k. Again, there are two subcases for the element Ax. If −kx2k tan2θ < −kx2k <
x1tan θ < kx2k, then we have Ax /∈ Kn∪ (−Kn). Thus, it follows that dist(Ax, Kn) = kAx − (Ax)+k
= r1
2(x1tan θ − kx2k)2−+ 1
2(x1tan θ + kx2k)2−
= r1
2(x1tan θ − kx2k)2. This together with θ ∈ (π4,π2) yields
dist(Ax, Kn) > dist(x, Lθ).
Moreover, by the expressions of dist(Ax, Kn) and dist(x, Lθ), it is easy to verify dist(Ax, Kn) =
r 2
1 + tan2θdist(x, Lθ).
On the other hand, if −kx2k tan2θ < x1tan θ ≤ −kx2k < kx2k, then we have Ax ∈ −Kn, which implies
dist(Ax, Kn) = kAxk = q
x21tan2θ + kx2k2. Therefore, it follows that
dist(x, Lθ) = s
cot2θ
1 + cot2θ(x1tan θ − kx2k)2
<
r1
2(x1tan θ − kx2k)2 ≤ q
x21tan2θ + kx2k2
= dist(Ax, Kn).
Since
tan2θ(x1tan θ − kx2k)2− (1 + tan2θ)(x21tan2θ + kx2k2)
= −2x1kx2k tan3θ − x21tan2θ − kx2k2
= −x1tan θ · (kx2k tan2θ + x1tan θ) + kx2k(−x1tan3θ − kx2k)
≥ kx2k(−x1tan3θ − kx2k)
≥ 0,
where the first inequality holds due to −kx2k tan2θ < x1tan θ, and the second inequality holds due to x1tan θ < −kx2k and θ ∈ (π4,π2), we have
dist(Ax, Kn) ≤ tan θ · dist(x, Lθ).
From all the above analysis and the fact that max{tan θ,q
2
1+tan2θ} = tan θ for θ ∈ (π4,π2), we can conclude that
dist(x, Lθ) ≤ dist(Ax, Kn) ≤ tan θ · dist(x, Lθ).
Thus, the proof is complete. 2
Remark 3.2. We point out that Theorem 2.3 is equivalent to the results in [14, Theorem 2.2], that is, for any x, z ∈ IRn, we have
kAk−1 dist(Az, Kn) ≤ dist(z, Lθ) ≤ kA−1k dist(Az, Kn) (7) and
kA−1k−1 dist(A−1x, Lθ) ≤ dist(x, Kn) ≤ kAk dist(A−1x, Lθ). (8) However, the above inequalities depends on the factors kAk and kA−1k. Here, we provide a more concrete and simple expression for the inequality. What is the benefit of such a new expression? Indeed, the new approach provides the situation where the equality holds, which is helpful for further studying in the relationship between the second-order cone programming problems and the circular cone programming problems. In particular, from the proof of Theorem 2.3, it is clear that dist(Ax, Kn) = tan θ · dist(x, Lθ) holds only under the cases of x = (x1, x2) ∈ Lθ or x2 = 0; otherwise we would have the strict inequality dist(Ax, Kn) < tan θ · dist(x, Lθ). To the contrast, it takes tedious algebraic manipulations to obtain such situations by using (7) and (8).
The following example elaborates more about why dist(Ax, Kn) = tan θ · dist(x, Lθ) holds only under the cases of x = (x1, x2) ∈ Lθ or x2 = 0.
Example 3.1. Let x = (x1, x2) ∈ IR × IRn−1 and A = tan θ 0 0 In−1
. When x ∈ Lθ, we have Ax ∈ Kn. It is clear to see that dist(Ax, Kn) = tan θ · dist(x, Lθ) = 0. When x2 = 0, i.e., x = (x1, 0) ∈ IR × IRn−1, it follows that Ax = (x1tan θ, 0). If x1 ≥ 0, we have x ∈ Lθ and Ax ∈ Kn, which implies that dist(Ax, Kn) = tan θ · dist(x, Lθ) = 0. In other case, if x1 < 0, we get that x ∈ −L∗θ and Ax ∈ −Kn. All the above gives that dist(Ax, Kn) = kAxk = |x1| tan θ = tan θ · dist(x, Lθ).
Competing interests The authors declare that none of the authors have any competing interests in the manuscript.
Authors contributions All authors participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Acknowledgments We would like to thank the Editor and the anonymous referees for their careful reading and constructive comments which have helped us to significantly improve the presentation of the paper. The first author’s work is supported by Na- tional Natural Science Foundation of China (No. 11471241). The third author’s work is supported by Ministry of Science and Technology, Taiwan.
References
[1] Chen, JS: A new merit function and its related properties for the second-order cone complementarity problem. Pac. J. Optim. 2, 167–179 (2006).
[2] Chen, JS, Chen, X, Tseng, P: Analysis of nonsmooth vector-valued functions associ- ated with second-order cones. Math. Program. 101, 95–117 (2004).
[3] Pznto Da Costa, A, Seeger, A: Numerical resolution of cone-constrained eigenvalue problems. Comput. Appl. Math. 28, 37–61 (2009).
[4] Faraut, U, Kor´anyi, A: Analysis on symmetric cones. Oxford Mathematical Mono- graphs, Oxford University Press, New York (1994).
[5] Fukushima, M, Luo, ZQ, Tseng, P: Smoothing functions for second-order cone com- plementarity problems. SIAM J. Optim. 12, 436–460 (2001).
[6] Facchinei, F, Pang, J: Finite-dimensional Variational Inequalities and Complemen- tarity Problems. Springer, New York (2003).
[7] Hayashi, S, Yamashita, N, Fukushima, M: A combined smoothing and regularization method for monotone second-order cone complementarity problems. SIAM J. Optim.
15, 593–615 (2005).
[8] Jeyakumar, V, Yang, XQ: characterizing the solution sets of pseudolinear programs.
J. Optim. Theory Appl. 87, 747–755 (1995).
[9] Kong, LC, Sun, J, Xiu, NH: A regularized smoothing Newton method for symmetric cone complementarity problems. SIAM J. Optim. 19, 1028–1047 (2008).
[10] Mangasarian, OL: A simple characterization of solution sets of convex programs.
Oper. Res. Lett. 7, 21–26 (1988).
[11] Miao, XH, Chen, JS: Characterizations of solution sets of cone-constrained convex programming problems. Optim. Lett. 9, 1433–1445 (2015).
[12] Miao, XH, Guo, SJ, Qi, N, Chen, JS: Constructions of complementarity functions and merit functions for circular cone complementarity problem. Comput. Optim.
Appl. 63, 495–522 (2016).
[13] Nesterov, Y: Towards nonsymmetric conic optimization. Center for Operations Re- search and Econometrics. 28, (2006).
[14] Zhou, JC, Chen, JS: Properties of circular cone and spectral factorization associated with circular cone. J. Nonlinear convex A. 14, 807–816 (2013).
[15] Zhou, JC, Chen, JS, Hung, HF: Circular cone convexity and some inequalities asso- ciated with circular cones. J. Inequal. Appl. 2013, Article ID 571, 17 pages (2013).
[16] Zhang, Y, Huang, ZH: A nonmonotone smoothing-type algorithm for solving a sys- tem of equalities and inequalities. J. Comput. Appl. Math. 233, 2312–2321 (2010).