Next Article in Journal
Trees in Positive Entropy Subshifts
Next Article in Special Issue
Mathematical Approach for System Repair Rate Analysis Used in Maintenance Decision Making
Previous Article in Journal
A Comparison of a Priori Estimates of the Solutions of a Linear Fractional System with Distributed Delays and Application to the Stability Analysis
Previous Article in Special Issue
Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximations of an Equilibrium Problem without Prior Knowledge of Lipschitz Constants in Hilbert Spaces with Applications

by
Chainarong Khanpanuk
1,
Nuttapol Pakkaranang
2,
Nopparat Wairojjana
3,* and
Nattawut Pholasa
4,*
1
Department of Mathematics, Faculty of Science and Technology, Phetchabun Rajabhat University, Phetchabun 67000, Thailand
2
Department of Mathematics, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), Bangkok 10140, Thailand
3
Applied Mathematics Program, Faculty of Science and Technology, Valaya Alongkorn Rajabhat University under the Royal Patronage (VRU), 1 Moo 20 Phaholyothin Road, Klong Neung, Klong Luang, Pathumthani 13180, Thailand
4
School of Science, University of Phayao, Phayao 56000, Thailand
*
Authors to whom correspondence should be addressed.
Axioms 2021, 10(2), 76; https://doi.org/10.3390/axioms10020076
Submission received: 21 March 2021 / Revised: 20 April 2021 / Accepted: 23 April 2021 / Published: 27 April 2021
(This article belongs to the Special Issue Numerical Analysis and Computational Mathematics)

Abstract

:
The objective of this paper is to introduce an iterative method with the addition of an inertial term to solve equilibrium problems in a real Hilbert space. The proposed iterative scheme is based on the Mann-type iterative scheme and the extragradient method. By imposing certain mild conditions on a bifunction, the corresponding theorem of strong convergence in real Hilbert space is well-established. The proposed method has the advantage of requiring no knowledge of Lipschitz-type constants. The applications of our results to solve particular classes of equilibrium problems is presented. Numerical results are established to validate the proposed method’s efficiency and to compare it to other methods in the literature.

1. Introduction

Suppose that C is a nonempty closed and convex subset of a real Hilbert space H . The inner product and induced norm are denoted by · , · and · , respectively. Let f : H × H R be a bifunction and f ( y , y ) = 0 , for all y C . The equilibrium problem (EP) [1,2] for a bifunction f on C is defined in the following way:
Find u * C such that f ( u * , y ) 0 , y C .
The equilibrium problem is a general mathematical problem in the sense that it unifies various mathematical problems, i.e., fixed-point problems, vector and scalar minimization problems, problems of variational inequality, complementarity problems, Nash equilibrium problems in noncooperative games, saddle point problems, and inverse optimization problems [2,3,4]. The equilibrium problem is also known as the well-known Ky Fan inequality due to the result [1]. Many authors established and generalized several results on the existence and nature of the solution of the equilibrium problems (see for more detail [1,4,5]). Due to the importance of this problem (EP) in both pure and applied sciences, many researchers studied it in recent years [6,7,8,9,10,11,12,13,14,15,16,17] and other in [18,19,20,21,22].
Tran et al. in [23] introduced iterative sequence { u n } in the following way:
u 0 C , y n = arg min z C { χ f ( u n , z ) + 1 2 u n z 2 } , u n + 1 = arg min z C { χ f ( y n , z ) + 1 2 u n z 2 } ,
where 0 < χ < min 1 2 c 1 , 1 2 c 2 . This method is also known as the extragradient method in [23] due to the previous contribution of Korpelevich [24] to solve the saddle-point problems. The iterative sequence generated by the above-mentioned method is weakly convergent to the solution with prior knowledge of Lipschitz-type constants. These Lipschitz-like constants are often not known or are difficult to compute. Recently, Hieu et al. [25] introduced an extension of the method (1) for solving the equilibrium problem. Let us consider that [ p ] + : = max { p , 0 } and choose u 0 C , μ ( 0 , 1 ) with χ 0 > 0 such that
y n = arg min z C { χ n f ( u n , z ) + 1 2 u n z 2 } , u n + 1 = arg min z C { χ n f ( y n , z ) + 1 2 u n z 2 } ,
where { χ n } is updated in the following manner:
χ n + 1 = min χ n , μ ( u n y n 2 + u n + 1 y n 2 ) 2 [ f ( u n , u n + 1 ) f ( u n , y n ) f ( y n , u n + 1 ) ] + .
Inertial-like methods are well-known two-step iterative methods in which the next iteration is derived from the previous two iterations (see [26,27] for more details). To speed up the iterative sequence convergence rate, an inertial extrapolation term is used. Numerical examples show that inertial effects improve numerical performance in terms of execution time and the expected number of iterations. Recently, many existing methods were established for the case of equilibrium problems (see [28,29,30,31] for more details).
In this paper, inspired by the methods in [23,25,26,32], we introduce a general inertial Mann-type subgradient extragradient method to evaluate the approximate solution of the equilibrium problems involving pseudomonotone bifunction. A strong convergence result corresponding to the proposed algorithm is well-established by assuming certain mild conditions. Some of the applications for our main results are considered to solve the fixed-point problems. Lastly, computational results show that the new method is more successful than existing ones [23,33,34].

2. Preliminaries

A metric projection P C ( u ) of u H onto a closed and convex subset C of H is defined by
P C ( u ) = arg min y C { y u } .
In this study, the equilibrium problem under the following conditions:
(c1).
A bifunction f : H × H R is said to be pseudomonotone [3,35] on C if
f ( y 1 , y 2 ) 0 f ( y 2 , y 1 ) 0 , y 1 , y 2 C .
(c2).
A bifunction f : H × H R is said to be Lipschitz-type continuous [36] on C if there exist constants c 1 , c 2 > 0 such that
f ( y 1 , y 3 ) f ( y 1 , y 2 ) + f ( y 2 , y 3 ) + c 1 y 1 y 2 2 + c 2 y 2 y 3 2 , y 1 , y 2 , y 3 C .
(c3).
lim sup n f ( y n , y ) f ( q * , y ) for all y C and { y n } C satisfy y n q * .
(c4).
f ( u , · ) is convex and subdifferentiable on H for each u H .
A cone on C at u C is defined by
N C ( u ) = { t H : t , y u 0 , y C } .
Let a convex function : C R and subdifferential of ℸ at u C is defined by
( u ) = { t H : ( y ) ( u ) t , y u , y C } .
Lemma 1.
 [37] Let : C R be a subdifferentiable, lower semicontinuous, and convex function on C . Then, u C is said to be a minimizer of ℸ if and only if 0 ( u ) + N C ( u ) , where ( u ) stands for the subdifferential of ℸ at u C and N C ( u ) is a normal cone of C on u .
Lemma 2.
 [38] Assume that P C : H C be a metric projection such that
