Fedoryuk values and stability of global H\"{o}lderian error bounds for polynomial functions
Huy-Vui H\`a, Phi-D\~ung Ho\`ang

TL;DR
This paper investigates the stability of global H"olderian error bounds for polynomial functions, providing criteria, explicit formulas involving Fedoryuk values, and classifications of stability types under perturbations.
Contribution
It introduces criteria and explicit formulas for the existence and stability of global H"olderian error bounds for polynomial sublevel sets, using Fedoryuk values.
Findings
Criteria for existence of global H"olderian error bounds.
Explicit formulas for the set of thresholds with error bounds.
Classification of stability types for error bounds.
Abstract
Let be a polynomial function of variables. In this paper, we study stability of global H\"{o}lderian error bound for a nonempty sublevel set under a perturbation of . In this paper, we give: * Criteria for the existence of a global H\"{o}lderian error bound of ; * Formulas for computing explicitly the set via some Fedoryuk values of and definition of threshold for the existence of global H\"{o}lderian error bound of ; * Definition of all types of stability of global H\"{o}lderian error bound of .
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsPolynomial and algebraic computation · Complexity and Algorithms in Graphs · Advanced Optimization Algorithms Research
Fedoryuk values and stability of global Hölderian error bounds for polynomial functions
HUY-VUI HÀ*†*
*†*Thang Long Institute of Mathematics and Applied Sciences,
Nghiem Xuan Yem Road,
Hoang Mai, District, Hanoi, Vietnam
and
PHI-DŨNG HOÀNG*‡*
*‡*Department of Mathematics - Faculty of Fundamental Sciences,
Laboratory of Applied Mathematics and Computing,
Posts and Telecommunications Institute of Technology,
Km10 Nguyen Trai Road, Ha Dong District, Hanoi, Vietnam
Abstract.
In this paper we study the stability of a global Hölderian error bound of the sublevel set under perturbation of , where is a polynomial function in real variables. Firstly, we give two formulas which compute the set
[TABLE]
via some special Fedoryuk values of . Then, based on these formulas, we can determine the stability type of a global Hölderian error bound of for any value .
Key words and phrases:
Error bounds, Stability, Polynomial Optimization
2010 Mathematics Subject Classification:
49K40, 14P10, 90C26
1. Introduction
Let be a polynomial function. For , put
[TABLE]
and .
Definition 1.1**.**
[Ha] We say that the nonempty set has a global Hölderian error bound (GHEB for short) if there exist such that
[TABLE]
Note that, if , then (1) becomes a global Lipschitzian error bound for .
The existence of error bounds have many important applications, including sensitivity analysis, convergence analysis in optimization problems, variational inequalities… After the earliest work by Hoffman ([Hoff]) and extended paper of Robinson ([Ro]), the study of error bounds has received rising awareness in many papers of mathematical programming in recent years, see [LL, WP, LS, Y, LiG1, LiG2, Ha, Ng, LMP, DHP] (for the case of polynomial functions) and [Hoff, Ro, M, AC, LiW, K, KL, P, LP, Luo, Jo, NZ, CM, LTW, I, BNPS, DL] (for non-polynomial cases). The reader is referred to survey papers [LP, P, Az, I] and the references therein for the theory and applications of error bounds.
Studying the stability of error bounds under perturbation is fundamental and hard problem. It has been investigated recently in the works of Daniel, Luo-Tseng, Deng, Ngai-Kruger-Théra, Kruger-Ngai-Théra, Kruger-López-Théra,… (see [Da, LT, D, NKT, KNT, KLT]).
In this paper, we study stability of a global Hölderian error bound for the set under a perturbation of , i.e. the perturbation of by a constant term. The following questions arise
Suppose that has a GHEB, when does there exist an open interval , such that for any , has also a GHEB? 2. 2.
Suppose that does not have GHEB, when does there exist an open interval , such that for any , also does not have GHEB? 3. 3.
Are there other types of stability which are different from types in questions 1 and 2?
To classify the stability types of GHEB, our idea is computing the set
[TABLE]
It turns out that the set can be determined via some speacial values of the Fedoryuk set of .
According [KOS], the Fedoryuk set of a polynomial is defined by
[TABLE]
We will show that there exists a value , which will be called the threshold of global Hölderian error bounds of and a subset of , such that
[TABLE]
Since is a semialgebraic subset of , this formula allows us answer the questions 1 and 2. Moreover, we can discover some other types of stability which are different from the types in questions 1-2 and give the list of all possible types of stability.
The paper is organized as follows. In Section 2, we give two different formulas for computing the set . The first formula is based on criterion for the existence of GHEB for , given in [Ha]. The second formula follows from a new criterion for the existence of global Hölderian error bounds. In Section 3, the relationship between and the set of Fedoryuk values of will be established. In Section 4, we use the formulas of and relationship between and to study our problems. It turns out that is a semialgebraic subset of , hence is either empty, or a finite set or a disjoint of finite number of points and intervals. Therefore, it is convenient to consider each of these cases separately.
In Subsection 4.1, we consider the case . In this case, or (Theorem 4.1). Therefore, there are two stability types of GHEB if . Namely, any point of is y-stable, by this we mean that and there exists an open interval such that . Besides, is y-right stable, by this we mean that and there exists such that and . Note that, for almost every polynomial , . Hence, or if is generic (Remark 4.1).
In Subsection 4.2, we consider the case when is a non-empty finite set. In this case, we show that
- •
(Proposition 4.1);
- •
Beside of y-stable type and y-right stable, there are at most 4 other stability types of GHEB. We have
**Case A: **
If , then there are 2 types
- **(i): **
is y-stable. 2. **(ii): **
is a n-isolated point: and for sufficiently small, .
**Case B: **
If is a finite value, then there are 5 types for all
- **1.: **
is y-stable; 2. **2.: **
is y-right stable; 3. **1’.: **
is n-stable: and there exists an open interval such that ; 4. **2’.: **
is n-right stable: and there exists such that and ; 5. **3’.: **
is n-left stable: and there exists such that and ; 6. **4’.: **
is a n-isolated point;
Note that:
- **–: **
If is y-right stable or is n-left stable, then it is necessarily that ;
- **–: **
If is n-right stable, then it is necessarily that and .
- •
We can determine the type of stability of any (Theorem 4.3);
- •
We give an estimation of the number of connected components of (Theorem 4.4);
In Subsection 4.3, we consider the case when . In this case
- •
Any value of belongs to one of the following types
is y-stable; 2. 2.
is y-right stable; 3. 3.
is y-left stable: and there exists such that and ; 4. 4.
is an y-isolated point: and for sufficiently small, ; 5. 1’.
is n-stable; 6. 2’.
is n-right stable; 7. 3’.
is n-left stable; 8. 4’.
is an n-isolated point.
- •
We can determine the type of stability of any (Theorem 4.5).
We conclude with some examples which illustrates some types of stability.
2. The set
2.1. The first formula of
Let be a polynomial function and .
Definition 2.1** ([DHN, Ha]).**
We say that
- (i)
A sequence is the first type of if
[TABLE] 2. (ii)
A sequence is the second type of if
[TABLE]
Theorem 2.1** ([Ha]).**
The following statements are equivalent:
- (i)
There are no sequences of the first or second types of . 2. (ii)
* has a GHEB, i.e. there exist such that*
[TABLE]
Put
[TABLE]
Definition 2.2**.**
Put
[TABLE]
We call the threshold of global Hölderian error bounds of .
Theorem 2.2** (The first formula of ).**
We have
- (i)
If , then ; 2. (ii)
If , then ; 3. (iii)
If , then ; 4. (iv)
If , then .
Proof.
Clearly, if and , then . Hence,
[TABLE]
Therefore, Theorem 2.2 follows from Theorem 2.1. ∎
2.2. A new criterion of the existence of a GHEB of and the second formula of
Let be the degree of a polynomial . By a linear change of coordinates, we can put in the form
[TABLE]
where and are polynomials in , where degrees .
Put .
Definition 2.3**.**
We say that
- (i)
A sequence is of the first type of w.r.t if
[TABLE] 2. (ii)
A sequence is of the second type of w.r.t if
[TABLE]
Let be of the form . Put
[TABLE]
Theorem 2.3**.**
Let be of the form . Then the following statements are equivalent
- (i)
There are no sequences of the first or second types of w.r.t ; 2. (ii)
* such that*
[TABLE]
for all ; 3. (iii)
* such that*
[TABLE]
for all ; 4. (iv)
* has a global Hölderian error bound.*
Proof.
We will prove that .
Proof of
For , put
[TABLE]
By (i), is well defined on . Moreover, it follows from Tarski-Seidenberg theorem (see, for example, [BCR, C, HP]), is a semialgebraic function.
To prove (ii), it is important to know the behavior of , as or . We distinguish 4 possibilities
- (a)
for sufficiently small and for sufficiently large; 2. (b)
for sufficiently small and for sufficiently large; 3. (c)
for sufficiently small and for sufficiently large; 4. (d)
both for sufficiently small and sufficiently large.
We will prove (i) (ii) for the case (d) because the proofs of other cases are similar.
In this case, since is semialgebraic and for any , we have
[TABLE]
and
[TABLE]
Clearly, . It follows from (2) that there exists such that
[TABLE]
for .
It follows from (3) that there exists sufficiently large, such that for any . We have
[TABLE]
if and
[TABLE]
if .
Since, by (i), there are no sequences of the second type, the function is bounded on the set
[TABLE]
This fact, together with (4), (5) and (6), give the proof of (i) (ii).
Proof of (ii) (iii):
The proof is based on the following classical result
Lemma** (van der Corput, [G]).**
Let be a real valued -function, , that satisfies for all . Then the following estimate is valid for all :
[TABLE]
Suppose that we have (ii). Then
- •
If , then and (iii) holds automatically.
- •
If , then (iii) follows from (ii).
Assume that .
Clearly
- •
(ii) holds if and only if there exists such that
[TABLE]
for all .
- •
(iii) holds if and only if there exists
[TABLE]
for all .
Let . We put
[TABLE]
and
[TABLE]
Since , it follows from the van der Corput Lemma that there exists a constant , independent of such that
[TABLE]
Clearly, and . Since is a closed semi-algebraic subset of , we have
[TABLE]
where , and
[TABLE]
Firstly, we see that . In fact, since is an isolated point of , is a local extremum of . Hence,
[TABLE]
or i.e. , while by assumption, . Thus, .
Without loss of generality, we may assume that . Since , we distinguish two cases
- •
If , then there exists such that , which means that or . Hence
[TABLE]
Then, by (9), (iii) holds.
- •
If , then, by Rolle’s Theorem, there exists such that
[TABLE]
which means that . Applying (7), there exists such that
[TABLE]
Moreover, since , we have
[TABLE]
Let be the point of such that
[TABLE]
We have
[TABLE]
Now:
- –
If , then
[TABLE]
- –
If , then
[TABLE]
Then (iii) follows from (10).
Hence, the implication (ii) (iii) is proved.
Proof of (iii) (iv):
Clearly, if (iii) holds, then there are no sequences of the first or second types of . Hence, by Theorem 2.1, (iv) holds.
The proof of (iv) (i) is straightforward. ∎
Proposition 2.1**.**
Let be a polynomial function and be a linear isomorphism. Then we have
[TABLE]
Proof.
Let and put .
Firstly, we prove that .
We have . This implies that
[TABLE]
Since , then there exists such that
[TABLE]
Suppose that , where or . Since and is a linear isomorphism, we have and there exists such that
[TABLE]
It follows that
[TABLE]
Combining (11), (12) and above fact, we have
[TABLE]
i.e., . The claim is proved similarly. ∎
We have the following theorem
Theorem 2.4** (The second formula of ).**
Let be a polynomial of the form . Then we have
- (i)
; 2. (ii)
If , then ; 3. (iii)
If , then ; 4. (iv)
If and , then ; 5. (v)
If and , then .
3. The relationship between and Fedoryuk values
The relationship between Fedoryuk values and the existence of global Hölderian error bounds is well-known and has been explored in many previous works, see, for example, [Az, CM, LP, Ha, I]. In this section, we will establish this relationship by proving that and . We recall
Definition 3.1**.**
Let be a polynomial function. The set of Fedoryuk values of is defined by
[TABLE]
Moreover, we have
Lemma 3.1**.**
* is a semialgebraic subset of .*
Remark 3.1**.**
It follows from Lemma 3.1 that either is empty or is finite set or is a union of finitely many points and intervals.
Note that can be an infinite set, for example (see [Par]), if , then and (see also [KOS] and [Sch]).
To prove the lemma, it is more convenient to use the logical formulation of the Tarski-Seidenberg Theorem. Let us to recall it.
A first-order formula is obtained as follows recursively (see, for example, [BCR, C, HP])
- (1)
If , then and are first-order formulas (with free variables ) and and are respectively the subsets of such that the formulas and hold. 2. (2)
If and are first-order formulas, then (conjunction), (disjunction) and (negation) are also first-order formulas. 3. (3)
If is a formula and is a variable ranging over , then and are first-order formulas.
Theorem** (Logical formulation of the Tarski–Seidenberg Theorem [BCR, C, HP]).**
If is a first-order formula, then the set
[TABLE]
is semialgebraic.
Proof of Lemma 3.1.
We have
[TABLE]
It follows from above that the set can be determined by a first-order formula, hence by the Tarski-Seidenberg Theorem, it is a semialgebraic subset of . ∎
The following proposition is contained implicitly in [Ha, Proof of Theorem B].
Proposition 3.1**.**
.
Proof.
Put . By the metric induced from that of , is a complete metric space and the function is bounded from below. Let and be a sequence of the first type of :
[TABLE]
Let . Then and as . Set . By the Ekeland’s Variational Principle ([E]), there exists a sequence such that
[TABLE]
and for any , we have
[TABLE]
Since and , the ball is contained in . Then, inequality (13) implies that
[TABLE]
holds true for every and . This gives us
[TABLE]
Putting , we get .
Clearly . Therefore . ∎
Proposition 3.2**.**
If there is a sequence of the second type of :
[TABLE]
then there exists a sequence of the second type of :
[TABLE]
with additional conditions
[TABLE]
In particular, the segment contains at least one point of F(f).
Proof.
Put and .
As in the proof of Proposition 3.1, we can find a sequence such that
[TABLE]
Since
[TABLE]
we have . The proposition is proved. ∎
Proposition 3.3**.**
If and , then .
Proof.
Assume that . By contradiction, suppose that . Hence, either or is a non-empty finite set.
By definition of , has a sequence of second type. Hence, it follows from Proposition 3.2, . Thus, is a non-empty finite set. Then, for any sufficiently small, we have and .
Let be a sequence of the second type of :
[TABLE]
By Proposition 3.2, we may assume that and there exists and .
Let such that and . Since , we can assume that for all . Let be the point of such that . Clearly, .
Claim: is a sequence of second type of .
Proof of Claim.
Since . Hence, for some , we have for all .
Let be the point of such that . We have
[TABLE]
This shows that and the claim is proved. ∎
Since is a sequence of the second type of and , by Proposition 3.2, there exists . Choose such that and . Similarly as in the proof of Claim, we can find a sequence of the second type of such that and such that .
Making this process iteratively, we see that the interval contains a infinite number of points in , which is a contradiction. ∎
4. Types of stability of global Hölderian error bounds
We will distinguish 3 cases.
4.1. Case 1 -
Theorem 4.1**.**
If then or .
Proof.
Assume that . Then by Proposition 3.1, . Moreover, it follows from Proposition 3.2 that is also empty.
Hence, by Theorem 2.1, or . ∎
Definition 4.1**.**
Let .
is called y-stable if and there exists an open interval such that ; 2. 2.
is called y-right stable if and there exists such that and .
Corollary 4.1**.**
If , then we have two cases
If , then there is only one type of stability of GHEB. Namely, for all , is y-stable. 2. 2.
If , then then there are two stability types of GHEB. Namely, for all , is y-stable and for , is y-right stable.
Remark 4.1**.**
We recall here results of [Ha] about the role that Newton polyhedron plays in studying GHEB’s.
Let be a polynomial in variables. Put and denote the convex hull in of the set . Following [Kou] we call the Newton polyhedron at infinity of .
Let be a face (of any dimension) of , set:
[TABLE]
Definition** ([Kou]).**
We say that a polynomial is nondegenerate with respect to its Newton boundary at infinity (nondegenerate for short), if for every face of not containing the origin, the system
[TABLE]
has no solution in .**
Definition**.**
A polynomial in variables is said to be convenient if for every , there exists a monomial of of the form , with a non-zero coefficient.**
Theorem 4.2** ([Ha]).**
If is convenient and nondegenerate w.r.t. its Newton polyhedron at infinity, then there exist such that
[TABLE]
In particular, .
Let denote the ring of polynomials in variables over .
For , as before, denotes the Newton polyhedron at infinity of . Let be a convenient polynomial.
Put and
[TABLE]
The set can be identified to the space , where is the number of integer points of .
Put \mathcal{B}_{\Gamma}=\{h\in\mathcal{A}_{\Gamma}:\Gamma_{h}=\Gamma\ \text{and h is nondegenerate}\}. According to [Kou], is an open and dense subset of . Hence, Theorem 4.1 and 4.2 show that if is a generic polynomial, then or . By Corollary 4.1, any value is y-stable and is y-right stable where
4.2. Case 2 - is non-empty finite set
Proposition 4.1**.**
If , then .
Proof.
By contradiction, assume that . Since , we have (Proposition 3.1). Then, it follows from the first formula that if and only if but the later is impossible, since we have
Claim: If , then .
Proof of Claim.
Take , since , has a sequence of the second type. By Proposition 3.2, there exists and . Take such that , then has a sequence of the second type. Hence, there exists and such that . Continuing this way, we find an infinite sequence of . Therefore, . ∎
∎
Now, we classify the stability types of GHEB in the case when is a non-empty finite set.
Definition 4.2**.**
Let .
Recall that is called y-stable if and there exists an open interval such that ; 2. 2.
Recall that is called y-right stable if and there exists such that and ; 3. 1’.
is called n-stable if and there exists an open interval such that ; 4. 2’.
is called n-right stable if and there exists such that and ; 5. 3’.
is called n-left stable if and there exists such that and ; 6. 4’.
is called n-isolated if and for sufficiently small, .
It follows from the first formula that
Theorem 4.3**.**
Let be a non-empty finite set and . Then, is one of the following types
**Case A: **
If , then
- **(i): **
* is y-stable if and only if .* 2. **(ii): **
* is a n-isolated point if and only if .*
**Case B: **
If is a finite value, then
- **1.: **
* is y-stable if and only if and ;* 2. **2.: **
* is y-right stable if and only if and ;* 3. **1’.: **
* is n-stable if and only if ;* 4. **2’.: **
* is n-right stable if and only if and ;* 5. **3’.: **
* is n-left stable if and only if and ;* 6. **4’.: **
* is a n-isolated point if and only if and .*
Remark 4.2**.**
Here, if we have item 2, then we does not have item 3’ and vice versa.
Now, to complete this subsection, we add an estimation of the number of connected components of for the case .
Let us denote the number of connected components of , we have the following result
Theorem 4.4**.**
Let be an any polynomial of degree . Then, if , we have
[TABLE]
Proof.
Put
[TABLE]
Since , we have . Then, according to Theorem 1.1 of [Je], we have
[TABLE]
Hence, it follows from the first formula that
[TABLE]
∎
4.3. Case 3 - is an infinite set
In this case, the following lemma tells us that the set has still very simple structure
Lemma 4.1**.**
* is a semialgebraic subset of .*
Using the first formula for (Theorem 2.2), it is enough to show that is semialgebraic.
Proof of Lemma 4.1.
We have
[TABLE]
[TABLE]
It follows from (a) and (b) that the set can be determined by a first-order formula, hence it is a semialgebraic subset of . ∎
Since is a semialgebraic subset of , we have
Corollary 4.2**.**
If and , then it is a union of finitely many points and intervals.
By Corollary 4.2, we have to consider three cases
- (a)
; 2. (b)
; 3. (c)
is a non-empty proper semialgebraic subset of .
- •
In the case (a), we have only one stable type: is y-stable for all ;
- •
In the case (b), we have only one stable type: is n-stable for all ;
- •
In the case (c), is a disjoint union of the sets of the following types:
[TABLE]
Where
- (1)
or ; 2. (2)
or ; 3. (3)
or ; 4. (4)
or ; 5. (5)
or , where are isolated points; 6. (6)
or or , where ; 7. (7)
or or , where .
Similarly, is a disjoint union of the sets of the following types:
[TABLE]
We have the following definition
Definition 4.3**.**
Let .
Recall that is said to be y-stable if and there exists an open interval such that ; 2. 2.
Recall that is said to be y-right stable if and there exists such that and ; 3. 3.
is said to be y-left stable if and there exists such that and ; 4. 4.
is said to be y-isolated if and for sufficiently small, ; 5. 1’.
Recall that is called n-stable if and there exists an open interval such that ; 6. 2’.
Recall that is called n-right stable if and there exists such that and ; 7. 3’.
Recall that is called n-left stable if and there exists such that and ; 8. 4’.
Recall that is called n-isolated if and for sufficiently small, .
Using the first formula of , we have
Theorem 4.5**.**
Let be of the form (c) and . Then we have
* is y-stable if and only if is an interior point of the sets*
[TABLE] 2. 2.
* is y-right stable if and only if we have or or (where );* 3. 3.
* is y-left stable if and only if we have or or (where );* 4. 4.
* is an y-isolated point if and only if .* 5. 1’.
* is n-stable if and only if is an interior point of the set:*
[TABLE] 6. 2’.
* is n-right stable if and only if we have or or (where );* 7. 3’.
* is n-left stable if and only if we have or or (where );* 8. 4’.
* is an n-isolated point if and only if .*
Remark 4.3**.**
In the above list, we collect all types of stability that could theoretically exist. The problem of deciding when this or that type really appears, seems to be very difficult.
We finish our paper by considering the following simple example
Example 4.1**.**
Let ([HT]). Clearly, is of the form .
We have . Hence, the roots of are:
[TABLE]
We have
[TABLE]
Hence .
It is not difficult to show that
- •
, hence ;
- •
and .
Therefore, by the second formula, . In this example, for any :
- •
If , then is y-stable;
- •
If , then is y-right stable;
- •
If , then is n-stable;
- •
If , then is n-right stable.
Acknowledgments
This research was partially supported by National Foundation for Science and Technology Development (NAFOSTED), Vietnam; Grant numbers 101.04-2017.12 of the first author and 101.04-2019.302 of the second author.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[AC] A. Auslender and Crouzeix, Global regularity theorem , Math. Oper. Res., 13 (1988), 243-253.
- 2[Az] D. Azé, A survey on error bounds for lower semicontinuous functions , Proceedings of 2003 MODE-SMAI Conference of ESAIM Proceedings, EDP Sci., Les Ulis, vol. 13, (2003), 1-17.
- 3[BCR] J. Bochnak, M. Coste and M. F. Roy, Real algebraic geometry , Springer, 1998.
- 4[BNPS] J. Bolte, T. P. Nguyen, J. Peypouquet and B. W. Suter, From error bounds to the complexity of first-order descent methods for convex functions , Math. Program., Ser. A, vol. 165 (2017), 2, 471-507.
- 5[CM] J. N. Corvellec, V. V. Montreanu, Nonlinear error bounds for lower semi-continuous functions on metric spaces , Math. Progam., Ser. A, vol. 114 (2008), 2, 291-319.
- 6[C] M. Coste, An introduction to Semialgebraic Geometry , Dottorato di ricerca in matematica/ Universita di Pisa, Dipartimento di Matematica, 2002.
- 7[Da] J. W. Daniel, On perturbations in systems of linear inequalities , SIAM J. Numer. Anal., 10 (1973), pp. 299-307.
- 8[D] S. Deng, Perturbation analysis of a condition number for convex inequality systems and global error bounds for analytic systems , Math. Program., vol. 83 (1998), 263-276.
