Home An extragradient inertial algorithm for solving split fixed-point problems of demicontractive mappings, with equilibrium and variational inequality problems
Article Open Access

An extragradient inertial algorithm for solving split fixed-point problems of demicontractive mappings, with equilibrium and variational inequality problems

  • Chibueze C. Okeke , Godwin C. Ugwunnadi and Lateef O. Jolaoso EMAIL logo
Published/Copyright: September 22, 2022
Become an author with De Gruyter Brill

Abstract

The purpose of this article is to study and analyse a new extragradient-type algorithm with an inertial extrapolation step for solving split fixed-point problems for demicontractive mapping, equilibrium problem, and pseudomonotone variational inequality problem in real Hilbert spaces. One of the advantages of the proposed algorithm is that a strong convergence result is achieved without a prior estimate of the Lipschitz constant of the cost operator, which is very difficult to find. In addition, the stepsize is generated at each iteration by some simple computations, which allows it to be easily implemented without the prior knowledge of the Lipschitz constant of the cost operator. Some numerical experiments are reported to show the performance and behaviour of the sequence generated by our algorithm. The obtained results in this article extend and improve many related recent results in this direction in the literature.

MSC 2010: 47H09; 47H10; 49J20; 49J40

1 Introduction

Let H denote a real Hilbert space with inner product , and induced norm . Consider a nonempty, closed, and convex subset C of H and F : C H as a continuous mapping. The variational inequality problem (for short, V I ( C , F ) ) is defined as: find x C such that

(1) F x , y x 0 , y C .

Several applications of V I ( C , F ) (1) can be found in [1,2, 3,4,5, 6,7].

A point x C is a solution of V I ( C , F ) (1) if and only if (see [3] for the details)

x = P C ( x γ F x ) , γ > 0 , and r λ ( x ) x P C ( x γ F x ) = 0 .

This leads to the introduction of the fixed-point approach for solving V I ( C , F ) (1) (see, e.g., [6,8,9]).

If F is η -strongly monotone and L -Lipschitz-continuous, then the sequence generated by the gradient-projection method is given as follows:

(2) x n + 1 = P C ( x n λ F x n )

converges to a solution of V I ( C , F ) (1), provided the step-size λ satisfies λ 0 , 2 η L 2 . However, the gradient-projection method (2) fails to converge to any solution of V I ( C , F ) if F is monotone. For example, take C = R 2 and A a rotation with π 2 angle. Then, F is monotone and L -Lipschitz-continuous, and ( 0 , 0 ) is the unique solution of V I ( C , F ) (1). However, { x n } generated by gradient-projection methods (2) satisfies the property x n + 1 > x n for all n .

In order to address this setback, Korpelevich [10] introduced the extragradient method defined as follows:

(3) x 1 C , y n = P C ( x n λ n F x n ) x n + 1 = P C ( x n λ n F y n ) , n 1 ,

where λ n 0 , 1 L . It is shown in [10] that { x n } converges to a solution of V I ( C , F ) (1) when F is monotone and L -Lipschitz-continuous. In recent years, the extragradient method has been further extended to infinite-dimensional spaces in various ways (see, e.g., [11,12, 13,14,15, 16,17]).

In fact, in the extragradient method, one needs to calculate two projections onto the closed convex set C in each iteration, which seriously affects the efficiency of the algorithm. Censor et al. [18] (see also [19,20]) presented the following subgradient extragradient method (SEGM), which replaces the second projection onto C with a projection onto a specific constructible half-space:

(4) x 1 H , λ > 0 y n = P C ( x n λ F x n ) , T n = { w H x n λ F x n y n , w y n 0 } , x n + 1 = P T n ( x n λ F y n ) , n 1 .

It was shown that (4) converges weakly to a solution of the V I P ( C , F ) if the stepsize satisfies the condition λ 0 , 1 L .

Recently, Gibali [21] introduced a self-adaptive subgradient extragradient method (SA-SEGM) by adopting Armijo-like searches and obtained a weak convergence result for V I P ( C , F ) (1) in R n under the pseudo-monotonicity and continuity of F . Gibali remarked that the Armijo-like searches he adopted can be viewed as a local approximation of the Lipschitz constant of F . Very recently, Shehu and Iyiola [22] proposed an algorithm that combines the viscosity method and the subgradient extragradient method (VSEGM). The VSEGM is of the following form:

Algorithm 1: VSEGM
Initialization : Given l ( 0 , 1 ) , μ ( 0 , 1 ) . Let x 0 H be arbitrary.
Iterative steps: Calculate x n + 1 as follows:
- Step 1. Compute
y n = P C ( x n λ n F x n ) ,
where λ n = l m n where m n is the smallest nonnegative integer of m such that
λ n F x n F y n μ x n y n .
If y n = x n , then stop, and x n is the solution of VIP. Otherwise,
- Step 2. Construct T n = { x H : x n λ n F x n y n , x y n 0 } and compute
z n = P T n ( x n λ n F y n ) .
- Step 3. Compute
x n + 1 = α n f ( x n ) + ( 1 α n ) z n ,
where f : H H is a contraction mapping with constant ρ [ 0 , 1 ) .
Set n n + 1 and back to Step 1.

Under appropriate conditions, they proved that the sequence { x n } generated by Algorithm 1 converges strongly to q V I ( C , F ) , where q = P V I ( C , F ) f ( q ) .

Since the operator F in this article is pseudomonotone, it is worthwhile to provide a detailed example of a map F which is pseudomonotone but not monotone.

Example 1.1

Let H be the Hilbert space L 2 ( [ 0 , 1 ] ) of square-integrable, measurable functions with the inner product x , y = 0 1 x ( t ) y ( t ) d t and the induced norm as follows:

x = x , x .

Take two positive real numbers R and r such that R > r > k k + 1 R for some k > 1 . Consider the operator F : H H defined by F ( x ) = ( R x ) x for each x H , and set Γ = { x H : x r } . Since R > r > k k + 1 R and k > 1 , we obtain immediately that R k + 1 < r k < r . Thus, there exists x 0 Γ such that R k + 1 < x 0 < r k . Take y 0 = k x 0 . Then y 0 = k x 0 < k . r k = r , which implies that y 0 also belongs to Γ . By simple computation, we obtain

F ( x 0 ) F ( y 0 ) , x 0 y 0 = ( 1 k ) 2 x 0 2 ( R ( k + 1 ) x 0 ) .

Thus, by R k + 1 < x 0 , we obtain F ( x 0 ) F ( y 0 ) , x 0 y 0 < 0 . This implies that the operator F is not monotone on Γ . Indeed, F ( x ) , y x 0 for some x , y Γ , i.e. ( R x ) x , y x 0 . Then x , y x 0 because R x > r x 0 . Thus,

F ( y ) , y x = ( R y ) y , y x ( R y ) ( y , y x x , y x ) = ( R y ) y x 2 0

because y Γ and R y > r y 0 . Hence, the operator F is pseudomonotone on Γ .

On the other hand, the split fixed-point problem (SFPP) for mappings T and S , first introduced by Censor and Segal [23], is to find

(5) v F ( S ) such that A v F ( T ) ,

where A : H 1 H 2 is a bounded linear operator, S : H 1 H 1 and T : H 2 H 2 are two mappings with F ( S ) and F ( T ) nonempty, where F ( S ) = { x H 1 : S x = x } and F ( T ) = { x H 2 : T x = x } are the fixed-point sets of T and S , respectively.

More so, let C be a nonempty subset of H and g : C × C R be a bifunction. The equilibrium problem is to find a point v C such that

(6) g ( v , y ) 0

for all y C . The set of all solutions of (6) is denoted by E P ( g ) . Equilibrium problems have been extensively studied by many authors because of their various important applications in nonlinear analysis and optimization, such as Nash equilibrium problems, variational inequalities problem, saddle point problems, and game theory (see [24,25, 26,27] and references therein).

In 2010, Moudafi [37] introduced the following iterative method for solving (5) for two demicontractive mappings:

(7) x 1 H 1 , choose arbitrarily, u n = x n + γ A ( T I ) A x n , x n + 1 = ( 1 α n ) u n + α n S u n , n N .

He proved that the generated sequence { x n } by (7) converges weakly to a solution of SFPP under some control conditions on the sequence { α n } while the step size γ depends on the spectral radius of operator A A .

Later, Shehu and Ogbuisi [28] introduced the following algorithm for solving SFP in real Hilbert spaces:

(8) x 1 H 1 choose arbitrarily, w n = ( 1 β n ) x n , u n = x n + γ A ( T I ) A w n , x n + 1 = ( 1 α n ) u n + α n S u n , n N ,