(i)
y 1 P C ( y 2 ) 2 + P C ( y 2 ) y 2 2 y 2 y 1 2 , y 1 C , y 2 H .  
(ii)
y 3 = P C ( y 1 ) if and only if y 1 y 3 , y 2 y 3 0 , y 2 C .  
(iii)
y 1 P C ( y 1 ) y 1 y 2 , y 2 C , y 1 H .
Lemma 3.
 [39] Assume that { n } ( 0 , + ) is a sequence satisfying, i.e., n + 1 ( 1 n ) n + n ð n , for all n N . Moreover, let { n } ( 0 , 1 ) and { ð n } R be two sequences, such that lim n n = 0 , n = 1 n = + a n d lim sup n ð n 0 . Then, lim n n = 0 .
Lemma 4.
 [40] Assume that { n } be a sequence of real numbers such that there exists a subsequence { n i } of { n } such that n i < n i + 1 for all i N . Then, there is a nondecreasing sequence m k N such that m k as k , and the following conditions are fullfiled by all (sufficiently large) numbers k N :
m k m k + 1 and k m k + 1 .
In fact, m k = max { j k : j j + 1 } .
Lemma 5.
[41] For all y 1 , y 2 H and ð R , the following inequalities hold.
(i)
ð y 1 + ( 1 ð ) y 2 2 = ð y 1 2 + ( 1 ð ) y 2 2 ð ( 1 ð ) y 1 y 2 2 .  
(ii)
y 1 + y 2 2 y 1 2 + 2 y 2 , y 1 + y 2 .

3. Main Results

We propose an iterative method for solving equilibrium problems involving a pseudomonotone that is based on Tran et al. in [23], and the Mann-type method [32] and the inertial scheme [26]. For clarity in the presentation, we use notation [ t ] + = max { 0 , t } and follow conventions 0 0 = + and a 0 = + ( a 0 ) .
Lemma 6.
A sequence { χ n } generated by (5) is monotonically decreasing, converges to χ > 0 , and has a lower bound min μ 2 max { c 1 , c 2 } , χ 0 .
Proof. 
Assume that f ( t n , z n ) f ( t n , y n ) f ( y n , z n ) > 0 such that
μ ( t n y n 2 + z n y n 2 ) 2 [ f ( t n , z n ) f ( t n , y n ) f ( y n , z n ) ] μ ( t n y n 2 + z n y n 2 ) 2 [ c 1 t n y n 2 + c 2 z n y n 2 ] μ 2 max { c 1 , c 2 } .
This implies that { χ n } has a lower bound min μ 2 max { c 1 , c 2 } , χ 0 . Moreover, there exists a fixed real number χ > 0 , such that lim n χ n = χ .  ☐
Lemma 7.
Suppose that Conditions (c1)(c4) are satisfied. Then, sequence { u n } generated by the Algorithm 1 is a bounded sequence.
Algorithm 1 (Explicit Accelerated Strong Convergence Iterative Scheme)
  • Step 0: Choose u 1 , u 0 C , ϕ > 0 , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) satisfies the following conditions:
    lim n ϱ n = 0 a n d n = 1 + ϱ n = + .
  • Step 1: Compute t n = u n + ϕ n ( u n u n 1 ) and choose ϕ n such that
    0 ϕ n ϕ n ^ and ϕ n ^ = min ϕ 2 , ς n u n u n 1 if u n u n 1 , ϕ 2 ϕ 2 otherwise ,
    where ς n = ( ϱ n ) , i.e., lim n ς n ϱ n = 0 .
  • Step 2: Compute
    y n = arg min y C { χ n f ( t n , y ) + 1 2 t n y 2 } .
    If t n = y n , then STOP the sequence. Else, go to Step 3.
  • Step 3: Construct a half-space H n = { z H : t n χ n ω n y n , z y n 0 } where ω n 2 f ( t n , y n ) and compute
    z n = arg min y H n { χ n f ( y n , y ) + 1 2 t n y 2 } .
  • Step 4: Compute u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n .
  • Step 5: Compute
    χ n + 1 = min χ n , μ t n y n 2 + μ z n y n 2 2 [ f ( t n , z n ) f ( t n , y n ) f ( y n , z n ) ] + .
    Set n : = n + 1 and go back to Step 1.
