Startseite Inertial shrinking projection algorithm with self-adaptive step size for split generalized equilibrium and fixed point problems for a countable family of nonexpansive multivalued mappings
Artikel Open Access

Inertial shrinking projection algorithm with self-adaptive step size for split generalized equilibrium and fixed point problems for a countable family of nonexpansive multivalued mappings

  • Musa A. Olona , Timilehin O. Alakoya , Abd-semii O.-E. Owolabi und Oluwatosin T. Mewomo EMAIL logo
Veröffentlicht/Copyright: 16. April 2021
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

In this paper, we introduce a shrinking projection method of an inertial type with self-adaptive step size for finding a common element of the set of solutions of a split generalized equilibrium problem and the set of common fixed points of a countable family of nonexpansive multivalued mappings in real Hilbert spaces. The self-adaptive step size incorporated helps to overcome the difficulty of having to compute the operator norm, while the inertial term accelerates the rate of convergence of the proposed algorithm. Under standard and mild conditions, we prove a strong convergence theorem for the problems under consideration and obtain some consequent results. Finally, we apply our result to solve split mixed variational inequality and split minimization problems, and we present numerical examples to illustrate the efficiency of our algorithm in comparison with other existing algorithms. Our results complement and generalize several other results in this direction in the current literature.

MSC 2010: 65K15; 47J25; 65J15; 90C33

1 Introduction

Let H be a real Hilbert space with inner product , and induced norm . Let C be a nonempty closed convex subset of H and ϕ : C × C R , F : C × C R be two bifunctions. The generalized equilibrium problem (GEP) is to find a point x C such that

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

The solution set of the GEP is denoted by GEP ( F , ϕ ) . In particular, if we set ϕ = 0 in (1), then the GEP reduces to the classical equilibrium problem (EP), which is to find a point x C such that F ( x , y ) 0 , y C . The solution set of EP is denoted by EP ( F ) .

The EP is a generalization of many mathematical models such as variational inequality problems (VIPs), fixed point problems (FPPs), certain optimization problems (OPs), Nash EPs, minimization problems (MPs), and others, see [1,2]. Many authors have studied and proposed several iterative algorithms for solving EPs and related OPs, see [318].

In 2013, Kazmi and Rizvi [19] introduced and studied the following split generalized equilibrium problem (SGEP): let C H 1 and Q H 2 , where H 1 and H 2 are real Hilbert spaces. Let F 1 , ϕ 1 : C × C R and F 2 , ϕ 2 : Q × Q R be nonlinear bifunctions, and A : H 1 H 2 be a bounded linear operator. The SGEP is defined as follows: find x C such that

(2) F 1 ( x , x ) + ϕ 1 ( x , x ) 0 , x C ,

and such that

(3) y = A x Q solves F 2 ( y , y ) + ϕ 2 ( y , y ) 0 , y Q .

We denote the solution set of SGEP (2)–(3) by

SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) { x C : x GEP ( F 1 , ϕ 1 ) and A x GEP ( F 2 , ϕ 2 ) } .

Furthermore, an iterative algorithm was also presented by the authors for approximating the solution of SGEP in a real Hilbert space. If ϕ 1 = 0 and ϕ 2 = 0 , then the SGEP reduces to split equilibrium problem (SEP), which is to find x C such that

(4) F 1 ( x , x ) 0 , x C ,

and such that

(5) y = A x Q solves F 2 ( y , y ) 0 , y Q .

Observe that (4) is the classical EP. Therefore, the inequalities (4) and (5) comprise a pair of EPs, which involves finding the image y = A x under a given bounded linear operator A , of the solution x of (4) in H 1 , which is the solution of (5) in H 2 . The solution set of SEP (4)–(5) is denoted by SEP ( F 1 , F 2 ) { z EP ( F 1 ) : A z EP ( F 2 ) } .

Another important problem in fixed point theory is the fixed point problem (FPP), which is defined as follows:

(6) Find a point x C such that S x = x ,

where S : C C is a nonlinear operator. If S is a multivalued mapping, i.e., S : C 2 C , then x C is called a fixed point of S if

(7) x S x .

We denote the set of fixed points of S by F ( S ) . The fixed point theory for multivalued mappings can be utilized in various areas such as game theory, control theory, and mathematical economics.

In this article, we are interested in studying the problem of finding a common solution for both the SGEP (2)–(3) and the common FPP for multivalued mappings. The motivation for studying such problems is in its potential application to mathematical models whose constraints can be expressed as FPP and SGEP. This occurs, in particular, in practical problems such as signal processing, network resource allocation, and image recovery. A scenario is in network bandwidth allocation problem for two services in heterogeneous wireless access networks in which the bandwidth of the services is mathematically related (see, for instance, [20,21] and references therein).

In 2016, Suantai et al. [22] introduced the following iterative scheme for solving SEP and FPP of nonspreading multi-valued mapping in Hilbert spaces:

(8) x 1 C arbitrarily , u n = T r n F 1 ( I γ A ( I T r n F 2 ) A ) x n , x n + 1 α n x n + ( 1 α n ) S u n ,

for all n 1 , where C is a nonempty closed convex subset of a real Hilbert space H , { α n } ( 0 , 1 ) , r n ( 0 , ) , S is a nonspreading multivalued mapping, and γ 0 , 1 L such that L is the spectral radius of A A and A is the adjoint of the bounded linear operator A . Under the following conditions on the control sequences:

  1. 0 < lim inf n α n lim sup n α n < 1 ; and

  2. lim inf n r n > 0 ,

the authors proved that the sequence { x n } defined by (8) converges weakly to p F ( S ) SEP ( F 1 , F 2 ) .

Bauschke and Combettes [23] pointed out that in solving OPs, strong convergence of iterative schemes is more desirable than their weak convergence counterparts. Hence, there is a the need to construct iterative schemes that generate a strong convergence sequence.

Takahashi et al. [24] introduced an iterative scheme known as the shrinking projection method for finding a fixed point of a nonexpansive single-valued mapping in Hilbert spaces. The shrinking projection method is a famous method, which plays a significant role in mastering strong convergence for finding fixed points of nonlinear mappings. The method has received much attention due to its applications, and it has been developed to solve many problems, such as, EPs, VIPs, and FPPs in Hilbert spaces (see, for example, [25]).

Very recently, Phuengrattana and Lerkchaiyaphum [26] introduced the following shrinking projection method for solving SGEP and FPP of a countable family of nonexpansive multivalued mappings: for x 1 C and C 1 = C , then

(9) u n = T r n ( F 1 , ϕ 1 ) ( I γ A ( I T r n ( F 2 , ϕ 2 ) ) A ) x n , z n = α n ( 0 ) x n + α n ( 1 ) y n ( 1 ) + + α n ( n ) y n ( n ) , y n ( i ) S i u n , C n + 1 = { p C n : z n p 2 x n p 2 } , x n + 1 = P C n + 1 x 1 , n N .

They proved that if

  1. lim inf n r n > 0 ,

  2. The limits lim n α n ( i ) ( 0 , 1 ) exist for all i 0 ,

then the sequence { x n } generated by (9) converges strongly to P Γ x 1 , where Γ = i = 1 F ( S i ) SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) and S i is a countable family of nonexpansive multivalued mappings.

It is important to point out at this point that the step size γ of the aforementioned algorithm plays an essential role in the convergence properties of iterative methods. The result obtained by the authors in [22,26] and several other related results in the literature involve step size that requires prior knowledge of the operator norm A . One of the drawbacks of such algorithms is that they are usually not easy to implement because they require computation of the operator norm A , which is very difficult if not impossible to calculate or even estimate. Moreover, the step size defined by such algorithms is often very small and deteriorates the convergence rate of the algorithm. In practice, a larger stepsize can often be used to yield better numerical results.

Based on the heavy ball methods of a two-order time dynamical system, Polyak [27] first proposed an inertial extrapolation as an acceleration process to solve the smooth convex minimization problem. The inertial algorithm is a two-step iteration where the next iterate is defined by making use of the previous two iterates. Recently, several researchers have constructed some fast iterative algorithms by using inertial extrapolation (see, e.g., [1,2832]).

Motivated by the above results and the ongoing research interest in this direction, in this paper, we present a new inertial shrinking projection algorithm, which does not require any prior knowledge of the operator norm for finding a common element of the set of solutions of SGEP and the set of common fixed points of a countable family of nonexpansive multivalued mappings in Hilbert spaces. We prove strong convergence theorem for the proposed algorithm and obtain some consequent results. Moreover, we apply our results to solve split mixed variational inequality problem (SMVIP) and split minimization problem (SMP), and we provide numerical examples to illustrate the efficiency of the proposed algorithm in comparison with the existing results in the current literature.

The remaining sections of the paper are organized as follows. In Section 2, we recall some basic definitions and results that will be employed in the convergence analysis of our proposed algorithm. Our new inertial shrinking projection algorithm is presented and analyzed in Section 3, and we also obtain some consequent results. In Section 4, we apply our result to solve SMVIP and SMP. In Section 5, we present some numerical experiments to demonstrate the validity and efficiency of our proposed method in comparison with some recent results in the literature. Finally, in Section 6, we give the concluding remarks.

2 Preliminaries

Let C be a nonempty, closed, and convex subset of a real Hilbert space H with an inner product , and norm . The nearest point projection of H onto C is denoted by P C , that is, x P C x x y for all x H and y C . P C is called the metric projection of H onto C . It is known that P C is firmly nonexpansive, i.e.,

(10) P C x P C y 2 P C x P C y , x y ,

for all x , y H . Moreover, x P C x , y P C x 0 holds for all x H and y C , see [33,34]. We denote the strong convergence and weak convergence of a sequence { x n } to a point x in a Hilbert space H by x n x and x n x , respectively. It is well known [35] that a Hilbert space H satisfies Opial condition, that is, for any sequence { x n } with x n x , the inequality

(11) lim sup n x n x < lim sup n x n y

holds for every y H with y x .

Definition 2.1

A single-valued mapping S : C C is said to be

  • nonexpansive, if and only if

    S x S y x y , x , y C ;

  • δ -inverse strongly monotone [36], if there exists a positive real number δ such that

    x y , S x S y δ S x S y 2 , x , y C ;

  • monotone, if and only if

    y x , S y S x 0 , x , y C .

If S is δ -inverse strongly monotone, for each γ ( 0 , 2 δ ] , it is known [26] that I γ S is a nonexpansive single-valued mapping.

A subset K of H is called proximal if for each x H , there exists y K such that

x y = d ( x , K ) .

We denote by C B ( C ) , C C ( C ) , K ( C ) , and P ( C ) the families of all nonempty closed bounded subsets of C , nonempty closed convex subset of C , nonempty compact subsets of C , and nonempty proximal bounded subsets of C , respectively. The Pompeiu-Hausdorff metric on C B ( C ) is defined by

H ( A , B ) max { sup x A d ( x , B ) , sup y B d ( y , A ) } ,

for all A , B C B ( C ) . Let S : C 2 C be a multivalued mapping. We say that S satisfies the endpoint condition if S p = { p } for all p F ( S ) . For multivalued mappings S i : C 2 C ( i N ) with i = 1 F ( S i ) , we say S i satisfies the common endpoint condition if S i ( p ) = { p } for all i N , p i = 1 F ( S i ) . We recall some basic and useful definitions on multivalued mappings.

Definition 2.2

A multivalued mapping S : C C B ( C ) is said to be nonexpansive if

H ( S x , S y ) x y , x , y C .

The class of nonexpansive multivalued mappings contains the class of nonexpansive single-valued mappings. If S is a nonexpansive single-valued mapping on a closed convex subset of a Hilbert space, then F ( S ) is closed and convex. The closedness of F ( S ) can easily be extended to the multivalued case. However, the convexity of F ( S ) cannot be extended (see, e.g., [37]). But, if S is a nonexpansive multivalued mapping which satisfies the endpoint condition, then F ( S ) is always closed and convex as shown by the following result:

Lemma 2.3

[38] Let C be a nonempty closed convex subset of a real Hilbert space H . Let S : C C B ( C ) be a nonexpansive multivalued mapping with F ( S ) and S p = { p } for each p F ( S ) . Then, F ( S ) is a closed and convex subset of C .

The best approximation operator P S for a multivalued mapping S : C P ( C ) is defined by

P S ( x ) { y S x : x y = d ( x , S x ) } .

It is known that F ( S ) = F ( P S ) and P S satisfies the endpoint condition. Song and Cho [39] gave an example of a best approximation operator P S which is nonexpansive, but where S is not necessarily nonexpansive.

The following results will be needed in the sequel:

Lemma 2.4

[40] In a real Hilbert space H , the following inequalities hold for all x , y H :

  1. x + y 2 x 2 + 2 y , x + y ;

  2. x + y 2 = x 2 + 2 x , y + y 2 ;

  3. x y 2 = x 2 2 x , y + y 2 .

Lemma 2.5

[41] Let H be a Hilbert space, { x n } be a sequence in H , and α 1 , α 2 , , α N be real numbers such that i = 1 N α i = 1 . Then,

(12) i = 1 N α i x i 2 = i = 1 N α i x i 2 1 i , j N α i α j x i x j 2 .

Lemma 2.6

[42] Let H be a Hilbert space, and let { x n } be a sequence in H . Let u , v H be such that lim n x n u and lim n x n v exist. If { x n k } and { x m k } are subsequences of { x n } that converge weakly to u and v respectively, then u = v .

Lemma 2.7

[43] Let C be a nonempty closed convex subset of a real Hilbert space H . Given x , y , z H and a real number α , the set { u C : y u 2 x u 2 + z , u + α } is closed and convex.

Lemma 2.8

[44,45] Let C be a nonempty closed convex subset of a real Hilbert space H , and let P C : H C be the metric projection. Then,

y P C x 2 + x P C x 2 x y 2 , x H , y C .

Assumption 2.9

Let C be a nonempty closed and convex subset of a Hilbert space H 1 . Let F 1 : C × C R and ϕ 1 : C × C R be two bifunctions satisfying the following conditions:

  1. F 1 ( x , x ) = 0 for all x C ,

  2. F 1 is monotone, that is, F 1 ( x , y ) + F 1 ( y , x ) 0 for all x , y C ,

  3. F 1 is upper hemicontinuous, that is, for all x , y , z C , lim t 0 F 1 ( t z + ( 1 t ) x , y ) F 1 ( x , y ) ,

  4. for each x C , y F 1 ( x , y ) is convex and lower semicontinuous,

  5. ϕ 1 ( x , x ) 0 for all x C ,

  6. for each y C , x ϕ 1 ( x , y ) is upper semicontinuous,

  7. for each x C ϕ 1 ( x , y ) is convex and lower semicontinuous,

and assume that for fixed r > 0 and z C , there exists a nonempty compact convex subset K of H 1 and x C K such that

F 1 ( y , x ) + ϕ 1 ( y , x ) + 1 r y x , x z < 0 , y C \ K .

Lemma 2.10

[46] Let C be a nonempty closed and convex subset of a Hilbert space H 1 . Let F 1 : C × C R and ϕ 1 : C × C R be two bifunctions satisfying Assumption 2.9. Assume that ϕ 1 is monotone. For r > 0 and x H 1 , define a mapping T r ( F 1 , ϕ 1 ) : H 1 C as follows:

(13) T r ( F 1 , ϕ 1 ) ( x ) = { z C : F 1 ( z , y ) + ϕ 1 ( z , y ) + 1 r y z , z x 0 , y C } ,

for all x H 1 , Then,

  1. for each x H 1 , T r ( F 1 , ϕ 1 ) ,

  2. T r ( F 1 , ϕ 1 ) is single-valued,

  3. T r ( F 1 , ϕ 1 ) is firmly nonexpansive, that is, for any x , y H 1 ,

    T r ( F 1 , ϕ 1 ) x T r ( F 1 , ϕ 1 ) y 2 T r ( F 1 , ϕ 1 ) x T r ( F 1 , ϕ 1 ) y , x y ,

  4. F ( T r ( F 1 , ϕ 1 ) ) = GEP ( F 1 , ϕ 1 ) ,

  5. GEP ( F 1 , ϕ 1 ) is compact and convex.

Furthermore, assume that F 2 : Q × Q R and ϕ 2 : Q × Q R satisfy Assumption 2.9, where Q is a nonempty closed and convex subset of a Hilbert space H 2 . For all s > 0 and w H 2 , define the mapping T s ( F 2 , ϕ 2 ) : H 2 Q by

(14) T s ( F 2 , ϕ 2 ) ( v ) = w Q : F 2 ( w , d ) + ϕ 2 ( w , d ) + 1 s d w , w v 0 , d Q .

Then, we have

  1. for each v H 2 , T s ( F 2 , ϕ 2 ) ,

  2. T s ( F 2 , ϕ 2 ) is single-valued,

  3. T s ( F 2 , ϕ 2 ) is firmly nonexpansive,

  4. F ( T s ( F 2 , ϕ 2 ) ) = GEP ( F 2 , ϕ 2 ) ,

  5. GEP ( F 2 , ϕ 2 ) is closed and convex,

where GEP ( F 2 , ϕ 2 ) is the solution set of the following GEP: find y Q such that

F 2 ( y , y ) + ϕ 2 ( y , y ) 0 y Q .

Moreover, SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) is a closed and convex set.

Lemma 2.11

[47] Let C be a nonempty closed and convex subset of a Hilbert space H 1 . Let F 1 : C × C R and ϕ 1 : C × C R be two bifunctions satisfying Assumption 2.9, and let T r ( F 1 , ϕ 1 ) be defined as in Lemma 2.10 for r > 0 . Let x , y H 1 and r 1 , r 2 > 0 . Then,

T r 2 ( F 1 , ϕ 1 ) y T r 1 ( F 1 , ϕ 1 ) x y x + r 2 r 1 r 2 T r 2 ( F 1 , ϕ 1 ) y y .

3 Main results

In this section, we state and prove our strong convergence theorem for finding a common element of the set of solutions of SGEP and the set of common fixed points of a countable family of nonexpansive multivalued mappings in real Hilbert spaces.

Theorem 3.1

Let C and Q be nonempty closed convex subsets of real Hilbert spaces H 1 and H 2 , respectively. Let A : H 1 H 2 be a bounded linear operator, and let { S i } be a countable family of nonexpansive multivalued mappings of C into C B ( C ) . Let F 1 , ϕ 1 : C × C R , F 2 , ϕ 2 : Q × Q R be bifunctions satisfying Assumption 2.9. Let ϕ 1 , ϕ 2 be monotone, ϕ 1 be upper hemicontinuous, and F 2 and ϕ 2 be upper semicontinuous in the first argument. Assume that Ω = i = 1 F ( S i ) SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) and S i p = { p } for each p i = 1 F ( S i ) . Let x 0 , x 1 C with C 1 = C , and let { x n } be a sequence generated as follows:

(15) w n = x n + θ n ( x n x n 1 ) , u n = T r n ( F 1 , ϕ 1 ) ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n , z n = α n , 0 u n + i = 1 n α n , i y n , i , y n , i S i u n , C n + 1 = { p C n : z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 } , x n + 1 = P C n + 1 x 1 , n N ,

γ n = τ n ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 if A w n T r n ( F 2 , ϕ 2 ) A w n , γ otherwise ( γ being any nonnegative real number ) ,

where 0 < a τ n b < 1 , { θ n } R , { α n , i } ( 0 , 1 ) , such that i = 0 n α n , i = 1 , and { r n } ( 0 , ) . Suppose that the following conditions hold:

  1. lim inf n r n > 0 ,

  2. the limits lim n α n , i ( 0 , 1 ) exist for all i 0 .

Then, the sequence { x n } generated by (15), converges strongly to P Ω x 1 .

Proof

We divide the proof into several steps as follows:

Step 1: First, we show that { x n } is well-defined for every n N .

By Lemmas 2.3 and 2.10, we have that SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) and i = 1 F ( S i ) are closed and convex subsets of C . Therefore, the solution set Ω is a closed and convex subset of C . By Lemma 2.7, it then follows that C n + 1 is closed and convex for each n N . Let p Ω , then we have p = T r n F 1 , ϕ 1 p and A p = T r n ( F 2 , ϕ 2 ) ( A p ) . Since T r n ( F 1 , ϕ 1 ) is nonexpansive, by Lemma 2.4, we have

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

By the firmly nonexpansivity of I T r n ( F 2 , ϕ 2 ) , we get

(17) w n p , A ( I T r n ( F 2 , ϕ 2 ) ) A w n = A w n A p , ( I T r n ( F 2 , ϕ 2 ) ) A w n = A w n A p , ( I T r n ( f 2 , ϕ 2 ) ) A w n ( I T r n ( F 2 , ϕ 2 ) ) A p ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 .

By substituting (17) into (16), applying the definition of γ n and the condition on τ n , we obtain

(18) u n p 2 w n p 2 + γ n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 2 γ n ( I T r n ( F 1 , ϕ 1 ) ) A w n 2 = w n p 2 γ n [ 2 ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 γ n A ( I T r n ( F 1 , ϕ 1 ) ) A w n 2 ] = w n p 2 γ n ( 2 τ n ) ( I T r n ( F 2 , ϕ 2 ) ) A w n 2

(19) w n p 2 .

Applying Lemma 2.5 and using (19), we have

(20) z n p 2 = α n , 0 u n + i = 1 n α n , i y n , i p 2 α n , 0 u n p 2 + i = 1 n α n , i y n , i p 2 α n , 0 i = 1 n α n , i u n y n , 1 2 = α n , 0 u n p 2 + i = 1 n α n , i d ( y n , i , S i p ) 2 α n , 0 i = 1 n α n , i u n y n , i 2 α n , 0 u n p 2 + i = 1 n α n , i H ( S i u n , S i p ) 2 α n , 0 i = 1 n α n , i u n y n , i 2 α n , 0 u n p 2 + i = 1 n α n , i u n p 2 α n , 0 i = 1 n α n , i u n y n , i 2 u n p 2 α n , 0 i = 1 n α n , i u n y n , i 2

(21) u n p 2 .

Also, by applying Lemma 2.4(iii), we get

(22) w n p 2 = ( x n p θ n ( x n 1 x 1 ) ) 2 = x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 .

By using (19) and (22) in (21), we have

(23) z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 .

This shows that p C n + 1 , and it follows that Ω C n + 1 C n . Therefore, P C n + 1 x 1 is well-defined for every x 1 C and the sequence { x n } is well defined.

Step 2: Next, we show that lim n x n = q for some q C .

We know that Ω is a nonempty closed convex subset of H 1 , then there exists a unique w Ω such that w = P Ω x 1 . Since x n = P C n x 1 and x n + 1 C n + 1 C n for all n N , we have

(24) x n x 1 x n + 1 x 1 , n N .

Similarly, since Ω C n , we have

(25) x n x 1 w x 1 , n N .

Therefore, { x n x 1 } is bounded, and it follows that { x n } is bounded. Consequently, { w n } , { u n } , { z n } , and { y n , i } are bounded. Hence, lim n x n x 1 exists. From the construction of C n , it is clear that x m = P C m x 1 C m C n for m > n 1 . By Lemma 2.8, we have that

(26) x m x n 2 x m x 1 2 x n x 1 2 0 as m , n .

Since lim n x n x 1 exists, then it follows that { x n } is a Cauchy sequence. By the completeness of H 1 and the closedness of C , we have that there exists an element q C such that lim n x n = q .

Step 3: We next show that lim n y n , i u n = 0 for all i N .

From (26), we get

(27) lim n x n + 1 x n = 0 .

Since x n + 1 C n + 1 , we have

z n x n + 1 2 x n x n + 1 2 2 θ n x n x n + 1 , x n 1 x n + θ n 2 x n 1 x n 2 .

By (27), we obtain

(28) lim n z n x n + 1 = 0 .

By applying (27) and (28), we get

(29) z n x n z n x n + 1 + x n + 1 x n 0 , n .

Hence, lim n z n = q .

By the triangle inequality, we have that

w n x n = x n + θ n ( x n x n 1 ) x n x n x n + θ n x n x n 1 .

By (27), we obtain

(30) lim n w n x n = 0 .

Applying (29) and (30), we get

(31) z n w n z n x n + x n w n 0 , n .

From (19) and (20), we obtain

z n p 2 w n p 2 α n , 0 i = 1 n α n , i u n y n , i 2 ,

which implies that

α n , 0 α n , i u n y n , i 2 α n , 0 i = 1 n α n , i u n y n , i 2 w n p 2 z n p 2 w n z n ( w n p + z n p ) .

By the conditions on { α n , i } and using (31), we get

(32) lim n u n y n , i = 0 , i N .

Step 4: We show that u n x n = 0 .

Substituting (18) into (21), we have

(33) z n p 2 w n p 2 γ n ( 2 τ n ) ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 .

From this, we obtain

γ n ( 2 τ n ) ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 w n p 2 z n p 2 w n z n ( w n p + z n p ) .

By the definition of γ n , condition on τ n and (31), we get

τ n ( 2 τ n ) ( I T r n ( F 2 , ϕ 2 ) ) A w n 4 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 0 , n ,

which implies that

( I T r n ( F 2 , ϕ 2 ) ) A w n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 0 , n .

Since A ( I T r n ( F 2 , ϕ 2 ) ) A w n is bounded, it follows that

(34) ( I T r n ( F 2 , ϕ 2 ) ) A w n 0 , n .

From this, we obtain

(35) A ( I T r n ( F 2 , ϕ 2 ) ) A w n A ( I T r n ( F 2 , ϕ 2 ) ) A w n = A ( I T r n ( F 2 , ϕ 2 ) ) A w n 0 , n .

Since T r n ( F 1 , ϕ 1 ) is firmly nonexpansive and I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A is non expansive by invoking Lemma 2.4(ii), we obtain

u n p 2 = T r n ( F 1 , ϕ 1 ) ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n T r n ( F 1 , ϕ 1 ) p 2 T r n ( F 1 , ϕ 1 ) ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n T r n ( F 1 , ϕ 1 ) p , ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n p = u n p , ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n p = 1 2 [ u n p 2 + ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n p 2 u n w n + γ n A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 ] 1 2 [ u n p 2 + w n p 2 ( u n w n 2 + γ n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 + 2 γ n u n w n , A ( I T r n ( F 2 , ϕ 2 ) ) A w n ) ] ,

which implies that

(36) u n p 2 w n p 2 u n w n 2 γ n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n + 2 γ n w n u n , A ( I T r n ( F 2 , ϕ 2 ) ) A w n w n p 2 u n w n 2 + 2 γ n w n u n A ( I T r n ( F 2 , ϕ 2 ) ) A w n .

Substituting (36) into (20), we have

z n p 2 w n p 2 u n w n 2 + 2 γ n w n u n A ( I T r n ( F 2 , ϕ 2 ) ) A w n .

From this, we get

(37) u n w n 2 w n p 2 z n p 2 + 2 γ n w n u n A ( I T r n ( F 2 , ϕ 2 ) ) A w n w n p 2 z n p 2 + 2 γ n M A ( I T r n ( F 2 , ϕ 2 ) ) A w n w n z n ( w n p + z n p ) + 2 γ n M A ( I T r n ( F 2 , ϕ 2 ) ) A w n ,

where M = sup { w n u n : n N } .

By applying (31) and (35) in (37), we get

(38) lim n u n w n = 0 .

Combining this together with (30) and (31), we have

(39) u n z n u n w n + w n z n 0 , n

and

(40) u n x n u n w n + w n x n 0 , n .

Step 5: Next, we show that q i = 1 F ( S i ) .

By (32), for all i N , we get that

(41) lim n d ( u n , S i u n ) lim n u n y n , i = 0 .

For each i N , we have

d ( q , S i q ) q u n + u n y n , i + d ( y n , i , S i q ) q u n + d ( u n , S i u n ) + H ( S i u n , S i q ) 2 q u n + d ( u n , S i u n ) .

By (40), we have that lim n u n = q . Then, it follows from (41) that

d ( q , S i q ) = 0 i N .

This shows that q S i q for all i N , which implies that q i = 1 F ( S i ) .

Step 6: Next, we show that q GEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) .

First, we will show that q GEP ( F 1 , ϕ 1 ) . Since u n = T r n ( F 1 , ϕ 1 ) ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n , then by Lemma 2.10, we obtain

F 1 ( u n , y ) + ϕ 1 ( u n , y ) + 1 r n y u n , u n w n γ n A ( I T r n ( F 2 , ϕ 2 ) ) A w n 0 , y C ,

which implies that

F 1 ( u n , y ) + ϕ 1 ( u n , y ) + 1 r n y u n , u n w n 1 r n y u n , γ n A ( I T r n ( F 2 , ϕ 2 ) ) A w n 0 , y C .

Since F 1 and ϕ 1 are monotone, we have

1 r n y u n , u n w n 1 r n y u n , γ n A ( I T r n ( F 2 , ϕ 2 ) ) A w n F 1 ( y , u n ) + ϕ 1 ( y , u n ) , y C .

By (30) and (38), and lim n x n = q , we obtain lim n u n = q . Then, by Condition (C1), (34), (38), Assumption 2.9, (A4) and (A7), it follows that

0 F 1 ( y , q ) + ϕ 1 ( y , q ) y C .

Let y t = t y + ( 1 t ) q for all t ( 0 , 1 ] and y C . Then, y t C , and thus, F 1 ( y t , q ) + ϕ 1 ( y t , q ) 0 . Therefore, by Assumption 2.9, (A1)–(A7), we obtain

0 F 1 ( y t , y t ) + ϕ 1 ( y t , y t ) t ( F 1 ( y t , y ) + ϕ 1 ( y t , y ) ) + ( 1 t ) ( F 1 ( y t , q ) + ϕ 1 ( y t , q ) ) t ( F 1 ( y t , y ) + ϕ 1 ( y t , y ) ) + ( 1 t ) ( F 1 ( q , y t ) + ϕ 1 ( q , y t ) ) F 1 ( y t , y ) + ϕ 1 ( y t , y ) .

This implies that

F 1 ( y t , y ) + ϕ 1 ( y t , y ) 0 , y C .

Letting t 0 , and by using assumption together with the upper hemicontinuity of ϕ 1 , we obtain

F 1 ( q , y ) + ϕ 1 ( q , y ) 0 , y C .

This implies that q GEP ( F 1 , ϕ 1 ) .

We next show that A q GEP ( F 2 , ϕ 2 ) . Since A is a bounded linear operator, A w n A q . Thus, from (34) we have

(42) T r n ( F 2 , ϕ 2 ) A w n A q .

By the definition of T r n ( F 2 , ϕ 2 ) A w n , we have

F 2 ( T r n ( F 2 , ϕ 2 ) A w n , y ) + ϕ 2 ( T r n ( F 2 , ϕ 2 ) A w n , y ) + 1 r n y T r n ( F 2 , ϕ 2 ) A w n , T r n ( F 2 , ϕ 2 ) A w n A w n 0 , y Q .

Since F 2 and ϕ 2 are upper semicontinuous in the first argument, it follows from (42) that,

F 2 ( A q , y ) + ϕ 2 ( A q , y ) 0 , y Q .

This implies that A q GEP ( F 2 , ϕ 2 ) . Hence, q SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) .

Step 7: Finally, we show that q = P Ω x i .

We know that x n = P c n x 1 and Ω C n , then it follows that x 1 x n , x n p 0 for all p Ω . Hence, we have x 1 q , q p 0 for all p Ω . This implies that q = P Ω x 1 .

Consequently, we can conclude by steps 1–8 that { x n } converges strongly to q = P Ω x 1 as required.□

If ϕ 1 = ϕ 2 = 0 in (2)–(3), then the SGEP reduces to the SEP. Hence, from Theorem 3.1, we obtain the following consequent result for approximating a common element of the set of solutions of SEP and the set of common fixed points of a countable family of nonexpansive multivalued mappings.

Corollary 3.2

Let C and Q be nonempty closed convex subsets of real Hilbert spaces H 1 and H 2 , respectively. Let A : H 1 H 2 be a bounded linear operator, and let { S i } be a countable family of nonexpansive multivalued mappings of C into C B ( C ) . Let F 1 : C × C R , F 2 : Q × Q R be bifunctions satisfying Assumption 2.9. Let F 2 be upper semicontinuous in the first argument. Assume that Ω = i = 1 F ( S i ) SEP ( F 1 , F 2 ) and S i p = { p } for each p i = 1 F ( S i ) . Let x 0 , x 1 C with C 1 = C , and let { x n } be a sequence generated as follows:

(43) w n = x n + θ n ( x n x n 1 ) , u n = T r n F 1 ( I γ n A ( I T r n F 2 ) A ) w n , z n = α n , 0 u n + i = 1 n α n , i y n , i , y n , i S i u n , C n + 1 = { p C n : z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 } , x n + 1 = P C n + 1 x 1 , n N ,

γ n = τ n ( I T r n F 2 ) A w n 2 A ( I T r n F 2 ) A w n 2 if A w n T r n F 2 A w n , γ otherwise ( γ being any nonnegative real number ) ,

where 0 < a τ n b < 1 , { θ n } R , { α n , i } ( 0 , 1 ) , such that i = 0 n α n , i = 1 , and { r n } ( 0 , ) . Suppose that the following conditions hold:

  1. lim inf n r n > 0 ,

  2. the limits lim n α n , i ( 0 , 1 ) exist for all i 0 .

Then, the sequence { x n } generated by (43), converges strongly to P Ω x 1 .

By the properties of the best approximation operator, we obtain the following consequent result.

Corollary 3.3

Let C and Q be nonempty closed convex subsets of real Hilbert spaces H 1 and H 2 , respectively. Let A : H 1 H 2 be a bounded linear operator, and let { S i } be a countable family of multivalued mappings of C into P ( C ) such that P S i is nonexpansive. Let F 1 , ϕ 1 : C × C R , F 2 , ϕ 2 : Q × Q R be bifunctions satisfying Assumption 2.9. Let ϕ 1 , ϕ 2 be monotone, ϕ 1 be upper hemicontinuous, and F 2 and ϕ 2 be upper semicontinuous in the first argument. Assume that Ω = i = 1 F ( S i ) SGEP ( F 1 , ϕ 1 , F 2 , ϕ 2 ) . Let x 0 , x 1 C with C 1 = C , and let { x n } be a sequence generated as follows:

(44) w n = x n + θ n ( x n x n 1 ) , u n = T r n ( F 1 , ϕ 1 ) ( I γ n A ( I T r n ( F 2 , ϕ 2 ) ) A ) w n , z n = α n , 0 u n + i = 1 n α n , i y n , i , y n , i P S i u n , C n + 1 = { p C n : z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 } , x n + 1 = P C n + 1 x 1 , n N ,

γ n = τ n ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 if A w n T r n ( F 2 , ϕ 2 ) A w n , γ otherwise ( γ being any nonnegative real number ) ,

where 0 < a τ n b < 1 , { θ n } R , { α n , i } ( 0 , 1 ) , such that i = 0 n α n , i = 1 , and { r n } ( 0 , ) . Suppose that the following conditions hold:

  1. lim inf n r n > 0 ,

  2. the limits lim n α n , i ( 0 , 1 ) exist for all i 0 .

Then the sequence { x n } generated by (44), converges strongly to P Ω x 1 .

Proof

Since P S i satisfies the common endpoint condition and F ( S i ) = F ( P S i ) for each i N , then the result follows from Theorem 3.1.□

4 Applications

In this section, we apply our results to approximate solutions of some important optimization problems.

4.1 Split mixed variational inequality and fixed point problems

Let H be a real Hilbert space and C be a nonempty closed convex subset of H . Let B : H H be a single-valued mapping and ϕ : C × C R be a bifunction. The mixed variational inequality problem (MVIP) is defined as follows:

(45) Find x C such that y x , B x + ϕ ( x , y ) 0 , y C .

We denote the set of solutions of MVIP by MVI ( C , B , ϕ ) . If we take ϕ = 0 in (45), then the MVIP reduces to the VIP, which is to find a point x C such that y x , B x 0 , y C . The solution set of the VIP is denoted by VI ( C , B ) . Variational inequality was first introduced independently by Fichera [48] and Stampacchia [49]. The VIP is a useful mathematical model that unifies many important concepts in applied mathematics, such as necessary optimality conditions, complementarity problems, network EPs, and systems of nonlinear equations (see [3,50,51]). Several methods have been proposed and analyzed for solving VIP and related OPs, see [5,37,52,53] and references therein.

Here, we apply our result to study the following SMVIP:

(46) Find x i = 1 F ( S i ) such that x x , B 1 x + ϕ 1 ( x , x ) 0 , x C

and such that

(47) y = A x Q solves y y , B 2 y + ϕ 2 ( y , y ) 0 , y Q ,

where C and Q are nonempty closed and convex subsets of real Hilbert spaces H 1 and H 2 , respectively, { S i } is a countable family of nonexpansive multivalued mappings of C into C B ( C ) , A : H 1 H 2 is a bounded linear operator, B 1 : C H 1 , B 2 : Q H 2 are monotone mappings, and ϕ 1 : C × C R , ϕ 2 : Q × Q R are bifunctions satisfying Assumptions (A5)–(A7). Moreover, ϕ 1 , ϕ 2 are monotone with ϕ 1 being upper hemicontinuous and ϕ 2 upper semicontinuous in the first argument. We denote the solution set of problems (46)–(47) by Ω and assume that Ω . By taking F j ( x , y ) y x , B j x , j = 1 , 2 , then the SMVIP (46)–(47) becomes the problem of finding a solution of the SGEP (2)–(3), which is also a solution of the countable family of nonexpansive multivalued mappings { S i } . In addition, all the conditions of Theorem 3.1 are satisfied. Hence, Theorem 3.1 provides a strong convergence theorem for approximating a common solution of SMVIP and fixed point of a countable family of nonexpansive multivalued mappings.

4.2 Split minimization and fixed point problems

Let C and Q be nonempty closed convex subsets of real Hilbert spaces H 1 and H 2 , respectively. Let f : C R , g : Q R be two operators and A : H 1 H 2 be a bounded linear operator, then the SMP is defined as follows:

(48) Find x C such that f ( x ) f ( x ) , x C

and such that

(49) y = A x Q solves g ( y ) g ( y ) , y Q .

We denote the solution set of SMP (48)–(49) by Φ and assume that Φ . For some recent results on iterative algorithms for solving MP, see [54,55] and references therein. Let F 1 ( x , y ) f ( y ) f ( x ) for all x , y C and F 2 ( u , v ) f ( v ) f ( u ) for all u , v Q , and taking ϕ 1 = ϕ 2 = 0 in the SGEP (2)–(3). Then, F 1 ( x , y ) and F 2 ( u , v ) satisfy Assumptions (A1)–(A4) provided f and g are convex and lower semi-continuous on C and Q , respectively. Clearly, ϕ 1 and ϕ 2 satisfy Assumptions (A5)–(A7). Therefore, from Theorem 3.1, we obtain a strong convergence theorem for approximating a common solution of SMP and fixed point problem for a countable family of nonexpansive multivalued mappings in real Hilbert spaces.

5 Numerical experiments

In this section, we present some numerical experiments to illustrate the performance of our algorithm as well as comparing it with Algorithm 9 in the literature. All numerical computations were carried out using Matlab version R2019(b).

We define the sequences { α n , i } as follows for each i N { 0 } and n N :

(50) α n , i = 1 b i + 1 n n + 1 , n > i , 1 n n + 1 k = 1 n 1 b k , n = i , 0 , n < i ,

where b > 1 .

Example 5.1

Let H 1 = H 2 = R and C = Q = [ 0 , 10 ] . Let A : H 1 H 2 be defined by A x = x 3 for all x H 1 . Then, we have that A y = y 3 for all y H 2 . For x C , i N , we define the multivalued mappings S i : C C B ( C ) as follows:

(51) S i ( x ) = 0 , x 10 i , i N .

It can easily be checked that S i is nonexpansive for all i N , S i ( 0 ) = { 0 } , and i = 1 F ( S i ) = { 0 } . We define the bifunctions F 1 , ϕ 1 : C × C R by F 1 ( x , y ) = y 2 + 3 x y 4 x 2 and ϕ 1 ( x , y ) = y 2 x 2 for x , y C , and F 2 , ϕ 2 : Q × Q R by F 2 ( w , v ) = 2 v 2 + w v 3 w 2 and ϕ 2 ( w , v ) = w v for w , v Q . Choose r n = n 3 n + 2 , θ n = 0.8 , and τ n = 0.7 . It can easily be verified that all the conditions of Theorem 3.1 are satisfied with Ω = { 0 } . Now, we compute T r ( F 1 , ϕ 1 ) ( x ) . We find u C such that for all z C

0 F 1 ( u , z ) + ϕ 1 ( u , z ) + 1 r z u , u x = 2 z 2 + 3 u z 5 u 2 + 1 r z u , u x 0 2 r z 2 + 3 r u z 5 r u 2 + ( z u ) ( u x ) = 2 r z 2 + 3 r u z 5 r u 2 + u z x z u 2 + u x = 2 r z 2 + ( 3 r u + u x ) z + ( 5 r u 2 u 2 + u x ) .

Let h ( z ) = 2 r z 2 + ( 3 r u + u x ) z + ( 5 r u 2 u 2 + u x ) . Then, h ( z ) is a quadratic function of z with coefficients a = 2 r , b = 3 r u + u x , and c = 5 r u 2 u 2 + u x . We determine the discriminant Δ of h ( z ) as follows:

(52) Δ = ( 3 r u + u x ) 2 4 ( 2 r ) ( 5 r u 2 u 2 + u x ) = 49 r 2 u 2 + 14 r u 2 14 r u x + u 2 2 u x + x 2 = ( ( 7 r + 1 ) u x ) 2 .

By Lemma 2.10, T r ( F 1 , ϕ 1 ) is single-valued. Hence, it follows that h ( z ) has at most one solution in R . Therefore, from (52), we have that u = x 7 r + 1 . This implies that T r ( F 1 , ϕ 1 ) ( x ) = x 7 r + 1 . Similarly, we compute T r ( F 2 , ϕ 2 ) ( v ) . Find w Q such that for all d Q

T s ( F 2 , ϕ 2 ) ( v ) = w Q : F 2 ( w , d ) + ϕ 2 ( w , d ) + 1 s d w , w v 0 , d Q .

By following similar procedure as above, we obtain w = v + s 5 s + 1 . This implies that T s ( F 2 , ϕ 2 ) ( v ) = v + s 5 s + 1 . We take y n , i = u n 10 i for all i N . Then, Algorithm (15) becomes

w n = x n + θ n ( x n x n 1 ) , u n = w n 7 r n + 1 γ n 15 w n r n + 2 w n 3 r n 9 ( 7 r n + 1 ) ( 5 r n + 1 ) , z n = α n , 0 u n + i = 1 n α n , i u n 10 i , C n + 1 = { p C n : z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 } , x n + 1 = P C n + 1 x 1 , n N ,

where

γ n = τ n ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 if A w n T r n ( F 2 , ϕ 2 ) A w n , γ otherwise ( γ being any nonnegative real number ) .

In this example, we set the parameter b on { α n , i } in (50) to be b = 50 , and we choose different initial values as follows:

Case Ia: x 0 = 11 2 , x 1 = 2 5 ;

Case Ib: x 0 = 8 , x 1 = 1 ;

Case Ic: x 0 = 5 , x 1 = 7 10 ;

Case Id: x 0 = 6 , x 1 = 4 5 .

We compare the performance of our Algorithm (15) with Algorithm (9). The stopping criterion used for our computation is x n + 1 x n < 1 0 4 . We plot the graphs of errors against the number of iterations in each case. The numerical results are reported in Figure 1 and Table 1.

Figure 1 
               Top left: Case Ia; top right: Case Ib; bottom left: Case Ic; and bottom right: Case Id.
Figure 1

Top left: Case Ia; top right: Case Ib; bottom left: Case Ic; and bottom right: Case Id.

Table 1

Numerical results for Example 5.1

Alg. 9 Alg. 15
Case Ia CPU time (s) 2.1794 0.1722
No of iter. 13 3
Case Ib CPU time (s) 2.2136 0.1514
No. of iter. 14 3
Case Ic CPU time (s) 2.2338 0.1517
No of iter. 14 3
Case Id CPU time (s) 2.1757 0.1495
No of iter. 14 3

Example 5.2

Let H 1 = H 2 = L 2 ( [ 0 , 1 ] ) with the inner product defined as

x , y = 0 1 x ( t ) y ( t ) d t , x , y L 2 ( [ 0 , 1 ] ) .

Let

C { x H 1 : a , x = d } ,

where a = 2 t 2 and d 0 . Here, we have

P C ( x ) = x + d a , x a 2 a .

Also, let

Q { x H 2 : c , x e } ,

where c = t 3 and e = 1 , we get

P Q ( x ) = x + max 0 , e c , x c 2 c .

We define F 1 : C × C R and F 2 : Q × Q R by F 1 ( x , y ) = L 1 x , y x and F 2 ( x , y ) = L 2 x , y x , where L 1 x ( t ) = x ( t ) 2 and L 2 x ( t ) = x ( t ) 5 . It can easily be verified that F 1 and F 2 satisfy Conditions (A1)–(A4). Also, take ϕ 1 = ϕ 2 = 0 . Moreover, let A : L 2 ( [ 0 , 1 ] ) L 2 ( [ 0 , 1 ] ) be defined by A x ( t ) = x ( t ) 2 and A y ( t ) = y ( t ) 2 . Then, A is a bounded linear operator. We consider the case for which the countable family of nonexpansive multivalued mappings { S i } are singled-valued. Define a countable family of nonexpansive mappings S i : L 2 ( [ 0 , 1 ] ) L 2 ( [ 0 , 1 ] ) by

( S i x ) ( t ) = 0 1 t i x ( s ) d s for all t [ 0 , 1 ] .

Observe that S i is nonexpansive for each i N . Choose θ n = 0.9 , τ n = 0.8 , r n = n n + 1 . It can easily be checked that all the conditions on the control sequences in Theorem 3.1 are satisfied. Next, we compute T r ( F 1 , ϕ 1 ) ( x ) . We find z C such that for all y C

(53) F 1 ( z , y ) + ϕ 1 ( z , y ) + 1 r y z , z x 0 z 2 , y z + 1 r y z , z x 0 z 2 ( y z ) + 1 r ( y z ) ( z x ) 0 ( y z ) [ r z + 2 ( z x ) ] 0 ( y z ) [ ( r + 2 ) z 2 x ] 0 .

According to Lemma 2.10,

T r ( F 1 , ϕ 1 ) ( x ) = z C : F 1 ( z , y ) + ϕ 1 ( z , y ) + 1 r y z , z x 0 , y C

is single-valued for all x H 1 . Hence, from (53) we have that z = 2 x r + 2 . This implies that T r ( F 1 , ϕ 1 ) ( x ) = 2 x r + 2 . Similarly, we compute T r ( F 2 , ϕ 2 ) ( v ) . We find w Q such that for all d Q

T s ( F 2 , ϕ 2 ) ( v ) = w Q : F 2 ( w , d ) + ϕ 2 ( w , d ) + 1 s d w , w v 0 , d Q .

Following similar procedure as above, we obtain w = 5 v s + 5 . This implies that T s ( F 2 , ϕ 2 ) ( v ) = 5 v s + 5 . Then, Algorithm (15) becomes

w n = x n + θ n ( x n x n 1 ) , u n = 2 w n r n + 2 γ n 2 r n + 5 2 ( r n + 5 ) ( r n + 2 ) w n , z n = α n , 0 u n + i = 1 n α n , i S i u n , C n + 1 = { p C n : z n p 2 x n p 2 2 θ n x n p , x n 1 x n + θ n 2 x n 1 x n 2 } , x n + 1 = P C n + 1 x 1 , n N ,

where

γ n = τ n ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 A ( I T r n ( F 2 , ϕ 2 ) ) A w n 2 if A w n T r n ( F 2 , ϕ 2 ) A w n , γ otherwise ( γ being any nonnegative real number ) .

Here, we set the parameter b on { α n , i } in (50) to be b = 2 , and we choose different initial values as follows:

Case Ia: x 0 = t 3 , x 1 = t 2 + t 4 ;

Case Ib: x 0 = t 2 + t 6 + t 8 , x 1 = t 3 ;

Case Ic: x 0 = t 5 + t 9 + t 11 , x 1 = t 5 ;

Case Id: x 0 = t + t 2 + t 4 + t 6 , x 1 = t 2 + t 7 .

We compare the performance of our Algorithm (15) with Algorithm (9). The stopping criterion used for our computation is x n + 1 x n < 1 0 4 . We plot the graphs of errors against the number of iterations in each case. The numerical results are reported in Figure 2 and Table 2.

Figure 2 
               Top left: Case Ia; top right: Case Ib; bottom left: Case Ic; and bottom right: Case Id.
Figure 2

Top left: Case Ia; top right: Case Ib; bottom left: Case Ic; and bottom right: Case Id.

Table 2

Numerical results for Example 5.2

Alg. 9 Alg. 15
Case Ia CPU time (s) 2.2241 1.3724
No. of iter. 23 19
Case Ib CPU time (s) 2.2247 1.2772
No. of iter. 23 18
Case Ic CPU time (s) 2.1359 1.3056
No of iter. 22 18
Case Id CPU time (s) 2.3458 1.4506
No of iter. 25 20

6 Conclusion

In this article, we proposed a new inertial shrinking projection algorithm with self-adaptive step size for approximating a common solution of SGMEP and FPP for a countable family of nonexpansive multivalued mappings. We proved strong convergence results for the considered problems without a prior knowledge of the operator norm. Finally, we applied our results to solve some other important OPs and presented some numerical experiments to demonstrate the efficiency of our proposed method in comparison with other existing methods. Our results extend and improve several existing results in this direction in the current literature.



Acknowledgments

The authors sincerely thank the anonymous reviewers for their careful reading and constructive comments.

  1. Funding information: Oluwatosin T. Mewomo was supported in part by the National Research Foundation (NRF) of South Africa Incentive Funding for Rated Researchers (Grant Number 119903). Opinions expressed and conclusions arrived are those of the authors and are not necessarily to be attributed to the NRF.

  2. Conflict of interest: The authors declare that they have no competing interests.

References

[1] L. O. Jolaoso , T. O. Alakoya , A. Taiwo , and O. T. Mewomo , Inertial extragradient method via viscosity approximation approach for solving equilibrium problem in Hilbert space, Optimization 70 (2021), no. 2, 387–412, https://doi.org/10.1080/02331934.2020.1716752 . 10.1080/02331934.2020.1716752Suche in Google Scholar

[2] S. Suantai and P. Cholamjiak , Algorithms for solving generalize equilibrium problems and fixed points of nonexpansive semigroups in Hilbert spaces, Optimization 63 (2014), no. 5, 799–815, https://doi.org/10.1080/02331934.2012.684355 . 10.1080/02331934.2012.684355Suche in Google Scholar

[3] T. O. Alakoya , A. Taiwo , O. T. Mewomo , and Y. J. Cho , An iterative algorithm for solving variational inequality, generalized mixed equilibrium, convex minimization and zeros problems for a class of nonexpansive-type mappings, Ann. Univ. Ferrara Sez. VII Sci. Mat. (2021), https://doi.org/10.1007/s11565-020-00354-2 .10.1007/s11565-020-00354-2Suche in Google Scholar

[4] T. O. Alakoya , L. O. Jolaoso , and O. T. Mewomo , A general iterative method for finding common fixed point of finite family of demicontractive mappings with accretive variational inequality problems in Banach spaces, Nonlinear Stud. 27 (2020), no. 1, 1–24. Suche in Google Scholar

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

[6] P. Cholamjiak and S. Suantai , Iterative methods for solving equilibrium problems, variational inequalities and fixed points of nonexpansive semigroups, J. Glob. Optim. 57 (2013), 1277–1297. 10.1007/s10898-012-0029-7Suche in Google Scholar

[7] D. V. Hieu , L. D. Muu , and P. K. Anh , Parallel hybrid extragradient methods for pseudomotone equilibrium problems and nonexpansive mappings, Numer. Algorithms 73 (2016), 197–217. 10.1007/s11075-015-0092-5Suche in Google Scholar

[8] C. Izuchukwu , G. N. Ogwo , and O. T. Mewomo , An inertial method for solving generalized split feasibility problems over the solution set of monotone variational inclusions, Optimization (2020), https://doi.org/10.1080/02331934.2020.1808648 .10.1080/02331934.2020.1808648Suche in Google Scholar

[9] T. O. Alakoya , L. O. Jolaoso , and O. T. Mewomo , Strong convergence theorems for finite families of pseudomonotone equilibrium and fixed point problems in Banach spaces, Afr. Mat. (2021), https://doi.org/10.1007/s13370-020-00869-z .10.1007/s13370-020-00869-zSuche in Google Scholar

[10] L. O. Jolaoso , A. Taiwo , T. O. Alakoya , and O. T. Mewomo , A unified algorithm for solving variational inequality and fixed point problems with application to the split equality problem, Comput. Appl. Math. 39 (2020), 38, https://doi.org/10.1007/s40314-019-1014-2 .10.1007/s40314-019-1014-2Suche in Google Scholar

[11] L. O. Jolaoso , A. Taiwo , T. O. Alakoya , and O. T. Mewomo , Strong convergence theorem for solving pseudo-monotone variational inequality problem using projection method in a reflexive Banach space, J. Optim. Theory Appl. 185 (2020), no. 3, 744–766. 10.1007/s10957-020-01672-3Suche in Google Scholar

[12] G. N. Ogwo , C. Izuchukwu , K. O. Aremu , and O. T. Mewomo , A viscosity iterative algorithm for a family of monotone inclusion problems in an Hadamard space, Bull. Belg. Math. Soc. Simon Stevin 27 (2020), 127–152. 10.36045/bbms/1590199308Suche in Google Scholar

[13] A. O.-E. Owolabi , T. O. Alakoya , A. Taiwo , and O. T. Mewomo , A new inertial-projection algorithm for approximating common solution of variational inequality and fixed point problems of multivalued mappings, Numer. Algebra Control Optim. (2021), https://doi.org/10.3934/naco.2021004 .10.3934/naco.2021004Suche in Google Scholar

[14] A. Taiwo , A. O.-E. Owolabi , L. O. Jolaoso , O. T. Mewomo , and A. Gibali , A new approximation scheme for solving various split inverse problems, Afr. Mat. (2020), https://doi.org/10.1007/s13370-020-00832-y .10.1007/s13370-020-00832-ySuche in Google Scholar

[15] A. Taiwo , T. O. Alakoya , and O. T. Mewomo , Halpern-type iterative process for solving split common fixed point and monotone variational inclusion problem between Banach spaces, Numer. Algorithms 86 (2020), 1359–1389, https://doi.org/10.1007/s11075-020-00937-2 .10.1007/s11075-020-00937-2Suche in Google Scholar

[16] A. Taiwo , L. O. Jolaoso , and O. T. Mewomo , Inertial-type algorithm for solving split common fixed-point problem in Banach spaces, J. Sci. Comput. 86 (2020), 12, https://doi.org/10.1007/s10915-020-01385-9 .10.1007/s10915-020-01385-9Suche in Google Scholar

[17] A. Taiwo , L. O. Jolaoso , O. T. Mewomo , and A. Gibali , On generalized mixed equilibrium problem with alpha-beta-eta bifunction and mu-tau monotone mapping, J. Nonlinear Convex Anal. 21 (2020), no. 6, 1381–1401. Suche in Google Scholar

[18] A. Taiwo , T. O. Alakoya , and O. T. Mewomo , Strong convergence theorem for solving equilibrium problem and fixed point of relatively nonexpansive multi-valued mappings in a Banach space with applications, Asian-Eur. J. Math. (2020), https://doi.org/10.1142/S1793557121501370 .10.1142/S1793557121501370Suche in Google Scholar

[19] K. R. Kazmi and S. H. Rizvi , Iterative approximation of a common solution of a split generalized equilibrium problem and a fixed point problem for nonexpansive semigroup, Math. Sci. (Springer) 7 (2013), 1. 10.1186/2251-7456-7-1Suche in Google Scholar

[20] H. Iiduka , Fixed point optimization algorithm and its application to network bandwidth allocation, J. Comp. App. Math. 236 (2012), 1733–1742. 10.1016/j.cam.2011.10.004Suche in Google Scholar

[21] C. Luo , H. Ji , and Y. Li , Utility-based multi-service bandwidth allocation in the 4G heterogeneous wireless networks , IEEE Wireless Communication and Networking Conference , 2009, https://doi.org/10.1109/WCNC.2009.4918017 .10.1109/WCNC.2009.4918017Suche in Google Scholar

[22] S. Suantai , P. Cholamjiak , Y. J. Cho , and W. Cholamjiak , On solving split equilibrium problems and fixed point problems of nonspreading multi-valued mappings in Hilbert spaces, Fixed Point Theory Appl. 2016 (2016), 35. 10.1186/s13663-016-0509-4Suche in Google Scholar

[23] H. H. Bauschke and P. L. Combettes , A weak-to-strong convergence principle for Fejer-monotone methods in Hilbert spaces, Math. Oper. Res. 26 (2001), no. 2, 248–264. 10.1287/moor.26.2.248.10558Suche in Google Scholar

[24] W. Takahashi , Y. Takeuchi , and R. Kubota , Strong convergence theorems by hybrid methods for families of nonexpansive mappings in Hilbert spaces, J. Math. Anal. Appl. 341 (2008), 276–286. 10.1016/j.jmaa.2007.09.062Suche in Google Scholar

[25] Y. Kimura , Convergence of a sequence of sets in a Hadamard space and the shrinking projection method for a real Hilbert ball, Abstr. Appl. Anal. 2010 (2010), 582475, https://doi.org/10.1155/2010/582475 .10.1155/2010/582475Suche in Google Scholar

[26] W. Phuengrattana and K. Lerkchaiyaphum , On solving the split generalized equilibrium problem and the fixed point problem for a countable family of nonexpansive multivalued mappings, Fixed Point Theory Appl. 2018 (2018), 6, https://doi.org/10.1186/s13663-018-0631-6 .10.1186/s13663-018-0631-6Suche in Google Scholar

[27] B. T. Polyak , Some methods of speeding up the convergence of iterative methods, Zh. Vychisl. Mat. Mat. Fiz. 4 (1964), 1–17. 10.1016/0041-5553(64)90137-5Suche in Google Scholar

[28] T. O. Alakoya , L. O. Jolaoso , and O. T. Mewomo , Modified inertial subgradient extragradient method with self adaptive stepsize for solving monotone variational inequality and fixed point problems, Optimization 70 (2020), no. 3, 545–574, https://doi.org/10.1080/02331934.2020.1723586 .10.1080/02331934.2020.1723586Suche in Google Scholar

[29] T. O. Alakoya , L. O. Jolaoso , and O. T. Mewomo , Two modifications of the inertial Tseng extragradient method with self-adaptive step size for solving monotone variational inequality problems, Demonstr. Math. 53 (2020), 208–224, https://doi.org/10.1515/dema-2020-0013 .10.1515/dema-2020-0013Suche in Google Scholar

[30] P. Cholamjiak and Y. Shehu , Inertial forward-backward splitting method in Banach spaces with application to compressed sensing, Appl. Math. 64 (2019), 409–435. 10.21136/AM.2019.0323-18Suche in Google Scholar

[31] Q. Dong , D. Jiang , P. Cholamjiak , and Y. Shehu , A strong convergence result involving an inertial forward-backward algorithm for monotone inclusions, J. Fixed Point Theory Appl. 19 (2017), 3097–3118, https://doi.org/10.1007/s11784-017-0472-7 .10.1007/s11784-017-0472-7Suche in Google Scholar

[32] A. Gibali , L. O. Jolaoso , O. T. Mewomo , and A. Taiwo , Fast and simple Bregman projection methods for solving variational inequalities and related problems in Banach spaces, Results Math. 75 (2020), 179, https://doi.org/10.1007/s00025-020-01306-0 .10.1007/s00025-020-01306-0Suche in Google Scholar

[33] E. C Godwin , C. Izuchukwu , and O. T. Mewomo , An inertial extrapolation method for solving generalized split feasibility problems in real Hilbert spaces, Boll. Unione Mat. Ital. (2021), https://doi.org/10.1007/s40574-020-00272-3 .10.1007/s40574-020-00272-3Suche in Google Scholar

[34] C. Izuchukwu , A. A. Mebawondu , and O. T. Mewomo , A new method for solving split variational inequality problems without co-coerciveness, J. Fixed Point Theory Appl. 22 (2020), 98, https://doi.org/10.1007/s11784-020-00834-0 .10.1007/s11784-020-00834-0Suche in Google Scholar

[35] Z. Opial , Weak convergence of the sequence of successive approximation for nonexpansive mappings, Bull. Amer. Math. Soc. 73 (1967), 591–597. 10.1090/S0002-9904-1967-11761-0Suche in Google Scholar

[36] H. Iiduka and W. Takahashi , Weak convergence theorem by Cesáro means for nonexpansive mappings and inverse-strongly monotone mappings, J. Nonlinear Convex Anal. 7 (2006), 105–113. Suche in Google Scholar

[37] A. R. Khan , Properties of fixed point set of a multivalued map, J. Appl. Math. Stoch. Anal. 3 (2005), 323–331. 10.1155/JAMSA.2005.323Suche in Google Scholar

[38] W. Cholamjiak and S. Suantai , A hybrid method for a countable family of multivalued maps, equilibrium problems, and variational inequality problems, Discrete Dyn. Nat. Soc. 2010 (2010), 349158, https://doi.org/10.1155/2010/349158 .10.1155/2010/349158Suche in Google Scholar

[39] Y. Song and Y. J. Cho , Some note on Ishikawa iteration for multivalued mappings, Bull. Korean Math. Soc. 48 (2011), no. 3, 575–584. 10.4134/BKMS.2011.48.3.575Suche in Google Scholar

[40] A. Taiwo , L. O. Jolaoso , and O. T. Mewomo , Viscosity approximation method for solving the multiple-set split equality common fixed-point problems for quasi-pseudocontractive mappings in Hilbert Spaces, J. Ind. Manag. Optim. (2020), https://doi.org/10.3934/jimo.2020092 .10.3934/jimo.2020092Suche in Google Scholar

[41] L. O. Jolaoso , T. O. Alakoya , A. Taiwo , and O. T. Mewomo , A parallel combination extragradient method with Armijo line searching for finding common solution of finite families of equilibrium and fixed point problems, Rend. Circ. Mat. Palermo (2) 69 (2020), no. 3, 711–735, https://doi.org/10.1007/s12215-019-00431-2 .10.1007/s12215-019-00431-2Suche in Google Scholar

[42] S. Suantai , Weak and strong convergence criteria of Noor iterations for asymptotically nonexpansive mappings, J. Math. Anal. Appl. 311 (2005), 506–517. 10.1016/j.jmaa.2005.03.002Suche in Google Scholar

[43] C. Martinez-Yanesa and H. K. Xu , Strong convergence of the CQ method for fixed point iteration processes, Nonlinear Anal. 64 (2006), 2400–2411. 10.1016/j.na.2005.08.018Suche in Google Scholar

[44] K. Goebel and S. Reich , Uniform Convexity, Hyperbolic Geometry, and Nonexpansive Mappings, Dekker, New York, 1984. Suche in Google Scholar

[45] K. Nakajo and W. Takahashi , Strong convergence theorems for nonexpansive mappings and nonexpansive semigroups, J. Math. Anal. Appl. 279 (2003), 372–379. 10.1016/S0022-247X(02)00458-4Suche in Google Scholar

[46] Z. Ma , L. Wang , S. S. Chang , and W. Duan , Convergence theorems for split equality mixed equilibrium problems with applications, Fixed Point Theory Appl. 2015 (2015), 31, https://doi.org/10.1186/s13663-015-0281-x .10.1186/s13663-015-0281-xSuche in Google Scholar

[47] F. Cianciaruso , G. Marino , L. Muglia , and Y. Yao , A hybrid projection algorithm for finding solutions of mixed equilibrium problem and variational inequality problem, Fixed Point Theory Appl. 2010 (2009), 383740, https://doi.org/10.1155/2010/383740 .10.1155/2010/383740Suche in Google Scholar

[48] G. Fichera , Sul problema elastostatico di Signorini con ambigue condizioni al contorno, Atti Accad. Naz. Lincei VIII. Ser. Rend. Cl. Sci. Fis. Mat. Nat. 34 (1963), 138–142. Suche in Google Scholar

[49] G. Stampacchia , Formes bilineaires coercitives sur les ensembles convexes, C. R. Acad. Sci. Paris 258 (1964), 4413–4416. Suche in Google Scholar

[50] A. Gibali , S. Reich , and R. Zalas , Outer approximation methods for solving variational inequalities in Hilbert space, Optimization 66 (2017), no. 3, 417–437, https://doi.org/10.1080/02331934.2016.1271800 .10.1080/02331934.2016.1271800Suche in Google Scholar

[51] G. Kassay , S. Reich , and S. Sabach , Iterative methods for solving systems of variational inequalities in reflexive Banach spaces, SIAM J. Optim. 21 (2011), 1319–1344. 10.1137/110820002Suche in Google Scholar

[52] T. O. Alakoya , L. O. Jolaoso , and O. T. Mewomo , A self adaptive inertial algorithm for solving split variational inclusion and fixed point problems with applications, J. Ind. Manag. Optim. (2020), https://doi.org/10.3934/jimo.2020152 .10.3934/jimo.2020152Suche in Google Scholar

[53] S. H. Khan , T. O. Alakoya , and O. T. Mewomo , Relaxed projection methods with self-adaptive step size for solving variational inequality and fixed point problems for an infinite family of multivalued relatively nonexpansive mappings in Banach spaces, Math. Comput. Appl. 25 (2020), no. 3, 54, https://doi.org/10.3390/mca25030054 .10.3390/mca25030054Suche in Google Scholar

[54] K. O. Aremu , H. A. Abass , C. Izuchukwu , and O. T. Mewomo , A viscosity-type algorithm for an infinitely countable family of (f,g) -generalized k-strictly pseudononspreading mappings in CAT(0) spaces, Analysis 40 (2020), no. 1, 19–37, https://doi.org/10.1515/anly-2018-0078 .10.1515/anly-2018-0078Suche in Google Scholar

[55] K. O. Aremu , C. Izuchukwu , G. N. Ogwo , and O. T. Mewomo , Multi-step Iterative algorithm for minimization and fixed point problems in p-uniformly convex metric spaces, J. Ind. Manag. Optim. (2020), https://doi.org/10.3934/jimo.2020063 .10.3934/jimo.2020063Suche in Google Scholar

Received: 2020-09-03
Revised: 2021-02-21
Accepted: 2021-03-02
Published Online: 2021-04-16

© 2021 Musa A. Olona et al., published by De Gruyter

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

Artikel in diesem Heft

  1. Regular Articles
  2. Graded I-second submodules
  3. Corrigendum to the paper “Equivalence of the existence of best proximity points and best proximity pairs for cyclic and noncyclic nonexpansive mappings”
  4. Solving two-dimensional nonlinear fuzzy Volterra integral equations by homotopy analysis method
  5. Chandrasekhar quadratic and cubic integral equations via Volterra-Stieltjes quadratic integral equation
  6. On q-analogue of Janowski-type starlike functions with respect to symmetric points
  7. Inertial shrinking projection algorithm with self-adaptive step size for split generalized equilibrium and fixed point problems for a countable family of nonexpansive multivalued mappings
  8. On new stability results for composite functional equations in quasi-β-normed spaces
  9. Sampling and interpolation of cumulative distribution functions of Cantor sets in [0, 1]
  10. Meromorphic solutions of the (2 + 1)- and the (3 + 1)-dimensional BLMP equations and the (2 + 1)-dimensional KMN equation
  11. On the equivalence between weak BMO and the space of derivatives of the Zygmund class
  12. On some fixed point theorems for multivalued F-contractions in partial metric spaces
  13. On graded Jgr-classical 2-absorbing submodules of graded modules over graded commutative rings
  14. On almost e-ℐ-continuous functions
  15. Analytical properties of the two-variables Jacobi matrix polynomials with applications
  16. New soft separation axioms and fixed soft points with respect to total belong and total non-belong relations
  17. Pythagorean harmonic summability of Fourier series
  18. More on μ-semi-Lindelöf sets in μ-spaces
  19. Range-Kernel orthogonality and elementary operators on certain Banach spaces
  20. A Cauchy-type generalization of Flett's theorem
  21. A self-adaptive Tseng extragradient method for solving monotone variational inequality and fixed point problems in Banach spaces
  22. Robust numerical method for singularly perturbed differential equations with large delay
  23. Special Issue on Equilibrium Problems: Fixed-Point and Best Proximity-Point Approaches
  24. Strong convergence inertial projection algorithm with self-adaptive step size rule for pseudomonotone variational inequalities in Hilbert spaces
  25. Two strongly convergent self-adaptive iterative schemes for solving pseudo-monotone equilibrium problems with applications
  26. Some aspects of generalized Zbăganu and James constant in Banach spaces
  27. An iterative approximation of common solutions of split generalized vector mixed equilibrium problem and some certain optimization problems
  28. Generalized split null point of sum of monotone operators in Hilbert spaces
  29. Comparison of modified ADM and classical finite difference method for some third-order and fifth-order KdV equations
  30. Solving system of linear equations via bicomplex valued metric space
  31. Special Issue on Computational and Theoretical Studies of free Boundary Problems and their Applications
  32. Dynamical study of Lyapunov exponents for Hide’s coupled dynamo model
  33. A statistical study of COVID-19 pandemic in Egypt
  34. Global existence and dynamic structure of solutions for damped wave equation involving the fractional Laplacian
  35. New class of operators where the distance between the identity operator and the generalized Jordan ∗-derivation range is maximal
  36. Some results on generalized finite operators and range kernel orthogonality in Hilbert spaces
  37. Structures of spinors fiber bundles with special relativity of Dirac operator using the Clifford algebra
  38. A new iteration method for the solution of third-order BVP via Green's function
  39. Numerical treatment of the generalized time-fractional Huxley-Burgers’ equation and its stability examination
  40. L -error estimates of a finite element method for Hamilton-Jacobi-Bellman equations with nonlinear source terms with mixed boundary condition
  41. On shrinkage estimators improving the positive part of James-Stein estimator
  42. A revised model for the effect of nanoparticle mass flux on the thermal instability of a nanofluid layer
  43. On convergence of explicit finite volume scheme for one-dimensional three-component two-phase flow model in porous media
  44. An adjusted Grubbs' and generalized extreme studentized deviation
  45. Existence and uniqueness of the weak solution for Keller-Segel model coupled with Boussinesq equations
  46. Special Issue on Advanced Numerical Methods and Algorithms in Computational Physics
  47. Stability analysis of fractional order SEIR model for malaria disease in Khyber Pakhtunkhwa
Heruntergeladen am 22.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/dema-2021-0006/html?lang=de
Button zum nach oben scrollen