where T and S are demicontractive mappings. They obtained that the sequence { x n } converges strongly to a solution of the SFP under some control conditions on the sequences { β n } , { α n } , and the stepsize which depends on the norm of A .

It is noted that the control conditions on the step size of algorithms (7) and (8) depend on the spectral radius and the norm of the operators. In many practical problems, computation of the norm of a bounded linear operator is not easy. So it is natural to ask how can we choose a step size that is independent of the norm of operator A ?

Motivated and inspired by the aforementioned works, Hanjing and Suantai [29] modified the algorithms (7) and (8) by applying the inertial term θ n ( x n x n 1 ) to the algorithms for improving the rate of convergence and using a new technique for choosing the control conditions and the step sizes, which do not require any prior knowledge of the operator norm. They proved the weak and strong convergence of the proposed algorithms for approximating solutions of the split fixed-point problem of demicontractive mappings and equilibrium problem in a real Hilbert space. In particular, they proved the following theorem:

Theorem 1.2

[29] Let H 1 and H 2 be real Hilbert spaces, C be a closed convex subset H 1 , A : H 1 H 2 be a bounded linear operator with its adjoint operator A . Let S : H 1 H 1 and T : H 2 H 2 k 2 be k 1 -demicontractive mappings such that S I and T I are demiclosed at zero, g be a bifunction of C × C into R satisfying the conditions ( A 1 ) ( A 4 ) in Lemma 2.3. Suppose F ( S ) A 1 F ( T ) E P ( g ) . Let { x n } be a sequence generated by

  1. Choose x 0 , x 1 H 1 and n = 1 .

  2. Given x n compute u n and y n using

    y n = ( 1 β n ) ( x n + θ n ( x n x n 1 ) ) , u n = T σ n g ( y n + γ n A ( T I ) A y n ) ,

    where { γ n } is bounded and γ n is chosen in such a way that { γ n } is bounded and

    γ n ϵ , ( 1 k ) ( T I ) A w n 2 A ( T I ) A w n 2 ϵ i f n Γ = { κ : ( T I ) A w k 0 } , ϵ > 0 ,

    otherwise γ n = γ ( γ is any nonnegative real number).

  3. Compute

    x n + 1 = ( 1 α n ) u n + α n S u n .

Set n n + 1 and go to Step 2.

The control sequences { θ n } , { β n } , { α n } ( 0 , 1 ) , and { σ n } are real sequences in [ σ , ) ( 0 , ) satisfying the following conditions:

  1. 0 < a α n b < 1 k 2 , for all n N ;

  2. lim n β n = 0 and n = 1 β n = ;

  3. lim n θ n β n x n x n 1 = 0 .

Then { x n } strongly converges to x F ( S ) A 1 F ( T ) E P ( g ) .

By utilizing the idea in [29,30] (for developing a simple self-adaptive method) and coupled with the inertial extrapolation technique known for its sufficiency in improving the convergence rate of algorithms, in this article, we introduced a new inertial extragradient-type algorithm for solving split fixed-point problems for demicontractive mapping, equilibrium problem, and variational inequality with a uniformly continuous pseudomonotone operator.

The remainder of this article is organized as follows. Section 2 recalls some definitions and preliminary results for further use. Section 3 deals with analysing the convergence of the proposed algorithm. In Section 4, we present a numerical experiment that illustrates the performance of our algorithm.

2 Preliminaries

Let H be a real Hilbert space and C be a nonempty closed convex subset of H . The weak convergence of { x n } n = 1 to x is denoted by x n x as n , while the strong convergence of { x n } n = 1 to x is written as x n x as n . For each x , y H and α ( 0 , 1 ) , we have

(9) x + y 2 x 2 + 2 y , x + y

(10) α x + ( 1 α ) y 2 = α x 2 + ( 1 α ) y 2 α ( 1 α ) x y 2 .

Definition 2.1

Let C be a nonempty closed convex subset of of a real Hilbert space H . A mapping T : C C is said to be

  1. L -Lipschitz-continuous with L > 0 if

    T x T y L x y , x , y C ;

    if L = 1 , then the operator T is called nonexpansive, and if L ( 0 , 1 ) , then T is called contraction.

  2. firmly nonexpansive if

    T x T y 2 x y 2 ( I T ) x ( I T ) y 2 , x , y C ;

  3. k -demicontractive if F ( T ) and there exists k [ 0 , 1 ) such that

    T x y 2 x y 2 + k x T x 2 , x C , y F ( T ) .

  4. monotone if

    T x T y , x y 0 x , y C .

  5. pseudomonotone if

    T x , y x 0 T y , y x 0 x , y C .

  6. α -strongly monotone if there exists a constant α > 0 such that

    T x T y , x y α x y 2 x , y C .

Lemma 2.2

Let C be a nonempty, closed, and convex set in H and x H . Then

  1. P C x x , y P C x 0 , y C .

  2. y P C x 2 x y 2 P C x x 2 , y C .

Lemma 2.3

[25] Let C be a nonempty, closed, and convex subset of H and g be a bifunction of C × C into R satisfying the following:

  1. g ( x , x ) = 0 for all x C ;

  2. g is monotone, i.e. g ( x , y ) + g ( y , x ) 0 for all x , y C ;

  3. for each x , y , z C ,

    limsup t 0 g ( t z + ( 1 t ) x , y ) g ( x , y ) ;

  4. g ( x , ) is a convex and lower semicontinuous for all x C .

If g : C × C R is a bifunction satisfying conditions ( A 1 ) ( A 4 ) and let r > 0 and x X . Then, there exists a unique z C such that

g ( z , y ) + 1 r y z , z x 0 , f o r a l l y C .

Lemma 2.4

[31] Let C be a nonempty, closed, and convex subset of H and g be a bifunction of C × C into R satisfying conditions ( A 1 ) ( A 4 ) . For r > 0 , define a mapping T r g : X C of g by

T r g x = z C : g ( z , y ) + 1 r y z , z x 0 , y C , x H .

Then, the following hold:

  1. T r g is single-valued;

  2. T r g is firmly nonexpansive;

  3. Fix ( T r g ) = E P ( g ) ;

  4. E P ( g ) is closed and convex.

Lemma 2.5

[32] Let λ ( 0 , 1 ] , T : C H be a nonexpansive mapping, and the mapping T l : C H defined by T l x T x l μ G ( T x ) x C , where G : H H is κ -Lipschitzian and η -strongly monotone. Then T λ is a contraction provided 0 < μ < 2 η / κ 2 , i.e. T l x T l y ( 1 l τ ) x y for all x , y C , where τ = 1 1 μ ( 2 η μ κ 2 ) ( 0 , 1 ] .

Lemma 2.6

[33] Let H 1 and H 2 be real Hilbert spaces and A : H 1 H 2 be a bounded linear operator with adjoint operator A . If T : H 2 H 2 is a κ -demicontractive mapping with κ < 1 and A 1 Fix ( T ) , then

( T I ) A x = 0 A ( T I ) A x = 0 x H 1 .

Lemma 2.7

[34] Let C be a nonempty, closed, and convex subset of a real Hilbert space H, and let F : C H be pseudomonotone and continuous. Then, x is a solution of V I ( C , F ) if and only if

F x , x x 0 , x C .

Lemma 2.8

[35] Assume { a n } is a sequence of nonnegative real numbers such that

a n + 1 ( 1 γ n ) a n + γ n δ n , n 0 ,

where { γ n } is a sequence in ( 0 , 1 ) and { δ n } is a sequence in R such that

  1. n = 0 γ n = ,

  2. limsup n δ n 0 .

  3. Then lim n a n = 0 .

Lemma 2.9

[36] Let { a n } be a sequence of real numbers such that there exists a subsequence { n j } of { n } such that a n j < a n j + 1 for all j N . Then, there exists a nondecreasing sequence { m k } N such that m k and the following properties are satisfied by all (sufficiently large) number k N :

a m k a m k + 1 and a k a m k + 1 .

In fact, m k = max { j k : a j < a j + 1 } .

3 Main results

We start with the following important lemma:

Lemma 3.1

[29] Let H 1 and H 2 be real Hilbert spaces and C be nonempty, closed, and convex subset of H 1 . Let A : H 1 H 2 be a bounded linear operator with its adjoint operator A , T : H 2 H 2 be a k-demicontractive mapping, and g be a bifunction of C × C into R satisfying conditions ( A 1 ) ( A 4 ) . Let { w n } be a sequence in H 1 and let

(11) u n = T σ n g ( w n + γ n A ( T I ) A w n ) , n N ,

where { γ n } is a real sequence in ( 0 , ) and { σ n } is a real sequence in [ σ , ) for some σ > 0 . Then, we have

u n p 2 w n p 2 γ n [ ( 1 k ) ( T I ) A w n 2 γ n A ( T I ) A w n 2 ]

for all n N and all p E P ( g ) such that A p Fix ( T ) .

Proof

The proof is included for completeness. Let p E P ( g ) be such that A p Fix ( T ) . Since T σ n g ( p ) = p and T σ n g is firmly nonexpansive, we have

(12) u n p 2 w n + γ n A ( T I ) A w n p 2 = w n p 2 + 2 γ n w n p , A ( T I ) A w n + γ n 2 A ( T I ) A w n 2 .

Since A is a bounded linear operator with its adjoint A and T is a k -demicontractive mapping, by Lemma 4(ii), we deduce that

(13) w n p , A ( T I ) A w n = A w n A p , ( T I ) A w n = T A w n A p , T A w n A w n ( T I ) A w n 2 = 1 2 [ T A w n A p 2 + T A w n A w n 2 A w n A p 2 ] ( T I ) A w n 2 1 2 [ A w n A p 2 + k T A w n A w n 2 + T A w n A p 2 A w n A p 2 ] ( T I ) A w n 2 = k 1 2 ( T I ) A w n 2 .

From (12) and (13), we obtain

u n p 2 w n p 2 γ n [ ( 1 k ) ( T I ) A w n 2 γ n A ( T I ) A w n 2 ] .

Now, we introduce our algorithm for approximating a solution of the split fixed-point problem for demicontractive mapping, equilibrium problem, and variational inequality problem.

Assumption 3.2

Let H 1 and H 2 be real Hilbert spaces. Suppose the following conditions hold:

  1. The set C is a nonempty, closed and convex subset of H 1 .

  2. F : H 1 H 1 is pseudomonotone, Lipschitz continuous on H 1 , and satisfy the following property:

    whenever { u n } C , u n z one has F z liminf F ( u n ) .

  3. The solution set ϒ { z V I ( C , F ) E P ( g ) : A z F ( T ) } is nonempty.

  4. T : H 2 H 2 is k -demicontractive mapping such that T I is demiclosed at zero and A : H 1 H 2 be a bounded linear operator with its adjoint operator A .

  5. We assume that f : C C is a contractive mapping with coefficient δ [ 0 , 1 ) and g : C × C R is a bifunction satisfying conditions ( A 1 ) ( A 4 ) in Lemma 2.3.

  6. ν δ < τ 1 1 ρ ( 2 η ρ κ 2 ) for ν 0 and ρ 0 , 2 η κ 2 .

Assumption 3.3

Suppose that { θ n } { l n } , { ϵ n } { α n } , and { σ n } are positive sequences satisfying the following conditions:

  1. { α n } ( 0 , 1 ) with lim n α n = 0 , n = 1 α n = .

  2. lim n lim ϵ n α n = 0 and sup n 1 ( θ n / α n ) < .

  3. 0 < liminf n l n liminf n l n < 1 and α n + l n 1 1 .

  4. { σ n } is a real sequence in [ σ , ) ( 0 , ) .

Algorithm 3.4

Step 0. Choose sequence { θ n }   { l n } , { ϵ n }   { α n } , and { σ n } such that the condition from Assumption 3.3 holds and let λ 1 > 0 , μ ( 0 , 1 ) , α 3 , and x 0 , x 1 H 1 be given arbitrarily. Set n 1 .
Step 1. Given the iterates x n 1 and x n ( n 1 ), choose θ n such that 0 θ n θ n ¯ , where
(14) θ n ¯ min n 1 n + θ 1 , ϵ n x n x n 1 if x n x n 1 , n 1 n + θ 1 otherwise .
Step 2. Set
w n = x n + θ n ( x n x n 1 ) .
Then, compute
(15) u n = T σ n g ( w n + γ n A ( T I ) A w n ) , y n = P C ( u n λ n F u n ) ,
where γ n is chosen in such a way that { γ n } is bounded and
(16) γ n ϵ , ( 1 k ) ( T I ) A w n 2 A ( T I ) A w n 2 ϵ if n Γ = { κ : ( T I ) A w k 0 } , ϵ > 0 ,
otherwise γ n = γ ( γ is any nonnegative real number).
Step 3. Construct T n = { x H u n λ n F ( y n ) y n , x y n 0 } and compute
(17) z n = P T n ( u n λ n F ( y n ) ) , x n + 1 = α n ν f ( x n ) + l n x n + [ ( 1 l n ) I α n ρ G ] z n .
Step 4. Compute
(18) λ n + 1 = min μ ( u n y n 2 + z n y n 2 ) 2 F ( u n ) F ( y n ) , z n y n , λ n , F ( u n ) F ( y n ) , z n y n > 0 λ n , otherwise .
Update n n + 1 and return to Step 1.

Remark 3.5

  1. The stepsize given by (18) are generated at each iteration by some simple computations, which allows it to be easily implemented without the prior knowledge of the Lipschitz contact of the operator F .

  2. Step 1 Algorithm 3.4 is also easily implemented since the value of x n x n 1 is a priori known before choosing α n .

  3. It is obvious that { λ n } is a monotonically decreasing sequence since F is Lipschitz-continuous mapping with constant L > 0 , and the case of F ( u n ) F ( y n ) , z n y n 0 , we have

    μ ( u n y n 2 + z n y n 2 ) 2 F ( u n ) F ( y n ) , z n y n 2 μ u n y n z n y n 2 F ( u n ) F ( y n ) z n y n μ u n y n L u n y n = μ L .

    Clearly, the sequence { λ n } has the lower bound min { μ / L , λ n } .

  4. The limit of { λ n } exists and we denote λ = lim n λ n as it is obvious that λ > 0 .

Lemma 3.6

Assume that conditions C 1 C 3 , hold. Let { u n } be generated by Algorithm 3.4 if there exists a subsequence { u n k } of { u n } such that { u n k } converges weakly to z C and lim k u n k y n k = 0 then z V I ( C , F ) .

Proof

Since y n k = P C ( u n k λ n k F u n k ) , we have

u n k λ n k F u n k y n k , x y n k 0 , x C .

Equivalently

u n k y n k , x y n k λ n k F u n k , x y n k x C .

This implies that

(19) u n k y n k λ n k , x y n k + F u n k , y n k u n k F u n k , x u n k x C .

Taking k in (19), since u n k y n k 0 and { F u n k } is bounded, we obtain

(20) liminf k F u n k , x u n k 0 .

Let us choose a sequence { ϵ k } of positive numbers decreasing and tending to 0. For each { k } , we denote by N k the smallest positive integer such that

(21) F u n j , x u n j + ϵ k 0 j N k ,

where the existence of N k follows from (20). Since { ϵ k } is decreasing, it is clear that the sequence { N k } is increasing. Furthermore, for each k , F u N k 0 and setting

v N k = F u N k F u N k 2 ,

we have F u N k , v N k = 1 for each k . Now, we can deduce from (21) that for each k

F u N k , x + ϵ k v N k u N k 0 ,

and, since F is pseudomonotone on H , then

(22) F ( x + ϵ k v N k ) , x + ϵ k v N k u N k 0 .

This implies that

F x , x u N k F x F ( x + ϵ k v N k ) , x + ϵ k v N k u N k ϵ k F x , v N k .

Next, we show that lim k ϵ k v N k = 0 . Indeed, we have u n k z C as k . Since F satisfies condition C 2 , we have

0 < F z liminf k F u n k .

Since { u N k } { u n k } and ϵ k 0 as k , we obtain

0 limsup k ϵ k v N k = limsup k ϵ k F u n k limsup k ϵ k limsup k F u n k 0 F z = 0 ,

which implies that lim k ϵ k v N k = 0 . F is uniformly continuous, the sequences { u N k } and { v N k } are bounded, and lim k v N k = 0 . Thus, we obtain that

liminf k F x , x u N k 0 .

Hence, for all x C , we have

F x , x z = lim k F x , x u N k = liminf k F x , x u N k 0 .

By Lemma 2.7, we obtain that z V I ( C , F ) .□

Lemma 3.7

Let { x n } be a sequence generated by Algorithm (3.4). Then the sequence { x n } is bounded.

Proof

Let p ϒ . Noting that y n C , we have F ( p ) , y n p 0 , n 0 . Since F is pseudo-monotone, it implies that F ( y n ) , y n p 0 , n 0 . That is

F ( y n ) , y n z n + z n p 0 n 0 .

Thus,

(23) F ( y n ) , p z n F ( y n ) , y n z n , n 0 .

By the definition of T n , we have u n λ n F ( u n ) y n , z n y n 0 . Then

(24) u n λ n F ( y n ) y n , z n y n = u n λ n F ( u n ) y n , z n y n + λ n F ( u n ) F ( y n ) , z n y n = λ n F ( u n ) F ( y n ) , z n y n .

By Lemma 2.2 ( i i ) and (23), we obtain the following:

(25) z n p 2 u n λ n F ( y n ) p 2 u n λ n F ( y n ) z n 2 = u n p 2 u n z n 2 + 2 λ n F ( y n ) , p z n u n p 2 u n z n 2 + 2 λ n F ( y n ) , y n z n = u n p 2 u n y n + y n z n 2 + 2 λ n F ( y n ) , y n z n = u n p 2 u n y n 2 y n z n 2 + 2 u n λ n F ( y n ) y n , z n y n .

It follows from (24) and (25) that

(26) z n p 2 u n p 2 u n y n 2 y n z n 2 + 2 u n λ n F ( y n ) y n , z n y n u n p 2 u n y n 2 y n z n 2 + 2 λ n F ( u n ) F ( y n ) , z n y n .

By the definition of λ n + 1 and (26), we obtain the following:

(27) z n p 2 u n p 2 u n y n 2 y n z n 2 + 2 λ n F ( u n ) F ( y n ) , z n y n = u n p 2 u n y n 2 y n z n 2 + 2 λ n λ n + 1 λ n + 1 F ( u n ) F ( y n ) , z n y n u n p 2 u n y n 2 y n z n 2 + λ n λ n + 1 μ ( u n y n 2 + z n y n 2 ) .

Consider the limit

lim n λ n μ λ n + 1 = μ , 0 < μ < 1 .

That is, there exists N 0 such that n N , 0 < λ n μ λ n + 1 < 1 . It implies that from (27), we obtain

(28) z n p u n p , n N .

Recall from Lemma 3.1, we obtain that

(29) u n p 2 w n p 2 γ n [ ( 1 k ) ( T I ) A w n 2 γ n A ( T I ) A w n 2 ] .

Hence, we obtain from (16) and (29) that

(30) u n p w n p .

Now, from Step 1 and Assumption 3.3, we have that θ n x n x n 1 ϵ n n N , which implies that

θ n α n x n x n 1 ϵ n α n 0 , as n .

Hence, there exists M 1 > 0 such that

θ n α n x n x n 1 M 1 , n N .

Thus, we obtain from Step 2 that

(31) w n p x n p + θ n x n x n 1 = x n p + α n θ n α n x n x n 1 x n p + α n M 1 , n n 0 .

By combining (28), (29), and (31), we obtain

(32) z n p w n p x n p + α n M 1 n 1 .

Thus, from α n + l n < 1 , Lemma 2.5, and (32), also, noting that for all n n 0 , α n ν δ + γ n + ( 1 l n α n τ ) 1 α n ( τ ν δ ) , we obtain

x n + 1 p l n ( x n p ) + α n ( ν f ( x n ) ρ G p ) + ( ( 1 l n ) I α n ρ G ) z n ( ( 1 l n ) I α n ρ G ) p α n ν δ x n p + α n ( ν f ρ G ) p + l n x n p + ( ( 1 l n ) I α n ρ G ) z n ( ( 1 l n ) I α n ρ G ) p α n v δ x n p + α n ( ν f ρ G ) p + l n x n p + ( 1 l n ) I α n 1 l n ρ G z n I α n 1 l n ρ G p α n ν δ x n p + α n ( ν f ρ G ) p + l n x n p + ( 1 l n ) I α n 1 l n τ z n p α n ν δ x n p + α n ( ν f ρ G ) p + l n x n p + ( 1 l n α n τ ) ( x n p + α n M 1 ) = [ α n ν δ + l n + ( 1 l n α n τ ) ] x n p + ( 1 l n α n τ ) α n M 1 + α n ( ν f ρ G ) p [ 1 α n ( τ ν δ ) ] x n p + α n ( τ ν δ ) M 1 + ( ν + ρ G ) p τ ν δ max M 1 + ( ν f ρ G ) p τ ν δ , x n p .

By induction, we conclude that

x n p max M 1 + ( ν f ρ G ) p τ ν δ , x n 0 p , n n 0 .

Therefore, we obtain that { x n } is bounded, and so, the sequences { w n } , { y n } , { z n } , { u n } , and { F ( u n ) } are bounded too.□

Theorem 3.8

Assume that Assumptions 3.2 and 3.3hold. Then, any sequence { x n } generated by Algorithm 3.4 converges strongly to an element x ϒ , where x = P ϒ f ( x ) and x Γ is a unique solution to VIP:

( ν f ρ G ) x , p x 0 p ϒ .

Proof

We first show that P ϒ ( ν f + I ρ G ) is a contractive map. Indeed, by Lemma 2.5, we have

P ϒ ( ν f + I ρ G ) x P Γ ( ν f + I ρ G ) y f ( x ) f ( y ) + ( 1 ρ G ) x ( I ρ G ) y δ x y + ( 1 τ ) ( x y ) = [ 1 ( τ δ ) ] x y x , y H 1 ,

which implies that P ϒ ( ν f + I ρ G ) is a contraction. Banach’s Contraction Mapping Principle guarantees that P ϒ ( ν f + I ρ G ) has a unique fixed-point. Say x H 1 , that is x = P ϒ ( ν f + I ρ G ) x . Thus, there exists a unique solution to x ϒ to the VIP

(33) ( ν f ρ G ) x , p x 0 p ϒ .

For simplicity, we divide the rest of the proofs into claims.

Claim 1. we show that

( 1 l n α n τ ) 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] x n p 2 x n + 1 p 2 + α n M 4 ,

for some M 4 > 0 . Indeed by (27) and Lemmas 9 and 2.5, we obtain

(34) x n + 1 p 2 = α n ν ( f ( x n ) f ( p ) ) + l n ( x n p ) + ( ( 1 l n ) I α n ρ G ) z n ( ( 1 l n ) I α n ρ G ) p + α n ( ν f ρ G ) p 2 α n ν ( f ( x n ) f ( p ) ) + l n ( x n p ) + ( ( 1 l n ) I α n ρ G ) z n ( ( 1 l n ) I α n ρ G ) p 2 + 2 α n ( ν f ρ G ) p , x n + 1 p = α n ν ( f ( x n ) f ( p ) ) + l n ( x n p ) + ( 1 l n ) I α n 1 l n ρ G z n I α n 1 l n ρ G p 2 + 2 α n ( ν f ρ G ) p , x n + 1 p α n ν δ x n p + l n x n p + ( 1 l n ) I α n 1 l n τ z n p 2 + 2 α n ( ν f ρ G ) p , x n + 1 p α n ν δ x n p 2 + l n x n p 2 + ( 1 l n α n τ ) z n p 2 + 2 α n ( ν f ρ G ) p , x n + 1 p α n ν δ x n p 2 + l n x n p 2 + α n M 2 + ( 1 l n α n τ ) u n p 2 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] ,

where sup n 1 ( 2 ( ν f ρ G ) x n + 1 p ) M 2 for some M 2 > 0 . It follows from (31) that

(35) w n p 2 ( x n p 2 + α n M 1 ) 2 = x n p 2 + α n ( 2 M 1 x n p + α n M 1 2 ) x n p 2 + α n M 3 , for some M 3 > 0 ,

where sup n 1 ( 2 M 1 x n p + α n M 1 2 ) M 3 for some M 3 > 0 . Note that α n ν δ + l n + ( 1 l n α n τ ) 1 α n ( τ ν δ ) n n 0 . Hence, from (30), (35), and the aforementioned inequality, we obtain

x n + 1 p 2 α n ν δ x n p 2 + l n x n p 2 + α n M 2 + ( 1 l n α n τ ) [ x n p 2 + α n M 3 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] ]

[ α n ν δ + l n + ( 1 l n α n τ ) ] x n p 2 + α n M 2 + α n M 3 ( 1 l n α n τ ) 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] [ 1 α n ( τ ν δ ) ] x n p 2 ( 1 l n α n τ ) 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] + α n M 2 + α n M 3 x n p 2 ( 1 l n α n τ ) 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] + α n M 2 + α n M 3 .

Hence, we obtain

( 1 l n α n τ ) 1 λ n λ n + 1 μ [ u n y n 2 + y n z n 2 ] x n p 2 x n + 1 p 2 + α n M 4 .

Claim 2.

γ n ( 1 l n α n τ ) ( 1 k ) ( T I ) A w n 2 x n p 2 x n + 1 p 2 + γ n 2 ( 1 l n α n τ ) A ( T I ) A w n 2 + α n M 4 ,

where M 4 = M 2 + M 3 . If Γ is finite, then ( T I ) A w n = 0 for all n N Γ . It follows from Lemma 2.6 that A ( T I ) A w n = 0 . Suppose Γ is infinite. It is noted that if n Γ , then lim n A ( T I ) A w n = 0 . For n Γ , again from (28), (29), (34), and condition of γ n , we have

(36) x n + 1 p 2 α n ν δ x n p 2 + l n x n p 2 + α n M 2 + ( 1 l n α n τ ) [ w n p 2 γ n ( ( 1 k ) ( T I ) A w n 2 γ n A ( T I ) A w n 2 ) ] [ 1 α n ( τ ν δ ) ] x n p 2 + α n M 2 + α n M 3 γ n ( 1 l n α n τ ) [ ( 1 k ) ( T I ) A w n 2 γ n A ( T I ) A w n 2 ] x n p 2 + ( 1 l n α n τ ) [ γ n ( 1 k ) ( T I ) A w n 2 + γ n 2 A ( T I ) A w n 2 ] + α n M 4 x n p 2 γ n ϵ ( 1 l n α n τ ) A ( T I ) A w n 2 + α n M 4 .

Hence, we obtain

(37) γ n ϵ ( 1 l n α n τ ) A ( T I ) A w n 2 x n p 2 x n + 1 p 2 + α n M 4 .

We can also obtain from (36) that

(38) γ n ( 1 l n α n τ ) ( 1 k ) ( T I ) A w n 2 x n p 2 x n + 1 p 2 + γ n 2 ( 1 l n α n τ ) A ( T I ) A w n 2 + α n M 4 .

Claim 3. We show that

( 1 l n α n τ ) u n w n 2 x n p 2 x n + 1 p 2 + 2 γ n ( 1 l n α n τ ) u n p A ( T I ) A w n + α n M 4 .

Using the fact that T σ n g is firmly nonexpansive, we have

u n p 2 = T σ n g ( w n + γ n A ( T I ) A w n ) p 2 = T σ n g ( w n + γ n A ( T I ) A w n ) T σ n g ( p ) 2 u n p , w n + γ n A ( T I ) A w n p = 1 2 u n p 2 + 1 2 w n + γ n A ( T I ) A w n p 2 1 2 u n p w n γ n A ( T I ) A w n + p 2 = 1 2 u n p 2 + 1 2 w n + γ n A ( T I ) A w n p 2 1 2 u n w n γ n A ( T I ) A w n 2

= 1 2 u n p 2 + 1 2 w n p 2 + 1 2 γ n 2 A ( T I ) A w n 2 + w n p , γ n A ( T I ) A w n 1 2 u n w n 2 1 2 γ n 2 A ( T I ) A w n 2 + u n w n , γ n A ( T I ) A w n = 1 2 u n p 2 + 1 2 w n p 2 1 2 u n w n 2 + u n p , γ n A ( T I ) A w n .

This implies that

(39) u n p 2 w n p 2 u n w n 2 + 2 u n p , γ n A ( T I ) A w n w n p 2 u n w n 2 + 2 γ n u n p A ( T I ) A w n .

By combining (34) and (39), we obtain

(40) x n + 1 p 2 α n ν δ x n p 2 + l n x n p 2 + ( 1 l n α n τ ) z n p 2 + 2 α n ( ν f ρ G ) p , x n + 1 p α n ν δ x n p 2 + l n x n p 2 + ( 1 l n α n τ ) [ w n p 2 w n u n 2 + 2 γ n u n p A ( T I ) A w n ] + α n M 2 [ 1 α n ( τ ν δ ) ] x n p 2 ( 1 l n α n τ ) u n w n 2 + 2 γ n ( 1 l n α n τ ) u n p A ( T I ) A w n + α n M 4 x n p 2 ( 1 l n α n τ ) u n w n 2 + 2 γ n ( 1 l n α n τ ) u n p A ( T I ) A w n + α n M 4 .

Hence, we obtain

( 1 l n α n τ ) u n w n 2 x n p 2 x n + 1 p 2 + 2 γ n ( 1 l n α n τ ) u n p A ( T I ) A w n + α n M 4 .

Claim 4. We show that

x n + 1 p 2 [ 1 α n ( τ ν δ ) ] x n p 2 + α n ( τ ν δ ) 2 ( ν f ρ G ) p , x n + 1 p τ ν δ + θ n α n . x n x n 1 M τ ν δ .

For some M > 0 , indeed, we have

(41) w n p 2 = x n + θ n ( x n x n 1 ) p 2 = x n p 2 + 2 θ n x n p , x n x n 1 + θ n 2 x n x n 1 2 x n p 2 + 2 θ n x n p x n x n 1 + θ n 2 x n x n 1 2 = x n p 2 + θ n x n x n 1 [ 2 x n p + θ n x n x n 1 ] x n p 2 + θ n x n x n 1 M ,

where sup n 1 [ 2 x n p + θ n x n x n 1 ] M . Note that α n ν δ + l n + ( 1 l n α n τ ) ( 1 α n ( τ ν δ ) ) for all n n 0 . Thus, combining (34) and (41), we have that for all n n 0 ,

(42) x n + 1 p 2 α n ν δ x n p 2 + l n x n p 2 + ( 1 l n α n τ ) [ x n p 2 + θ n x n x n 1 M ] + 2 α n ( ν f ρ G ) p , x n + 1 p = [ α n ν δ + l n + ( 1 l n α n τ ) ] x n p 2 + ( 1 l n α n τ ) θ n x n x n 1 M + 2 α n ( ν f ρ G ) p , x n + 1 p [ 1 α n ( τ ν δ ) ] x n p 2 + α n ( τ ν δ ) 2 ( ν f ρ G ) p , x n + 1 p τ ν δ + θ n α n . x n x n 1 M τ ν δ .

Claim 5. The sequence { x n p 2 } converges to zero by considering two possible cases in the sequence { x n p 2 } :

Case 1: There exists N N such that x n + 1 p 2 x n p 2 for all n N . This implies that lim n x n p 2 exists. It follows from Claim 1 that

(43) lim n u n y n = 0 and lim n z n y n = 0 .

In addition, by (37), we have that

(44) lim n A ( T I ) A w n 2 = 0 .

We also obtain from (38) and (44) that

(45) lim n ( T I ) A w n = 0 .

We obtain from Claim 2. that

(46) lim n u n w n = 0 ,

(47) w n x n = θ n x n x n 1 = θ n α n x n x n 1 α n 0 as n .

On the other hand, we obtain that

u n z n u n y n + y n z n 0 , as n

and

z n x n z n u n + u n x n 0 as n .

Thus,

(48) x n + 1 x n ( 1 l n ) z n x n + α n ν f ( x n ) ρ G z n 0 as n .