Proof. 
From the value of z n , we have
0 2 χ n f ( y n , y ) + 1 2 t n y 2 ( z n ) + N H n ( z n ) .
For ω f ( y n , z n ) there exists ω ¯ N H n ( z n ) such that
χ n ω + z n t n + ω ¯ = 0 .
This implies that
t n z n , y z n = χ n ω , y z n + ω ¯ , y z n , y H n .
Due to ω ¯ N H n ( z n ) , it implies that ω ¯ , y z n 0 for each y H n . Thus, we have
t n z n , y z n χ n ω , y z n , y H n .
Moreover, ω f ( y n , z n ) and owing to the subdifferential, we have
f ( y n , y ) f ( y n , z n ) ω , y z n , y H .
From Expressions (6) and (7), we obtain
χ n f ( y n , y ) χ n f ( y n , z n ) t n z n , y z n , y H n .
Due to the definition of H n , we have
χ n ω n , z n y n t n y n , z n y n .
Now, using ω n f ( t n , y n ) , we obtain
f ( t n , y ) f ( t n , y n ) ω n , y y n , y H .
By letting y = z n , we obtain
f ( t n , z n ) f ( t n , y n ) ω n , z n y n , y H .
Combining Expressions (9) and (10), we obtain
χ n f ( t n , z n ) f ( t n , y n ) t n y n , z n y n .
By substituting y = u * in Expression (8), we obtain
χ n f ( y n , u * ) χ n f ( y n , z n ) t n z n , u * z n .
Since u * E p ( f , C ) , we have f ( u * , y n ) 0 . From the pseudomonotonicity of bifunction f, we achieve f ( y n , u * ) 0 . It follows from Expression (12) that
t n z n , z n u * χ n f ( y n , z n ) .
From the description of χ n + 1 , we obtain
f ( t n , z n ) f ( t n , y n ) f ( y n , z n ) μ t n y n 2 + μ z n y n 2 2 χ n + 1
From (13) and (14), we obtain
t n z n , z n u * χ n { f ( t n , z n ) f ( t n , y n ) } μ χ n 2 χ n + 1 t n y n 2 μ χ n 2 χ n + 1 z n y n 2 .
Combining Expressions (11) and (15), we have
t n z n , z n u * t n y n , z n y n μ χ n 2 χ n + 1 t n y n 2 μ χ n 2 χ n + 1 z n y n 2 .
We have the given formula in place:
2 t n z n , z n u * = t n u * 2 + z n t n 2 + z n u * 2 .
2 y n t n , y n z n = t n y n 2 + z n y n 2 t n z n 2 .
Combining (16)–(18), we obtain
z n u * 2 t n u * 2 1 μ χ n χ n + 1 t n y n 2 1 μ χ n χ n + 1 z n y n 2 .
Since χ n χ , then there is number ( 0 , 1 μ ) that
lim n 1 μ χ n χ n + 1 = 1 μ > > 0 .
Thus, there exists a finite number n 1 N , such that
1 μ χ n χ n + 1 > > 0 , n n 1 .
From Expression (19), we obtain
u n + 1 u * 2 t n u * 2 , n n 1 .
From Expression (4), we have ϕ n u n u n 1 ς n , for all n N and lim n ς n ϱ n = 0 implies that
lim n ϕ n ϱ n u n u n 1 lim n ς n ϱ n = 0 .
From Expression (21) and { t n } , we have
z n u * t n u * = u n + ϕ n ( u n u n 1 ) u * u n u * + ϕ n u n u n 1 u n u * + ϱ n ϕ n ϱ n u n u n 1 u n u * + ϱ n 1 ,
where for some fixed 1 > 0 and
ϕ n ϱ n u n u n 1 1 , n 1 .
It is given that u * E p ( f , C ) and by definition of { u n + 1 } , we have
u n + 1 u * = ( 1 ρ n ϱ n ) u n + ρ n z n u * = ( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) ϱ n u * ( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) + ϱ n u * .
Next, we compute
( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) 2 = ( 1 ρ n ϱ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + 2 ( 1 ρ n ϱ n ) ( u n u * ) , ρ n ( z n u * ) ( 1 ρ n ϱ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + 2 ρ n ( 1 ρ n ϱ n ) u n u * z n u * ( 1 ρ n ϱ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + ρ n ( 1 ρ n ϱ n ) u n u * 2 + ρ n ( 1 ρ n ϱ n ) z n u * 2
( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ρ n ( 1 ϱ n ) z n u * 2
( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ρ n ( 1 ϱ n ) ( u n u * + ϱ n 1 ) 2
( 1 ϱ n ) 2 u n u * 2 + ϱ n 2 1 2 + 2 ϱ n 1 ( 1 ϱ n ) u n u * 2 .
The above expression implies that
( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) ( 1 ϱ n ) u n u * + ϱ n 1 .
Combining Expressions (25) and (28), we obtain
u n + 1 u * ( 1 ϱ n ) u n u * + ϱ n 1 + ϱ n u * max u n u * , 1 + u * } max u 0 u * , 1 + u * } .
Therefore, we conclude that { u n } is bounded sequence. ☐
Theorem 1.
Let { u n } be a sequence generated by Algorithm 1, and Conditions (c1)(c4) are satisfied. Then, { u n } strongly converges to u * = P E p ( f , C ) ( 0 ) .
Proof. 
By using definition of { u n + 1 } , we have
u n + 1 u * 2 = ( 1 ρ n ϱ n ) u n + ρ n z n u * 2 = ( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) ϱ n u * 2 = ( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) 2 + ϱ n 2 u * 2 2 ( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) , ϱ n u * .
From Expression (26), we have
( 1 ρ n ϱ n ) ( u n u * ) + ρ n ( z n u * ) 2 ( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ρ n ( 1 ϱ n ) z n u * 2 .
Combining Expressions (30) and (31) (for some 2 > 0 ), we obtain
u n + 1 u * 2 ( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ρ n ( 1 ϱ n ) z n u * 2 + ϱ n 2 ( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ϱ n 2 + ρ n ( 1 ϱ n ) t n u * 2 1 μ χ n χ n + 1 t n y n 2 1 μ χ n χ n + 1 z n y n 2 .
From Expression (23), we have
t n u * 2 u n u * 2 + ϱ n 3 ,
for some 3 > 0 . Substituting (33) into (32), we obtain
u n + 1 u * 2 ( 1 ρ n ϱ n ) ( 1 ϱ n ) u n u * 2 + ϱ n 2 + ρ n ( 1 ϱ n ) u n u * 2 + ϱ n 3 1 μ χ n χ n + 1 t n y n 2 1 μ χ n χ n + 1 z n y n 2 = ( 1 ϱ n ) 2 u n u * 2 + ϱ n 2 + ρ n ( 1 ϱ n ) ϱ n 3 ρ n ( 1 ϱ n ) 1 μ χ n χ n + 1 t n y n 2 + 1 μ χ n χ n + 1 z n y n 2 u n u * 2 + ϱ n 4 ρ n ( 1 ϱ n ) 1 μ χ n χ n + 1 t n y n 2 + 1 μ χ n χ n + 1 z n y n 2 ,
for some 4 > 0 . It is given that u * = P E p ( f , C ) ( 0 ) and by using Lemma 2 (ii) ( E p ( f , C ) is a convex and closed set ([23,34])), we obtain
u * , u * y 0 , y E p ( f , C ) .
The remainder of the proof shall be taken into account in the following two parts:
Case 1: Assume that there is a fixed number n 2 N ( n 2 n 1 ) such as
u n + 1 u * u n u * , n n 2 .
It implies that lim n u n u * exists, and due to (34), we obtain
ρ n ( 1 ϱ n ) 1 μ χ n χ n + 1 t n y n 2 + 1 μ χ n χ n + 1 z n y n 2 u n u * 2 + ϱ n 4 u n + 1 u * 2 .
Due to the existence of lim n u n u * , ϱ n 0 and χ n χ , we infer that
lim n y n t n = lim n y n z n = 0 .
We can calculate that
lim n z n t n lim n t n y n + lim n y n z n = 0 .
It follows that
u n + 1 u n = ( 1 ρ n ϱ n ) u n + ρ n z n u n = u n ϱ n u n + ρ n z n ρ n u n u n ρ n z n u n + ϱ n u n .
The term is referred to above that
lim n u n + 1 u n = 0 .
Thus, this implies that { y n } and { z n } are bounded. The reflexivity of H and the boundedness of { u n } guarantee that there is a subsequence { u n k } , such that { u n k } x ^ H as k . Next, our aim to prove that x ^ E p ( f , C ) . Using (8), due to χ n + 1 and (11), we write
χ n k f ( y n k , y ) χ n k f ( y n k , z n k ) + t n k z n k , y z n k χ n k f ( t n k , z n k ) χ n k f ( t n k , y n k ) μ χ n k 2 χ n k + 1 t n k y n k 2 μ χ n k 2 χ n k + 1 y n k z n k 2 + t n k z n k , y z n k t n k y n k , z n k y n k μ χ n k 2 χ n k + 1 t n k y n k 2 μ χ n k 2 χ n k + 1 y n k z n k 2 + t n k z n k , y z n k ,
while y is an any arbitrary member in H n . It continues from (38) and (39) that the right-hand side approaches to zero. From χ > 0 , Condition (c3) and y n k x ^ , we have
0 lim sup k f ( y n k , y ) f ( x ^ , y ) , y H n .
The following is that f ( x ^ , y ) 0 , y C ; thus x ^ E p ( f , C ) . It continues from that
lim sup n u * , u * u n = lim sup k u * , u * u n k = u * , u * x ^ 0 .
Due to lim n u n + 1 u n = 0 , we can deduce that
lim sup n u * , u * u n + 1 lim sup k u * , u * u n + lim sup k u * , u n u n + 1 0 .
Next, consider the following value
t n u * 2 = u n + ϕ n ( u n u n 1 ) u * 2 = u n u * + ϕ n ( u n u n 1 ) 2 = u n u * 2 + ϕ n 2 u n u n 1 2 + 2 u n u * , ϕ n ( u n u n 1 ) u n u * 2 + ϕ n 2 u n u n 1 2 + 2 ϕ n u n u * u n u n 1 = u n u * 2 + ϕ n u n u n 1 [ 2 u n u * + ϕ n u n u n 1 u n u * 2 + ϕ n u n u n 1 5 ,
Substituting q n = ( 1 ρ n ) u n + ρ n z n , we have
u n + 1 = q n ϱ n u n = ( 1 ϱ n ) q n ϱ n ( u n q n ) = ( 1 ϱ n ) q n ϱ n ρ n ( u n z n ) .
where u n q n = u n ( 1 ρ n ) u n ρ n z n = ρ n ( u n z n ) . Consider that
q n u * 2 = ( 1 ρ n ) u n + ρ n z n u * 2 = ( 1 ρ n ) ( u n u * ) + ρ n ( z n u * ) 2 = ( 1 ρ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + 2 ( 1 ρ n ) ( u n u * ) , ρ n ( z n u * ) ( 1 ρ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + 2 ρ n ( 1 ρ n ) u n u * z n u * ( 1 ρ n ) 2 u n u * 2 + ρ n 2 z n u * 2 + ρ n ( 1 ρ n ) u n u * 2 + ρ n ( 1 ρ n ) z n u * 2 = ( 1 ρ n ) u n u * 2 + ρ n z n u * 2 ( 1 ρ n ) u n u * 2 + ρ n t n u * 2 ( 1 ρ n ) u n u * 2 + ρ n [ u n u * 2 + ϕ n u n u n 1 5 u n u * 2 + ϕ n u n u n 1 5 .
Next, consider that
u n + 1 u * 2 = ( 1 ϱ n ) q n + ρ n ϱ n ( z n u n ) u * 2 = ( 1 ϱ n ) ( q n u * ) + ρ n ϱ n ( z n u n ) ϱ n u * 2 ( 1 ϱ n ) 2 q n u * 2 + 2 ρ n ϱ n ( z n u n ) ϱ n u * , ( 1 ϱ n ) ( q n u * ) + ρ n ϱ n ( z n u n ) ϱ n u * = ( 1 ϱ n ) 2 q n u * 2 + 2 ρ n ϱ n ( z n u n ) ϱ n u * , q n ϱ n q n ϱ n ( u n q n ) u * = ( 1 ϱ n ) q n u * 2 + 2 ρ n ϱ n z n u n , u n + 1 u * + 2 ϱ n u * , u * u n + 1 ( 1 ϱ n ) q n u * 2 + 2 ρ n ϱ n z n u n u n + 1 u * + 2 ϱ n u * , u * u n + 1
for some 5 > 0 . Combining Expressions (46), (48), and (49), we obtain
u n + 1 u * 2 ( 1 ϱ n ) u n u * 2 + ( 1 ϱ n ) ϕ n u n u n 1 5 + 2 ρ n ϱ n z n u n u n + 1 u * + 2 ϱ n u * , u * u n + 1 ( 1 ϱ n ) u n u * 2 + ϱ n [ ϕ n ϱ n ( 1 ϱ n ) u n u n 1 5 + 2 ρ n z n u n u n + 1 u * + 2 u * , u * u n + 1 .
Due to (45), (50), and the implemented Lemma 3, we conclude that u n u * 0 as n .
Case 2: Assume there is a subsequence { n i } of { n } that
u n i u * u n i + 1 u * , i N .
Using Lemma 4, there is a { m k } N sequence, such as { m k } ,
u m k u * u m k + 1 u * and u k u * u m k + 1 u * , for all k N .
Similar to Case 1, Relation (37) gives that
ρ m k ( 1 ϱ m k ) 1 μ χ m k χ m k + 1 t m k y m k 2 + 1 μ χ m k χ m k + 1 z m k y m k 2 u m k u * 2 + ϱ m k 4 u m k + 1 u * 2 .
Due to ϱ m k 0 and χ m k χ , we deduce the following:
lim n t m k y m k = lim n z m k y m k = 0 .
It continues on from that
u m k + 1 u m k = ( 1 ρ m k ϱ m k ) u m k + ρ m k z m k u m k = u m k ϱ m k u m k + ρ m k z m k ρ m k u m k u m k ρ m k z m k u m k + ϱ m k u m k 0 .
We use the same reasoning as that in Case 1:
lim sup k u * , u * u m k + 1 0 .
Now, using Expressions (50) and (51), we have
u m k + 1 u * 2 ( 1 ϱ m k ) u m k u * 2 + ϱ m k [ ϕ m k ϱ m k ( 1 ϱ m k ) u m k u m k 1 5 + 2 ρ m k z m k u m k u m k + 1 u * + 2 u * , u * u m k + 1 . ( 1 ϱ m k ) u m k + 1 u * 2 + ϱ m k [ ϕ m k ϱ m k ( 1 ϱ m k ) u m k u m k 1 5 + 2 ρ m k z m k u m k u m k + 1 u * + 2 u * , u * u m k + 1 .
It implies that
u m k + 1 u * 2 [ ϕ m k ϱ m k ( 1 ϱ m k ) u m k u m k 1 5 + 2 ρ m k z m k u m k u m k + 1 u * + 2 u * , u * u m k + 1 .
Since ϱ m k 0 , and u m k u * is bounded. Thus, with Expressions (55) and (57), we have
u m k + 1 u * 2 0 , as k .
The above implies that
lim n u k u * 2 lim n u m k + 1 u * 2 0 .
As a result, u n u * . This completes the proof of the theorem. ☐
By letting ϕ n = 0 , we obtain a strong convergence of the result in [25].
Corollary 1.
Let f : C × C R be a bifunction satisfying Conditions (c1)(c4). Choosing u 0 C , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) satisfies the following conditions:
lim n ϱ n = 0 and n ϱ n = + .
Let { u n } be a sequence that is generated in the following manner:
y n = arg min y C { χ n f ( u n , y ) + 1 2 u n y 2 } , z n = arg min y H n { χ n f ( y n , y ) + 1 2 u n y 2 } , u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n ,
where H n = { z H : u n χ n ω n y n , z y n 0 } and ω n 2 f ( u n , y n ) . The step size is updated in the following way:
χ n + 1 = min χ n , μ u n y n 2 + μ z n y n 2 2 [ f ( u n , z n ) f ( u n , y n ) f ( y n , z n ) ] + .
Then, sequence { u n } converges strongly to u * E p ( f , C ) .

4. Applications to Solve Fixed-Point Problems

We propose our results to focus on fixed-point problems regarding κ -strict pseudocontraction mapping. The fixed-point problem (FPP) for S : H H is defined in the following manner:
Find u * C such that S ( u * ) = u * .
We assume that the following conditions were met:
(c1*)
A mapping S : C C is said to be κ-strict pseudocontraction [42] on C if
T y 1 T y 2 2 y 1 y 2 2 + κ ( y 1 T y 1 ) ( y 2 T y 2 ) 2 , y 1 , y 2 C ;
(c2*)
A mapping that is weakly sequentially continuous on C if
S ( y n ) S ( q * ) for any sequence in C satisfying y n q * .
If we consider that mapping S is weakly continuous and a κ -strict pseudocontraction, then f ( u , y ) = u S u , y u satisfies the conditions (c1)–(c4) (see [43]) and 2 c 1 = 2 c 2 = 3 2 κ 1 κ . The values of y n and z n in Algorithm 1 can be written as follows:
y n = arg min y C { χ n f ( t n , y ) + 1 2 t n y 2 } = P C t n χ n ( t n S ( t n ) ) , z n = arg min y H n { χ n f ( y n , y ) + 1 2 t n y 2 } = P H n t n χ n ( y n S ( y n ) ) .
Corollary 2.
Suppose C is a nonempty, convex, and closed subset of a Hilbert space H and S : C C is weakly continuous and κ-strict pseudocontraction with solution set F i x ( S ) . Let u 1 , u 0 C , ϕ > 0 , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) fulfill the items, i.e., lim n ϱ n = 0 a n d n = 1 ϱ n = + . Moreover, choose ϕ n satisfying 0 ϕ n ϕ n ^ such that
ϕ n ^ = min ϕ 2 , ς n u n u n 1 if u n u n 1 , ϕ 2 else ,
where ς n = ( ϱ n ) , i.e., lim n ς n ϱ n = 0 . Assume that { u n } is the sequence generated in the following manner:
t n = u n + ϕ n ( u n u n 1 ) , y n = P C t n χ n ( t n S ( t n ) ) , z n = P H n t n χ n ( y n S ( y n ) ) , u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n ,
where H n = { z H : ( 1 χ n ) t n + χ n S ( t n ) y n , z y n 0 } . Compute
χ n + 1 = min χ n , μ t n y n 2 + μ z n y n 2 2 ( t n y n ) [ T ( t n ) T ( y n ) ] , z n y n +
Then, { u n } strongly converges to u * F i x ( S , C ) .
Corollary 3.
Suppose C to be a convex and closed subset of a Hilbert space H and S : C C is weakly continuous and κ-strict pseudocontraction with solution set F i x ( S ) . Let u 0 C , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) fulfills the requirement, i.e., lim n ϱ n = 0 a n d n = 1 ϱ n = + . Assume that { u n } is the sequence formed as follows:
y n = P C u n χ n ( u n S ( u n ) ) , z n = P H n u n χ n ( y n S ( y n ) ) , u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n ,
where H n = { z H : ( 1 χ n ) u n + χ n S ( u n ) y n , z y n 0 } . Compute
χ n + 1 = min χ n , μ u n y n 2 + μ z n y n 2 2 ( u n y n ) [ T ( u n ) T ( y n ) ] , z n y n +
Then, sequence { u n } converges strongly to u * F i x ( S , C ) .

5. Applications to Solve Variational-Inequality Problems

Next, we consider the application of our results in the problem of classical variational inequalities [44,45]. The variational-inequality problem (VIP) for an operator L : H H is stated in the following manner:
Find u * C such that L ( u * ) , y u * 0 , y C .
We assume that the following conditions were met:
( L 1)
The solution set of problem (VIP) denoted by V I ( L , C ) is nonempty.
( L 2)
An operator L : H H is said to be pseudomonotone if
L ( y 1 ) , y 2 y 1 0 L ( y 2 ) , y 1 y 2 0 , y 1 , y 2 C .
( L 3)
An operator L : H H is said to be Lipschitz continuous through L > 0 , such that
L ( y 1 ) L ( y 2 ) L y 1 y 2 , y 1 , y 2 C ;
( L 4)
lim sup n L ( y n ) , y y n L ( q * ) , y q * for all y C and { y n } C satisfy y n q * .
If we define f ( u , y ) : = L ( u ) , y u for all u , y C . Then, problem (EP) becomes the problem of variational inequalities described above where L = 2 c 1 = 2 c 2 . From the above value of the bifunction f, we have
y n = arg min y C { χ n f ( t n , y ) + 1 2 t n y 2 } = P C ( t n χ n L ( t n ) ) , z n = arg min y H n { χ n f ( y n , y ) + 1 2 t n y 2 } = P H n ( t n χ n L ( y n ) ) .
Corollary 4.
Suppose that L : C H is a function satisfying the assumptions ( L 1)( L 4). Let u 1 , u 0 C , ϕ > 0 , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) satisfies the items, i.e., lim n ϱ n = 0 a n d n = 1 ϱ n = + . Moreover, choose ϕ n satisfying 0 ϕ n ϕ n ^ , such that
ϕ n ^ = min ϕ 2 , ς n u n u n 1 if u n u n 1 , ϕ 2 else ,
where ς n = ( ϱ n ) , i.e., lim n ς n ϱ n = 0 . Assume that { u n } is the sequence generated in the following manner:
t n = u n + ϕ n ( u n u n 1 ) , y n = P C ( t n χ n L ( t n ) ) , z n = P H n ( t n χ n L ( y n ) ) , u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n ,
where H n = { z H : t n χ n L ( t n ) y n , z y n 0 } . Compute
χ n + 1 = min χ n , μ t n y n 2 + μ z n y n 2 2 L ( t n ) L ( y n ) , z n y n + .
Then, sequences { u n } converge strongly to u * V I ( L , C ) .
Corollary 5.
Suppose that L : C H is a function meeting conditions ( L 1)( L 4). Let u 0 C , χ 0 > 0 , { ρ n } ( a , b ) ( 0 , 1 ϱ n ) and { ϱ n } ( 0 , 1 ) satisfies the conditions, i.e., lim n ϱ n = 0 a n d n = 1 ϱ n = + . Assume that { u n } is the sequence generated in the following manner:
y n = P C ( u n χ n L ( u n ) ) , z n = P H n ( u n χ n L ( y n ) ) , u n + 1 = ( 1 ρ n ϱ n ) u n + ρ n z n ,
where H n = { z H : u n χ n L ( u n ) y n , z y n 0 } .
Compute
χ n + 1 = min χ n , μ u n y n 2 + μ z n y n 2 2 L ( u n ) L ( y n ) , z n y n + .
Then, sequences { u n } converge strongly to u * V I ( L , C ) .
Remark 1.
Condition ( L 4) could be exempted when L is monotone. Indeed, this condition, which is a particular case of Condition (c3), is only used to prove (43). Without Condition ( L 4), inequality (42) can be obtained by imposing monotonocity on L . In that case,
L ( y ) , y y n L ( y n ) , y y n , y C .
By allowing f ( u , y ) = L ( u ) , y u in (42), we have
lim sup k L ( y n k ) , y y n k 0 , y H n .
Combining (65) with (66), we conclude that
lim sup k L ( y ) , y y n k 0 , y C .
Let y t = ( 1 t ) z + t y , for every t [ 0 , 1 ] . By using the convexity of set C , y t C for every t ( 0 , 1 ) . Since y n k z C and L ( y ) , y z 0 for every y C , we have
0 L ( y t ) , y t z = t L ( y t ) , y z .
Therefore, L ( y t ) , y z 0 , t ( 0 , 1 ) . Since y t z as t 0 and due to L continuity, we have L ( z ) , y z 0 , for each y C , which provides z V I ( L , C ) .
Remark 2.
From Remark 1, it can be concluded that Corollaries 4 and 5 still hold, even if we remove Condition ( L 4) in the case of monotone operators.

6. Numerical Illustrations

Numerical results are presented in this section to demonstrate the efficiency of our proposed method. The MATLAB codes were run in MATLAB version 9.5 (R2018b) on an Intel(R) Core(TM)i5-6200 CPU PC @ 2.30 GHz 2.40 GHz, RAM 4.00 GB.
Example 1.
Let there be m companies that manufacture the same product. Assume vector u of each item u i represents the quantity of the material produced by a company i. We consider that cost function P to be a declining affine function that relies on μ = i = 1 m u i , i.e., P i ( μ ) = ϕ i ψ i S , where ϕ i > 0 , ψ i > 0 . The formula for profit of every company i is taken as F i ( u ) = P i ( S ) u i q i ( u i ) , where q i ( u i ) is the tax value and cost for developing item u i . Moreover, consider that C i = [ u i min , u i max ] is the set of actions related to each company i , and the plan to figure out the model as C : = C 1 × C 2 × × C m . In addition, each member wants to achieve its peak turnover by a good level of production on the basis that the performance of other firms is an input parameter. The commonly used modelling methodology is based on the famous Nash equilibrium principle. A point u * C = C 1 × C 2 × × C m is the level of equilibrium of the model if
F i ( u * ) F i ( u * [ u i ] ) , u i C i , i = 1 , 2 , , m ,
wile u * [ u i ] is obtain from u * by letting ζ i * with u i . Furthermore, we consider f ( u , y ) : = Δ ( u , y ) Δ ( u , u ) while Δ ( u , y ) : = i = 1 m F i ( u [ y i ] ) . An equilibrium level of the model is defined by
F i n d u * C : f ( u * , z ) 0 , z C .
Bifunction f converts into the following form (see [23]):
f ( u , y ) = P u + Q y + c , y u
where c R m and P, Q matrices of order m . Matrix P is positive semidefinite, and matrix Q P is negative semidefinite with Lipschitz-type constants c 1 = c 2 = 1 2 P Q (see [23]) for details. P , Q are taken randomly. (Two diagonal matrices randomly A 1 and A 2 take elements from [ 0 , 2 ] and [ 2 , 0 ] respectively. Randomly O 1 = R a n d O r t h M a t ( m ) and O 2 = R a n d O r t h M a t ( m ) orthogonal matrices are generated. Then, a positive semidefinite matrix B 1 = O 1 A 1 O 1 T and a negative semidefinite matrix B 2 = O 2 A 2 O 2 T are achieved. Lastly, set Q = B 1 + B 1 T , S = B 2 + B 2 T and P = Q S . ). The constraint set C R m be defined by
C : = { u R m : 10 u i 10 } .
Numerical explanations for the first 200 iterations of three methods are considered in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 and Table 1 by letting initial points u 0 = u 1 = ( 1 , 1 , , 1 , 1 ) T . For Algorithm 3.2 (mAlg2) in [34]: χ = 1 4 c 1 and ρ n = 1 100 ( n + 2 ) ; For Algorithm (mAlg3) in (60): χ 0 = 0.20 , μ = 0.70 , ϱ n = 1 100 ( n + 2 ) , ρ n = 0.5 ( 1 ϱ n ) ; For Algorithm 1 (mAlg1): χ 0 = 0.20 , μ = 0.70 , ϕ = 0.60 , ς n = 1 ( n + 1 ) 2 , ϱ n = 1 100 ( n + 2 ) and ρ n = 0.5 ( 1 ϱ n ) .
Example 2.
Assume that set C L 2 ( [ 0 , 1 ] is defined by
C : = { u L 2 ( [ 0 , 1 ] ) : u 1 } .
Let us define an operator L : C H , such that
L ( u ) ( t ) = 0 1 u ( t ) H ( t , s ) f ( u ( s ) ) d s + g ( t ) ,
where H ( t , s ) = 2 t s e ( t + s ) e e 2 1 , f ( u ) = cos ( u ) and g ( t ) = 2 t e t e e 2 1 . In the above H = L 2 ( [ 0 , 1 ] ) is a Hilbert space with inner product u , y = 0 1 u ( t ) y ( t ) d t , u , y H and induced norm is u = 0 1 | u ( t ) | 2 d t . Numerical explanations for the first 200 iterations of three methods are considered in Figure 7, Figure 8, Figure 9 and Figure 10 by letting initial points u 0 = u 1 = ( 1 , 1 , , 1 , 1 ) T . For Algorithm 3.2 (mAlg2) in [34]: χ = 1 3 c 1 and ρ n = 1 100 ( n + 2 ) ; For Algorithm (mAlg3) in (60): χ 0 = 0.50 , μ = 0.50 , ϱ n = 1 100 ( n + 2 ) , ρ n = 0.7 ( 1 ϱ n ) ; For Algorithm 1 (mAlg1): χ 0 = 0.50 , μ = 0.50 , ϕ = 0.70 , ς n = 1 ( n + 1 ) 2 , ϱ n = 1 100 ( n + 2 ) and ρ n = 0.7 ( 1 ϱ n ) .

7. Conclusions

We studied a Mann-type extragradient-like scheme for determining the numerical solution of equilibrium problem involving pseudomonotone function and also prove a strong convergent theorem. Computational conclusions were established to illustrate the computational performance of our algorithms relative to other approaches. Such computational experiments showed that the inertial effect increases the efficacy of the iterative method in this sense.

Author Contributions

Formal analysis, C.K.; funding acquisition, N.P. (Nuttapol Pakkaranang), N.P. (Nattawut Pholasa) and C.K.; investigation, N.W., N.P. (Nuttapol Pakkaranang) and C.K.; methodology, C.K.; project administration, C.K., N.P. (Nattawut Pholasa) and C.K.; resources, N.P. (Nattawut Pholasa) and C.K.; software, N.P. (Nuttapol Pakkaranang); supervision, N.P. (Nattawut Pholasa); Writing—original draft, N.W. and N.P. (Nuttapol Pakkaranang); Writing—review and editing, N.P. (Nuttapol Pakkaranang). All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by University of Phayao.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Chainarong Khanpanuk would like to thank Phetchabun Rajabhat University. Nattawut Pholasa was financial supported by University of Phayao.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, K. A Minimax Inequality and Applications, Inequalities III; Shisha, O., Ed.; Academic Press: New York, NY, USA, 1972. [Google Scholar]
  2. Muu, L.; Oettli, W. Convergence of an adaptive penalty scheme for finding constrained equilibria. Nonlinear Anal. Theory Methods Appl. 1992, 18, 1159–1166. [Google Scholar] [CrossRef]
  3. Blum, E. From optimization and variational inequalities to equilibrium problems. Math. Stud. 1994, 63, 123–145. [Google Scholar]
  4. Bigi, G.; Castellani, M.; Pappalardo, M.; Passacantando, M. Existence and solution methods for equilibria. Eur. J. Oper. Res. 2013, 227, 1–11. [Google Scholar] [CrossRef] [Green Version]
  5. Antipin, A. Equilibrium programming: Proximal methods. Comput. Math. Math. Phys. 1997, 37, 1285–1296. [Google Scholar]
  6. Giannessi, F.; Maugeri, A.; Pardalos, P.M. Equilibrium Problems: Nonsmooth Optimization and Variational Inequality Models; Kluwer Academic Publisher: Dordrecht, The Netherlands, 2006; Volume 58. [Google Scholar]
  7. Dafermos, S. Traffic Equilibrium and Variational Inequalities. Transp. Sci. 1980, 14, 42–54. [Google Scholar] [CrossRef] [Green Version]
  8. ur Rehman, H.; Kumam, P.; Cho, Y.J.; Yordsorn, P. Weak convergence of explicit extragradient algorithms for solving equilibirum problems. J. Inequalities Appl. 2019, 2019. [Google Scholar] [CrossRef]
  9. Ur Rehman, H.; Pakkaranang, N.; Kumam, P.; Cho, Y.J. Modified subgradient extragradient method for a family of pseudomonotone equilibrium problems in real a Hilbert space. J. Nonlinear Convex Anal. 2020, 21, 2011–2025. [Google Scholar]
  10. Ur Rehman, H.; Kumam, P.; Dong, Q.L.; Cho, Y.J. A modified self-adaptive extragradient method for pseudomonotone equilibrium problem in a real Hilbert space with applications. Math. Methods Appl. Sci. 2020, 44, 3527–3547. [Google Scholar] [CrossRef]
  11. Ur Rehman, H.; Kumam, P.; Sitthithakerngkiet, K. Viscosity-type method for solving pseudomonotone equilibrium problems in a real Hilbert space with applications. AIMS Math. 2021, 6, 1538–1560. [Google Scholar] [CrossRef]
  12. Ur Rehman, H.; Kumam, P.; Gibali, A.; Kumam, W. Convergence analysis of a general inertial projection-type method for solving pseudomonotone equilibrium problems with applications. J. Inequalities Appl. 2021, 2021. [Google Scholar] [CrossRef]
  13. Rehman, H.U.; Kumam, P.; Dong, Q.L.; Peng, Y.; Deebani, W. A new Popov’s subgradient extragradient method for two classes of equilibrium programming in a real Hilbert space. Optimization 2020, 1–36. [Google Scholar] [CrossRef]
  14. Ur Rehman, H.; Kumam, P.; Shutaywi, M.; Alreshidi, N.A.; Kumam, W. Inertial Optimization Based Two-Step Methods for Solving Equilibrium Problems with Applications in Variational Inequality Problems and Growth Control Equilibrium Models. Energies 2020, 13, 3292. [Google Scholar] [CrossRef]
  15. Ferris, M.C.; Pang, J.S. Engineering and Economic Applications of Complementarity Problems. Siam Rev. 1997, 39, 669–713. [Google Scholar] [CrossRef] [Green Version]
  16. Ur Rehman, H.; Kumam, P.; Cho, Y.J.; Suleiman, Y.I.; Kumam, W. Modified Popov’s explicit iterative algorithms for solving pseudomonotone equilibrium problems. Optim. Methods Softw. 2020, 36, 82–113. [Google Scholar] [CrossRef]
  17. Ur Rehman, H.; Kumam, P.; Abubakar, A.B.; Cho, Y.J. The extragradient algorithm with inertial effects extended to equilibrium problems. Comput. Appl. Math. 2020, 39. [Google Scholar] [CrossRef]
  18. Ur Rehman, H.; Pakkaranang, N.; Hussain, A.; Wairojjana, N. A modified extra-gradient method for a family of strongly pseudomonotone equilibrium problems in real Hilbert spaces. J. Math. Comput. Sci. 2020, 22, 38–48. [Google Scholar] [CrossRef]
  19. Wairojjana, N.; ur Rehman, H.; Argyros, I.K.; Pakkaranang, N. An Accelerated Extragradient Method for Solving Pseudomonotone Equilibrium Problems with Applications. Axioms 2020, 9, 99. [Google Scholar] [CrossRef]
  20. Wairojjana, N.; ur Rehman, H.; la Sen, M.D.; Pakkaranang, N. A General Inertial Projection-Type Algorithm for Solving Equilibrium Problem in Hilbert Spaces with Applications in Fixed-Point Problems. Axioms 2020, 9, 101. [Google Scholar] [CrossRef]
  21. Wairojjana, N.; Pakkaranang, N.; Ur Rehman, H.; Pholasa, N.; Khanpanuk, T. Strong Convergence of Extragradient-Type Method to Solve Pseudomonotone Variational Inequalities Problems. Axioms 2020, 9, 115. [Google Scholar] [CrossRef]
  22. Wairojjana, N.; Younis, M.; Ur Rehman, H.; Pakkaranang, N.; Pholasa, N. Modified Viscosity Subgradient Extragradient-Like Algorithms for Solving Monotone Variational Inequalities Problems. Axioms 2020, 9, 118. [Google Scholar] [CrossRef]
  23. Tran, D.Q.; Dung, M.L.; Nguyen, V.H. Extragradient algorithms extended to equilibrium problems. Optimization 2008, 57, 749–776. [Google Scholar] [CrossRef]
  24. Korpelevich, G. The extragradient method for finding saddle points and other problems. Matecon 1976, 12, 747–756. [Google Scholar]
  25. Hieu, D.V.; Quy, P.K.; Vy, L.V. Explicit iterative algorithms for solving equilibrium problems. Calcolo 2019, 56. [Google Scholar] [CrossRef]
  26. Polyak, B. Some methods of speeding up the convergence of iteration methods. USSR Comput. Math. Math. Phys. 1964, 4, 1–17. [Google Scholar] [CrossRef]
  27. Beck, A.; Teboulle, M. A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems. SIAM J. Imaging Sci. 2009, 2, 183–202. [Google Scholar] [CrossRef] [Green Version]
  28. Hung, P.G.; Muu, L.D. The Tikhonov regularization extended to equilibrium problems involving pseudomonotone bifunctions. Nonlinear Anal. Theory Methods Appl. 2011, 74, 6121–6129. [Google Scholar] [CrossRef]
  29. Konnov, I. Application of the Proximal Point Method to Nonmonotone Equilibrium Problems. J. Optim. Theory Appl. 2003, 119, 317–333. [Google Scholar] [CrossRef]
  30. Moudafi, A. Proximal point algorithm extended to equilibrium problems. J. Nat. Geom. 1999, 15, 91–100. [Google Scholar]
  31. Oliveira, P.; Santos, P.; Silva, A. A Tikhonov-type regularization for equilibrium problems in Hilbert spaces. J. Math. Anal. Appl. 2013, 401, 336–342. [Google Scholar] [CrossRef]
  32. Mann, W.R. Mean value methods in iteration. Proc. Am. Math. Soc. 1953, 4, 506. [Google Scholar] [CrossRef]
  33. Censor, Y.; Gibali, A.; Reich, S. The Subgradient Extragradient Method for Solving Variational Inequalities in Hilbert Space. J. Optim. Theory Appl. 2010, 148, 318–335. [Google Scholar] [CrossRef] [Green Version]
  34. Hieu, D.V. Halpern subgradient extragradient method extended to equilibrium problems. Rev. Real Acad. Cienc. Exactas FÍSicas Nat. Ser. Matemíticas 2016, 111, 823–840. [Google Scholar] [CrossRef]
  35. Bianchi, M.; Schaible, S. Generalized monotone bifunctions and equilibrium problems. J. Optim. Theory Appl. 1996, 90, 31–43. [Google Scholar] [CrossRef]
  36. Mastroeni, G. On Auxiliary Principle for Equilibrium Problems. In Nonconvex Optimization and Its Applications; Springer: New York, NY, USA, 2003; pp. 289–298. [Google Scholar] [CrossRef]
  37. Tiel, J.V. Convex Analysis: An Introductory Text, 1st ed.; Wiley: New York, NY, USA, 1984. [Google Scholar]
  38. Kreyszig, E. Introductory Functional Analysis with Applications, 1st ed.; Wiley: New York, NY, USA, 1989. [Google Scholar]
  39. Xu, H.K. Another control condition in an iterative method for nonexpansive mappings. Bull. Aust. Math. Soc. 2002, 65, 109–113. [Google Scholar] [CrossRef] [Green Version]
  40. Maingé, P.E. Strong Convergence of Projected Subgradient Methods for Nonsmooth and Nonstrictly Convex Minimization. Set Valued Anal. 2008, 16, 899–912. [Google Scholar] [CrossRef]
  41. Bauschke, H.H.; Combettes, P.L. Convex Analysis and Monotone Operator Theory in Hilbert Spaces, 2nd ed.; CMS Books in Mathematics; Springer International Publishing: New York, NY, USA, 2017. [Google Scholar]
  42. Browder, F.; Petryshyn, W. Construction of fixed points of nonlinear mappings in Hilbert space. J. Math. Anal. Appl. 1967, 20, 197–228. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, S.; Zhang, Y.; Ping, P.; Cho, Y.; Guo, H. New extragradient methods with non-convex combination for pseudomonotone equilibrium problems with applications in Hilbert spaces. Filomat 2019, 33, 1677–1693. [Google Scholar] [CrossRef] [Green Version]
  44. Stampacchia, G. Formes bilinéaires coercitives sur les ensembles convexes. Comptes Rendus Hebd. Seances Acad. Sci. 1964, 258, 4413. [Google Scholar]
  45. Konnov, I.V. On systems of variational inequalities. Russ. Math. C/C Izv. Vyss. Uchebnye Zaved. Mat. 1997, 41, 77–86. [Google Scholar]
Figure 1. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 5 .
Figure 1. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 5 .
Axioms 10 00076 g001
Figure 2. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 10 .
Figure 2. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 10 .
Axioms 10 00076 g002
Figure 3. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 20 .
Figure 3. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 20 .
Axioms 10 00076 g003
Figure 4. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 50 .
Figure 4. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 50 .
Axioms 10 00076 g004
Figure 5. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 100 .
Figure 5. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 100 .
Axioms 10 00076 g005
Figure 6. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 200 .
Figure 6. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for R 200 .
Axioms 10 00076 g006
Figure 7. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + t + 2 t 2 .
Figure 7. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + t + 2 t 2 .
Axioms 10 00076 g007
Figure 8. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 2 t + 3 e t .
Figure 8. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 2 t + 3 e t .
Axioms 10 00076 g008
Figure 9. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 2 t + sin ( t ) .
Figure 9. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 2 t + sin ( t ) .
Axioms 10 00076 g009
Figure 10. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 3 t 2 + cos ( t ) .
Figure 10. Algorithm 1 compared to Algorithm (60) and Algorithm 3.2 in [34] for u 0 = 1 + 3 t 2 + cos ( t ) .
Axioms 10 00076 g010
Table 1. Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 execution time required for first 200 iterations.
Table 1. Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 execution time required for first 200 iterations.
Execution Time in Seconds
m mAlg2mAlg3mAlg1
52.558468122.736222482.923849848
102.898231332.998536853.341848537
203.238472543.518352123.332562246
503.936450464.054621574.084188882
1004.578374365.328735485.723835682
2005.862418366.281947136.825465869
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khanpanuk, C.; Pakkaranang, N.; Wairojjana, N.; Pholasa, N. Approximations of an Equilibrium Problem without Prior Knowledge of Lipschitz Constants in Hilbert Spaces with Applications. Axioms 2021, 10, 76. https://doi.org/10.3390/axioms10020076

AMA Style

Khanpanuk C, Pakkaranang N, Wairojjana N, Pholasa N. Approximations of an Equilibrium Problem without Prior Knowledge of Lipschitz Constants in Hilbert Spaces with Applications. Axioms. 2021; 10(2):76. https://doi.org/10.3390/axioms10020076

Chicago/Turabian Style

Khanpanuk, Chainarong, Nuttapol Pakkaranang, Nopparat Wairojjana, and Nattawut Pholasa. 2021. "Approximations of an Equilibrium Problem without Prior Knowledge of Lipschitz Constants in Hilbert Spaces with Applications" Axioms 10, no. 2: 76. https://doi.org/10.3390/axioms10020076

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop