**1. Introduction**

The classical variational inequality problem (VIP) is defined in a real Hilbert space setting as: find *x* ∈ *D* such that

$$
\langle Tx, y - x \rangle \ge 0 \,\,\forall \, y \in D,\tag{1}
$$

where *T* is a nonlinear operator defined on *D* and *D* is a nonempty subset of the Hilbert space. The theory of VIP combines concepts of nonlinear operators and convex analysis in such a way that it generalizes both and is used to model nonlinear problems of physical phenomena in economics, sciences and engineering (see [1] for details). The VIP (1) was first introduced in finite dimensional spaces by Stampacchia [2], and since then researchers have devoted a lot of attention to VIP in finite and infinite dimensional spaces (see [3–8] and other references therein). Another form of the VIP widely studied in real Hilbert space settings (see [9,10] and the references therein) is defined as: find *x* ∈ *D* such that

$$<\langle \mathbf{x} - T\mathbf{x}, y - \mathbf{x} \rangle \ge 0 \,\forall \, y \in D. \tag{2}$$

Several algorithms have been developed for solving VIP and related optimization problems in Hilbert and Banach spaces (see [3,7,11–16] and other references therein). It is well known that many of the problems in practical applications of optimization are constrained optimization problems, where the constraints are nonlinear, non-convex and non-smooth. Hence, it is pertinent to extend the study of these optimization problems to the nonlinear space settings, due to its ability to see non-convex and non-smooth constrained optimization problems as convex, smooth and unconstrained problems. For this reason, Németh [17] introduced and generalized the existence and uniqueness results of the classical VIP from Euclidean spaces to complete Riemannian manifolds. This development led to increasing interest from researchers in the study of VIPs and their generalizations in nonlinear spaces (see [18–27] and other references therein). Despite the increasing attention of researchers in this direction, little attention has been given to other more general nonlinear spaces apart from the Riemannian manifolds. In 2015, Khatibzadeh and Ranjbar [28] extended the study of VIP (2) to the framework of complete CAT(0) spaces. They formulated the VIP as follows:

$$\text{Find } \mathbf{x} \in D \text{ such that } \langle \overrightarrow{Txx}, \overrightarrow{xy} \rangle \ge 0 \,\forall \, y \in D,\tag{3}$$

where *D* is a nonempty, closed and convex subset of an Hadamard space *X* and *T* is a nonexpansive mapping. They established the existence of solutions for the VIP (3) and also employed an inexact proximal point algorithm to approximate a fixed point of the nonexpansive mapping which is also a solution of (3). They obtained convergence result for the algorithm under suitable conditions on the control sequences. Very recently, Alizadeh-Dehghan-Moradlou [29] introduced the notion of inverse strongly monotone mappings in metric spaces as follows: Let *D* be a nonempty subset of a metric space *X* and *T* : *D* → *X* be a mapping. *T* is called *α*-inverse strongly monotone if there exists *α* > 0 such that

$$d^2(\mathbf{x}, \mathbf{y}) - \langle \overrightarrow{\operatorname{Tx}T\tilde{y}}, \overrightarrow{\operatorname{xy}} \rangle \le a \Phi\_T(\mathbf{x}, \mathbf{y}), \forall \ \mathbf{x}, \mathbf{y} \in D,\tag{4}$$

where <sup>Φ</sup>*T*(*x*, *<sup>y</sup>*) = *<sup>d</sup>*2(*x*, *<sup>y</sup>*) + *<sup>d</sup>*2(*Tx*, *Ty*) <sup>−</sup> <sup>2</sup> −−→ *TxTy*, −*xy* →.

Additionally, in [29], the authors studied the VIP (3) in an Hadamard space, where *T* is an inverse strongly monotone mapping. They established the existence of solutions for the VIP (3) associated with an inverse strongly monotone mapping. Furthermore, they introduced the following iterative algorithm to solve the VIP (3): for arbitrary *x*<sup>1</sup> ∈ *D*, the sequence {*xn*} is generated by

$$\begin{cases} y\_n = P\_D[\beta\_n \mathbf{x}\_n \oplus (1 - \beta\_n) T \mathbf{x}\_n]\_\prime \\ \mathbf{x}\_{n+1} = P\_D[a\_n \mathbf{x}\_n \oplus (1 - a\_n) S \mathbf{x}\_n]\_\prime, n \ge 1, \end{cases} \tag{5}$$

where {*βn*} and {*αn*} ∈ (0, 1), *PD* is a metric projection, *T* is inverse strongly monotone and *S* is nonexpansive mapping. They obtained that Algorithm (5) Δ-converges to a solution of the VIP (3), which is also a fixed point of the nonexpansive mapping *S*.

Very recently, Osisiogu et al. [30] proposed and studied the following Halpern-type algorithm in Hadamard spaces for approximating a common solution of a finite family of the VIP (3):

$$\begin{cases} y\_n = \bigoplus\_{i=1}^N \beta\_i P\_D T\_{\lambda\_i} x\_n \\ x\_{n+1} = a\_n u \oplus \beta\_n x\_n \oplus \gamma\_n Sy\_n \; \forall \; n \ge 1, \end{cases} \tag{6}$$

where *Tλ<sup>i</sup>* = (1 − *λi*)*x* ⊕ *λiTx*, 0 < *λ<sup>i</sup>* < 2*αi*, for each *i* = 1, 2, ··· , *N*, {*βn*}, {*αn*} and {*γn*} ∈ (0, 1), *T* is an inverse strongly monotone mapping and *S* is a nonexpansive mapping. They obtained a strong convergence result of Algorithm (6) under some suitable conditions.

Motivated by the results of Khatibzadeh and Ranjbar [28], Alizadeh-Dehghan-Moradlou [29] and Osisiogu et al. [30], we propose and study a viscosity iterative algorithm (from the fact that viscosity-type algorithms converge faster than Halpern-type algorithms and also Halpern-type algorithms are particular cases of viscosity-type algorithms, see [31,32]) that comprises of a finite family of inverse strongly monotone mappings (3) and a finite family of Lipchitz demicontractive mappings in an Hadamard space. Additionally, we establish that the proposed algorithm converges strongly to a common solution of a finite family of VIPs, which is also a common fixed point of a finite family of Lipchitz demicontractive mappings in the framework of Hadamard spaces. Furthermore, we provide a

numerical experiment to demonstrate the applicability of our results. Our result generalizes the works of Alizadeh-Dehghan-Moradlou [29] and Osisiogu et al. [30] and other similar works in literature.

#### **2. Preliminaries**

Let (*X*, *d*) be a metric space, *x*, *y* ∈ *X* and *I* = [0, *d*(*x*, *y*)] be an interval. A curve *c* (or simply a geodesic path) joining *x* to *y* is an isometry *c* : *I* → *X* such that *c*(0) = *x*, *c*(*d*(*x*, *y*)) = *y* and *d*(*c*(*t*), *c*(*t* )) = |*t* − *t* | for all *t*, *t* ∈ *I*. The image of a geodesic path is called the geodesic segment, which is denoted by [*x*, *y*] whenever it is unique. We say that a metric space *X* is a geodesic space if for every pair of points *x*, *y* ∈ *X*, there is a minimal geodesic from *x* to *y*. A geodesic triangle Δ(*x*1, *x*2, *x*3) in a geodesic metric space (*X*, *d*) consists of three vertices (points in *X*) with geodesic segments between each pair of vertices. For any geodesic triangle, there is a comparison (Alexandrov) triangle <sup>Δ</sup>¯ <sup>⊂</sup> <sup>R</sup><sup>2</sup> such that *<sup>d</sup>*(*xi*, *xj*) = *<sup>d</sup>*R<sup>2</sup> (*x*¯*i*, *<sup>x</sup>*¯*j*) for *<sup>i</sup>*, *<sup>j</sup>* ∈ {1, 2, 3}. A geodesic space *<sup>X</sup>* is a CAT(0) space if the distance between arbitrary pair of points on a geodesic triangle Δ does not exceed the distance between its pair of corresponding points on its comparison triangle Δ¯ . If Δ is a geodesic triangle and Δ¯ is its comparison triangle in *X*, then Δ is said to satisfy the CAT(0) inequality for all points *x*, *y* of Δ and *x*¯, *y*¯ of Δ¯ , if

$$d(\mathfrak{x}, \mathfrak{y}) = d\_{\mathbb{R}^2}(\mathfrak{x}, \mathfrak{y}). \tag{7}$$

Let *x*, *y*, *z* be points in *X* and *y*<sup>0</sup> be the midpoint of the segment [*y*, *z*]; then the CAT(0) inequality implies

$$d^2(\mathbf{x}, y\_0) \le \frac{1}{2} d^2(\mathbf{x}, y) + \frac{1}{2} d^2(\mathbf{x}, z) - \frac{1}{4} d(y, z). \tag{8}$$

Inequality (8) is known as CN inequality of Bruhat and Tits [33]. A geodesic space *X* is said to be a CAT(0) space if all geodesic triangles satisfy the CN inequality. Equivalently, *X* is called a CAT(0) space if and only if it satisfies the CN inequality. Examples of CAT(0) spaces includes Hadamard manifold, R-trees [34], pre-Hilbert spaces [35], hyperbolic metric [36] and Hilbert balls [37].

Let *D* be a nonempty subset of a metric space (*X*, *d*). A point *x* ∈ *X* is called a fixed point of a nonlinear mapping *T* : *D* → *X*, if *x* = *Tx*. We denote by *F*(*T*) the set of fixed points of *T*. The mapping *T* is said to be:

1. *L*-Lipschitz, if there exists *L* > 0 such that

$$d(Tx, Ty) \le L d(x, y), \forall \ x, y \in X; y \in X$$

if *L* = 1, then *T* is called nonexpansive;

2. Firmly nonexpansive (see [38]), if

$$d^2(Tx, Ty) \le \langle \overrightarrow{Tx}\overrightarrow{Ty}, \overrightarrow{xy} \rangle \,\,\forall \, x, y \in X; y$$

3. Quasi-nonexpansive, if *F*(*T*) = ∅ and

$$d(Tx, p) \le d(x, p) \,\forall \, x \in X \text{ and } p \in F(T);$$

4. *k*-demicontractive, if *F*(*T*) = ∅ and there exists *k* ∈ [0, 1) such that

$$d^2(T\mathbf{x}, p) \le d^2(\mathbf{x}, p) + kd^2(\mathbf{x}, T\mathbf{x}) \,\forall \, \mathbf{x}, \boldsymbol{y} \in X, \text{ and } p \in F(T).$$

Obviously, the class of quasi-nonexpansive are *k*-demicontractive mappings. However, the converse is not true (see [39] Example 1.1).

**Definition 1.** *[40] Let a pair* (*a*, *<sup>b</sup>*) <sup>∈</sup> *<sup>X</sup>* <sup>×</sup> *X, denoted by* −→*ab*, *be called a vector in <sup>X</sup>* <sup>×</sup> *X. The quasilinearization map* ., . : (*X* × *X*) × (*X* × *X*) → R *is defined by*

$$
\langle \overrightarrow{ab}, \overrightarrow{cd} \rangle = \frac{1}{2} (d^2(a,d) + d^2(b,c) - d^2(a,c) - d^2(b,d)), \forall \, a, b, c, d \in X. \tag{9}
$$

*It is easy to see that* −→*ba*, −→*cd* <sup>=</sup> −−→*ab*, −→*cd*, −→*ab*, −→*cd* <sup>=</sup> −→*ae*, −→*cd* <sup>+</sup> −→*eb* , −→*cd and* −→*ab*, −→*cd* <sup>=</sup> −→*cd*, −→*ab for all a*, *b*, *c*, *d*,*e* ∈ *X. Furthermore, a geodesic space X is said to satisfy the Cauchy–Schwarz inequality if*

$$\langle \overrightarrow{ab}, \overrightarrow{cd} \rangle \le d(a, b)d(c, d), \forall \, a, b, c, d \in X.$$

It is known from [41] that a geodesically connected metric space is a CAT(0) space if and only if it satisfies the Cauchy–Schwarz inequality.

We state some known and useful results which will be needed in the proof of our main results. In the sequel, we denote strong and Δ-convergence by "→ " and " " respectively.

Let {*xn*} be a bounded sequence in *X* and *r*(., {*xn*}) : *X* → [0, ∞) be a continuous mapping defined by *r*(*x*, {*xn*}) = lim sup *n*→∞ *d*(*x*, *xn*). The asymptotic radius of {*xn*} is given by *r*({*xn*}) := inf{*r*(*x*, {*xn*}) : *x* ∈ *X*} while the asymptotic center of {*xn*} is the set *A*({*xn*}) = {*x* ∈ *X* : *r*(*x*, {*xn*}) = *r*({*xn*})}. It is known that in an Hadamard space *X*, *A*({*xn*}) consists of exactly one point. A sequence {*xn*} in *X* is said to be Δ-convergent to a point *x* ∈ *X*, if *A*({*xnk*}) = {*x*} for every subsequence {*xnk*} of {*xn*}. In this case, we write <sup>Δ</sup>- lim*n*→∞*xn* <sup>=</sup> *<sup>x</sup>* (see [42,43]).

**Definition 2.** *Let D be a nonempty, closed and convex subset of an Hadamard space X*. *The metric projection is a mapping PD* : *X* → *D which assigns to each x* ∈ *X*, *the unique point PDx* ∈ *D such that*

$$d(\mathbf{x}, P\_D \mathbf{x}) = \inf \{ d(\mathbf{x}, y) : y \in D \}.$$

**Lemma 1.** *[44] Let D be a nonempty, closed convex subset of an Hadamard space X*, *x* ∈ *X and u* ∈ *D*. *Then u* = *PDx if and only if* −*xu* →, −*uy* → ≥ 0, *for all y* <sup>∈</sup> *<sup>D</sup>*.

**Lemma 2.** *[29] Let D be a nonempty convex subset of an Hadamard space X and T* : *D* → *X be an α-inverse strongly monotone mapping. Assume λ* ∈ [0, 1] *and define T<sup>λ</sup>* : *D* → *X by Tλx* = (1 − *λ*)*x* ⊕ *λTx*. *If* 0 < *λ* < 2*α*, *then T<sup>λ</sup> is nonexpansive and F*(*Tλ*) = *F*(*T*).

**Lemma 3.** *[29] Let D be a nonempty convex subset of an Hadamard space X and T* : *D* → *X be a mapping. Then*

$$VI(D,T) = VI(D,T\_\lambda)\_{\prime\prime}$$

*where λ* ∈ (0, 1] *and T<sup>λ</sup>* : *D* → *X is a mapping defined by Tλx* = (1 − *λ*)*x* ⊕ *λTx*, *for all x* ∈ *D*.

**Remark 1.** *Observe from Lemma 2 that*

$$F(P\_D T) = VI(D, T) = VI(D, T\_\lambda) = F(P\_D T\_\lambda).$$

**Lemma 4.** *[29] Let D be a nonempty bounded closed convex subset of an Hadamard space X and T* : *D* → *X be an α-inverse-strongly monotone mapping. Then V I*(*D*, *T*) *is nonempty, closed and convex.*

**Lemma 5.** *[41,45] Let X be an Hadamard space. Then for all x*, *y*, *z* ∈ *X and all t*,*s* ∈ [0, 1], *we have*


**Lemma 6.** *[46] Let X be a CAT(0) space and z* ∈ *X. Let x*1, ··· , *xN* ∈ *X and γ*1, ··· , *γ<sup>N</sup> be real numbers in* [0, 1] *such that <sup>N</sup>* ∑ *i*=1 *γ<sup>i</sup>* = 1*. Then the following inequality holds:*

$$d^2\left(\sum\_{i=1}^N \oplus \gamma\_i \mathbf{x}\_i, z\right) \le \sum\_{i=1}^N \gamma\_i d^2\left(\mathbf{x}\_i, z\right) - \sum\_{i,j=1, i \ne j}^N \gamma\_i \gamma\_j d^2\left(\mathbf{x}\_i, \mathbf{x}\_j\right).$$

**Lemma 7.** *[47] Every bounded sequence in an Hadamard space has a* Δ*-convergent subsequence.*

**Lemma 8.** *[48] Let X be an Hadamard space,* {*xn*} *be a sequence in X and x* ∈ *X*. *Then* {*xn*} Δ*-converges to x if and only if*

$$\limsup\_{n \to \infty} \langle \overrightarrow{\mathfrak{x}\_n} \mathfrak{X}, \overrightarrow{y} \widehat{\mathfrak{X}} \rangle \le 0, \quad \forall \ y \in X.$$

**Definition 3.** *Let D be a nonempty, closed and convex subset of an Hadamard space X*. *A mapping T* : *<sup>D</sup>* <sup>→</sup> *<sup>D</sup> is said to be* <sup>Δ</sup>*-demiclosed, if for any bounded sequence* {*xn*} *in <sup>X</sup>*, *such that* <sup>Δ</sup> <sup>−</sup> lim*n*→<sup>∞</sup> *xn* <sup>=</sup> *<sup>x</sup> and* lim*n*→<sup>∞</sup> *<sup>d</sup>*(*xn*, *Txn*) = 0, *then x* <sup>=</sup> *Tx.*

**Lemma 9.** *[49] Let X be an Hadamard space and T* : *X* → *X be a nonexpansive mapping. Then T is* Δ*-demiclosed.*

**Lemma 10.** *[50,51] Let* {*an*} *be a sequence of non-negative real numbers satisfying*

$$a\_{n+1} \le (1 - \alpha\_n)a\_n + \delta\_{n\prime} \quad n \ge 0\_{\prime\prime}$$

*where* {*αn*} *and* {*δn*} *satisfy the following conditions:*


**Lemma 11.** *[52] Let* {*an*} *be a sequence of non-negative real numbers such that there exists a subsequence* {*nj*} *of* {*n*} *with anj* < *anj*<sup>+</sup><sup>1</sup> *for all j* ∈ N*. Then there exists a nondecreasing sequence* {*mk*} ⊂ N *such that mk* → ∞ *and the following properties are satisfied by all (sufficiently large) numbers k* ∈ N*:*

*amk* ≤ *amk*<sup>+</sup><sup>1</sup> *and ak* ≤ *amk*+1.

*In fact, mk* = max{*i* ≤ *k* : *ai* < *ai*+1}*.*

#### **3. Main Results**

In this section, we present our strong convergence results.

**Theorem 1.** *Let X be an Hadamard space and D be a nonempty, closed and convex subset of X*. *Let Si* : *D* → *D be a finite family of Li-Lipschitz ki-demicontractive mappings and* Δ*-demiclosed such that L* > 0, *L* = max{*Li*, *i* = 1, 2, ··· , *N*}, *k* ∈ [0, 1), *k* = max{*ki*, *i* = 1, 2, ··· , *N*}, *ki* ∈ [0, 1), *i* = 1, 2, ··· , *N*. *Let Ti* : *D* → *X be a finite family of αi-inverse strongly monotone mappings and f be a contraction on D with coefficient <sup>θ</sup>* <sup>∈</sup> (0, 1). *Suppose that* <sup>Γ</sup> :<sup>=</sup> <sup>∩</sup>*<sup>N</sup> <sup>i</sup>*=1*V I*(*D*, *Ti*) ∩ ∩*<sup>N</sup> <sup>i</sup>*=1*F*(*Si*) = ∅. *For arbitrary x*<sup>1</sup> ∈ *D*, *let the sequence* {*xn*} *be generated by*

$$\begin{cases} w\_n = \gamma\_n f(\mathbf{x}\_n) \oplus (1 - \gamma\_n) \mathbf{x}\_{n\prime} \\ y\_n = \beta\_{n,0} w\_n \oplus \sum\_{i=1}^N \oplus \beta\_{n,i} P\_D T\_{\lambda\_i} w\_{n\prime} \\ \mathbf{x}\_{n+1} = \mathbf{a}\_{n,0} y\_n \oplus \sum\_{i=1}^N \oplus \mathbf{a}\_{n,i} \mathbf{S}\_i y\_{n\prime} \quad \forall \ n \ge 1 \,, \end{cases} \tag{10}$$

*where Tλ<sup>i</sup>* = (1 − *λi*)*x* ⊕ *λiTix*, 0 < *λ<sup>i</sup>* < 2*αi*, *for each i* = 1, 2, ··· , *N*, {*γn*}, {*βn*,*i*} *and* {*αn*,*i*} ∈ (0, 1) *such that the following conditions are satisfied:*

$$\begin{aligned} &(A1) \lim\_{\substack{n \to \infty \\ \infty \\ \end{n}} \gamma\_{\mathbb{N}} = 0; \\ &(A2) \sum\_{n=1}^{\infty} \gamma\_{\mathbb{N}} = \infty; \\ &(A3) \ 0 < a \le \beta\_{\mathbb{N}, i}, a\_{\mathbb{N}, i} \le b < 1; \sum\_{i=0}^{N} a\_{\mathbb{N}, i} = 1 \text{ and } \sum\_{i=0}^{N} \beta\_{\mathbb{N}, i} = 1; \\ &(A4) \ 0 < c \le a\_{\mathbb{N}, 0} - k \le d < 1. \end{aligned}$$

*Then, the sequence* {*xn*} *converges strongly to an element z*¯ = *P*<sup>Γ</sup> *f*(*z*¯), *where P*<sup>Γ</sup> *is the metric projection of X onto* Γ.

**Proof.** Let *p* ∈ Γ; then by (10), condition (A3), Lemma 6 and the fact that *PDT<sup>λ</sup>* is nonexpansive, we have

$$\begin{split} d^2(y\_n, p) &= d^2(\beta\_{n,0} w\_n \sum\_{i=1}^N \oplus \beta\_{n,i} P\_D T\_{\lambda\_i} w\_{n'} p) \\ &\leq \beta\_{n,0} d^2(w\_n, p) + \sum\_{i=1}^N \beta\_{n,i} d^2(P\_D T\_{\lambda\_i} w\_{n'} p) - \sum\_{i=1}^N \beta\_{n,0} \beta\_{n,i} d^2(P\_D T\_{\lambda\_i} w\_{n'} w\_n) \\ &\leq d^2(w\_n, p) - \sum\_{i=1}^N \beta\_{n,0} \beta\_{n,i} d^2(P\_D T\_{\lambda\_i} w\_{n'} w\_n) .\end{split} \tag{11}$$

Additionally, since *Si* is demicontractive, we have from (10) and Lemma 6 that

$$\begin{split} d^2(\mathbf{x}\_{n+1}, p) &\leq a\_{n,0} d^2(y\_n, p) + \sum\_{i=1}^N a\_{n,i} d^2(S\_i y\_n, p) - \sum\_{i=1}^N a\_{n,0} a\_{n,i} d^2(y\_n, S\_i y\_n) \\ &\leq a\_{n,0} d^2(y\_n, p) + \sum\_{i=1}^N a\_{n,i} \left[ d^2(y\_n, p) + k d^2(S\_i y\_n, y\_n) \right] - \sum\_{i=1}^N a\_{n,0} a\_{n,i} d^2(y\_n, S\_i y\_n) \\ &= d^2(y\_n, p) - \sum\_{i=1}^N a\_{n,i} (a\_{n,0} - k) d^2(y\_n, S\_i y\_n) .\end{split} \tag{12}$$

From (11), (12) and condition (A4), we have

$$\begin{split}d(\mathbf{x}\_{n+1\prime}p) &\leq \gamma\_{\rm n}d(f(\mathbf{x}\_{n}),p) + (1-\gamma\_{\rm n})d(\mathbf{x}\_{n\prime}p) \\ &\leq \gamma\_{\rm n}\theta d(\mathbf{x}\_{n\prime}p) + \gamma\_{\rm n}d(f(p),p) + (1-\gamma\_{\rm n})d(\mathbf{x}\_{n\prime}p) \\ &\leq (1-\gamma\_{\rm n}(1-\theta))d(\mathbf{x}\_{n\prime}p) + \gamma\_{\rm n}d(f(p),p). \end{split}$$
 
$$\leq \max\{d(\mathbf{x}\_{n\prime}p), \frac{1}{1-\theta}d(f(p),p)\}$$
 
$$\vdots$$
 
$$\leq \max\{d(\mathbf{x}\_{1\prime}p), \frac{1}{1-\theta}d(f(p),p)\}.$$

Thus, {*xn*} is bounded. Consequently, {*wn*}, {*yn*}, {*Syn*} and {*PDT<sup>λ</sup><sup>i</sup> wn*} are also bounded. Now we divide the rest of the proof into two cases:

*Axioms* **2020**, *9*, 143

**Case 1:** Assume that {*d*2(*xn*, *<sup>p</sup>*)} is a monotonically non-increasing sequence. Then, {*d*2(*xn*, *<sup>p</sup>*)} is convergent and

$$d^2(\mathfrak{x}\_{n\prime}p) - d^2(\mathfrak{x}\_{n+1}, p) \to 0, \text{ as } n \to \infty.$$

Hence, from (11) and (12), we have

$$\begin{split} d^2(x\_{n+1}, p) &\leq d^2(w\_n, p) - \sum\_{i=1}^N \beta\_{n,0} \beta\_{n,i} d^2(P\_{\mathcal{D}} T\_{\lambda\_i} w\_n, w\_n) - \sum\_{i=1}^N a\_{n,i} (a\_{n,0} - k) d^2(y\_n, S\_i y\_n) \\ &\leq \gamma\_n d^2(f(x\_n), p) + (1 - \gamma\_n) d^2(x\_n, p) - \sum\_{i=1}^N \beta\_{n,0} \beta\_{n,i} d^2(P\_{\mathcal{D}} T\_{\lambda\_i} w\_n, w\_n) - \sum\_{i=1}^N a\_{n,i} (a\_{n,0} - k) d^2(y\_n, S\_i y\_n), \end{split} \tag{13}$$

which implies

$$\sum\_{i=1}^{N} \beta\_{n,0} \beta\_{n,i} d^2(P\_D T\_{\lambda\_i} w\_{n\prime} w\_n) \le \gamma\_n (d^2(f(\mathbf{x}\_n), p) - d^2(\mathbf{x}\_{n\prime} p)) + d^2(\mathbf{x}\_{n\prime} p) - d^2(\mathbf{x}\_{n+1\prime} p).$$

By conditions (A1) and (A3), we obtain that

$$\lim\_{n \to \infty} d(P\_D T\_{\lambda\_i} w\_{n\prime} w\_n) = 0, \; i = 0, 1, 2, \cdots \; , N. \tag{14}$$

Similarly, from (13) and condition (A1), we have

$$\lim\_{n \to \infty} d(S\_i y\_n, y\_n) = 0, \; i = 0, 1, 2, \cdots \; \text{,} \; N. \tag{15}$$

Additionally,

$$d(w\_{\hbar}, \mathbf{x}\_{\hbar}) \le \gamma\_{\hbar} d(f(\mathbf{x}\_{\hbar}), \mathbf{x}\_{\hbar}) \to 0 \text{ as } \mathbf{n} \to \infty. \tag{16}$$

Again from (10) and Lemma 6, we have

$$\begin{aligned} d^2(y\_n, \mathbf{x}\_n) &\leq \beta\_{n,0} d^2(w\_n, \mathbf{x}\_n) + \sum\_{i=1}^N \beta\_{n,i} d^2(P\_D T\_{\lambda\_i} w\_n, \mathbf{x}\_n) \\ &\leq \beta\_{n,0} d^2(w\_n, \mathbf{x}\_n) + \sum\_{i=1}^N \beta\_{n,i} \left[ d(P\_D T\_{\lambda\_i} w\_n, w\_n) + d(w\_n, \mathbf{x}\_n) \right]^2. \end{aligned}$$

Hence, from (14) and (16), we obtain

$$d(y\_n, \mathbf{x}\_n) \to 0 \text{ as } n \to \infty. \tag{17}$$

Additionally, from (15) and (17), we have

$$d(S\_i y\_n, \mathbf{x}\_n) \le d(S\_i y\_n, y\_n) + d(y\_n, \mathbf{x}\_n) \to 0, \text{ as } n \to \infty, \text{ i } = 0, 1, 2, \dots, N. \tag{18}$$

Since *Si* is Lipschitz, then from (17) and (18), we have that

$$\begin{split}d(S\_i \mathbf{x}\_n, \mathbf{x}\_n) &\leq d(S\_i \mathbf{x}\_n, S\_i y\_n) + d(S\_i y\_n, \mathbf{x}\_n) \\ &\leq \text{Ld}(\mathbf{x}\_n, y\_n) + d(S\_i y\_n, \mathbf{x}\_n) \rightarrow 0, \text{ as } n \rightarrow \infty, \text{ i } i = 0, 1, 2, \cdots \text{ }, N. \end{split} \tag{19}$$

Hence, from (14), (16) and Lemma 2, we obtain that

$$\begin{split}d\left(P\_{\mathcal{D}}T\_{\lambda\_{i}}\mathbf{x}\_{n},\mathbf{x}\_{n}\right) &\leq d\left(P\_{\mathcal{D}}T\_{\lambda\_{i}}\mathbf{x}\_{n},P\_{\mathcal{D}}T\_{\lambda\_{i}}\mathbf{w}\_{n}\right) + d\left(P\_{\mathcal{D}}T\_{\lambda\_{i}}\mathbf{w}\_{n},\mathbf{w}\_{n}\right) + d\left(\mathbf{w}\_{n},\mathbf{x}\_{n}\right) \\ &\leq d\left(\mathbf{x}\_{n},\mathbf{w}\_{n}\right) + d\left(P\_{\mathcal{D}}T\_{\lambda\_{i}}\mathbf{w}\_{n},\mathbf{w}\_{n}\right) + d\left(\mathbf{w}\_{n},\mathbf{x}\_{n}\right) \rightarrow 0,\text{ as }n \rightarrow \infty,\text{ }i = 0,1,2,\cdots \; ,N.\end{split}\tag{20}$$

Since {*xn*} is bounded, by Lemma 7 there exists a subsequence {*xnk*} of {*xn*} such that Δ − lim *<sup>k</sup>*→<sup>∞</sup> *xnk* <sup>=</sup> *<sup>z</sup>* for some *<sup>z</sup>* <sup>∈</sup> *<sup>D</sup>*. Then, it follows from (16) that there exists a subsequence {*wnk*} of {*wn*}, such that Δ- lim *<sup>k</sup>*→<sup>∞</sup> *wnk* <sup>=</sup> *<sup>z</sup>*. Additionally, from (17), we have that <sup>Δ</sup>- lim *<sup>k</sup>*→<sup>∞</sup> *ynk* <sup>=</sup> *<sup>z</sup>*. Since *Si* is <sup>Δ</sup>-demiclosed for each *<sup>i</sup>* <sup>=</sup> 1, 2, ··· , *<sup>N</sup>*, it follows from (19) that *<sup>z</sup>* ∈ ∩*<sup>N</sup> <sup>i</sup>*=0*F*(*Si*). Additionally, *PDT<sup>λ</sup><sup>i</sup>* is nonexpansive (by Lemma 2) for each *i* = 1, 2, ··· , *N*, thus we obtain from (20) and Remark 1 that *<sup>z</sup>* ∈ ∩*<sup>N</sup> <sup>i</sup>*=0*F*(*PDT<sup>λ</sup><sup>i</sup>* ) = <sup>∩</sup>*<sup>N</sup> <sup>i</sup>*=0*F*(*PDTi*). Hence, *z* ∈ Γ.

Next we show that {*xn*} converges strongly to *z*¯ = *P*<sup>Γ</sup> *f*(*z*¯). Since {*xn*} is bounded, we may choose without loss of generality, a subsequence {*xnk*} of {*xn*} such that {*xnk*} Δ-converges to *z* and

$$\limsup\_{n \to \infty} \langle \overline{f(\overline{z})} \overline{\overline{z}}, \overline{x\_n} \overline{\overline{z}} \rangle = \lim\_{k \to \infty} \langle \overline{f(\overline{z})} \overline{\overline{z}}, \overline{x\_{n\_k}} \overline{\overline{z}} \rangle. \tag{21}$$

Thus, by (21) and Lemma 1, we obtain that

$$\limsup\_{k \to \infty} \langle \overline{f(\overline{z})} \overline{z}, \overline{x\_n} \overline{\overline{z}} \rangle = \langle \overline{f(\overline{z})} \overline{z}, \overline{z} \overline{\overline{z}} \rangle \le 0. \tag{22}$$

From (12), Lemma 5 (iii) and quasilinearization properties in Definition 1, we have that

*<sup>d</sup>*2(*xn*+1, *<sup>z</sup>*¯) <sup>≤</sup> *<sup>γ</sup>*<sup>2</sup> *nd*<sup>2</sup>(*f*(*xn*), *<sup>z</sup>*¯)+(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*)2*d*2(*xn*, *<sup>z</sup>*¯) + <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*) −−−→ *f*(*xn*)*z*¯, −→*xnz*¯ <sup>≤</sup> *<sup>γ</sup>*<sup>2</sup> *nd*<sup>2</sup>(*f*(*xn*), *<sup>z</sup>*¯)+(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*)2*d*2(*xn*, *<sup>z</sup>*¯) + <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*) −−−−−−→ *<sup>f</sup>*(*xn*)*f*(*z*¯), −→*xnz*¯ <sup>+</sup> −−−→*f*(*z*¯)*z*¯, −→*xnz*¯ <sup>≤</sup> *<sup>γ</sup>*<sup>2</sup> *nd*<sup>2</sup>(*f*(*xn*), *<sup>z</sup>*¯)+(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*)2*d*2(*xn*, *<sup>z</sup>*¯) + <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*) *<sup>θ</sup>d*2(*xn*, *<sup>z</sup>*¯) + −−−→*f*(*z*¯)*z*¯, −→*xnz*¯ <sup>≤</sup> (<sup>1</sup> <sup>−</sup> <sup>2</sup>*γ<sup>n</sup>* <sup>+</sup> <sup>2</sup>*γnθ*)*d*2(*xn*, *<sup>z</sup>*¯) + <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>γ</sup>n*) −−−→*f*(*z*¯)*z*¯, −→*xnz*¯ <sup>+</sup> *<sup>γ</sup>*<sup>2</sup> *n d*2(*f*(*xn*), *z*¯) + *d*2(*xn*, *z*¯) <sup>=</sup> <sup>1</sup> <sup>−</sup> <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>θ</sup>*)*d*2(*xn*, *<sup>z</sup>*¯) + <sup>2</sup>*γn*(<sup>1</sup> <sup>−</sup> *<sup>θ</sup>*) 1 − *γ<sup>n</sup>* <sup>1</sup> <sup>−</sup> *<sup>θ</sup>* −−−→*f*(*z*¯)*z*¯, −→*xnz*¯ <sup>+</sup> *<sup>γ</sup><sup>n</sup>* 2(1 − *θ*) *d*2(*f*(*xn*), *z*¯) + *d*2(*xn*, *z*¯) .

That is,

$$d^2(\mathbf{x}\_{n+1}, \boldsymbol{\bar{z}}) \le 1 - 2\gamma\_n (1 - \theta) d^2(\mathbf{x}\_n, \boldsymbol{\bar{z}}) + 2\gamma\_n (1 - \theta) M\_{n\nu} \tag{23}$$

.

where

$$M\_{\rm ll} = \left[\frac{1-\gamma\_n}{1-\theta} \langle \overrightarrow{f(\vec{z})}\overrightarrow{\vec{z}}, \overrightarrow{\vec{x}\_n}\overrightarrow{\vec{z}}\rangle + \frac{\gamma\_n}{2(1-\theta)} \langle d^2(f(\vec{x}\_n), \vec{z}) + d^2(\mathfrak{x}\_n, \vec{z})\rangle\right]^2$$

Thus from (22), (23) and condition (A1), we conclude by Lemma 10 that {*xn*} converges strongly to *z*¯ = *P*<sup>Γ</sup> *f*(*z*¯).

**Case 2:** Suppose there exists a subsequence {*nk*} of {*n*} such that *<sup>d</sup>*2(*xnk* , *<sup>p</sup>*) <sup>≤</sup> *<sup>d</sup>*2(*xk*+1, *<sup>p</sup>*) for all *<sup>k</sup>* <sup>∈</sup> <sup>N</sup>. Then by Lemma 11, there exists a nondecreasing sequence {*mk*} ⊂ N such that *mk* → ∞ :

$$d(\mathbf{x}\_{\mathfrak{m}\_k}, p) < d(\mathbf{x}\_{\mathfrak{m}\_{k+1}\prime}, p), \text{ and } d(\mathbf{x}\_k, p) < d(\mathbf{x}\_{k+1}, p) \,\,\forall \, k \in \mathbb{N}.\tag{24}$$

Therefore

$$\begin{split} 0 &\leq \lim\_{k \to \infty} \left( d(\mathbf{x}\_{m\_{k+1}}, p) - d(\mathbf{x}\_{m\_{k}}, p) \right) \\ &\leq \limsup\_{n \to \infty} \left( d(\mathbf{x}\_{n+1}, p) - d(\mathbf{x}\_{n}, p) \right) \\ &\leq \limsup\_{n \to \infty} \left( \gamma\_{n} d(f(\mathbf{x}\_{n}), p) + (1 - \gamma\_{n}) d(\mathbf{x}\_{n}, p) - d(\mathbf{x}\_{n}, p) \right) \\ &= \limsup\_{n \to \infty} \left( \gamma\_{n} \left( d(f(\mathbf{x}\_{n}), p) - d(\mathbf{x}\_{n}, p) \right) \right) = 0. \end{split}$$

This implies that

$$\lim\_{k \to \infty} \left( d(\mathbf{x}\_{m\_{k+1}\prime}, p) - d(\mathbf{x}\_{m\_{k\prime}\prime} p) \right) = 0. \tag{25}$$

Following the arguments as in Case 1, we get

$$\lim\_{k \to \infty} \langle \overrightarrow{f(\vec{z})} \overrightarrow{\vec{z}}, \overrightarrow{x\_{m\_k}} \overrightarrow{\vec{z}} \rangle \le 0. \tag{26}$$

Hence, from (23), we obtain that