Since the sequence { x n } is bounded, then there exists a subsequence { x n k } such that x n k z , from (43), there exists a subsequence { u n k } of { u n } such that u n k z , by Lemma 3.6 implies that z V I ( C , F ) . By (47), we have w n k z H 1 . Since A is a bounded linear operator, we have A w n k A z . By (45) and demiclosedness of T I , we obtain A z F ( T ) . From u n = T σ n g ( I + λ n A ( T I A ) w n , we have

g ( u n , y ) + 1 σ n y u n , u n w n λ n A ( T I ) A w n 0 y C ,

which implies,

g ( u n , y ) + 1 σ n y u n , u n w n 1 σ n y u n , λ n A ( T I ) A w n 0 y C .

By condition ( A 2 ) , we have

1 σ n k y u n k , u n k w n k 1 σ n k y u n k , λ n k A ( T I ) A w n k g ( y , u n k ) y C .

This, together with (45) and (46), gives

(49) g ( y , z ) 0 , y C .

For any t ( 0 , 1 ] and y C , let y t = t y + ( 1 t ) z . We see that y t C . From (49), we have g ( y t , z ) 0 . By conditions ( A 1 ) and ( A 4 ) , we have

0 = g ( y t , y t ) t g ( y t , y ) + ( 1 t ) g ( y t , z ) t g ( y t , y ) ,

hence g ( y t , y ) 0 . By condition ( A 3 ) , we obtain that g ( z , y ) 0 , that is z E P ( g ) . Hence, z ϒ .

On the other hand, from the boundedness of { x n } , it follows that { x n k } { x n } such that

(50) limsup n ( ν f ρ G ) z , x n z = lim k ( ν f ρ G ) z , x n k z .

Utilizing the reflexivity of H 1 and the boundedness of { x n } , one may suppose that x n k x . Therefore, one obtains from (50) that

(51) limsup n ( ν f ρ G ) z , x n z = ( ν f ρ G ) z , x z 0 ,

which together with x n x n + 1 0 implies that

limsup n ( ν f ρ G ) z , x n + 1 z = limsup n [ ( ν f ρ G ) z , x n + 1 x n + ( ν f ρ G ) z , x n z ] = ( ν f ρ G ) z , x z 0 .

Observe that α n ( τ ν δ ) [ 0 , 1 ] , n = 1 α n ( τ ν δ ) = , and

(52) limsup n 2 τ ν δ ( ν f ρ G ) z , x n + 1 z + θ n α n . 2 M τ ν δ x n x n 1 0 .

Consequently, by Lemma 2.8 and from Claim 4, we obtain x n z 0 , as n .

Case 2. There exists a subsequence { x n j p 2 } of { x n p 2 } such that x n j p 2 < x n j + 1 p 2 for all j N . In this case, it follows from Lemma 2.9 that there exists a non-decreasing sequence { m k } of N such that lim k m k = , and the following inequality holds for all k N .

(53) x m k p 2 x m k + 1 p 2 and x k p 2 x m k p 2 .

According to Claim 1, we have

( 1 l m k α m k τ ) 1 λ m k λ m k + 1 μ [ u m k y m k 2 + y m k z m k 2 ] x m k p 2 x m k + 1 p 2 + α m k M 4 α m k M 4 .

Therefore,

lim k u m k y m k = 0 and lim k y m k z m k = 0 .

According to Claim 2, we have

γ m k ( 1 l m k α m k τ ) ( 1 k ) ( T I ) A w m k 2 x m k p 2 x m k + 1 p 2 + γ m k 2 ( 1 l m k α m k τ ) ( 1 k ) A ( T I ) A w m k 2 + α m k M 4 α m k M 4 0 as k .

Hence, we obtain

lim k ( T I ) A w m k = 0 .

Using the same argument as in the proof of Case 1, we obtain

lim k x m k + 1 x m k = 0 , lim k u m k z m k = 0

and

limsup k f ( p ) p , x m k + 1 p 0 .

According to Claim 5, we have

(54) x m k + 1 p 2 ( 1 ( τ ν δ ) α m k ) x m k p 2 + ( τ ν δ ) α m k × 2 τ ν δ ( ν f ρ G ) p , x m k + 1 p + M τ ν ρ θ m k α m k x m k x m k 1 .

From (53) and (54), we obtain

x m k + 1 p 2 ( 1 ( τ ν δ ) α m k ) x m k + 1 p 2 + ( τ ν δ ) α m k × 2 τ ν δ ( ν f ρ G ) p , x m k + 1 p + M τ ν ρ θ m k α m k x m k x m k 1 .

Thus,

x m k + 1 p 2 ( ν f ρ G ) p , x m k + 1 p + M τ ν ρ θ m k α m k x m k x m k 1 .

Therefore,

(55) limsup k x m k + 1 p 0 .

Combining (53) and (55), we obtain limsup k x k p 0 , that is x k p .□

Applying Algorithm 3.4 with u n = w n , we obtain the following corollary.

Assumption 3.9

Let H 1 and H 2 be real Hilbert spaces. Suppose the following conditions hold:

  1. The set C is a nonempty, closed and convex subset of H 1 .

  2. F : H 1 H 1 is pseudomonotone, uniformly continuous on C , and satisfy the following property:

    whenever { u n } C , u n z one has F z liminf F ( u n ) .

  3. The solution set V I ( C , F ) is nonempty.

  4. We assume that f : C C is a contractive mapping with coefficient δ [ 0 , 1 ) .

  5. ν δ < τ 1 1 ρ ( 2 η ρ κ 2 ) for ν 0 and ρ 0 , 2 η κ 2 .

Assumption 3.10

Suppose that { θ n } { l n } , { ϵ n } , and { α n } are positive sequences satisfying the following conditions:

  1. { α n } ( 0 , 1 ) with lim n α n = 0 , n = 1 α n = .

  2. lim n ϵ n α n = 0 and sup n 1 ( θ n / α n ) < .

  3. 0 < liminf n l n limsup n l n < 1 and α n + l n 1 1 .

Algorithm 3.11

Step 0. Choose sequence { θ n }   { l n } , { ϵ n } , and { α n } such that the condition from Assumption 3.10 holds and let λ 1 > 0 , μ ( 0 , 1 ) , α 3 , and x 0 , x 1 H 1 be given arbitrarily. Set n 1 .
(56) θ n ¯ min n 1 n + θ 1 , ϵ n x n x n 1 if x n x n 1 , n 1 n + θ 1 otherwise .
Step 2. Set
u n = x n + θ n ( x n x n 1 )
Then, compute
(57) y n = P C ( u n λ n F u n ) .
Step 3. Construct T n = { x H u n λ n F ( y n ) y n , x y n 0 } and compute
(58) z n = P T n ( u n λ n F ( y n ) ) , x n + 1 = α n ν f ( x n ) + l n x n + [ ( 1 l n ) I α n ρ G ] z n .
Step 4. Compute
(59) λ n + 1 = min μ ( u n y n 2 + z n y n 2 ) 2 F ( u n ) F ( y n ) , z n y n , λ n , F ( u n ) F ( y n ) , z n y n > 0 λ n , otherwise .
Update n n + 1 and return to Step 1.

Corollary 3.12

Assume that Assumptions 3.9 and 3.10 hold. Then, any sequence { x n } generated by Algorithm 3.11 converges strongly to an element x V I ( C , F ) , where x = P V I ( C , F ) f ( x ) and x V I ( C , F ) is a unique solution to VIP:

( ν f ρ G ) x , p x 0 p ϒ .

4 Numerical example

In this section, we present a numerical example to illustrate the behaviour of the sequence generated by Algorithm 3.4.

Example 4.1

Let H 1 = H 2 = 2 ( R ) be the linear space whose elements are all 2-summable sequences { x i } i = 1 of scalars in R , that is

2 ( R ) x = ( x 1 , x 2 , , x 3 , ) , x i R and i = 1 x i 2 < ,

with the inner product , : 2 × 2 R be defined by x , y i = 1 x i y i and the norm : 2 R be defined by x 2 i = 1 x i 2 1 2 , where x = { x i } i = 1 and y = { y i } i = 1 . Let C { x 2 ( R ) : x 2 1 } , F : C 2 be defined by

F ( x ) = x 2 + 1 x 2 + 1 x x C ,

and g : C × C R be defined by g ( x , y ) = x 2 + x y 2 y 2 for all x , y C . It is easy to verify that F is pseudomonotone, not Lipschitz continuous and conditions (A1)–(A4) are satisfied. More so,

T σ n g z = z 3 σ n + 1 , z C .

Now let T : 2 ( R ) 2 ( R ) be defined by

(60) T ( x 1 , x 2 , x 3 , ) = ( α + 1 ) 2 ( x 1 , x 2 , x 3 , ) ,

where α > 1 . We then show that T is k -demicontractive mapping but not quasi-nonexpansive.

It is obvious that F ( T ) = { 0 } . Choose any x 2 ( R ) , then we have

(61) T x 0 2 2 = ( α + 1 ) 2 x 0 2 2 = α + 1 2 2 x 0 2 = x 0 2 2 + α 2 + 2 α 3 4 x 0 2 2 .

Thus, T is not quasi-nonexpansive since α 2 + 2 α 3 4 > 0 for any α > 1 . Moreover,

x T x 2 2 = x + α + 1 2 x 2 2 = α + 3 2 2 x 0 2 2 .

This implies that

(62) x 0 2 2 = 4 ( α + 3 ) 2 x T x 2 2 .

It follows from (61) and (62) that

T x 0 2 2 = x 0 2 2 + α 2 + 2 α 3 4 × 4 ( α + 3 ) 2 x T x 2 2 = x 0 2 2 + α 1 α + 3 x T x 2 2 .

Therefore, T is k -demicontractive mapping with k = α 1 α + 3 ( 0 , 1 ) for any α > 1 .

In this example, we show the numerical behaviour of the sequences generated by Algorithm 3.4 using different starting points. We choose A x = x , f ( x ) = x 2 , l = 0.05 , μ = 2 , λ = 0.34 , α n = 1 n + 1 , θ n = 1 ( n + 1 ) 2 , σ n = n n + 1 , and λ n = 2 n 5 n + 4 , and compare the performance of Algorithm 3.4 with the non-inertial type, i.e., by taking θ n = 0 . The stopping criterion used was x n + 1 x n < 1 0 5 , and we plot the graphs of Error ( x n + 1 x n ) against number of iteration in each case as follows:

  1. x 0 = ( 1 , 1 , 1 , 0 , ) and x 1 = ( 1 , 2 , 3 , 4 , ) ,

  2. x 0 = ( 2 , 2 , 2 , 0 , ) and x 1 = ( 5 , 5 , 5 , 0 , ) ,

  3. x 0 = ( 3 , 3 , 3 , 0 , ) and x 1 = ( 1 , 0 , 1 , 0 , ) ,

  4. x 0 = ( 2 , 0 , 2 , 0 , ) and x 1 = ( 4 , 4 , 4 , 0 , ) .

The computation results are shown in Table 1 and Figure 1.

Table 1

Computation result for Example 4.1

Algorithm 3.4 Non-inertial alg.
Case I No of Iter. 15 90
CPU time (s) 0.0038 0.0630
Case II No of Iter. 18 121
CPU time (s) 0.0047 0.0998
Case III No of Iter. 10 49
CPU time (s) 0.0036 0.0325
Case IV No of Iter. 9 37
CPU time (s) 0.0030 0.0111
Figure 1 
               Example 4.1, Top Left: Case I; Top Right: Case II; Bottom Left: Case III; Bottom Right: Case IV.
Figure 1

Example 4.1, Top Left: Case I; Top Right: Case II; Bottom Left: Case III; Bottom Right: Case IV.

5 Conclusion

In this article, we proposed an inertia extragradient-type algorithm for solving split fixed-point problems of demicontractive mapping, variational inequalities problems, and equilibrium problems in a real Hilbert spaces. The proposed algorithm shows the strong convergence property under pseudomonotonicity and non-Lipschitz continuity of the mapping F . The algorithms require the calculation of only two projections onto the feasibility set C per iteration. These two properties of pseudo-monotonicity and non-Lipschitz continuity of the mapping F emphasize the applicability and advantages of the mapping over several existing results in the literature.



Acknowledgment

The research was completed when L.O. Jolaoso was visiting the Federal University of Agriculture Abeokuta Nigeria. The authors acknowledge their respective institutions for making their facilities available for the research.

  1. Funding information: L.O. Jolaoso is supported by the Postdoctoral research grant from the Sefako Makgatho Health Sciences University, South Africa.

  2. Author contributions: All authors contribute equally to prepare the manuscript.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

References

[1] J. P. Aubin and I. Ekeland, Applied Nonlinear Analysis, Willey, New York, 1984. Search in Google Scholar

[2] C. Baiocchi and A. Capelo, Variational and Quasivariational Inequalities; Applications to Free Boundary Problems, Willey, New York, 1984. Search in Google Scholar

[3] R. Glowinski, J.-L. Lions, and R. Trémolières, Numerical Analysis of Variational Inequalities, Elsevier, North-Holland, Amsterdam, 1981. Search in Google Scholar

[4] E. N. Khobotov, Modification of the extra gradient method for solving variational inequalities and certain optimization problems, USSR Comput. Math. Phys. 27 (1989), 120–127. 10.1016/0041-5553(87)90058-9Search in Google Scholar

[5] D. Kinderlehrer and G. Stampacchia, An Introduction to Variational Inequalities and their Applications, Academic Press, New York, 1980. Search in Google Scholar

[6] I. V. Konnov, Combined Relaxation Methods for Variational Inequalities, Springer-Verlag, Berlin, 2001. 10.1007/978-3-642-56886-2Search in Google Scholar

[7] P. Marcotte, Applications of Khobotov’s algorithm to variational and network equilibrium problems, Inf. Syst. Oper. Res. 29 (1991), 258–270. 10.1080/03155986.1991.11732174Search in Google Scholar

[8] F. Facchinei and J.-S. Pang, Finite-Dimensional Variational Inequalities and Complementary Problems, Vol II, Springer Series in Operations Research, Springer, New York, 2003. 10.1007/b97544Search in Google Scholar

[9] A. Nagurney, Network Economics: A variational Inequalities Approach, Kluwer Academic Publishers, Dordrecht, 1999. 10.1007/978-1-4757-3005-0Search in Google Scholar

[10] G. M. Korpelevich, The extragradient method for finding saddle points and other problems, Ekonomika Mat. Metody 12 (1976), 747–756. Search in Google Scholar

[11] A. Gibali, Two simple relaxed perturbed extragradient methods for solving variational inequalities in Euclidean spaces, J. Nonlinear Var. Anal. 2 (2018), 49–61. 10.23952/jnva.2.2018.1.05Search in Google Scholar

[12] B. S. He, Z.-H. Yang, and X.-M. Yuan, An approximate proximal-extragradient-type method for monotone variational inequalities, J. Math. Anal. Appl. 300 (2004), 362–374. 10.1016/j.jmaa.2004.04.068Search in Google Scholar

[13] A. N. Iusem and B. F. Savaiter, A variant of Korpelevich’s method for variational inequalities with a new search strategy, Optimization 42 (1997), 309–321. 10.1080/02331939708844365Search in Google Scholar

[14] N. Nadezhkina and W. Takahashi, Weak convergence theorem by an extragradient method for nonexpansive mappings and monotone mappings, J. Optim. Theory Appl. 128 (2006), 191–201. 10.1007/s10957-005-7564-zSearch in Google Scholar

[15] M. V. Solodov and B. F. Svaiter, A new projection method for variational inequality problems, SIAM J. Control Optim. 37 (1999), 765–776. 10.1137/S0363012997317475Search in Google Scholar

[16] P. Tseng, A modified forward-backward splitting method for maximal monotone mappings, SIAM J. Control Optim. 38 (2000), 431–446. 10.1137/S0363012998338806Search in Google Scholar

[17] Y. Yao, G. Marino, and L. Muglia, A modified Korpelevich’s method convergent to minimum-norm solution of a variational inequality, Optimization 63 (2014), 559–569. 10.1080/02331934.2012.674947Search in Google Scholar

[18] Y. Censor, A. Gibali, and S. Reich, The subgradient extragradient method for solving variational inequalities in Hilbert space, J. Optim. Theory Appl. 148 (2011), 318–335. 10.1007/s10957-010-9757-3Search in Google Scholar

[19] Y. Censor, A. Gibali, and S. Reich, Strong convergence of subgradient extragradient methods for the variational inequality problem in Hilbert space, Optim. Meth. Softw. 26 (2011), 827–845. 10.1080/10556788.2010.551536Search in Google Scholar

[20] Y. Censor, A. Gibali, and S. Reich, Extensions of Korpelevichas extragradient method for the variational inequality problem in Euclidean space, Optimization 61 (2012), 1119–1132. 10.1080/02331934.2010.539689Search in Google Scholar

[21] A. Gibali, A new non-Lipschitzian projection method for solving variational inequalities in Euclidean spaces, J. Nonlinear Anal. Optim. 6 (2015), 41–51. Search in Google Scholar

[22] Y. Shehu and O. S. Iyiola, Strong convergence result for monotone variational inequalities, Numer. Algor. 76 (2017), 259–282, https://doi.org/10.1007/s11075-016-0253-1. Search in Google Scholar

[23] Y. Censor and Y. Segal, The split common fixed-point for directed operators, J. Convex Anal. 16 (2009), 587–600. 10.1111/j.1475-3995.2008.00684.xSearch in Google Scholar

[24] F. Alvarez and H. Attouch, An inertial proximal method for monotone operators via discretization of nonlinear oscillator with damping, Set-Valued Anal. 9 (2001), 3–11. 10.1023/A:1011253113155Search in Google Scholar

[25] E. Blum and W. Oettli, From optimization and variational inequalities to equilibrium problems, Math. Stud. 63 (1993), 1–23. Search in Google Scholar

[26] C. C. Okeke and O. T. Mewomo, On split equilibrium problem, variational inequality problem and fixed-point problem for multi-valued mappings, Ann. Acad. Rom. Sci. Ser. Math. Appl. 9 (2017), no. 2, 223–248. Search in Google Scholar

[27] C. C. Okeke and O. T. Mewomo, Regularized gradient-projection algorithm for solving one-parameter nonexpansive semigroup, constrained convex minimization and generalized equilibrium problems, Buletinul Academiei de Ştiinţe a Republicii Moldova Matematica 3 (2018), no. 88, 32–56. Search in Google Scholar

[28] Y. Shehu and F. U. Ogbuisi, An iterative algorithm for approximating a solution of split common fixed-point problem for demicontractive maps, Dynam. Cont. Dis. Ser. B 23 (2012), 225–240. Search in Google Scholar

[29] A. Hanjing and S. Suantai, Hybrid inertial accelerated algorithms for split fixed-point problems of demicontractive mappings and equilibrium problems, Numer. Algor. 85 (2020), 1051–1073, https://doi.org/10.1007/s11075-019-00855-y. Search in Google Scholar

[30] J. Yang, H. Liu, and Z. Liu, Modified subgradient extragradient algorithms for solving monotone variational inequalities, Optimization 67 (2018), no. 12, 2247–2258. 10.1080/02331934.2018.1523404Search in Google Scholar

[31] P. L. Combettes and A. Hirostoaga, Equilibrium programming in Hilbert spaces, J. Nonlinear Convex Anal. 6 (2005), 117–136. Search in Google Scholar

[32] I. Yamada, The hybrid steepest-descent method for variational inequalities problems over the intersection of the fixed-point sets of nonexpansive mappings, in: D. Butnariu, Y. Censor, S. Reich (eds), Inherently Parallel Algorithms in Feasibility and Optimization and Their Applications, The Netherlands, North-Holland: Amsterdam, 2001, pp. 473–504. 10.1016/S1570-579X(01)80028-8Search in Google Scholar

[33] H. Cui and F. Wang, Iterative methods for the split common fixed-point problem in a Hilbert spaces, Fixed Point Theory Appl. 2014 (2014), 78. 10.1186/1687-1812-2014-78Search in Google Scholar

[34] R. W. Cottle and J. C. Yao, Pseudomonotone of complementarity problems in Hilbert space, J. Optim Theory Appl. 75 (1992), 281–295. 10.1007/BF00941468Search in Google Scholar

[35] H. K. Xu, Iterative algorithms for nonlinear operators, J. Lond. Math. Soc. 66 (2002), 1–17. 10.1112/S0024610702003332Search in Google Scholar

[36] P. E. Maingé, Strong convergence of projected subgradient methods for nonsmooth and nonstrictly convex minimization, Set-Valued Anal. 16 (2008), 899–912. 10.1007/s11228-008-0102-zSearch in Google Scholar

[37] A. Moudafi, The split common fixed-point problem for demicontractive mappings, Inverse Probl. 26 (2010), 587–600. 10.1088/0266-5611/26/5/055007Search in Google Scholar PubMed PubMed Central

Received: 2020-12-28
Revised: 2021-08-09
Accepted: 2021-12-28
Published Online: 2022-09-22

© 2022 Chibueze C. Okeke et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. On some summation formulas
  3. A study of a meromorphic perturbation of the sine family
  4. Asymptotic behavior of even-order noncanonical neutral differential equations
  5. Unconditionally positive NSFD and classical finite difference schemes for biofilm formation on medical implant using Allen-Cahn equation
  6. Starlike and convexity properties of q-Bessel-Struve functions
  7. Mathematical modeling and optimal control of the impact of rumors on the banking crisis
  8. On linear chaos in function spaces
  9. Convergence of generalized sampling series in weighted spaces
  10. Persistence landscapes of affine fractals
  11. Inertial iterative method with self-adaptive step size for finite family of split monotone variational inclusion and fixed point problems in Banach spaces
  12. Various notions of module amenability on weighted semigroup algebras
  13. Regularity and normality in hereditary bi m-spaces
  14. On a first-order differential system with initial and nonlocal boundary conditions
  15. On solving pseudomonotone equilibrium problems via two new extragradient-type methods under convex constraints
  16. Local linear approach: Conditional density estimate for functional and censored data
  17. Some properties of graded generalized 2-absorbing submodules
  18. Eigenvalue inclusion sets for linear response eigenvalue problems
  19. Some integral inequalities for generalized left and right log convex interval-valued functions based upon the pseudo-order relation
  20. More properties of generalized open sets in generalized topological spaces
  21. An extragradient inertial algorithm for solving split fixed-point problems of demicontractive mappings, with equilibrium and variational inequality problems
  22. An accurate and efficient local one-dimensional method for the 3D acoustic wave equation
  23. On a weighted elliptic equation of N-Kirchhoff type with double exponential growth
  24. On split feasibility problem for finite families of equilibrium and fixed point problems in Banach spaces
  25. Entire and meromorphic solutions for systems of the differential difference equations
  26. Multiplication operators on the Banach algebra of bounded Φ-variation functions on compact subsets of ℂ
  27. Mannheim curves and their partner curves in Minkowski 3-space E13
  28. Characterizations of the group invertibility of a matrix revisited
  29. Iterates of q-Bernstein operators on triangular domain with all curved sides
  30. Data analysis-based time series forecast for managing household electricity consumption
  31. A robust study of the transmission dynamics of zoonotic infection through non-integer derivative
  32. A Dai-Liao-type projection method for monotone nonlinear equations and signal processing
  33. Review Article
  34. Remarks on some variants of minimal point theorem and Ekeland variational principle with applications
  35. Special Issue on Recent Methods in Approximation Theory - Part I
  36. Coupled fixed point theorems under new coupled implicit relation in Hilbert spaces
  37. Approximation of integrable functions by general linear matrix operators of their Fourier series
  38. Sharp sufficient condition for the convergence of greedy expansions with errors in coefficient computation
  39. Approximation of conic sections by weighted Lupaş post-quantum Bézier curves
  40. On the generalized growth and approximation of entire solutions of certain elliptic partial differential equation
  41. Existence results for ABC-fractional BVP via new fixed point results of F-Lipschitzian mappings
  42. Linear barycentric rational collocation method for solving biharmonic equation
  43. A note on the convergence of Phillips operators by the sequence of functions via q-calculus
  44. Taylor’s series expansions for real powers of two functions containing squares of inverse cosine function, closed-form formula for specific partial Bell polynomials, and series representations for real powers of Pi
  45. Special Issue on Recent Advances in Fractional Calculus and Nonlinear Fractional Evaluation Equations - Part I
  46. Positive solutions for fractional differential equation at resonance under integral boundary conditions
  47. Source term model for elasticity system with nonlinear dissipative term in a thin domain
  48. A numerical study of anomalous electro-diffusion cells in cable sense with a non-singular kernel
  49. On Opial-type inequality for a generalized fractional integral operator
  50. Special Issue on Advances in Integral Transforms and Analysis of Differential Equations with Applications
  51. Mathematical analysis of a MERS-Cov coronavirus model
  52. Rapid exponential stabilization of nonlinear continuous systems via event-triggered impulsive control
  53. Novel soliton solutions for the fractional three-wave resonant interaction equations
  54. The multistep Laplace optimized decomposition method for solving fractional-order coronavirus disease model (COVID-19) via the Caputo fractional approach
  55. Special Issue on Problems, Methods and Applications of Nonlinear Analysis
  56. Some recent results on singular p-Laplacian equations
  57. Infinitely many solutions for quasilinear Schrödinger equations with sign-changing nonlinearity without the aid of 4-superlinear at infinity
  58. Special Issue on Recent Advances for Computational and Mathematical Methods in Scientific Problems
  59. Existence of solutions for a nonlinear problem at resonance
  60. Asymptotic stability of solutions for a diffusive epidemic model
  61. Special Issue on Computational and Numerical Methods for Special Functions - Part I
  62. Fully degenerate Bernoulli numbers and polynomials
  63. Wigner-Ville distribution and ambiguity function associated with the quaternion offset linear canonical transform
  64. Some identities related to degenerate Stirling numbers of the second kind
  65. Two identities and closed-form formulas for the Bernoulli numbers in terms of central factorial numbers of the second kind
  66. λ-q-Sheffer sequence and its applications
  67. Special Issue on Fixed Point Theory and Applications to Various Differential/Integral Equations - Part I
  68. General decay for a nonlinear pseudo-parabolic equation with viscoelastic term
  69. Generalized common fixed point theorem for generalized hybrid mappings in Hilbert spaces
  70. Computation of solution of integral equations via fixed point results
  71. Characterizations of quasi-metric and G-metric completeness involving w-distances and fixed points
  72. Notes on continuity result for conformable diffusion equation on the sphere: The linear case
Downloaded on 13.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/dema-2020-0120/html
Scroll to top button