Next Article in Journal
Existence and Asymptotic Stability of the Solution for the Timoshenko Transmission System with Distributed Delay
Next Article in Special Issue
Common Fixed Point Results on a Double-Controlled Metric Space for Generalized Rational-Type Contractions with Application
Previous Article in Journal
The Stability Analysis of Linear Systems with Cauchy—Polynomial-Vandermonde Matrices
Previous Article in Special Issue
Positive Solutions for a System of Hadamard Fractional Boundary Value Problems on an Infinite Interval
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified Inertial-Type Subgradient Extragradient Methods for Variational Inequalities and Fixed Points of Finite Bregman Relatively Nonexpansive and Demicontractive Mappings

Department of Mathematics, Shanghai Normal University, Shanghai 200234, China
*
Author to whom correspondence should be addressed.
Submission received: 23 July 2023 / Revised: 17 August 2023 / Accepted: 21 August 2023 / Published: 28 August 2023
(This article belongs to the Special Issue Research on Fixed Point Theory and Application)

Abstract

:
In this paper, we design two inertial-type subgradient extragradient algorithms with the linear-search process for resolving the two pseudomonotone variational inequality problems (VIPs) of and the common fixed point problem (CFPP) of finite Bregman relatively nonexpansive operators and Bregman relatively demicontractive operators in Banach spaces of both p-uniform convexity and uniform smoothness, which are more general than Hilbert ones. By the aid of suitable restrictions, it is shown that the sequences fabricated by the suggested schemes converge weakly and strongly to a solution of a pair of VIPs with a CFPP constraint, respectively. Additionally, the illustrative instance is furnished to back up the practicability and implementability of the suggested methods. This paper reveals the competitive advantage of the proposed algorithms over the existing algorithms; that is, the existing hybrid projection method for a single VIP with an FPP constraint is extended to develop the modified inertial-type subgradient extragradient method for a pair of VIPs with an CFPP constraint.

1. Introduction

Suppose that the real Hilbert space H has both inner product · , · and induced norm · , and let P C be the metric projection of H onto a nonempty, convex and closed C H given a nonlinear operator S : C C . We denote by Fix ( S ) the fixed-point set of S. Also, the R , and → are used to represent the real-number set, the weak convergence and the strong convergence, respectively. A mapping S : C C is said to be strictly pseudocontractive (see [1]) if ξ [ 0 , 1 ) s.t. S x S y 2 x y 2 + ξ ( I S ) x ( I S ) y 2 x , y C . In particular, in case ξ = 0 , S reduces to a nonexpansive mapping. Moreover, S is said to be demicontractive if Fix ( S ) and ξ [ 0 , 1 ) s.t. S x y 2 x y 2 + ξ x S x 2 x C , y Fix ( S ) . In particular, in case ξ = 0 , S reduces to a quasi-nonexpansive mapping. During the past few decades, the fixed point theory has played a vital part in solving many problems arising in nonlinear analysis and optimization theory, such as differential hemivariational inequalities (see [2]), monotone bilevel equilibrium problems (see [3]), fractional set-valued projected dynamical systems (see [4]) and so on.
Let A : H H be a mapping. Consider the classical variational inequality problem (VIP) of finding u C s.t. A u , v u 0 v C . The solution set of the VIP is denoted by VI ( C , A ) . In 1976, to seek a point in VI ( C , A ) , via relative weak conditions, Korpelevich [5] put forward an extragradient approach below; i.e., for any initial x 0 C , the sequence { x n } is generated by
y n = P C ( x n τ A x n ) , x n + 1 = P C ( x n τ A y n ) n 0 ,
with τ ( 0 , 1 L ) . If VI ( C , A ) , then the sequence { x n } converges weakly to an element in VI ( C , A ) . To the best of our knowledge, the Korpelevich extragradient approach is one of the most effective methods for solving the VIP at present. The literature on the VIP is vast and the Korpelevich extragradient approach has attained wide attention paid by many scholars, who ameliorated it in various forms; see, e.g., [1,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26].
Furthermore, in 2018, Thong and Hieu [18] first put forward the inertial subgradient extragradient method; that is, for any initial x 0 , x 1 H , the sequence { x n } is generated by
g n = x n + α n ( x n x n 1 ) , y n = P C ( g n A g n ) , C n = { v H : g n A g n y n , v y n 0 } , x n + 1 = P C n ( g n A y n ) n 1 ,
with constant ( 0 , 1 L ) . Under suitable conditions, they proved the weak convergence of { x n } to an element of VI ( C , A ) . Subsequently, Ceng et al. [14] proposed the modification of an inertial subgradient extragradient approach for settling the VIP of pseudomonotonicity and common fixed-point problem (CFPP) of finite nonexpansive mappings. Let S i : H H be nonexpansive for i = 1 , , N , A : H H be of both L-Lipschitz continuity and pseudomonotonicity on H, and of sequentially weak continuity on C, s.t. Ω = j = 1 N Fix ( S j ) VI ( C , A ) . Let f be δ -contractive self-mapping on H for δ [ 0 , 1 ) and the self-mapping F on H be of both η -strong monotone and κ -Lipschitz continuity s.t. τ : = 1 1 ρ ( 2 η ρ κ 2 ) > δ with ρ ( 0 , 2 η / κ 2 ) . Presume that { β n } , { γ n } , { τ n } are positive sequences s.t. β n + γ n < 1 , n = 1 β n = , lim n β n = 0 , lim inf n γ n ( 1 γ n ) > 0 and τ n = o ( β n ) . Moreover, one writes S n : = S n mod N for integer n 1 with the mod function taking values in the set { 1 , , N } ; i.e., if n = j N + m for some integers j 0 and 0 m < N , then S n = S N if m = 0 and S n = S m if 0 < m < N .
Under appropriate conditions, they proved the strong convergence of { x n } to an element of Ω = i = 1 N Fix ( S i ) VI ( C , A ) . In addition, combining the subgradient extragradient method and Halpern’s iteration method, Kraikaew and Saejung [19] proposed the Halpern subgradient extragradient rule for solving the VIP in 2014. They proved the strong convergence of the proposed method to an element in VI ( C , A ) . In 2021, Reich et al. [23] invented two gradient-projection schemes for settling the VIP for uniformly continuous pseudomonotone mapping. In particular, they used a novel Armijo-type line search to acquire a hyperplane that strictly separates the current iterate from the solutions of the VIP under consideration. They proved that the sequences generated by two schemes converge weakly and strongly to a point in VI ( C , A ) for uniformly continuous pseudomonotone mapping A, respectively.
On the other hand, let 1 / p + 1 / q = 1 for p , q > 1 , and suppose that E is a Banach space of both p-uniform convexity and uniform smoothness and the nonempty C E is of both convexity and closedness. The dual space of E is denoted by E * . The norm and the duality pairing between E and E * are denoted by · and · , · , respectively. Let J E p and J E * q be the duality mappings of E and E * , respectively. Let f p ( u ) = u p / p u E , D f p be the Bregman distance with respect to (w.r.t) f p and the surjective Π C : E C be the Bregman’s projection w.r.t. f p , and presume that { α n } , { β n } ( 0 , 1 ) s.t. lim n α n = 0 , lim inf n β n ( 1 β n ) > 0 and n = 1 α n = . Assume that A : E E * is uniformly continuous and pseudo-monotone mapping and S : C C is Bregman relatively nonexpansive mapping. Very recently, inspired by the research outcomes in [23], Eskandani et al. [25] invented the hybrid projection method with linear search term in order to seek a solution of a VIP with an FPP constraint of S.
By the aid of mild restrictions, it was proven in [23] that the sequence { u n } converges strongly to Π Ω u . Motivated by the existing outcomes as above, we design two inertial-type subgradient extragradient algorithms with a linear-search process for resolving the two pseudomonotone VIPs and the CFPP of finite Bregman’s relative nonexpansivity operators and a Bregman’s relative demicontractivity operator in Banach spaces of both p-uniform convexity and uniform smoothness. With the help of appropriate assumptions, it is proven that the sequences fabricated by the suggested algorithms converge weakly and strongly to a solution of a pair of VIPs with a CFPP constraint, respectively. Additionally, an illustrative instance is furnished to back up the practicability and implementability of the proposed approaches.
The structure for the paper is built as follows: Section 2 releases certain terminologies and preliminary results. In Section 3, we discuss the convergent behavior of the sequences generated by the proposed approaches. In Section 4, the major outcomes are employed to deal with a pair of VIPs with a CFPP constraint in an illustrative instance. Our algorithms are of both advantage and flexibility over Algorithms 1 and 2 as above due to their solving a pair of VIPs with a CFPP constraint. Our outcomes are the improvement and extension of the existing ones in the literature; see, e.g., [14,23,25].
Algorithm 1: ([14], Algorithm 3)
  Inertial subgradient extragradient method. Initialization: Given λ 1 > 0 , α > 0 , μ ( 0 , 1 ) . Let x 0 , x 1 H be arbitrary.
Iterative steps: Calculate x n + 1 as follows:
Step 1. Given the iterates x n 1 and x n ( n 1 ) , choose α n s.t. 0 α n α n ¯ , where
α n ¯ = min { α , τ n x n x n 1 } if x n x n 1 , α otherwise .
Step 2. Compute w n = S n x n + α n ( S n x n S n x n 1 ) and y n = P C ( w n λ n A w n ) .
Step 3. Construct the half-space C n : = { z H : w n λ n A w n y n , z y n 0 } , and compute z n = P C n ( w n λ n A y n ) .
Step 4. Calculate x n + 1 = β n f ( x n ) + γ n x n + ( ( 1 γ n ) I β n ρ F ) z n , and update
λ n + 1 = min { μ w n y n 2 + z n y n 2 2 A w n A y n , z n y n , λ n } A w n A y n , z n y n > 0 , λ n otherwise .
Set n : = n + 1 and go to Step 1.
Algorithm 2: ([25])
  Hybrid projection method. Initial step: Let positive l < 1 , ν > 0 , λ ( 0 , 1 ν ) , and put u , u 1 C arbitrarily.
Iterations: Compute u n + 1 below:
Step 1. Calculate v n = Π C ( J E * q ( J E p u n λ A u n ) ) and r λ ( u n ) : = u n v n . If r λ ( u n ) = 0 and S u n = u n , then stop; u n Ω = Fix ( S ) VI ( C , A ) . Otherwise,
Step 2. Compute t n = u n τ n r λ ( u n ) , where τ n : = l k n and k n is the smallest nonnegative integer k satisfying A u n A ( u n l k r λ ( u n ) ) , r λ ( u n ) ν 2 D f p ( u n , v n ) .
Step 3. Compute w n = J E * q ( β n J E p u n + ( 1 β n ) J E p ( S Π C n u n ) ) and u n + 1 = Π C ( J E * q ( α n J E p u + ( 1 α n ) J E p w n ) ) , with C n : = { v C : h ¯ n ( v ) 0 } and h ¯ n ( v ) = A t n , v u n + τ n 2 λ D f p ( u n , v n ) .
Again set n : = n + 1 and go to Step 1.
In the end, it is worthy to mention that the existing method in [25] is most closely related to our proposed method; that is, the hybrid projection method for resolving a single VIP with an FPP constraint in [25] is extended to develop our modified inertial-type subgradient extragradient method for resolving a pair of VIPs with a CFPP constraint. Compared with the corresponding results in [25], our results improve, extend and develop them in the two aspects below: (i) the problem of finding a solution of a single VIP with an FPP constraint (involving a Bregman relatively nonexpansive mapping) in [25] is extended to develop our problem of finding a solution of a pair of VIPs with a CFPP constraint (involving finite Bregman relatively nonexpansive mappings and a Bregman relatively demicontractive mapping); (ii) the hybrid projection method with line-search process in [25] is extended to develop our modified inertial-type subgradient extragradient method with the line-search process.
Algorithm 3: The 1st modified inertial-type subgradient extragradient method
  The 1st modified inertial-type subgradient extragradient method. Initialization: Given x 0 , x 1 C arbitrarily and let ϵ > 0 , μ i > 0 , λ i ( 0 , 1 μ i ) , l i ( 0 , 1 ) for i = 1 , 2 . Choose { n } , { β n } ( 0 , 1 ) and { α n } ( ξ , 1 ) s.t. n = 1 n < , lim inf n β n ( 1 β n ) > 0 and lim inf n ( α n ξ ) ( 1 α n ) > 0 . Moreover, given the iterates x n 1 and x n ( n 1 ) , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where
ϵ n ¯ = min { ϵ , n J E p S n x n J E p ( 2 S n x n S n x n 1 ) } if S n x n S n x n 1 , ϵ otherwise .
Iterative steps: Calculate x n + 1 as follows:
Step 1. Calculate g n = J E * q ( ( 1 ϵ n ) J E p S n x n + ϵ n J E p ( 2 S n x n S n x n 1 ) ) and calculate u n = J E * q ( β n J E p x n + ( 1 β n ) J E p g n ) , y n = Π C ( J E * q ( J E p u n λ 1 A 1 u n ) ) , r λ 1 ( u n ) : = u n y n and s n = u n τ n r λ 1 ( u n ) , with τ n : = l 1 k n and k n being the smallest k 0 s.t.
A 1 u n A 1 ( u n l 1 k r λ 1 ( u n ) ) , u n y n μ 1 2 D f p ( u n , y n ) .
Step 2. Calculate w n = Π C n ( u n ) , with C n : = { y C : h n ( y ) 0 } and
h n ( y ) = A 1 s n , y u n + τ n 2 λ 1 D f p ( u n , y n ) .
Step 3. Calculate y ˜ n = Π C ( J E * q ( J E p w n λ 2 A 2 w n ) ) , r λ 2 ( w n ) : = w n y ˜ n and t n = w n τ ˜ n r λ 2 ( w n ) , with τ ˜ n : = l 2 j n and j n is the smallest j 0 s.t.
A 2 w n A 2 ( w n l 2 j r λ 2 ( w n ) ) , w n y ˜ n μ 2 2 D f p ( w n , y ˜ n ) .
Step 4. Calculate v n = J E * q ( α n J E p w n + ( 1 α n ) J E p ( S 0 w n ) ) and x n + 1 = Π C ˜ n Q n ( w n ) , with Q n : = { y C : D f p ( y , v n ) D f p ( y , w n ) } , C ˜ n : = { y C : h ˜ n ( y ) 0 } and
h ˜ n ( y ) = A 2 t n , y w n + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) .
Again set n : = n + 1 and go to Step 1.
In Section 4, we have provided a numerical example to show the competitive advantage of our proposed algorithms over the existing algorithms, e.g., the ones in [25]. In fact, we have provided the illustrative example of a pair of VIPs with a CFPP constraint in Section 4. Note that the existing method in [25] is only utilized for solving a single VIP with an FPP constraint. So, there is no way for this method to handle the numerical example in Section 4; that is, it is invalid for a pair of VIPs with a CFPP constraint. However, our suggested method can settle the illustrative example in Section 4. This ensures the competitive advantage of our proposed algorithms over the existing algorithms.

2. Preliminaries

Let the real Banach space ( E , · ) possess the dual space E * . The g n g (resp., g n g ) is used to stand for the strong (resp., weak) convergence of { g n } to g E . Moreover, the set of weak cluster points of { g n } is denoted by ω w ( g n ) , i.e., ω w ( g n ) = { g E : g n k g for some { g n k } { g n } } . Let { U } = { g E : g = 1 } and 1 < q 2 p with 1 p + 1 q = 1 . A Banach space E is referred to as being strictly convex if, for each y , g U with y g , one has y + g / 2 < 1 . E is referred to as being of uniform convexity if ϵ ( 0 , 2 ] , δ > 0 s.t. y , g { U } with y g ϵ , one has y + g / 2 1 δ . It is known that a uniformly convex Banach space is reflexive and strictly convex. The modulus of convexity of E is the function δ : [ 0 , 2 ] [ 0 , 1 ] defined by δ ( ϵ ) = inf { 1 y + g / 2 : y , g { U } with y g ϵ } . It is also known that E is uniformly convex if and only if δ ( ϵ ) > 0 ϵ ( 0 , 2 ] . Moreover, E is referred to as being p-uniformly convex if c > 0 s.t. δ ( ϵ ) c ϵ p ϵ [ 0 , 2 ] .
A function ρ E : [ 0 , ) [ 0 , ) is the modulus of smoothness iff it is written as ρ E ( τ ) = sup { ( y + τ g + y τ g ) / 2 1 : y , g { U } } . E is said to be uniformly smooth if lim τ 0 ρ E ( τ ) / τ = 0 , and q-uniformly smooth if C q > 0 s.t. ρ E ( τ ) C q τ q τ > 0 . It is known that E is p-uniformly convex if and only if E * is q-uniformly smooth. For example, see [27] for more details. Putting B ( 0 , r ) = { g E : g r } for each r > 0 , we say that f : E R is uniformly convex on bounded sets (see [25]) if ρ r ( t ) > 0 r , t > 0 , where ρ r ( t ) : [ 0 , ) [ 0 , ] is specified below
ρ r ( t ) = inf { [ α f ( y ) + ( 1 α ) f ( g ) f ( α y + ( 1 α ) g ) ] / α ( 1 α ) : α ( 0 , 1 ) and y , g B ( 0 , r ) with y g = t } ,
for all t 0 . ρ r is called the gauge function of f with uniform convexity. It is obvious that the function ρ r is nondecreasing.
Let f : E R be a convex function. If the limit lim t 0 + f ( y + t g ) f ( y ) t exists for each g E , then f is referred to as being of Gâteaux differentiability at y. In this case, the gradient of f at y is the linear function f ( y ) , which is defined by f ( y ) , g : = lim t 0 + f ( y + t g ) f ( y ) t for each g E . The function f is referred to as being of Gâteaux’s differentiability if it is of Gâteaux’s differentiability at each y E . Whenever lim t 0 + f ( y + t g ) f ( y ) t is achieved uniformly for any g { U , one says that f is of Fréchet’s differentiability at y. Furthermore, f is termed as being of uniform Fréchet differentiability on K E if lim t 0 + f ( y + t g ) f ( y ) t is achieved uniformly for ( y , g ) K × U . A Banach space E is called smooth if its norm is Gâteaux differentiable.
Let 1 / p + 1 / q = 1 for p , q > 1 . The duality mapping J E p : E E * is specified as follows
J E p ( y ) = { ϕ E * : ϕ , y = y p and ϕ = y p 1 } y E .
Recall that E is of smoothness iff the duality mapping J E p is single-valued. Also, E is of reflexivity iff J E p is of surjectivity, and E is of strict convexity iff J E p is an injection. So, it follows that, if E is smooth, strictly convex and reflexive Banach space, then J E p is a single-valued bijection, and, in this case, J E p = ( J E * q ) 1 . Furthermore, E is of uniform smoothness if and only if function f p ( y ) = y p / p is of uniform Fréchet differentiability on any bounded set if and only if the single-valued J E p is of uniform continuity on any bounded set. It is easy to see that E is of uniform convexity if and only if function f p is of uniform convexity (see [27]).
Let the convex function f : E R be of Gâteaux’s differentiability. Bregman’s distance w.r.t. f is specified below
D f ( y , g ) : = f ( y ) f ( g ) f ( g ) , y g y , g E .
It is worth mentioning that Bregman’s distance does not become a metric in the common sense of the terminology. Evidently, D f ( y , y ) = 0 , but D f ( y , g ) = 0 cannot yield y = g . Generally, D f is of no symmetry and does not fulfill the triangle inequality. But, D f fulfills the three point identity
D f ( y , g ) + D f ( g , w ) = D f ( y , w ) f ( g ) f ( w ) , y g .
See [28] for more details on Bregman functions and distances.
It is noteworthy that, if E is a smooth Banach space, then J E p is Gâteaux’s derivative of f p . Then, Bregman’s distance w.r.t. f p is specified below
D f p ( y , g ) = y p / p g p / p J E p ( g ) , y g = y p / p + g p / q J E p ( g ) , y = ( g p y p ) / q J E p ( g ) J E p ( y ) , y .
In the Banach space E of both smoothness and p-uniformly convexity for p 2 , there holds the following relationship between the metric and Bregman distance:
τ y g p D f p ( y , g ) J E p ( y ) J E p ( g ) , y g ,
where τ > 0 is some fixed number (see [29]). Using (5), one knows that, for any bounded { g n } E ,
g n g D f p ( g , g n ) 0 ( n ) .
Let the Banach space E be of reflexivity, smoothness and strict convexity. Let C E with C be of convexity and closedness. Bregman’s projections are formulated as minimizers of Bregman’s distances. The Bregman’s projection of y E onto C w.r.t. f p is only a point Π C y C s.t. D f p ( Π C y , y ) = min g C D f p ( g , y ) . In Hilbert spaces, the Bregman projection w.r.t. f 2 reduces to the metric projection.
Using Corollary 4.4 in [30] and Theorem 2.1 [31] in Banach spaces of uniform convexity, Bregman projections can be featured as the relation below:
J E p ( y ) J E p ( Π C y ) , g Π C y 0 g C .
Moreover, this inequality is equivalent to the descent property
D f p ( g , Π C y ) + D f p ( Π C y , y ) D f p ( g , y ) g C .
In case p = 2 , the duality mapping J E p reduces to the normalized duality mapping and is denoted by J. The function ϕ : E 2 R is formulated below
ϕ ( y , g ) = y 2 2 J g , y + g 2 y , g E ,
and Π C ( y ) = argmin g C ϕ ( g , y ) y E .
In terms of [25], the V f p : E × E * [ 0 , ) w.r.t. f p is specified below
V f p ( y , y * ) = y p / p y * , y + y * q / q ( y , y * ) E × E * .
So, V f p ( y , y * ) = D f p ( y , J E * q ( y * ) ) ( y , y * ) E × E * . Moreover, by the subdifferential inequality, we obtain
V f p ( y , y * ) + g * , J E * q ( y * ) y V f p ( y , y * + g * ) y E , y * , g * E * .
In addition, V f p is convex in the second variable. Thus, one has
D f p ( y , J E * q ( j = 1 n t j J E p ( g j ) ) ) j = 1 n t j D f p ( y , g j ) y E , { g j } j = 1 n E , j = 1 n t j = 1 for { t j } j = 1 n [ 0 , 1 ] .
Lemma 1
([31]). Suppose that the Banach space E is of uniform convexity and { y n } , { g n } are two sequences in E such that the first one is bounded. If lim n D f p ( g n , y n ) = 0 , then lim n g n y n = 0 .
Let S : C C be a mapping. We denote by Fix ( S ) the set of fixed points of S; that is, Fix ( S ) = { y C : y = S y } . A point g C is referred to as an asymptotic fixed point of S if { g n } C s.t. g n g and g n S g n 0 . We denote by Fix ^ ( S ) the set of asymptotic fixed points of S. The notion of asymptotic fixed points was invented in Reich [32]. A self-mapping S on C is termed as Bregman’s relatively ξ -demicontractive operator w.r.t. f p iff Fix ( S ) = Fix ^ ( S ) , and ξ [ 0 , 1 ) s.t. for each bounded { g n } C satisfying sup n 1 S g n < ; the following holds:
D f p ( y , S g n ) D f p ( y , g n ) + ξ ρ b * J E p g n J E p S g n y Fix ( S ) ,
with b = sup n 1 { g n p 1 , S g n p 1 } < . In particular, putting b x = max { x p 1 , S x p 1 } for each x C , one has
D f p ( y , S x ) D f p ( y , x ) + ξ ρ b x * J E p x J E p S x y Fix ( S ) .
In addition, if ξ = 0 , then S reduces to a mapping of Bregman’s relative nonexpansivity w.r.t. f p ; that is, S is named as a mapping of Bregman’s relative nonexpansivity w.r.t. f p if Fix ( S ) = Fix ^ ( S ) and D f p ( y , S g ) D f p ( y , g ) g C , y Fix ( S ) .
Definition 1.
Let C be a nonempty closed convex subset of E. A mapping A : C E * is known as being
(i) 
of monotonicity iff A y A g , y g 0 y , g C ;
(ii) 
of pseudo-monotonicity iff A y , g y 0 A g , g y 0 y , g C ;
(iii) 
ℓ-Lipschitz continuous or ℓ-Lipschitzian iff > 0 s.t. A y A g L y g y , g C ;
(iv) 
of weakly sequential continuity iff { g n } C ; the relation holds: g n g A g n A g .
Lemma 2
([25]). Given a constant r > 0 . If the function f : E R is of uniform convexity on any bounded subset of a Banach space E, then
f ( k = 1 n α k g k ) k = 1 n α k f ( g k ) α i α j ρ r ( g i g j ) ,
i , j { 1 , 2 , , n } , { g k } k = 1 n B ( 0 , r ) and { α k } k = 1 n ( 0 , 1 ) with k = 1 n α k = 1 , where ρ r is the gauge of uniform convexity of f.
Proof. 
It is easy to show the conclusion.    □
Lemma 3
([24]). Let E 1 and E 2 be two Banach spaces. Suppose that A : E 1 E 2 is uniformly continuous on bounded subsets of E 1 and D is a bounded subset of E 1 . Then, A ( D ) is bounded.
Lemma 4
([33]). Let C E with C being closed and convex in a Banach space E and let A : C E * be of both pseudo-monotonicity and continuity. Given x C . Then, A x , y x 0 y C A y , y x 0 y C .
Lemma 5.
Given p 2 . Suppose that the Banach space E is of both smoothness and p-uniform convexity s.t. the duality mapping J E p is of sequentially weak continuity. Let { q n } E and Ω E . If lim n { D f p ( g , q n ) } exists for each g Ω , and ω w ( q n ) Ω . Then, { q n } is weakly convergent to an element of Ω.
Proof. 
Using (5), we obtain τ y q n p D f p ( y , q n ) y Ω . This ensures that { q n } is of boundedness. Hence, from the reflexivity of E, we have ω w ( q n ) . Also, let us show the weak convergence of { q n } to a point in Ω . Indeed, let q ¯ , q ^ ω w ( q n ) with q ¯ q ^ . Then, { q n k } { q n } and { q m k } { q n } s.t. q n k q ¯ and q m k q ^ . From the sequentially weak continuity of J E p , we obtain J E p ( q n k ) J E p q ¯ and J E p ( q m k ) J E p q ^ . Note that D f p ( q ¯ , q ^ ) + D f p ( q ^ , q n ) = D f p ( q ¯ , q n ) J E p q ^ J E p q n , q ¯ q ^ . So, exploiting the convergence of the sequences { D f p ( q ¯ , q n ) } and { D f p ( q ^ , q n ) } , we deduce that
J E p q ^ J E p q ¯ , q ¯ q ^ = lim k [ J E p q ^ J E p q n k , q ¯ q ^ ] = lim n [ D f p ( q ¯ , q ^ ) + D f p ( q ^ , q n ) D f p ( q ¯ , q n ) ] = lim k [ J E p q ^ J E p q m k , q ¯ q ^ ] = J E p q ^ J E p q ^ , q ¯ q ^ = 0 ,
which, hence, yields J E p q ¯ J E p q ^ , q ¯ q ^ = 0 . From (5), we obtain 0 < τ q ¯ q ^ p D f p ( q ¯ , q ^ ) J E p q ¯ J E p q ^ , q ¯ q ^ = 0 . This arrives at a contradiction. Consequently, the sequence { q n } converges weakly to a point in Ω .    □
The lemma below was put forth in R m by [34]. It is easy to verify that the proof of the lemma in a Banach space E is actually the same as in R m . Here, we present the lemma but drop its demonstration.
Lemma 6.
Assume the nonempty C E with C being convex and closed in E. Suppose that K : = { y C : h ( y ) 0 } , where h : E R is real-valued. If K and h is Lipschitz continuous on C with modulus θ > 0 , then θ dist ( y , K ) max { h ( y ) , 0 } y C , where dist ( y , K ) stands for the distance of y to K.
Lemma 7
([35]). Let { Γ m } be a sequence of real numbers that does not decrease at infinity in the sense that { Γ m i } { Γ m } s.t. Γ m i < Γ m i + 1 i 1 . Let the sequence { ψ ( m ) } m m 0 of integers be defined as ψ ( m ) = max { i m : Γ i < Γ i + 1 } , with integer m 0 1 satisfying { i m 0 : Γ i < Γ i + 1 } . Then, the following holds:
(i) 
ψ ( m 0 ) ψ ( m 0 + 1 ) and ψ ( m ) ;
(ii) 
Γ ψ ( m ) Γ ψ ( m ) + 1 and Γ m Γ ψ ( m ) + 1 m m 0 .
Lemma 8
([36]). Let { a n } be a sequence in [ 0 , ) satisfying a n + 1 ( 1 μ n ) a n + μ n ν n n 1 , where { μ n } and { ν n } both are real sequences such that (i) { μ n } [ 0 , 1 ] and n = 1 μ n = , and (ii) lim sup n ν n 0 or n = 1 | μ n ν n | < . Then, lim n a n = 0 .
Lemma 9
([37]). Let { a n } , { b n } and { δ n } be sequences of nonnegative real numbers satisfying the inequality a n + 1 ( 1 + δ n ) a n + b n n 1 . If n = 1 δ n < and n = 1 b n < , and then lim n a n  exists.

3. Main Results

In this section, let E be a p-uniformly convex and uniformly smooth Banach space with 2 p < . Let C E with C be closed and convex in E. We are now in a position to state and analyze our iterative algorithms for settling a pair of VIPs with CFPP constraint, where the pair of VIPs implicates two mappings of both uniform continuity and pseudomonotonicity and the CFPP involves finite mappings of Bregman’s relative nonexpansivity and a mapping of Bregman’s relative demicontractivity in E. Assume always that the conditions hold below:
(C1)
For i = 1 , , N , the self-mapping S i on C is of both uniform continuity and Bregman’s relative nonexpansivity and self-mapping S 0 on C is of both uniform continuity and Bregman’s relative ξ -demicontractivity.
(C2)
{ S n } is defined as S n : = S n mod N for integer n 1 with the mod function taking values in the set { 1 , , N } ; i.e., if n = j N + m for some integers j 0 and 0 m < N , then S n = S N if m = 0 and S n = S m if 0 < m < N .
(C3)
For i = 1 , 2 , A i : E E * is pseudomonotone and uniformly continuous on C, s.t. A i x lim inf n A i x n { x n } C with x n x .
(C4)
Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 N Fix ( S i ) ) .
We will make use of Lemmas 10–13 below to derive our major outcomes in this paper.
Lemma 10.
Let { x n } be the constructed sequence in Algorithm 3. Then, the relations hold: 1 λ 1 D f p ( u n , y n ) A 1 u n , r λ 1 ( u n ) and 1 λ 2 D f p ( w n , y ˜ n ) A 2 w n , r λ 2 ( w n ) .
Proof. 
Observe that the last two relations are similar. Then, it suffices to show that the latter relation holds. In fact, using the definition of y ˜ n and properties of Π C , one has
J E p w n λ 2 A 2 w n J E p y ˜ n , y y ˜ n 0 y C .
Setting y = w n in the last inequality, from (5), we obtain
D f p ( w n , y ˜ n ) J E p w n J E p y ˜ n , w n y ˜ n λ 2 A 2 w n , w n y ˜ n .
This completes the proof.    □
Lemma 11.
The linear-search rules (1), (3) and the sequence { x n } constructed in Algorithm 3 are well defined.
Proof. 
Observe that the rules (1) and (3) are similar. Then, it suffices to show that the latter rule (3) is valid. Using the uniform continuity of A 2 on C, from l 2 ( 0 , 1 ) , one obtains lim j A 2 w n A 2 ( w n l 2 j r λ 2 ( w n ) ) , r λ 2 ( w n ) = 0 . In case r λ 2 ( w n ) = 0 , it is evident that j n = 0 . In case r λ 2 ( w n ) 0 , we know that j n 0 s.t. (3) holds.
It is evident to see that, for each n , Q ˜ n and C n are convex and closed. In what follows, we assert that Ω lies in Q n C ˜ n . Let z Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 N Fix ( S i ) ) . Using Lemma 2 and the Bregman relative ξ -demicontractivity of S 0 , from { α n } ( ξ , 1 ) , we obtain
D f p ( z , v n ) α n D f p ( z , w n ) + ( 1 α n ) D f p ( z , S 0 w n ) α n ( 1 α n ) ρ b w n * J E P w n J E P S 0 w n α n D f p ( z , w n ) + ( 1 α n ) [ D f p ( z , w n ) + ξ ρ b w n * J E P w n J E P S 0 w n ] α n ( 1 α n ) ρ b w n * J E P w n J E P S 0 w n = D f p ( z , w n ) ( α n ξ ) ( 1 α n ) ρ b w n * J E P w n J E P S 0 w n D f p ( z , w n ) ,
which, hence, leads to z Q n . Meanwhile, by Lemma 4, one obtains A 2 t n , z t n 0 . Therefore,
h ˜ n ( z ) = A 2 t n , z w n + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) = A 2 t n , w n t n A 2 t n , t n z + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) τ ˜ n A 2 t n , r λ 2 ( w n ) + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) .
So, it follows from (3) that
A 2 w n A 2 t n , r λ 2 ( w n ) μ 2 2 D f p ( w n , y ˜ n ) .
Using Lemma 10, we have
A 2 t n , r λ 2 ( w n ) A 2 w n , r λ 2 ( w n ) μ 2 2 D f p ( w n , y ˜ n ) ( 1 λ 2 μ 2 2 ) D f p ( w n , y ˜ n ) .
This together with (11) arrives at
h ˜ n ( z ) τ ˜ n 2 ( 1 λ 2 μ 2 ) D f p ( w n , y ˜ n ) 0 .
Consequently, Ω lies in Q n C ˜ n . So, { x n } is well defined.    □
Lemma 12.
Suppose that { y n } and { y ˜ n } are the sequences generated by Algorithm 3. If lim n u n y n = 0 and lim n w n y ˜ n = 0 , then ω w ( u n ) VI ( C , A 1 ) and ω w ( w n ) VI ( C , A 2 ) .
Proof. 
Observe that the last two inclusions are similar. Then, it suffices to show that the latter inclusion is valid. In fact, let z ω w ( w n ) . Whereby, we know that there exists a subsequence { w n k } of { w n } , satisfying both w n k z and w n k y ˜ n k 0 . Thus, one obtains y ˜ n k z . Noticing the convexity and closedness of C, from { y ˜ n } C and y ˜ n k z , we obtain z C . In what follows, one discusses two aspects. In case A 2 z = 0 , one has z VI ( C , A 2 ) because A 2 z , y z 0 y C . In case A 2 z 0 , by the condition on A 2 , one has lim inf k A 2 w n k A 2 z > 0 . So, we could assume that A 2 w n k 0 k 1 . From (6), we obtain
J E p w n k λ 2 A 2 w n k J E p y ˜ n k , x y ˜ n k 0 x C ,
and hence
1 λ 2 J E p w n k J E p y ˜ n k , x y ˜ n k + A 2 w n k , y ˜ n k w n k A 2 w n k , x w n k x C .
Note that A 2 is uniformly continuous. Then, we know that { A 2 w n k } is of boundedness by Lemma 3. Since { y ˜ n k } is of boundedness as well, using the uniform continuity of J E p on any bounded subset of E, from (12), we deduce that, for all x in C,
lim inf k A 2 w n k , x w n k 0 .
To show z VI ( C , A 2 ) , one picks a positive { ε ˜ k } ( 0 , 1 ) s.t. ε ˜ k 0 . For every k, l k is denoted as the smallest j 1 satisfying
A 2 w n j , y w n j + ε ˜ k 0 j l k .
Since { ε ˜ k } is decreasing, it is explicit that { l k } is increasing. For simplicity, { A 2 w n l k } is still written as { A 2 w l k } . Noticing A 2 w l k 0 k 1 (due to { A 2 w l k } { A 2 w n k } ), we put g ˜ l k = A 2 w l k A 2 w l k q q 1 and hence obtain A 2 w l k , J E * q g ˜ l k = 1 k 1 . Indeed, it is evident that A 2 w l k , J E * q g ˜ l k = A 2 w l k , ( 1 A 2 w l k q q 1 ) q 1 J E * q A 2 w l k = ( 1 A 2 w l k q q 1 ) q 1 A 2 w l k q = 1 k 1 . So, by (14), one has A 2 w l k , y + ε ˜ k J E * q g ˜ l k w l k 0 k 1 . Again, from the pseudomonotonicity of A 2 , one has that, for all y C ,
A 2 ( y + ε ˜ k J E * q g ˜ l k ) , y + ε ˜ k J E * q g ˜ l k w l k 0 .
Let us show lim k ε ˜ k J E * q g ˜ l k = 0 . Indeed, since { w l k } { w n k } and ε ˜ k 0 , we obtain that
0 lim sup k ε ˜ k J E * q g ˜ l k = lim sup k ε ˜ k A 2 w l k lim sup k ε ˜ k lim inf k A 2 w n k = 0 .
Hence, one obtains ε ˜ k J E * q g ˜ l k 0 as k . Thus, taking the limit as k in (15), by condition (C2), one has A 2 y , y z 0 y C . This implies that z VI ( C , A 2 ) .    □
Lemma 13.
Suppose that { y n } and { y ˜ n } are the sequences generated by Algorithm 3. Then, the following hold:
(i) 
lim n τ n D f p ( u n , y n ) = 0 lim n D f p ( u n , y n ) = 0 ;
(ii) 
lim n τ ˜ n D f p ( w n , y ˜ n ) = 0 lim n D f p ( w n , y ˜ n ) = 0 .
Proof. 
Observe that the assertions (i) and (ii) are similar. Then, it suffices to show that assertion (ii) is valid. To verify assertion (ii), we consider two cases. In case lim inf n τ ˜ n > 0 , we might presume that τ ˜ > 0 s.t. τ ˜ n τ ˜ > 0 n 1 , which, hence, arrives at
D f p ( w n , y ˜ n ) = 1 τ ˜ n τ ˜ n D f p ( w n , y ˜ n ) 1 τ ˜ · τ ˜ n D f p ( w n , y ˜ n ) .
Combining (16) and lim n τ ˜ n D f p ( w n , y ˜ n ) = 0 attains lim n D f p ( w n , y ˜ n ) = 0 .
In case lim inf n τ ˜ n = 0 , we presume lim sup n D f p ( w n , y ˜ n ) = a 2 > 0 . Whereby, one knows that { n k } { n } s.t.
lim k τ ˜ n k = 0 and lim k D f p ( w n k , y ˜ n k ) = a 2 > 0 .
One puts t n k ^ = 1 l 2 τ ˜ n k y ˜ n k + ( 1 1 l 2 τ ˜ n k ) w n k k . From (5) and lim k τ ˜ n k D f p ( w n k , y ˜ n k ) = 0 , we have lim k τ ˜ n k w n k y ˜ n k p = 0 and hence
lim k t n k ^ w n k p = lim k τ ˜ n k p 1 l 2 p · τ ˜ n k w n k y ˜ n k p = 0 .
Because A 2 is of uniform continuity on any bounded subset of C, one obtains
lim k A 2 w n k A 2 t n k ^ = 0 .
From the step size rule (3) and the definition of t n k ^ , it follows that
A 2 w n k A 2 t n k ^ , w n k y ˜ n k > μ 2 2 D f p ( w n k , y ˜ n k ) .
Now, taking the limit as k , from (18), we have lim k D f p ( w n k , y ˜ n k ) = 0 . This, however, reaches a contradiction. So, it follows that lim n D f p ( w n , y ˜ n ) is zero.    □
In what follows, we intend to demonstrate the first convergence result in this paper.
Theorem 1.
Suppose that the Banach space E is of both p-uniform convexity and uniform smoothness s.t. J E p is of sequentially weak continuity. If { x n } is the constructed sequence in Algorithm 3, then x n z Ω sup n 0 x n < .
Proof. 
It is clear that the necessity of Theorem 1 is valid. Next, it suffices to show that the sufficiency is valid. Assume that sup n 0 x n < . Let z Ω . It is clear that S n x n S n x n 1 J E p S n x n J E p ( 2 S n x n S n x n 1 ) . Using the definition of ϵ n , we obtain ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) n n 1 . From (5) and (10) and the three point identity of D f p , we obtain
D f p ( z , g n ) ( 1 ϵ n ) D f p ( z , S n x n ) + ϵ n D f p ( z , 2 S n x n S n x n 1 ) = D f p ( z , S n x n ) + ϵ n [ D f p ( z , 2 S n x n S n x n 1 ) D f p ( z , S n x n ) ] = D f p ( z , S n x n ) + ϵ n [ D f p ( S n x n , 2 S n x n S n x n 1 ) + J E p S n x n J E p ( 2 S n x n S n x n 1 ) , z S n x n ] D f p ( z , x n ) + ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) , z + S n x n 1 2 S n x n D f p ( z , x n ) + ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) z + S n x n 1 2 S n x n D f p ( z , x n ) + n M ,
where sup n 1 z + S n x n 1 2 S n x n M for some M > 0 . Using Lemma 2, we obtain
D f p ( z , u n ) = V f p ( z , β n J E p x n + ( 1 β n ) J E p g n ) 1 p z p β n J E p x n , z ( 1 β n ) J E p g n , z + β n q J E p x n q + ( 1 β n ) q J E p g n q β n ( 1 β n ) ρ b * J E p x n J E p g n = 1 p z p β n J E p x n , z ( 1 β n ) J E p g n , z + β n q x n p + ( 1 β n ) q g n p β n ( 1 β n ) ρ b * J E p x n J E p g n = β n D f p ( z , x n ) + ( 1 β n ) D f p ( z , g n ) β n ( 1 β n ) ρ b * J E p x n J E p g n D f p ( z , x n ) + n M β n ( 1 β n ) ρ b * J E p x n J E p g n .
Since w n = Π C n u n , by (5) and (7), we obtain
D f p ( z , w n ) D f p ( z , u n ) D f p ( w n , u n ) = D f p ( z , u n ) D f p ( Π C n u n , u n ) D f p ( z , u n ) τ Π C n u n u n p D f p ( z , u n ) τ P C n u n u n p = D f p ( z , u n ) τ [ dist ( C n , u n ) ] p .
Because x n + 1 = Π C ˜ n Q n w n , from (5) and (7), we obtain
D f p ( z , x n + 1 ) D f p ( z , w n ) D f p ( x n + 1 , w n ) D f p ( z , w n ) D f p ( Π C ˜ n w n , w n ) D f p ( z , w n ) τ P C ˜ n w n w n p = D f p ( z , w n ) τ [ dist ( C ˜ n , w n ) ] p .
This along with (19) arrives at
D f p ( z , x n + 1 ) D f p ( z , w n ) D f p ( x n + 1 , w n ) D f p ( z , u n ) D f p ( w n , u n ) D f p ( x n + 1 , w n ) D f p ( z , u n ) τ [ dist ( C n , u n ) ] p τ [ dist ( C ˜ n , w n ) ] p D f p ( z , x n ) + n M β n ( 1 β n ) ρ b * J E p x n J E p g n τ [ dist ( C n , u n ) ] p τ [ dist ( C ˜ n , w n ) ] p ,
which, hence, arrives at
D f p ( z , x n + 1 ) D f p ( z , x n ) + n M .
Since n = 1 n < , by Lemma 9, we deduce that lim n D f p ( z , x n ) exists. In addition, by the boundedness of { x n } , we conclude that { g n } , { u n } , { v n } , { w n } , { y n } , { y ˜ n } , { s n } , { t n } , { S n x n } and { S 0 w n } are also bounded. Using (20), we obtain
D f p ( w n , u n ) + D f p ( x n + 1 , w n ) D f p ( z , u n ) D f p ( z , x n + 1 ) D f p ( z , x n ) + n M β n ( 1 β n ) ρ b * J E p x n J E p g n D f p ( z , x n + 1 ) ,
which immediately yields
D f p ( w n , u n ) + D f p ( x n + 1 , w n ) + β n ( 1 β n ) ρ b * J E p x n J E p g n D f p ( z , x n ) D f p ( z , x n + 1 ) + n M .
Since lim n n = 0 , lim inf n β n ( 1 β n ) > 0 and lim n D f p ( z , x n ) exist, it follows that lim n D f p ( w n , u n ) = 0 , lim n D f p ( x n + 1 , w n ) = 0 and lim n ρ b * J E p x n J E p g n = 0 , which, hence, yields lim n J E p x n J E p g n = 0 . From u n = J E * q ( β n J E p x n + ( 1 β n ) J E p g n ) , it can be readily seen that lim n J E p u n J E p x n = 0 . Noticing g n = J E * q ( ( 1 ϵ n ) J E p S n x n + ϵ n J E p ( 2 S n x n S n x n 1 ) ) , we obtain from lim n n = 0 and the definition of ϵ n that
J E p g n J E p S n x n = ϵ n J E p ( 2 S n x n S n x n 1 ) J E p S n x n n 0 ( n ) .
Hence, using (5) and uniform continuity of J E p on bounded subsets of E, we conclude that lim n g n S n x n = 0 and
lim n w n u n = lim n x n + 1 w n = lim n x n S n x n = lim n u n x n = 0 .
Since { x n } is bounded and E is reflexive, then we know that ω w ( x n ) . In what follows, we claim that ω w ( x n ) lies in Ω . Let z ω w ( x n ) . Whereby, one knows that there exists a subsequence { x n k } of { x n } such that x n k converges weakly to z. From (21), one obtains w n k z . Since { A 1 s n } is of boundedness, one infers that L 1 > 0 such that A 1 s n L 1 . So, it follows that, for each y , g C n ,
| h n ( y ) h n ( g ) | = | A 1 s n , y g | A 1 s n y g L 1 y g ,
which implies that h n ( y ) is L 1 -Lipschitzian. Using Lemma 6, one obtains
dist ( C n , u n ) 1 L 1 h n ( u n ) = τ n 2 λ 1 L 1 D f p ( u n , y n ) .
Noticing x n + 1 Q n , from the definition of Q n and (20), we have
D f p ( x n + 1 , v n ) D f p ( x n + 1 , w n ) D f p ( z , w n ) D f p ( z , x n + 1 ) D f p ( z , u n ) D f p ( z , x n + 1 ) D f p ( z , x n ) D f p ( z , x n + 1 ) + n M .
Since lim n n = 0 and lim n D f p ( z , x n ) exists, we have lim n D f p ( x n + 1 , v n ) = 0 , which immediately yields x n + 1 v n 0 . Hence, from (21), we have
lim n w n v n = 0 .
Furthermore, by Lemma 2, we obtain that
D f p ( z , v n ) = V f p ( z , α n J E p w n + ( 1 α n ) J E p S 0 w n ) 1 p z p α n J E p w n , z ( 1 α n ) J E p S 0 w n , z + α n q J E p w n q + ( 1 α n ) q J E p S 0 w n q α n ( 1 α n ) ρ b * J E p w n J E p S 0 w n = 1 p z p α n J E p w n , z ( 1 α n ) J E p S 0 w n , z + α n q w n p + ( 1 α n ) q S 0 w n p α n ( 1 α n ) ρ b * J E p w n J E p S 0 w n = α n D f p ( z , w n ) + ( 1 α n ) D f p ( z , S 0 w n ) α n ( 1 α n ) ρ b * J E p w n J E p S 0 w n α n D f p ( z , w n ) + ( 1 α n ) [ D f p ( z , w n ) + ξ ρ b * J E p w n J E p S 0 w n ] α n ( 1 α n ) ρ b * J E p w n J E p S 0 w n = D f p ( z , w n ) ( α n ξ ) ( 1 α n ) ρ b * J E p w n J E p S 0 w n .
Therefore,
( α n ξ ) ( 1 α n ) ρ b * J E p w n J E p S 0 w n D f p ( z , w n ) D f p ( z , v n ) D f p ( z , w n ) D f p ( z , v n ) + D f p ( w n , v n ) = J E p v n J E p w n , z w n .
Taking the limit in the last inequality as n , and using uniform continuity of J E p on bounded subsets of E, (23) and lim inf n ( α n ξ ) ( 1 α n ) > 0 , we obtain lim n ρ b * J E p w n J E p S 0 w n = 0 and hence lim n J E p w n J E p S 0 w n = 0 . Since J E * q is uniformly continuous on any bounded subset of E * , we have
lim n w n S 0 w n = 0 .
Now, let us show z i = 1 2 VI ( C , A i ) . By the boundedness of { A 2 t n } , one knows that there exists L 2 > 0 , satisfying A 2 t n L 2 . So, it follows that, for each y , g C ˜ n ,
| h ˜ n ( y ) h ˜ n ( g ) | = | A 2 t n , y g | A 2 t n y g L 2 y g ,
which guarantees that h ˜ n ( y ) is L 2 -Lipschitzian. By Lemma 6, one obtains
dist ( C ˜ n , w n ) 1 L 2 h ˜ n ( w n ) = τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) .
Combining (20), (22) and (25), we obtain
D f p ( z , x n ) D f p ( z , x n + 1 ) + n M D f p ( z , u n ) D f p ( z , x n + 1 ) τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p .
Thus,
lim n τ ˜ n D f p ( w n , y ˜ n ) = lim n τ n D f p ( u n , y n ) = 0 .
By Lemma 13, we obtain
lim n w n y ˜ n = lim n u n y n = 0 .
Moreover, noticing x n k z and (21), we obtain that u n k z and w n k z . By Lemma 12, we deduce that z ω w ( u n ) VI ( C , A 1 ) and z ω w ( w n ) VI ( C , A 2 ) . Consequently,
z i = 1 2 VI ( C , A i ) .
Next, we claim that z i = 0 N Fix ( S i ) . Indeed, by (21), we immediately obtain
x n + 1 x n x n + 1 w n + w n u n + u n x n 0 ( n ) .
We first claim that lim n x n S r x n = 0 r { 1 , , N } . Indeed, by the formulation of S n , one obtains S n { S 1 , , S N } , which, hence, leads to S n + k { S 1 , , S N } n 1 , k { 1 , , N } . Note that
x n S n + k x n x n x n + k + x n + k S n + k x n + k + S n + k x n + k S n + k x n x n x n + k + x n + k S n + k x n + k + m = 1 N S m x n + k S m x n .
Utilizing the uniform continuity of each S m , one deduces from (21) and (27) that k , m { 1 , , N } , x n + k S n + k x n + k 0 and S m x n + k S m x n 0 . Hence, one obtains lim n x n S n + k x n = 0 . So, it follows that
lim n S r x n x n = 0 r { 1 , , N } .
This along with x n k z , leads to z Fix ^ ( S r ) = Fix ( S r ) r { 1 , , N } . Accordingly, z k = 1 N Fix ( S k ) . Also, using x n k z and (21), one has w n k z . As a result, from (24), we obtain z Fix ^ ( S 0 ) = Fix ( S 0 ) . Therefore, z k = 0 N Fix ( S k ) , and thus z Ω . Consequently, ω w ( x n ) Ω . Hence, by Lemma 5, one concludes that { x n } converges weakly to z.    □
On the other hand, let us show the strongly convergent result for a pair of VIPs with CFPP constraint, where the two VIPs implicate two mappings of both uniform continuity and pseudomonotonicity and the CFPP involves finite mappings of Bregman’s relative nonexpansivity and a mapping of Bregman’s relative demicontractivity.
Theorem 2.
Suppose that the conditions (C1)–(C3) hold. If { x n } is the constructed sequence in Algorithm 4, then x n Π Ω u sup n 0 x n < .
Proof. 
It is clear that the necessity of the theorem is true. Next, it suffices to show that the sufficiency is valid. Assume that sup n 0 x n < . In what follows, we divide our proof into four claims.
Claim 1. We show that
( 1 α n ) γ n ( 1 γ n ) ρ b * J E p x n J E p g n α n D f p ( u ^ , u ) + D f p ( z , x n ) D f p ( u ^ , x n + 1 ) + n M ,
for some M > 0 . Indeed, put u ^ = Π Ω u . Noticing w n = Π C n u n and z n = Π C ˜ n w n , we deduce from (5) and (7) that
D f p ( u ^ , w n ) D f p ( u ^ , u n ) D f p ( w n , u n ) D f p ( u ^ , u n ) τ [ dist ( C n , u n ) ] p ,
and
D f p ( u ^ , z n ) D f p ( u ^ , w n ) D f p ( z n , w n ) D f p ( u ^ , w n ) τ [ dist ( C ˜ n , w n ) ] p .
Using the same inferences as in the proof of Theorem 1, we know that
D f p ( u ^ , g n ) D f p ( u ^ , x n ) + ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) × u ^ + S n x n 1 2 S n x n D f p ( u ^ , x n ) + n M ,
where sup n 1 u ^ + S n x n 1 2 S n x n M for some M > 0 . This ensures that { g n } is bounded.
Algorithm 4: The 2nd modified inertial-type subgradient extragradient method
  The 2nd modified inertial-type subgradient extragradient method. Initialization: Given x 0 , x 1 C arbitrarily and let ϵ > 0 , l k ( 0 , 1 ) , μ k > 0 and λ k ( 0 , 1 μ k ) for k = 1 , 2 . Choose { n } , { γ n } , { α n } ( 0 , 1 ) and { β n } ( ξ , 1 ) s.t.
lim n n = 0 , n = 1 α n = , lim n α n = 0 , lim inf n ( β n ξ ) ( 1 β n ) > 0 and lim inf n γ n ( 1 γ n ) > 0 . Moreover, given the iterates x n 1 and x n ( n 1 ) , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where sup n 1 ϵ n α n < and
ϵ n ¯ = min { ϵ , n J E p S n x n J E p ( 2 S n x n S n x n 1 ) } if S n x n S n x n 1 , ϵ otherwise .
Iterations: Compute x n + 1 below:
Step 1. Put g n = J E * q ( ( 1 ϵ n ) J E p S n x n + ϵ n J E p ( 2 S n x n S n x n 1 ) ) , and calculate u n = J E * q ( γ n J E p x n + ( 1 γ n ) J E p g n ) , y n = Π C ( J E * q ( J E p u n λ 1 A 1 u n ) ) , r λ 1 ( u n ) : = u n y n and s n = u n τ n r λ 1 ( u n ) , with τ n : = l 1 k n and k n being the smallest k 0 s.t.
A 1 u n A 1 ( u n l 1 k r λ 1 ( u n ) ) , u n y n μ 1 2 D f p ( u n , y n ) .
Step 2. Calculate w n = Π C n ( u n ) , with C n : = { y C : h n ( y ) 0 } and
h n ( y ) = A 1 s n , y u n + τ n 2 λ 1 D f p ( u n , y n ) .
Step 3. Calculate y ˜ n = Π C ( J E * q ( J E p w n λ 2 A 2 w n ) ) , r λ 2 ( w n ) : = w n y ˜ n and t n = w n τ ˜ n r λ 2 ( w n ) , with τ ˜ n : = l 2 j n and j n being the smallest j 0 s.t.
A 2 w n A 2 ( w n l 2 j r λ 2 ( w n ) ) , w n y ˜ n μ 2 2 D f p ( w n , y ˜ n ) .
Step 4. Set z n = Π C ˜ n ( w n ) , and calculate v n = J E * q ( β n J E p z n + ( 1 β n ) J E p ( S 0 z n ) ) and x n + 1 = Π C ( J E * q ( α n J E p u + ( 1 α n ) J E p v n ) , with C ˜ n : = { y C : h ˜ n ( y ) 0 } and
h ˜ n ( y ) = A 2 t n , y w n + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) .
Again set n : = n + 1 and return to Step 1.
Using (10) and the last two inequalities, from { γ n } ( 0 , 1 ) and { β n } ( ξ , 1 ) , we obtain
D f p ( u ^ , x n + 1 ) α n D f p ( u ^ , u ) + ( 1 α n ) D f p ( u ^ , v n ) α n D f p ( u ^ , u ) + ( 1 α n ) [ β n D f p ( u ^ , z n ) + ( 1 β n ) D f p ( u ^ , S 0 z n ) β n ( 1 β n ) ρ b z n * z n S 0 z n ] α n D f p ( u ^ , u ) + ( 1 α n ) [ D f p ( u ^ , z n ) ( β n ξ ) ( 1 β n ) ρ b z n * z n S 0 z n ] α n D f p ( u ^ , u ) + ( 1 α n ) D f p ( u ^ , u n ) α n D f p ( u ^ , u ) + ( 1 α n ) [ γ n D f p ( u ^ , x n ) + ( 1 γ n ) D f p ( u ^ , g n ) γ n ( 1 γ n ) ρ b * J E p x n J E p g n ] α n D f p ( u ^ , u ) + D f p ( z , x n ) + n M ( 1 α n ) γ n ( 1 γ n ) ρ b * J E p x n J E p g n ,
which immediately arrives at the desired claim. In addition, it is easily known that { u n } , { v n } , { w n } , { y n } , { y ˜ n } , { z n } , { s n } , { t n } and { S 0 z n } are also bounded.
   Claim 2. We show that
D f p ( w n , u n ) + D f p ( z n , w n ) D f p ( u ^ , u n ) D f p ( u ^ , x n + 1 ) + α n J E p u J E p u ^ , ζ n u ^ .
Indeed, take b = sup n 1 { x n p 1 , g n p 1 , z n p 1 , S 0 z n p 1 } . Using Lemma 2, one obtains
D f p ( u ^ , u n ) = V f p ( u ^ , γ n J E p x n + ( 1 γ n ) J E p g n ) 1 p u ^ p γ n J E p x n , u ^ ( 1 γ n ) J E p g n , u ^ + γ n q J E p x n q + ( 1 γ n ) q J E p g n q γ n ( 1 γ n ) ρ b * J E p x n J E p g n = 1 p u ^ p γ n J E p x n , u ^ ( 1 γ n ) J E p g n , u ^ + γ n q x n p + ( 1 γ n ) q g n p γ n ( 1 γ n ) ρ b * J E p x n J E p g n = γ n D f p ( u ^ , x n ) + ( 1 γ n ) D f p ( u ^ , g n ) γ n ( 1 γ n ) ρ b * J E p x n J E p g n γ n D f p ( u ^ , x n ) + ( 1 γ n ) [ D f p ( u ^ , x n ) + n M ] γ n ( 1 γ n ) ρ b * J E p x n J E p g n D f p ( u ^ , x n ) + n M γ n ( 1 γ n ) ρ b * J E p x n J E p g n ,
and
D f p ( u ^ , v n ) = V f p ( u ^ , β n J E p z n + ( 1 β n ) J E p S 0 z n ) β n D f p ( u ^ , z n ) + ( 1 β n ) D f p ( u ^ , S 0 z n ) β n ( 1 β n ) ρ b * J E p z n J E p S 0 z n β n D f p ( u ^ , z n ) + ( 1 β n ) [ D f p ( u ^ , z n ) + ξ ρ b * J E p z n J E p S 0 z n ] β n ( 1 β n ) ρ b * J E p z n J E p S 0 z n = D f p ( u ^ , z n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n D f p ( u ^ , w n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n .
Set ζ n = J E * q ( α n J E p u + ( 1 α n ) J E p v n ) . From (9), we have
D f p ( u ^ , x n + 1 ) V f p ( u ^ , α n J E p u + ( 1 α n ) J E p v n ) V f p ( u ^ , α n J E p u + ( 1 α n ) J E p v n α n ( J E p u J E p u ^ ) ) + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) D f p ( u ^ , v n ) + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) [ D f p ( u ^ , w n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n ] + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) D f p ( u ^ , w n ) + α n J E p u J E p u ^ , ζ n u ^ .
Furthermore, from (29), one has
D f p ( u ^ , v n ) D f p ( u ^ , z n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n D f p ( u ^ , w n ) D f p ( z n , w n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n D f p ( u ^ , w n ) D f p ( z n , w n ) .
This, along with (30), leads to
D f p ( u ^ , x n + 1 ) ( 1 α n ) D f p ( u ^ , v n ) + α n J E p u J E p u ^ , ζ n u ^ D f p ( u ^ , w n ) D f p ( z n , w n ) + α n J E p u J E p u ^ , ζ n u ^ D f p ( u ^ , u n ) D f p ( w n , u n ) D f p ( z n , w n ) + α n J E p u J E p u ^ , ζ n u ^ .
Consequently,
D f p ( w n , u n ) + D f p ( z n , w n ) D f p ( u ^ , u n ) + α n J E p u J E p u ^ , ζ n u ^ D f p ( u ^ , x n + 1 ) .
Claim 3. Let us show
α n D f p ( u ^ , u ) + D f p ( u ^ , x n ) D f p ( u ^ , x n + 1 ) + n M ( 1 α n ) { τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } .
Indeed, by the analogous reasonings to these of (26), one obtains
D f p ( u ^ , z n ) D f p ( u ^ , w n ) τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p D f p ( u ^ , u n ) τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p .
Applying (28), (29) and (32), we have
D f p ( u ^ , x n + 1 ) α n D f p ( u ^ , u ) + ( 1 α n ) D f p ( u ^ , v n ) α n D f p ( u ^ , u ) + ( 1 α n ) [ D f p ( u ^ , z n ) ( β n ξ ) ( 1 β n ) ρ b * z n S 0 z n ] α n D f p ( u ^ , u ) + ( 1 α n ) { D f p ( u ^ , u n ) τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } α n D f p ( u ^ , u ) + D f p ( u ^ , u n ) ( 1 α n ) { τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } α n D f p ( u ^ , u ) + D f p ( u ^ , x n ) + n M γ n ( 1 γ n ) ρ b * J E p x n J E p g n ( 1 α n ) { τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } α n D f p ( u ^ , u ) + D f p ( u ^ , x n ) + n M ( 1 α n ) { τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } .
Claim 4. We show that x n u ^ as n . Indeed, since E is reflexive and { x n } is bounded, we know that ω w ( x n ) is nonempty. Let z ω w ( x n ) . Whereby, there exists a subsequence { x n k } of { x n } such that { x n k } converges weakly to z. One defines Γ n : = D f p ( u ^ , x n ) n . In what follows, let us demonstrate Γ n 0 ( n ) in both possible aspects.
    Aspect 1. Presume that there exists n 0 1 s.t. { Γ n } is non-increasing. Whereby, Γ n d < + and, hence, Γ n Γ n + 1 0 . From (28) and (31), we obtain
D f p ( w n , u n ) + D f p ( z n , w n ) D f p ( u ^ , u n ) + α n J E p u J E p u ^ , ζ n u ^ D f p ( u ^ , x n + 1 ) D f p ( u ^ , x n ) + n M γ n ( 1 γ n ) ρ b * J E p x n J E p g n D f p ( u ^ , x n + 1 ) + α n J E p u J E p u ^ , ζ n u ^ ,
which, hence, yields
D f p ( w n , u n ) + D f p ( z n , w n ) + γ n ( 1 γ n ) ρ b * J E p x n J E p g n Γ n Γ n + 1 + n M + α n J E p u J E p u ^ , ζ n u ^ .
Because n 0 , α n 0 , 0 < lim inf n γ n ( 1 γ n ) , Γ n Γ n + 1 0 and the sequence { ζ n } is of boundedness, one deduces that D f p ( w n , u n ) 0 , D f p ( z n , w n ) 0 and ρ b * J E p x n J E p g n 0 , which, hence, yields J E p x n J E p g n 0 . From u n = J E * q ( γ n J E p x n + ( 1 γ n ) J E p g n ) , it is easily known that lim n J E p u n   J E p x n = 0 . Noticing g n = J E * q ( ( 1 ϵ n ) J E p S n x n + ϵ n J E p ( 2 S n x n S n x n 1 ) ) , we deduce from lim n n = 0 and the definition of ϵ n that
J E p g n J E p S n x n = ϵ n J E p ( 2 S n x n S n x n 1 ) J E p S n x n n 0 ( n ) .
Hence, using (5) and uniform continuity of J E p on bounded subsets of E, we conclude that lim n g n x n = 0 and
lim n w n u n = lim n z n w n = lim n x n S n x n = lim n u n x n = 0 .
Furthermore, from (28) and (30), we have
( 1 α n ) ( β n ξ ) ( 1 β n ) ρ b * J E p z n J E p S 0 z n D f p ( u ^ , u n ) D f p ( u ^ , x n + 1 ) + α n J E p u J E p u ^ , ζ n u ^ D f p ( u ^ , x n ) D f p ( u ^ , x n + 1 ) + n M + α n J E p u J E p u ^ , ζ n u ^ .
According to the analogous reasonings, one obtains lim n J E p z n J E p S 0 z n = 0 , which, hence, leads to lim n J E p v n J E p z n = 0 (due to v n = J E * q ( β n J E p z n + ( 1 β n ) J E p S 0 z n ) ). Using uniform continuity of J E * q on any bounded subset of E * , we obtain
lim n z n S 0 z n = lim n v n z n = 0 .
This together with (34) implies that
v n x n v n z n + z n w n + w n u n + u n x n 0 ( n ) .
It is clear that
lim n z n x n = 0 .
Let us show that z i = 0 N Fix ( S i ) . Indeed, since ζ n = J E * q ( α n J E p u + ( 1 α n ) J E p v n ) , it can be readily seen that
lim n ζ n x n = 0 .
In addition, using (7), (28) and (29), we have
D f p ( u ^ , x n + 1 ) D f p ( u ^ , J E * q ( α n J E p u + ( 1 α n ) J E p v n ) D f p ( x n + 1 , ζ n ) α n D f p ( u ^ , u ) + D f p ( u ^ , w n ) D f p ( x n + 1 , ζ n ) α n D f p ( u ^ , u ) + D f p ( u ^ , u n ) D f p ( x n + 1 , ζ n ) α n D f p ( u ^ , u ) + D f p ( u ^ , x n ) + n M D f p ( x n + 1 , ζ n ) ,
which, hence, arrives at
D f p ( x n + 1 , ζ n ) α n D f p ( u ^ , u ) + Γ n Γ n + 1 + n M .
So, it follows that lim n D f p ( x n + 1 , ζ n ) = 0 and thus x n + 1 ζ n 0 . This along with (37) arrives at
x n + 1 x n x n + 1 ζ n + ζ n x n 0 ( n ) .
Observe that k { 1 , . . . , N } ,
x n S n + k x n x n x n + k + x n + k S n + k x n + k + S n + k x n + k S n + k x n x n x n + k + x n + k S n + k x n + k + m = 1 N S m x n + k S m x n .
Exploiting the uniform continuity of each S m , one deduces from (34) and (38) that k , m { 1 , , N } , x n + k S n + k x n + k 0 and S m x n + k S m x n 0 . Hence, one has that x n S n + k x n 0 k { 1 , , N } . So, it follows that x n S r x n 0 r { 1 , , N } . This along with x n k z , attains z Fix ^ ( S r ) = Fix ( S r ) r { 1 , . . . , N } . Consequently, z lies in k = 1 N Fix ( S k ) . Additionally, from (36) and x n k z , one has that z n k z . Thus, using (35), we obtain z Fix ^ ( S 0 ) = Fix ( S 0 ) . Consequently, z k = 0 N Fix ( S k ) ,
In what follows, we show that z i = 1 2 VI ( C , A i ) . From (33), we have
( 1 α n ) { τ [ τ n 2 λ 1 L 1 D f p ( u n , y n ) ] p + τ [ τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) ] p } α n D f p ( u ^ , u ) + Γ n Γ n + 1 + n M .
So, it follows that lim n τ n 2 λ 1 L 1 D f p ( u n , y n ) = lim n τ ˜ n 2 λ 2 L 2 D f p ( w n , y ˜ n ) = 0 , and, hence,
lim n τ n D f p ( u n , y n ) = lim n τ ˜ n D f p ( w n , y ˜ n ) = 0 .
By Lemma 13, one has that
lim n u n y n = lim n w n y ˜ n = 0 .
Applying (40) and Lemma 12, one obtains that z lies in i = 1 2 VI ( C , A i ) . Thus, we obtain ω w ( x n ) i = 1 2 VI ( C , A i ) . Consequently, ω w ( x n ) Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 N Fix ( S i ) ) . Lastly, we show that lim sup n J E p u J E p u ^ , ζ n u ^ 0 . We can pick a subsequence { x n j } of { x n } such that
lim sup n J E p u J E p u ^ , x n u ^ = lim j J E p u J E p u ^ , x n j u ^ .
Because E is reflexive and { x n } is bounded, we might assume that { x n j } converges weakly to z ˜ . This, along with z ˜ Ω and (6), arrives at
lim sup n J E p u J E p u ^ , x n u ^ = lim j J E p u J E p u ^ , x n j u ^ = J E p u J E p u ^ , z ˜ u ^ 0 .
In terms of (37), we have
lim sup n J E p u J E p u ^ , ζ n u ^ 0 .
From (28) and (30), we obtain
D f p ( u ^ , x n + 1 ) ( 1 α n ) D f p ( u ^ , w n ) + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) D f p ( u ^ , u n ) + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) [ D f p ( u ^ , x n ) + ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) × u ^ + S n x n 1 2 S n x n ] + α n J E p u J E p u ^ , ζ n u ^ ( 1 α n ) D f p ( u ^ , x n ) + ϵ n J E p S n x n J E p ( 2 S n x n S n x n 1 ) × u ^ + S n x n 1 2 S n x n + α n J E p u J E p u ^ , ζ n u ^ = ( 1 α n ) D f p ( u ^ , x n ) + α n { ϵ n α n J E p S n x n J E p ( 2 S n x n S n x n 1 ) × u ^ + S n x n 1 2 S n x n + J E p u J E p u ^ , ζ n u ^ } .
Using uniform continuity of each S i ( 1 i N ) on C and uniform continuity of J E p on bounded subsets of E, from (38) and the boundedness of { x n } , we obtain
lim n J E p S n x n J E p ( 2 S n x n S n x n 1 ) u ^ + S n x n 1 2 S n x n = 0 .
Noticing sup n 1 ϵ n α n < and lim sup n J E p u J E p u ^ , ζ n u ^ 0 , we infer that
lim sup n { ϵ n α n J E p S n x n J E p ( 2 S n x n S n x n 1 ) × u ^ + S n x n 1 2 S n x n + J E p u J E p u ^ , ζ n u ^ } 0 .
Because { α n } lies in ( 0 , 1 ) and n = 1 α n diverges, by the application of Lemma 8 to (42), one obtains D f p ( u ^ , x n ) 0 , and, hence, x n u ^ .
    Aspect 2. Presume that { Γ n k } { Γ n } s.t. Γ n k < Γ n k + 1 for all k 1 . Let the N indicate the natural-number set and ψ : N N be formulated below
ψ ( n ) : = max { k n : Γ k < Γ k + 1 } .
In terms of Lemma 7, one obtains
max { Γ ψ ( n ) , Γ n } Γ ψ ( n ) + 1 .
From (28) and (31), it follows that
D f p ( w ψ ( n ) , u ψ ( n ) ) + D f p ( z ψ ( n ) , w ψ ( n ) ) + γ ψ ( n ) ( 1 γ ψ ( n ) ) ρ b * J E p x ψ ( n ) J E p g ψ ( n ) Γ ψ ( n ) Γ ψ ( n ) + 1 + ψ ( n ) M + α ψ ( n ) J E p u J E p u ^ , ζ ψ ( n ) u ^ .
Noticing g ψ ( n ) = J E * q ( ( 1 ϵ ψ ( n ) ) J E p S ψ ( n ) x ψ ( n ) + ϵ ψ ( n ) J E p ( 2 S ψ ( n ) x ψ ( n ) S ψ ( n ) x ψ ( n ) 1 ) ) and u ψ ( n ) = J E * q ( γ ψ ( n ) J E p x ψ ( n ) + ( 1 γ ψ ( n ) ) J E p g ψ ( n ) ) ) , we obtain that lim n g ψ ( n ) x ψ ( n ) = 0 and
lim n w ψ ( n ) u ψ ( n ) = lim n z ψ ( n ) w ψ ( n ) = lim n x ψ ( n ) S ψ ( n ) x ψ ( n ) = lim n u ψ ( n ) x ψ ( n ) = 0 .
Also, from (28) and (30), one has
( 1 α ψ ( n ) ) ( β ψ ( n ) ξ ) ( 1 β ψ ( n ) ) ρ b * J E p z ψ ( n ) J E p S 0 z ψ ( n ) Γ ψ ( n ) Γ ψ ( n ) + 1 + ψ ( n ) M + α ψ ( n ) J E p u J E p u ^ , ζ ψ ( n ) u ^ .
Noticing v ψ ( n ) = J E * q ( β ψ ( n ) J E p z ψ ( n ) + ( 1 β ψ ( n ) ) J E p S 0 z ψ ( n ) ) and using the analogous reasonings to those in Aspect 1, one obtains
lim n z ψ ( n ) S 0 z ψ ( n ) = lim n v ψ ( n ) z ψ ( n ) = 0 .
This together with (44) implies that
lim n v ψ ( n ) x ψ ( n ) = lim n z ψ ( n ) x ψ ( n ) = 0 .
Noticing ζ ψ ( n ) = J E * q ( α ψ ( n ) J E p u + ( 1 α ψ ( n ) ) J E p v ψ ( n ) ) , from (45), we obtain
lim n ζ ψ ( n ) x ψ ( n ) = 0 .
Using the analogous reasonings to those in Aspect 1, one obtains x ψ ( n ) + 1 x ψ ( n ) 0 ,
lim n u ψ ( n ) y ψ ( n ) = lim n w ψ ( n ) y ˜ ψ ( n ) = 0 ,
and
lim sup n J E p u J E p u ^ , ζ ψ ( n ) u ^ 0 .
Using (42), we obtain
D f p ( u ^ , x ψ ( n ) + 1 ) ( 1 α ψ ( n ) ) D f p ( u ^ , x ψ ( n ) ) + α ψ ( n ) { ϵ ψ ( n ) α ψ ( n ) J E p S ψ ( n ) x ψ ( n ) J E p ( 2 S ψ ( n ) x ψ ( n ) S ψ ( n ) x ψ ( n ) 1 ) × u ^ + S ψ ( n ) x ψ ( n ) 1 2 S ψ ( n ) x ψ ( n ) + J E p u J E p u ^ , ζ ψ ( n ) u ^ } ,
which, together with (43), hence yields
Γ ψ ( n ) ϵ ψ ( n ) α ψ ( n ) J E p S ψ ( n ) x ψ ( n ) J E p ( 2 S ψ ( n ) x ψ ( n ) S ψ ( n ) x ψ ( n ) 1 ) × u ^ + S ψ ( n ) x ψ ( n ) 1 2 S ψ ( n ) x ψ ( n ) + J E p u J E p u ^ , ζ ψ ( n ) u ^ .
As a result, from (48), we deduce that
lim n Γ ψ ( n ) = 0 .
From (48)–(50), one has that
lim n Γ ψ ( n ) + 1 = 0 .
Again, using (43), one obtains Γ n 0 . Thus, x n u ^ . This completes the proof.    □
Remark 1.
It can be easily seen from the proof of Theorem 2 that, if the assumption that lim n n α n = 0 is used in place of the one that lim n n = 0 and sup n 1 ϵ n α n < , then Theorem 2 is still valid. It can be easily seen that the existing method in [25] is most closely related to our proposed method; that is, the hybrid projection method for resolving a single VIP with FPP constraint in [25] is extended to develop our modified inertial-type subgradient extragradient method for resolving a pair of VIPs with CFPP constraint. Compared with the corresponding results in [25], our results exhibit the novelty below:
First, the problem of finding a solution of a single VIP with FPP constraint (involving a Bregman relatively nonexpansive mapping) in [25] is extended to develop our problem of finding a solution of a pair of VIPs with CFPP constraint (involving finite Bregman relatively nonexpansive mappings and a Bregman relatively demicontractive mapping).
Second, the hybrid projection method with line-search process in [25] is extended to develop our modified inertial-type subgradient extragradient method with line-search process.
If we set A 2 = 0 , then Algorithm 3 reduces to the iterative algorithm for finding an element of Ω = VI ( C , A 1 ) ( i = 0 N Fix ( S i ) ) .
Corollary 1.
Assume the conditions (C1)–(C3) with A 2 = 0 , hold, and Ω = VI ( C , A 1 ) ( i = 0 N Fix ( S i ) ) . If { x n } is the fabricated sequence in Algorithm 5, and then x n z Ω sup n 0 x n < .
Algorithm 5: The 3rd modified inertial-type subgradient extragradient method
  The 3rd modified inertial-type subgradient extragradient method. Initialization: Given x 0 , x 1 C arbitrarily and let ϵ > 0 , μ 1 > 0 , λ 1 ( 0 , 1 μ 1 ) , l 1 ( 0 , 1 ) . Choose { n } , { β n } ( 0 , 1 ) and { α n } ( ξ , 1 ) s.t. n = 1 n < , lim inf n β n ( 1 β n ) > 0 and lim inf n ( α n ξ ) ( 1 α n ) > 0 . Moreover, given the iterates x n 1 and x n ( n 1 ) , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where
ϵ n ¯ = min { ϵ , n J E p S n x n J E p ( 2 S n x n S n x n 1 ) } if S n x n S n x n 1 , ϵ otherwise .
Iterative steps: Calculate x n + 1 as follows:
Step 1. Set g n = J E * q ( ( 1 ϵ n ) J E p S n x n + ϵ n J E p ( 2 S n x n S n x n 1 ) ) , and calculate u n = J E * q ( β n J E p x n + ( 1 β n ) J E p g n ) , y n = Π C ( J E * q ( J E p u n λ 1 A 1 u n ) ) , r λ 1 ( u n ) : = u n y n and s n = u n τ n r λ 1 ( u n ) , with τ n : = l 1 k n and k n being the smallest k 0 s.t.
A 1 u n A 1 ( u n l 1 k r λ 1 ( u n ) ) , u n y n μ 1 2 D f p ( u n , y n ) .
Step 2. Calculate w n = Π C n ( u n ) , with C n : = { y C : h n ( y ) 0 } and
h n ( y ) = A 1 s n , y u n + τ n 2 λ 1 D f p ( u n , y n ) .
Step 3. Calculate v n = J E * q ( α n J E p w n + ( 1 α n ) J E p ( S 0 w n ) ) and x n + 1 = Π Q n ( w n ) , with Q n : = { y C : D f p ( y , v n ) D f p ( y , w n ) } .
Again put n : = n + 1 and return to Step 1.
Next, let S 1 : E C be a Bregman relatively nonexpansive mapping and S i = S = I the identity mapping of E for i = 2 , , N . Then, we obtain Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 N Fix ( S i ) ) = ( i = 1 2 VI ( C , A i ) ) Fix ( S 1 ) . In this case, Algorithm 4 reduces to the following iterative scheme for solving a pair of VIPs and the FPP of S 1 . By Theorem 2, we obtain the following strong convergence result.
Corollary 2.
Suppose that the condition (C3) holds, and let Ω = ( i = 1 2 VI ( C , A i ) ) Fix ( S 1 ) . For initial x 0 , x 1 C , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where
ϵ n ¯ = min { ϵ , n J E p S 1 x n J E p ( 2 S 1 x n S 1 x n 1 ) } if S 1 x n S 1 x n 1 , ϵ otherwise .
Suppose that { x n } is the sequence constructed by
g n = J E * q ( ( 1 ϵ n ) J E p S 1 x n + ϵ n J E p ( 2 S 1 x n S 1 x n 1 ) ) , u n = J E * q ( γ n J E p x n + ( 1 γ n ) J E p g n ) , y n = Π C ( J E * q ( J E p u n λ 1 A 1 u n ) ) , s n = ( 1 τ n ) u n + τ n y n , w n = Π C n u n , y ˜ n = Π C ( J E * q ( J E p w n λ 2 A 2 w n ) ) , t n = ( 1 τ ˜ n ) w n + τ ˜ n y ˜ n , z n = Π C ˜ n w n , x n + 1 = Π C ( J E * q ( α n J E p u + ( 1 α n ) J E p z n ) n 1 ,
where τ n : = l 1 k n , τ ˜ n : = l 2 j n and k n , j n are the smallest nonnegative integers k and j satisfying, respectively,
A 1 u n A 1 ( u n l 1 k ( u n y n ) ) , u n y n μ 1 2 D f p ( u n , y n ) ,
A 2 w n A 2 ( w n l 2 j ( w n y ˜ n ) ) , w n y ˜ n μ 2 2 D f p ( w n , y ˜ n ) ,
and the sets C n , C ˜ n are constructed below:
(i) 
C n : = { y C : h n ( y ) 0 } and h n ( y ) = A 1 s n , y u n + τ n 2 λ 1 D f p ( u n , y n ) ;
(ii) 
C ˜ n : = { y C : h ˜ n ( y ) 0 } and h ˜ n ( y ) = A 2 t n , y w n + τ ˜ n 2 λ 2 D f p ( w n , y ˜ n ) .
Then, x n Π Ω u sup n 0 x n < .

4. Examples

In what follows, we furnish an illustrative instance to back up the practicability and implementability of the suggested approaches. Put ϵ = 1 3 , μ i = 1 and l i = λ i = 1 3 for i = 1 , 2 . We first provide an example of uniformly continuous and pseudomonotone mappings A i : E E * , i = 1 , 2 , Bregman relatively nonexpansive mapping S 1 : C C and Bregman relatively demicontractive mapping S 0 : C C with Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 1 Fix ( S i ) ) . Let C = [ 2 , 2 ] and E = H = R with the inner product a , b = a b and induced norm · = | · | . The initial points x 0 , x 1 are randomly chosen in C. For i = 1 , 2 , let A i : H H be defined as A 1 x : = 1 1 + | sin x | 1 1 + | x | and A 2 x : = x + sin x for all x H . Now, we first show that A 1 is Lipschitz continuous and pseudomonotone. In fact, for each y , g H , one has
A 1 y A 1 g = | 1 1 + sin y 1 1 + y 1 1 + sin g + 1 1 + g | | g y ( 1 + y ) ( 1 + g ) | + | sin g sin y ( 1 + sin y ) ( 1 + sin g ) | y g + sin y sin g 2 y g .
This ensures that A 1 is of Lipschitz continuity. Also, let us show that A 1 is of pseudomonotonicity. For each y , g H , it can be easily seen that
A 1 y , g y = ( 1 1 + | sin y | 1 1 + | y | ) ( g y ) 0 A 1 g , g y = ( 1 1 + | sin g | 1 1 + | g | ) ( g y ) 0 .
It is readily known that A 2 is Lipschitz continuous and monotone. Indeed, we deduce that A 2 x A 2 y x y + sin x sin y 2 x y and
A 2 x A 2 y , x y = x y 2 + sin x sin y , x y x y 2 x y 2 = 0 .
Now, let S 1 : C C and S 0 : C C be defined as S 1 x = sin x and S 0 x = 1 5 x + 3 5 sin x . It is easy to verify that Fix ( S 1 ) = Fix ( S 0 ) = { 0 } and S 1 : C C is Bregman relatively nonexpansive. Also, S 0 : C C is Bregman relatively ξ -demicontractive with ξ = 1 5 . Indeed, note that
S 0 x S 0 y 2 = 1 5 ( x y ) + 3 5 ( sin x sin y ) 2 x y 2 + 2 5 ( I S 0 ) x ( I S 0 ) y 2 .
Consequently,
Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 1 Fix ( S i ) ) = { 0 } .
In addition, putting β n = n + 2 2 ( n + 1 ) n 1 , we obtain
lim n ( β n ξ ) ( 1 β n ) = lim n ( n + 2 2 ( n + 1 ) 1 5 ) ( 1 n + 2 2 ( n + 1 ) ) = ( 1 2 1 5 ) ( 1 1 2 ) = 3 20 > 0
In this case, the conditions (C1)–(C3) are satisfied.
Example 1.
Let n = 1 2 ( n + 1 ) 2 and α n = β n = n + 2 2 ( n + 1 ) n 1 . Given the iterates x n 1 and x n ( n 1 ) , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where
ϵ n ¯ = min { ϵ , n S 1 x n S 1 x n 1 } if S 1 x n S 1 x n 1 , ϵ otherwise .
Algorithm 3 is rewritten as follows:
g n = S 1 x n + ϵ n ( S 1 x n S 1 x n 1 ) , u n = n + 2 2 ( n + 1 ) x n + n 2 ( n + 1 ) g n , y n = P C ( u n 1 3 A 1 u n ) , s n = ( 1 τ n ) u n + τ n y n , w n = P C n u n , y ˜ n = P C ( w n 1 3 A 2 w n ) , t n = τ ˜ n y ˜ n + ( 1 τ ˜ n ) w n , v n = n 2 ( n + 1 ) S 0 w n + n + 2 2 ( n + 1 ) w n , Q n = { y C : y v n y w n } , x n + 1 = P C ˜ n Q n w n n 1 ,
with the sets C n , C ˜ n and the step-sizes τ n , τ ˜ n being picked as in Algorithm 3 for each n. Then, by Theorem 1, we deduce that { x n } converges to 0 Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 1 Fix ( S i ) ) .
Example 2.
Let n = 1 2 ( n + 1 ) 2 , α n = 1 2 ( n + 1 ) and β n = γ n = n + 2 2 ( n + 1 ) n 1 . Given the iterates x n 1 and x n ( n 1 ) , choose ϵ n s.t. 0 ϵ n ϵ n ¯ , where
ϵ n ¯ = min { ϵ , n S 1 x n S 1 x n 1 } if S 1 x n S 1 x n 1 , ϵ otherwise .
Algorithm 4 is rewritten as follows:
g n = S 1 x n + ϵ n ( S 1 x n S 1 x n 1 ) , u n = n + 2 2 ( n + 1 ) x n + n 2 ( n + 1 ) g n , y n = P C ( u n 1 3 A 1 u n ) , s n = ( 1 τ n ) u n + τ n y n , w n = P C n u n , y ˜ n = P C ( w n 1 3 A 2 w n ) , t n = ( 1 τ ˜ n ) w n + τ ˜ n y ˜ n , z n = P C ˜ n w n , v n = n + 2 2 ( n + 1 ) z n + n 2 ( n + 1 ) S 0 z n , x n + 1 = P C ( 1 2 ( n + 1 ) u + 2 n + 1 2 ( n + 1 ) v n ) ,
with the sets C n , C ˜ n and the step-sizes τ n , τ ˜ n are picked as in Algorithm 4 for each n. Then, by Theorem 2, we deduce that { x n } converges to 0 Ω = ( i = 1 2 VI ( C , A i ) ) ( i = 0 1 Fix ( S i ) ) .
It is worthy to emphasize that the above numerical example shows the competitive advantage of our suggested algorithms over the existing algorithms, e.g., the ones in [25]. In fact, we have provided the illustrative example of a pair of VIPs with CFPP constraint as above. Note that the existing method in [25] is only utilized for solving a single VIP with an FPP constraint. Hence, there is no way for this method to handle the above numerical example; that is, it is invalid for a pair of VIPs with a CFPP constraint. However, our suggested method can settle the above illustrative example. This reveals the competitive advantage of our proposed algorithms over the existing algorithms in the literature.