$$d^2(\mathbf{x}\_{\mathfrak{m}\_{k+1}}, \bar{z}) \le 1 - 2\gamma\_{\mathfrak{m}\_k}(1-\theta)d^2(\mathbf{x}\_{\mathfrak{m}\_{k'}}, \bar{z}) + 2\gamma\_{\mathfrak{m}\_k}(1-\theta)M\_{\mathfrak{m}\_k}.$$

Additionally, from (24), we have that

$$d^2(\mathfrak{x}\_{m\_{k'}}z) \le \mathcal{M}\_{m\_{k'}}z$$

which implies that

$$\lim\_{k \to \infty} d^2(x\_{m\_k}, \bar{z}) = 0.$$

Thus, from cases 1 and 2, we conclude that {*xn*} converges to *z*¯ = *P*<sup>Γ</sup> *f*(*z*¯) which is an element of Γ.

We present some consequences of our main results.

Now, by setting *Si* to be a family of quasi-nonexpansive mappings in Theorem 1, we obtain the following result:

**Corollary 2.** *Let X be an Hadamard space and D be a nonempty, closed and convex subset of X*. *Let Si* : *D* → *D be a finite family of quasi-nonexpansive mappings, Ti* : *D* → *X be a finite family of αi-inverse strongly monotone mappings and <sup>f</sup> be a contraction on <sup>D</sup> with coefficient <sup>θ</sup>* <sup>∈</sup> (0, 1). *Suppose that* <sup>Γ</sup> :<sup>=</sup> <sup>∩</sup>*<sup>N</sup> <sup>i</sup>*=1*V I*(*D*, *Ti*) ∩ ∩*N <sup>i</sup>*=1*F*(*Si*) = ∅. *For arbitrary x*1, ∈ *D*, *let the sequence* {*xn*} *be generated by*

$$\begin{cases} w\_n = \gamma\_n f(\mathbf{x}\_n) \oplus (1 - \gamma\_n) \mathbf{x}\_{n\prime} \\ y\_n = \beta\_{n,0} w\_n \oplus \sum\_{i=1}^N \oplus \beta\_{n,i} P\_D T\_{\lambda\_i} w\_{n\prime} \\ \mathbf{x}\_{n+1} = \mathbf{a}\_{n,0} y\_n \oplus \sum\_{i=1}^N \oplus \mathbf{a}\_{n,i} \mathbf{S}\_i y\_{n\prime} \quad \forall \ n \ge 1 \,, \end{cases} \tag{27}$$

*where Tλ<sup>i</sup>* = (1 − *λi*)*x* ⊕ *λiTix*, 0 < *λ<sup>i</sup>* < 2*αi*, *for each i* = 1, 2, ··· , *N*, {*γn*}, {*βni* } *and* {*αni* } ∈ (0, 1) *such that conditions (A1)-(A3) of Theorem 1 are satisfied. Then, the sequence* {*xn*} *converges strongly to an element z*¯ = *P*<sup>Γ</sup> *f*(*z*¯), *where P*<sup>Γ</sup> *is the metric projection of X onto* Γ.

**Proof.** The proof follows from the proof of Theorem 1.

By setting *N* = 1 in Corollary 2, we obtain the following result:

**Corollary 3.** *Let X be an Hadamard space and D be a nonempty, closed and convex subset of X*. *Let S* : *D* → *D be a quasi-nonexpansive mapping, T* : *D* → *X be an α-inverse strongly monotone mapping and f be a* *contraction on D with coefficient θ* ∈ (0, 1). *Suppose that* Γ := *V I*(*D*, *T*) ∩ *F*(*S*) = ∅. *For arbitrary x*<sup>1</sup> ∈ *D*, *let the sequence* {*xn*} *be generated by*

$$\begin{cases} w\_n = \gamma\_n f(\mathbf{x}\_n) \oplus (1 - \gamma\_n) \mathbf{x}\_{n\prime} \\ y\_n = \beta\_{n,0} w\_n \oplus \beta\_{n,1} P\_D T\_\lambda w\_{n\prime} \\ x\_{n+1} = a\_{n,0} y\_n \oplus a\_{n,1} \mathbb{S} y\_{n\prime} \quad \forall \ n \ge 1, \end{cases} \tag{28}$$

*where T<sup>λ</sup>* = (1 − *λ*)*x* ⊕ *λTx*, 0 < *λ* < 2*α*, {*γn*} *and* {*αn*} ∈ (0, 1) *such that* 0 < *a* ≤ *βn*,*i*, *αn*,*<sup>i</sup>* ≤ *b* < 1, *for i* <sup>=</sup> 0, 1, <sup>1</sup> ∑ *i*=0 *<sup>α</sup>n*,*<sup>i</sup>* <sup>=</sup> 1, <sup>1</sup> ∑ *i*=0 *βn*,*<sup>i</sup>* = 1 *and conditions (A1)-(A2) of Theorem 1 are satisfied. Then, the sequence* {*xn*} *converges strongly to an element z*¯ = *P*<sup>Γ</sup> *f*(*z*¯), *where P*<sup>Γ</sup> *is the metric projection of X onto* Γ.

Suppose *u* = *f*(*x*) for arbitrary but fixed *u* ∈ *X* and for all *x* ∈ *X* in Theorem 1, we obtain the following result:

**Corollary 4.** *Let X be an Hadamard space and D be a nonempty, closed and convex subset of X*. *Let Si* : *D* → *D be a finite family of Li-Lipschitz demicontractive mappings and* Δ*-demiclosed such that L* > 0, *L* = max{*Li*, *i* = 1, 2, ··· , *N*}, *k* ∈ [0, 1), *k* = max{*ki*, *i* = 1, 2, ··· , *N*}, *ki* ∈ [0, 1), *i* = 1, 2, ··· , *N*. *Let Ti* : *<sup>D</sup>* <sup>→</sup> *<sup>X</sup> be a finite family of <sup>α</sup>i-inverse strongly monotone mappings and suppose that* <sup>Γ</sup> :<sup>=</sup> <sup>∩</sup>*<sup>N</sup> <sup>i</sup>*=1*V I*(*D*, *Ti*) ∩ ∩*N <sup>i</sup>*=1*F*(*Si*) = ∅. *For arbitrary x*1, *u* ∈ *D*, *let the sequence* {*xn*} *be generated by*

$$\begin{cases} w\_n = \gamma\_n u \oplus (1 - \gamma\_n) \mathbf{x}\_n \\ y\_n = \beta\_{n,0} w\_n \oplus \sum\_{i=1}^N \oplus \beta\_{n,i} P\_D T\_{\lambda\_i} w\_n \\ \mathbf{x}\_{n+1} = \mathbf{a}\_{n,0} y\_n \oplus \sum\_{i=1}^N \oplus \mathbf{a}\_{n,i} \mathbf{S}\_i y\_n \quad \forall \ n \ge 1 \end{cases} \tag{29}$$

*where Tλ<sup>i</sup>* = (1 − *λi*)*x* ⊕ *λiTix*, 0 < *λ<sup>i</sup>* < 2*αi*, *for each i* = 1, 2, ··· , *N*, {*γn*}, {*βni* } *and* {*αni* } ∈ (0, 1) *such that conditions (A1)-(A4) of Theorem 1 are satisfied. Then, the sequence* {*xn*} *converges strongly to an element z*¯ ∈ Γ *which is the nearest point to u*.

By setting *Si* ≡ *I* for all *i* = 1, 2, ··· , *N* in Theorem 1, we obtain the following result:

**Corollary 5.** *Let X be an Hadamard space and D be a nonempty, closed and convex subset of X*. *Let Ti* : *D* → *X be a finite family of αi-inverse strongly monotone mappings and f be a contraction on D with coefficient <sup>θ</sup>* <sup>∈</sup> (0, 1). *Suppose that* <sup>Γ</sup> :<sup>=</sup> <sup>∩</sup>*<sup>N</sup> <sup>i</sup>*=1*V I*(*D*, *Ti*) = ∅. *For arbitrary x*<sup>1</sup> ∈ *D*, *let the sequence* {*xn*} *be generated by*

$$\begin{cases} w\_n = \gamma\_n f(\mathbf{x}\_n) \oplus (1 - \gamma\_n) \mathbf{x}\_{n\prime} \\ y\_n = \beta\_{n,0} w\_n \oplus \sum\_{i=1}^N \oplus \beta\_{n,i} P\_D T\_{\lambda\_i} w\_n \\ x\_{n+1} = a\_{n,0} y\_n \oplus \sum\_{i=1}^N \oplus a\_{n,i} y\_n \quad \forall \ n \ge 1, \end{cases} \tag{30}$$

*where Tλ<sup>i</sup>* = (1 − *λi*)*x* ⊕ *λiTix*, 0 < *λ<sup>i</sup>* < 2*αi*, *for each i* = 1, 2, ··· , *N*, {*γn*}, {*βni* } *and* {*αni* } ∈ (0, 1) *such that conditions (A1)–(A3) of Theorem 1 are satisfied. Then, the sequence* {*xn*} *converges strongly to an element z*¯ = *P*<sup>Γ</sup> *f*(*z*¯), *where P*<sup>Γ</sup> *is the metric projection of X onto* Γ.

#### **4. Numerical Example**

In this section, we give a numerical experiment to show the applicability of our result.

**Example 1.** *[29] Let X* = R<sup>2</sup> *be an R-tree with radial metric dr*, *where dr*(*x*, *y*) = *d*(*x*, *y*) *if x and y are situated on a Euclidean straight line passing through the origin and dr*(*x*, *y*) = *d*(*x*, **0**) + *d*(*y*, **0**) *otherwise. We put p* = (0, 1), *q* = (1, 0) *and D* = *A* ∪ *B* ∪ *C*, *where*

$$A = \{(0, t) : t \in [2/3, 1] \}, \quad B = \{(t, 0) : t \in [2/3, 1] \}, \quad \mathbb{C} = \{(t, s) : t + s = 1, t \in (0, 1)\}.$$

*Define T* : *D* → *X by*

$$Tx = \begin{cases} q \text{ if } x \in A, \\ p \text{ if } x \in B, \\ x \text{ if } x \in C, \end{cases} \tag{31}$$

*then T is* <sup>1</sup> <sup>4</sup> *-inverse strongly monotone in* (*X*, *dr*).

*Now, define <sup>S</sup>* : *<sup>D</sup>* <sup>→</sup> *<sup>D</sup> by Sx* <sup>=</sup> <sup>5</sup> <sup>8</sup> *<sup>x</sup>*. *We make the following choices of parameters: <sup>λ</sup>* <sup>=</sup> <sup>1</sup> <sup>4</sup> , *<sup>γ</sup><sup>n</sup>* <sup>=</sup> <sup>1</sup> *<sup>n</sup>*+<sup>1</sup> , *<sup>α</sup>n*,0 <sup>=</sup> <sup>2</sup>*<sup>n</sup>* <sup>5</sup>*n*+<sup>7</sup> , *<sup>α</sup>n*,1 <sup>=</sup> <sup>3</sup>*n*+<sup>7</sup> <sup>5</sup>*n*+<sup>7</sup> , *<sup>β</sup>n*,0 <sup>=</sup> *<sup>n</sup>* <sup>3</sup>*n*+<sup>2</sup> , *<sup>β</sup>n*,1 <sup>=</sup> <sup>2</sup>*n*+<sup>2</sup> <sup>3</sup>*n*+<sup>2</sup> <sup>∀</sup> *<sup>n</sup>* <sup>≥</sup> <sup>1</sup> *and <sup>f</sup>*(*x*) = <sup>1</sup> <sup>2</sup> *x* ∀ *x* ∈ *X*; *then the conditions (A1)–(A3) of Theorem 2 are satisfied. Therefore, for x*<sup>1</sup> ∈ *X*, *Algorithm (28) becomes*

$$\begin{cases} w\_{\boldsymbol{n}} = \frac{1}{\boldsymbol{n}+1} f(\mathbf{x}\_{\boldsymbol{n}}) \oplus (1 - \frac{1}{\boldsymbol{n}+1}) \mathbf{x}\_{\boldsymbol{n}}, \\ y\_{\boldsymbol{n}} = \beta\_{\boldsymbol{n},0} w\_{\boldsymbol{n}} \oplus \beta\_{\boldsymbol{n},1} P\_{\boldsymbol{D}} T\_{\lambda} w\_{\boldsymbol{n}}, \\ x\_{\boldsymbol{n}+1} = a\_{\boldsymbol{n},0} y\_{\boldsymbol{n}} \oplus a\_{\boldsymbol{n},1} \mathbb{S} y\_{\boldsymbol{n}} \quad \forall \, \boldsymbol{n} \geq 1. \end{cases} \tag{32}$$

We now consider the following 3 cases for our numerical experiments given in Figure 1 above.

\*\*Case 1:  $\mathbf{x}\_1 = (-0.5, 0.5)^T$ .
\*\*Case 2:  $\mathbf{x}\_1 = (0.5, -0.5)^T$ .
\*\*Case 2:  $\mathbf{x}\_1 = (1, 2)^T$ .

**Figure 1.** *Cont.*

**Figure 1.** Errors vs. iteration numbers (n) for Example 1: case 1 (**a**); case 2 (**b**); case 3 (**c**).

**Author Contributions:** Conceptualization of the article was given by K.O.A., C.I. and O.T.M., methodology by H.A.A., K.O.A., and C.I., software by C.I., validation by O.T.M. and C.I., formal analysis, investigation, data curation, and writing–original draft preparation by K.O.A., H.A.A. and C.I. resources by K.O.A., C.I., H.A.A. and O.T.M. writing–review and editing by K.O.A., C.I. and O.T.M., visualization by C.I. and O.T.M., project administration by O.T.M., Funding acquisition by O.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** C.I. is funded by National Research Foundation (NRF) South Africa (S&F-DSI/NRF Free Standing Postdoctoral Fellowship; grant number: 120784); H.A.A. is funded by Department of Science and Technology and National Research Foundation, Republic of South Africa Center of Excellence in Mathematical and Statistical Sciences (DST-NRF COE-MaSS) and O.T.M. is funded by National Research Foundation (NRF) of South Africa Incentive Funding for Rated Researchers (grant number 119903).

**Acknowledgments:** The authors sincerely thank the anonymous reviewers for their careful reading, constructive comments and fruitful suggestions that substantially improved the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