5. Conclusions

Let 1 < q 2 p < with 1 p + 1 q = 1 and let E be a p-uniformly convex and uniformly smooth Banach space. Then, its dual space E * is q-uniformly smooth Banach space with 1 < q 2 . Utilizing the geometric properties of E and E * , we design two inertial-type subgradient extragradient algorithms with the line-search process for solving the pseudomonotone variational inequality problems (VIPs) and common fixed-point problem (CFPP), where the geometric properties involve the properties of the generalized duality mappings J E p , J E * q and Bregman projection operator Π C . Here, the CFPP indicates the common fixed-point problem of finite Bregman’s relative nonexpansivity mappings and a Bregman’s relative demicontractivity mapping in E. Under the properties of the generalized duality mappings J E p , J E * q and Bregman projection operator Π C , we have proved that the sequences generated by the suggested algorithms converge weakly and strongly to a solution of a pair of VIPs with a CFPP constraint, respectively. Additionally, an illustrated example is furnished to demonstrate the feasibility and implementability of our proposed approaches. Compared with the corresponding results in [25], our results reveal the novelty as follows: (a) the problem of finding a solution of a single VIP with an FPP constraint (involving Bregman relatively nonexpansive mapping) in [25] is extended to develop our problem of finding a solution of a pair of VIPs with a CFPP constraint (involving finite Bregman relatively nonexpansive mappings and a Bregman relatively demicontractive mapping); (b) the hybrid projection method with a line-search process in [25] is extended to develop our modified inertial-type subgradient extragradient method with the line-search process. In the end, it is worthy to point out that part of our future research is aimed at attaining the weak and strong convergence results for the modifications of our proposed approaches with Nesterov double inertial-type extrapolation steps (see [26]) and adaptive stepsizes.

Author Contributions

Conceptualization, B.L., H.-Y.H. and Y.-S.L.; Data curation, B.L.; Formal analysis, C.-S.W., B.L., S.-L.C. and Y.-S.L.; Funding acquisition, L.-C.C.; Investigation, C.-S.W., L.-C.C., B.L., S.-L.C., H.-Y.H. and Y.-S.L.; Methodology, L.-C.C., H.-Y.H. and Y.-S.L.; Project administration, L.-C.C.; Resources, C.-S.W.; Software, B.L. and S.-L.C.; Supervision, L.-C.C.; Validation, Y.-S.L.; Writing—original draft, C.-S.W., L.-C.C. and B.L.; Writing—review & editing, L.-C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the 2020 Shanghai Leading Talents Program of the Shanghai Municipal Human Resources and Social Security Bureau (20LJ2006100), the Innovation Program of Shanghai Municipal Education Commission (15ZZ068) and the Program for Outstanding Academic Leaders in Shanghai City (15XD1503100).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, S. Hybrid extragradient methods for asymptotically strict pseudo-contractions in the intermediate sense and variational inequality problems. Optimization 2011, 60, 739–754. [Google Scholar] [CrossRef]
  2. Migórski, S.; Zeng, S.D. A class of differential hemivariational inequalities in Banach spaces. J. Glob. Optim. 2018, 72, 761–779. [Google Scholar] [CrossRef]
  3. Anh, P.N.; An, L.T.H. New subgradient extragradient methods for solving monotone bilevel equilibrium problems. Optimization 2019, 68, 2097–2122. [Google Scholar] [CrossRef]
  4. Wu, X.K.; Wu, Z.B.; Zou, Y.Z. Sensitivity of the set of solutions for a class of fractional set-valued projected dynamical systems. Nonlinear Anal. Forum 2015, 20, 95–105. [Google Scholar]
  5. Korpelevich, G.M. The extragradient method for finding saddle points and other problems. Ekon. Mat. Metod. 1976, 12, 747–756. [Google Scholar]
  6. Jolaoso, L.O.; Shehu, Y.; Yao, J.C. Inertial extragradient type method for mixed variational inequalities without monotonicity. Math. Comput. Simul. 2022, 192, 353–369. [Google Scholar] [CrossRef]
  7. Cui, Y.L.; Ceng, L.C.; Zhang, F.F.; He, L.; Yin, J.; Wang, C.-S.; Hu, H.-Y. Mann hybrid deepest-descent extragradient method with line-search process for hierarchical variational inequalities for countable nonexpansive mappings. J. Math. 2023, 2023, 6177912. [Google Scholar] [CrossRef]
  8. Chen, J.; Liu, S. Extragradient-like method for pseudomontone equilibrium problems on Hadamard manifolds. J. Inequal. Appl. 2020, 2020, 205. [Google Scholar] [CrossRef]
  9. Censor, Y.; Gibali, A.; Reich, S. The subgradient extragradient method for solving variational inequalities in Hilbert space. J. Optim. Theory Appl. 2011, 148, 318–335. [Google Scholar] [CrossRef]
  10. Yao, Y.; Shahzad, N.; Yao, J.C. Convergence of Tseng-type self-adaptive algorithms for variational inequalities and fixed point problems. Carpathian J. Math. 2021, 37, 541–550. [Google Scholar] [CrossRef]
  11. Ceng, L.C.; Petrusel, A.; Qin, X.; Yao, J.C. On inertial subgradient extragradient rule for monotone bilevel equilibrium problems. Fixed Point Theory 2023, 24, 101–126. [Google Scholar] [CrossRef]
  12. Ceng, L.C.; Zhu, L.J.; Yin, T.C. Modified subgradient extragradient algorithms for systems of generalized equilibria with constraints. AIMS Math. 2023, 8, 2961–2994. [Google Scholar] [CrossRef]
  13. Yang, J.; Liu, H.; Liu, Z. Modified subgradient extragradient algorithms for solving monotone variational inequalities. Optimization 2018, 67, 2247–2258. [Google Scholar] [CrossRef]
  14. Ceng, L.C.; Petrusel, A.; Qin, X.; Yao, J.C. A modified inertial subgradient extragradient method for solving pseudomonotone variational inequalities and common fixed point problems. Fixed Point Theory 2020, 21, 93–108. [Google Scholar] [CrossRef]
  15. Li, H.; Wu, Y.; Wang, F. Inertial algorithms of common solutions for variational inequality, fixed point and equilibrium problems. J. Nonlinear Convex Anal. 2022, 23, 859–877. [Google Scholar]
  16. Ceng, L.C.; Huang, S.; Petrusel, A. Weak convergence theorem by a modified extragradient method for nonexpansive mappings and monotone mappings. Taiwan. J. Math. 2009, 13, 225–238. [Google Scholar] [CrossRef]
  17. Vuong, P.T.; Shehu, Y. Convergence of an extragradient-type method for variational inequality with applications to optimal control problems. Numer. Algorithms 2019, 81, 269–291. [Google Scholar] [CrossRef]
  18. Thong, D.V.; Hieu, D.V. Modified subgradient extragradient method for variational inequality problems. Numer. Algorithms 2018, 79, 597–610. [Google Scholar] [CrossRef]
  19. Kraikaew, R.; Saejung, S. Strong convergence of the Halpern subgradient extragradient method for solving variational inequalities in Hilbert spaces. J. Optim. Theory Appl. 2014, 163, 399–412. [Google Scholar] [CrossRef]
  20. Ceng, L.C.; Wen, C.F. Systems of variational inequalities with hierarchical variational inequality constraints for asymptotically nonexpansive and pseudocontractive mappings. Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Mat. 2019, 113, 2431–2447. [Google Scholar] [CrossRef]
  21. Hu, S.; Wang, Y.; Tan, B.; Wang, F. Inertial iterative method for solving variational inequality problems of pseudo-monotone operators and fixed point problems of nonexpansive mappings in Hilbert spaces. J. Ind. Manag. Optim. 2023, 19, 2655–2675. [Google Scholar] [CrossRef]
  22. Ceng, L.C.; Huang, S. A hybrid extragradient method for asymptotically strict pseudo-contractions in the intermediate sense. Fixed Point Theory 2010, 11, 189–201. [Google Scholar]
  23. Reich, S.; Thong, D.V.; Dong, Q.L.; Li, X.H.; Dung, V.T. New algorithms and convergence theorems for solving variational inequalities with non-Lipschitz mappings. Numer. Algorithms 2021, 87, 527–549. [Google Scholar] [CrossRef]
  24. Iusem, A.N.; Nasri, M. Korpelevich’s method for variational inequality problems in Banach spaces. J. Glob. Optim. 2011, 50, 59–76. [Google Scholar] [CrossRef]
  25. Eskandani, G.Z.; Lotfikar, R.; Raeisi, M. Hybrid projection methods for solving pseudo-monotone variational inequalities in Banach spaces. Fixed Point Theory, 2023; submitted. [Google Scholar]
  26. Yao, Y.; Iyiola, O.S.; Shehu, Y. Subgradient extragradient method with double inertial steps for variational inequalities. J. Sci. Comput. 2022, 90, 71. [Google Scholar] [CrossRef]
  27. Takahashi, Y.; Hashimoto, K.; Kato, M. On sharp uniform convexity, smoothness, and strong type, cotype inequalities. J. Nonlinear Convex Anal. 2002, 3, 267–281. [Google Scholar]
  28. Reem, D.; Reich, S.; Pierro, A.D. Re-examination of Bregman functions and new properties of their divergences. Optimization 2019, 68, 279–348. [Google Scholar] [CrossRef]
  29. Schöpfer, F.; Schuster, T.; Louis, A.K. An iterative regularization method for the solution of the split feasibility problem in Banach spaces. Inverse Probl. 2008, 24, 055008. [Google Scholar] [CrossRef]
  30. Butnariu, D.; Resmerita, E. Bregman distances, totally convex functions, and a method for solving operator equations in Banach spaces. Abstr. Appl. Anal. 2006, 2006, 084919. [Google Scholar] [CrossRef]
  31. Butnariu, D.; Iusem, A.N.; Resmerita, E. Total convexity for powers of the norm in uniformly convex Banach spaces. J. Convex Anal. 2000, 7, 319–334. [Google Scholar]
  32. Reich, S. A weak convergence theorem for the alternating method with Bregman distances. In Theory and Applications of Nonlinear Operators; Marcel Dekker: New York, NY, USA, 1996; pp. 313–318. [Google Scholar]
  33. Cottle, R.W.; Yao, J.C. Pseudo-monotone complementarity problems in Hilbert space. J. Optim. Theory Appl. 1992, 75, 281–295. [Google Scholar] [CrossRef]
  34. He, Y.R. A new double projection algorithm for variational inequalities. J. Comput. Appl. Math. 2006, 185, 166–173. [Google Scholar] [CrossRef]
  35. 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]
  36. Xu, H.K.; Kim, T.H. Convergence of hybrid steepest-descent methods for variational inequalities. J. Optim. Theory Appl. 2003, 119, 185–201. [Google Scholar] [CrossRef]
  37. Osilike, M.O.; Aniagbosor, S.C.; Akuchu, B.G. Fixed points of asymptotically demicontractive mappings in arbitrary Banach spaces. PanAm. Math. J. 2002, 12, 77–88. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.-S.; Ceng, L.-C.; Li, B.; Cao, S.-L.; Hu, H.-Y.; Liang, Y.-S. Modified Inertial-Type Subgradient Extragradient Methods for Variational Inequalities and Fixed Points of Finite Bregman Relatively Nonexpansive and Demicontractive Mappings. Axioms 2023, 12, 832. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12090832

AMA Style

Wang C-S, Ceng L-C, Li B, Cao S-L, Hu H-Y, Liang Y-S. Modified Inertial-Type Subgradient Extragradient Methods for Variational Inequalities and Fixed Points of Finite Bregman Relatively Nonexpansive and Demicontractive Mappings. Axioms. 2023; 12(9):832. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12090832

Chicago/Turabian Style

Wang, Cong-Shan, Lu-Chuan Ceng, Bing Li, Sheng-Long Cao, Hui-Ying Hu, and Yun-Shui Liang. 2023. "Modified Inertial-Type Subgradient Extragradient Methods for Variational Inequalities and Fixed Points of Finite Bregman Relatively Nonexpansive and Demicontractive Mappings" Axioms 12, no. 9: 832. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12090832

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