Convergence analysis of Suzuki’s generalized nonexpansive mappings using the Picard–Abbas iteration process
Bashir Nawaz, Krzysztof Gdawiec, Kifayat Ullah, Maha Noorwali, Maggie Aphane

TL;DR
This paper studies how a new iteration method works with certain types of mathematical mappings, showing it can effectively find fixed points and offering visual comparisons.
Contribution
The paper introduces convergence results for Suzuki’s nonexpansive mappings using the Picard–Abbas iteration and provides visual insights via polynomiographs.
Findings
Weak and strong convergence results are established for Suzuki’s generalized nonexpansive mappings.
A numerical example demonstrates the effectiveness of the Picard–Abbas iteration process.
Polynomiographs reveal visual advantages of the proposed method over existing ones.
Abstract
This manuscript investigates the convergence behavior of Suzuki’s generalized nonexpansive mappings using the recently introduced Picard–Abbas iteration process. We establish both weak and strong convergence results for the associated fixed-point approximations. To demonstrate the effectiveness of our approach, a numerical example is provided. Furthermore, we generate polynomiographs based on the proposed iteration process and compare them with those produced by existing methods, highlighting the advantages and visual insights offered by our scheme.
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Fig 1
Fig 2
Fig 3
Fig 4Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsOptimization and Variational Analysis · Fixed Point Theorems Analysis · Advanced Optimization Algorithms Research
1 Introduction
Fixed point theory is a versatile and powerful mathematical tool that plays a crucial role in various scientific and engineering disciplines. It is particularly effective for addressing complex nonlinear problems, where conventional analytical methods often prove inefficient or infeasible. The theory has broad applications, including population dynamics in biology [1], market equilibrium models in economics [2], stable strategy profiles in game theory [3], chemical equilibrium analysis in chemistry [4], stability analysis in engineering [4], and algorithm development in artificial intelligence [5]. By leveraging fixed point results, researchers can obtain optimal solutions while minimizing computational costs.
Given the complexity of these applications, standard analytical techniques are often either computationally expensive or incapable of providing exact solutions. Fixed point theory offers a powerful alternative by proving the existence of solutions and furnishing constructive methods to approximate them. A fundamental result in this field is Banach’s Contraction Principle (BCP) [6], which asserts that any contraction operator on a closed subset of a Banach space has a unique fixed point. Moreover, this fixed point can be effectively approximated using the Picard iteration method. This result forms a cornerstone for establishing the existence and approximation of solutions in a wide range of applied problems.
To formally define a contraction mapping, let be a nonempty subset of a Banach space . A self-mapping is said to be a contraction mapping if there exists a constant such that for all , the following inequality holds:
When , the mapping is said to be nonexpansive. Furthermore, a point is called a fixed point of if . Throughout this paper, will denote the set of all fixed points of . The mapping is said to be quasi-nonexpansive if:
for all , and .
Over time, various generalizations of contraction mappings have been proposed. One such extension is the class of nonexpansive mappings, introduced independently by Browder [7], Gohde [8], and Kirk [9]. To establish fixed point results for nonexpansive mappings, certain structural conditions such as closedness, boundedness, and uniform convexity are typically required [10]. Suzuki made a significant advancement in this direction [11], who proposed a generalization termed condition (C), characterizing a class of mappings now referred to as Suzuki’s generalized nonexpansive mappings. A mapping is said to satisfy condition (C) if, for all , the following holds:
Suzuki demonstrated that this class of mappings forms a broader category than quasi-nonexpansive mappings but is not as general as the class of nonexpansive mappings. Specifically, while every nonexpansive mapping satisfies condition (C), the converse does not necessarily hold. The following example illustrates this distinction.
Example 1.1 ([11]). Define a mapping by
In this example, satisfies condition (C) but is not a nonexpansive mapping.
Determining the fixed points of various classes of nonlinear mappings is a mathematically challenging task. This challenge is compounded by the failure of Picard iteration to converge for nonexpansive mappings in a complete metric space and by the inapplicability of the Banach Contraction Principle to such mappings. Consequently, numerous iterative procedures have been developed to approximate fixed points of these mappings. These methods have been extensively studied in the literature, notably in the works of Mann [12], Ishikawa [13], Noor [14], Abbas and Nazir [15], Sahu et al. [16], Thakur et al. [17], and Eke and Akewe [18], among others.
Let , , and be sequences in (0,1), where . The iteration scheme introduced by Noor [14] is recognized as the first three-step iteration process. This iteration process generates the sequence {un} as follows:
Abbas and Nazir proposed a faster iteration process than the Noor iteration, known as the Abbas iteration process [15], which generates the sequence {un} as follows:
Thakur et al. [17] introduced the following iteration process for approximating the fixed point of nonexpansive mappings:
Sahu et al. [16] proposed a new three-step iteration process to approximate fixed points of nonexpansive mappings, generating the sequence {un} as follows:
Eke and Akewe proposed a four-step iteration process, called the Picard–Noor iteration, which generates the sequence {un} as follows [18]:
A recent contribution by Manbhalang and Naveen [19] introduced the Picard–Abbas iteration process and established both weak and strong convergence results for contraction mappings. The Picard–Abbas iteration process is defined as follows:
In recent years, Suzuki’s generalized nonexpansive mappings have attracted considerable attention across various mathematical disciplines, leading to significant progress in fixed-point theory (see [20–23]). These mappings are particularly valuable for the development and analysis of iterative methods due to their rich structural properties and nuanced convergence behavior.
In this work, we investigate the convergence properties of Suzuki’s generalized nonexpansive mappings using the Picard–Abbas iteration process. Our study not only extends existing results but also offers a comparative perspective by analyzing the performance of several established iteration schemes, including those by Noor, Abbas, Thakur, Sahu, and the Picard–Noor iteration processes. To complement our theoretical findings, we present a new numerical example and employ polynomiography—a modern digital visualization technique—to depict the convergence behavior of the various iteration processes. This visual approach enhances the interpretability of the results and facilitates a deeper understanding of their dynamics.
The structure of this paper is as follows. Sect 2 introduces key definitions and fundamental lemmas. In Sect 3, we establish fixed-point results for the proposed iteration process. In Sect 4, we provide a numerical example to demonstrate the effectiveness of the scheme. Sect 5 illustrates the iteration process using visualizations generated through polynomiography. Sect 6 concludes the paper with final remarks.
2 Preliminaries
The following basic results are key to proving our main result.
Definition 2.1 ([10]). A Banach space is said to be a uniformly convex Banach space (UCBS) if, for all , there exists such that
Definition 2.2. Let be a convex and closed subset of a Banach space , and let {un} be a bounded sequence in . For any , the asymptotic radius of the sequence {un} with respect to is defined by
and the asymptotic center of {un} with respect to as
Definition 2.3 ([24]). A Banach space is said to have Opial’s property if, for every sequence {un} in that converges weakly to some (i.e., ), the following inequality holds:
for all with .
Proposition 2.4 ([11]). Let be a subset of a Banach space and let be a mapping:
(a) If is nonexpansive, then satisfies condition (C).**(b) Any mapping that satisfies condition (C) and has a fixed point is quasi-nonexpansive.**(c) If fulfills condition (C), then
Lemma 2.5 ([11]). Let be a subset of a Banach space equipped with Opial’s property. Let be a mapping satisfying condition (C). If {un} converges weakly to k and , then .
The concept of condition (I), originally introduced by Senter and Dotson [25], serves as an alternative approach for demonstrating the strong convergence of certain iterative processes in non-compact domains.
Definition 2.6. Let be a subset of a Banach space and let be a self-mapping defined on . The mapping is said to satisfy condition (I) if there exists a non-decreasing function with and g(u) > 0 for all u > 0, such that
where .
Lemma 2.7 ([11]). Let be a weakly compact convex subset of a UCBS , and let be a self-map on . Assume that satisfies condition (C), then has a fixed point.
Lemma 2.8([26]). Suppose is a UCBS and for all , we have . Suppose {un} and are two sequences in satisfying , and holds for some . Then, .
3 Main Results
This section presents convergence results for mappings that satisfy condition (C), utilizing the Picard–Abbas iteration process.
Lemma 3.1. Let be a closed and convex subset of a UCBS . Suppose that is a mapping satisfying condition (C) with . Let {un} be the sequence generated by the Picard–Abbas iteration process (10). Then, for any , the sequence {un} satisfies
Proof: Let and . Since satisfies condition (C), by Proposition 2.4(b), we have that is quasi non-expansive mapping, i.e.,
Using (10), we get
And,
Using (12), we have
Also
Using (12) and (13), we obtain
Similarly,
By using (14), we get
It follows from (12)–(15) that
Hence, the sequence is both bounded and non-increasing. Thus, we can conclude that exists for each . □
Next, we discuss the existence of a fixed point for mappings satisfying condition (C).
Theorem 3.2. Let be a UCBS, and let be a nonempty, closed and convex subset. Suppose that is a mapping satisfying condition (C), and let {un} be the sequence generated by the Picard–Abbas iteration process (10). Then, is nonempty if and only if the sequence {un} is bounded and .
Proof: Suppose that and let . By Lemma 3.1, we conclude that the sequence un is bounded and that the limit exists and it is finite. Define
From Lemma 3.1, we get
Thus,
As satisfies condition (C), and by Preposition 2.4(b), we get
Thus,
By Lemma 3.1, we get
Using (19) and (17), we have
From (17) and (20), we obtain that
From Lemma 3.1, one has
so
Using (21) and (22), we get
Using Lemma 2.8 with (16), (18) and (23), we have
Conversely, suppose that {un} is bounded and . Let , by Proposition 2.4(c), we obtain
which implies that . As is UCBS, therefore is a singleton, which means . Hence, . □
Next, we prove weak convergence using Opial’s property.
Theorem 3.3. Let be a closed and convex subset of a UCBS . Let be a mapping satisfying condition (C) with . Suppose that the space satisfies Opial’s condition. If {un} is the sequence generated by the Picard–Abbas iteration process (10), then {un} converges weakly to a fixed point of .
Proof: By Theorem 3.2, the sequence {un} is bounded and satisfies . Since is a uniformly convex Banach space, it is reflexive. Thus, there exists a subsequence of {un} that converges weakly to some . By Lemma 2.5, it follows that .
To show that {un} converges weakly to x1, assume, for the sake of contradiction, that it does not. Then there exists another subsequence of {un} such that converges weakly to with . Again, by Lemma 2.5, we have .
Now, by applying Opial’s condition together with Lemma 3.1, we obtain
This contradicts our supposition, so . Thus, {un} converges weakly to a point in . □
Now, we use the concept of compactness to prove strong convergence.
Theorem 3.4. Let be a mapping satisfying condition (C) defined on a nonempty, closed, and compact subset of a uniformly convex Banach space . Let {un} be a sequence generated by (10). Then, {un} converges strongly to a fixed point of .
Proof: By using Lemma 2.7, we obtain . Since , it follows from Lemma 3.2 that . As is given to be compact and closed, there exists a subsequence of {un} in such that it converges strongly to some , i.e., . Hence, using these facts together with Proposition 2.4(c), we obtain
Letting , we obtain . This implies which means that . Moreover, Lemma 3.1 implies that the limit exists. Hence, k is the strong limit of the sequence {un}. □
The following theorem proves strong convergence without requiring compactness.
Theorem 3.5. Let be a UCBS, and let be a nonempty, closed, and convex subset of . Suppose that is a mapping satisfying condition (C), and let {un} be a sequence generated by (10). Then, {un} converges to a point in if and only if
Proof: Suppose that the sequence {un} converges to some . Then, by the definition of convergence,
it follows that
Conversely, assume that . From Lemma 3.1, the limit exists, which gives
and this provides
Therefore constitutes a decreasing sequence that is bounded below by zero, so it may be obtained that exists. Since , so . We now show that {un} is a Cauchy sequence in .
Since , for any , there exists an integer such that for all ,
Especially,
Thus, we can choose some such that
For any , applying the triangle inequality, we obtain
Since both terms on the right-hand side are bounded by , it follows that
which implies that {un} is a Cauchy sequence in .
Since is a closed subset of the Banach space , then {un} is convereges in . Consider for any . Applying , one obtains
Thus, , hence . □
We now use condition (I) to prove the strong convergence of the Picard-Abbas iteration process. This condition imposes additional constraints that strengthen convergence, especially when generalized non-expansiveness alone is not enough. It ensures norm convergence by linking the iterative sequence to a fixed point and guaranteeing that the distance between successive iterates gradually decreases.
Theorem 3.6. Let be a closed and convex subset of UCBS . Suppose that is a mapping satisfying condition (I), and let {un} be a sequence generated by (10). Then, the sequence {un} converges strongly to a fixed point of .
Proof: From (25), one can get exist and by Theorem 3.2, we obtain
From condition (I) and (27), we have
Therefore, Since g is a nondecreasing function with , g(u) > 0 for each u > 0, therefore we have
Hence, all conditions of Theorem 3.5 are satisfied, therefore, {un} converges strongly to a fixed point of . □
4 Numerical example
This section introduces a novel numerical example to demonstrate the convergence properties of mappings satisfying condition (C), as analyzed through the Picard–Abbas iteration process.
Example 4.1. Define such that
First, we show that the given mapping is not a nonexpansive mapping. For and , we obtain
Now,
We can notice that
Therefore, the mapping given in (28) is not a nonexpansive mapping.
Now, we prove that given in (28) satisfies condition (C).
When , then , and
For , we must have . Thus,
So, and implies that
Now,
Thus,
which implies that satisfies condition (C).When then , and
For , we must have . Now, we have two cases:
- (i) x > p. In this case, we have . Thus,
Which implies that .Therefore, . So,
which means that implies , so satisfies condition (C).
- (ii) p>x. In this case, we have . Thus,
Which implies that . As , we get
Now, let and . As is previously discussed in (i), therefore, now, working for and , we have , . Thus,
We first suppose that and , then
Which implies . Now, let assume that and , then
which shows that
Thus, , which shows that satisfies condition (C).Hence, it is established that is a Suzuki generalized nonexpansive mapping.
To illustrate the faster convergence of the proposed Picard–Abbas iteration process (10), we compare it against the Noor, Abbas, Thakur, Sahu, and Picard–Noor iteration methods. The selected parameters are , , and , with the stopping criterion defined as and the initial point u0 = 0.1. The corresponding results are presented in Fig 1 and Table 1.
Convergence behavior of Picard–Abbas, Noor, Abbas, Thakur, Sahu and Picard–Noor iteration processes corresponding to Table 1.
Table 1: Iterates produced by various iteration processes for the mapping 𝒴 given in (28) and the starting point u0 = 0.1.
The findings indicate that, after the first iteration, the value obtained by the Picard–Abbas iteration process (0.98716480) is the closest to the fixed point 1 among all compared methods. As shown in Table 1, each iteration method converges at a different rate. The proposed Picard–Abbas iteration is the fastest, reaching the fixed point in 5 iterations. The Sahu and Picard–Noor iteration methods require 6 iterations to converge. The Thakur and Abbas processes exhibit similar convergence behavior. In contrast, the Noor iteration method shows the slowest convergence, taking 15 iterations to reach the fixed point.
5 Comparison via polynomiography
Mathematician and computer scientist Bahman Kalantari introduced polynomiography, a digital art form and visual analytic technique for exploring root-finding problems [27,28]. Although related concepts, such as basins of attraction, dynamical planes, and speed of convergence, had appeared earlier in the literature, Kalantari was the first to consolidate these ideas under a unified framework. He defined polynomiography as the art and science of visualizing the approximation of the zeros of complex polynomials through iterative functions, referring to the resulting images as polynomiographs. Various types of iteration processes have since been compared and analyzed using polynomiographic techniques (see [29–33]).
The general procedure for generating polynomiographs is outlined in Algorithm 1. Color assignment within this algorithm can follow various approaches; in this study, we adopt a method that integrates basins of attraction with convergence speed [34]. Each root of the polynomial is assigned a distinct non-black color, while points that do not converge are marked in black. For each initial point u0 in the region A, the iterative method I is applied for up to K iterations. If convergence occurs in fewer than K steps, we determine the root closest to the resulting point un and assign its corresponding color to u0. The brightness of the color reflects the speed of convergence: lighter shades indicate faster convergence, while darker shades represent slower convergence. If no convergence is achieved within K iterations, u0 is colored black. This scheme effectively visualizes both the destination of convergence (via color) and the convergence rate (via shading), providing intuitive insights into the behavior of the iterative process.
Algorithm 1. Creation of a polynomiograpgh.
One well-known root-finding algorithm is the Newton’s iteration method, also known as the Newton–Raphson method. Its definition is:
where is the starting point and is a polynomial with complex coefficients. We can write (29) in terms of a fixed point iteration process as follows:
where . Thus, this is the Picard iteration. If the iteration process (30) converges to any fixed point of , then one has
Thus, , which means that x is a root of . Finding the fixed points of is therefore equal to solving the problem of finding the roots of . This enables us to use various fixed point iteration processes for , such as the suggested Picard–Abbas iteration.
In the considered example, we use three sets of iterations’ parameters
, , ; , , ; , , .
For each of the three sets of iteration parameters, we generated polynomiographs of the polynomial , which has six roots: –1.0, –0.5–0.866025i, i, 0.5–0.866025i, i, and 1.0. The iteration schemes used include Picard–Abbas, Sahu, Abbas, Thakur, Noor, and Picard–Noor methods. The parameters for polynomiograph generation were: region A = [−2,2*]*^2^, maximum number of iterations K = 45, and . Additionally, for each polynomiograph, we computed the Average Number of Iterations (ANI) as proposed in [35].
The polynomiographs generated for the first set of parameter values are shown in Fig 2. Distinct convergence patterns are observed for the Picard–Abbas, Sahu, Abbas, Thakur, Noor, and Picard–Noor iteration processes. Visual inspection indicates that the proposed Picard–Abbas iteration exhibits the fastest convergence, followed by Abbas, Picard–Noor, Sahu, Thakur, and Noor. Notably, for the Noor iteration, no points within the considered region converged to any root, resulting in a completely black polynomiograph. ANI values in Table 2 corroborate these findings: Picard–Abbas (4.156), Abbas (6.110), Picard–Noor (8.392), Sahu (8.609), and Thakur (10.786).
Comparison of polynomiographs obtained from different iteration processes with parameters ρn = 0.06, ψn = 0.06, ηn = 0.06.
Table 2: ANI values of the polynomiographs given in Figs 2–4.
Comparison of polynomiographs obtained from different iteration processes with parameters ρn = 0.6, ψn = 0.6, ηn = 0.6.
Comparison of polynomiographs obtained from different iteration processes with parameters ρn = 0.9, ψn = 0.9, and ηn = 0.9.
The polynomiographs for the parameter settings , , and are shown in Fig 3. The results show that the Noor iteration exhibits the slowest convergence speed, with the highest ANI value of 12.068. Among the iterations studied, the Picard–Abbas method achieves the fastest convergence, yielding the lowest ANI value of 3.746. In terms of convergence speed, the Picard–Noor iteration ranks second with an ANI of 4.778, followed by the Sahu (5.216), Abbas (5.390), and Thakur (6.078) iterations.
The third configuration employs high values for the iteration parameters. Similar to the previous cases, the Noor iteration exhibits the slowest convergence, as shown in Fig 4. In contrast, the Picard–Abbas iteration once again achieves the fastest convergence. Interestingly, the high-parameter setting leads to faster convergence across all methods, requiring fewer iterations to reach the polynomial’s roots. The ANI values corresponding to this configuration are presented in Table 2. The Picard–Noor iteration yields the lowest ANI value of 3.720, followed closely by the Picard–Abbas iteration with an ANI of 3.817.
6 Conclusion
Our analysis of Suzuki mappings using the Picard–Abbas iteration process demonstrates its enhanced convergence performance. The numerical results in Table 1 confirm its efficiency relative to several established methods, including those of Noor, Sahu, Thakur, Abbas, and Picard–Noor. Furthermore, visualizations generated through polynomiography provide additional insight into the convergence behavior, highlighting the iteration’s faster convergence rate. Collectively, these findings suggest that the Picard–Abbas process is a robust and effective tool for solving fixed-point problems, with promising potential for broader applications in mathematical and computational contexts.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Amar AB, Jeribi A, Mnif M. Some fixed point theorems and application to biological model. Numerical Functional Analysis and Optimization. 2008;29(1–2):1–23. doi: 10.1080/01630560701749482 · doi ↗
- 2Danet RM, Popescu MV. Some applications of the fixed-point theory in economics. Creat Math Inform. 2008;17(3):392–8.
- 3Yuan A. Fixed point theorems and applications to game theory. University of Chicago; 2017.
- 4Constanda C, Ahues M, Largillier A. Integral methods in science and engineering: analytic and numerical techniques. Boston: Birkhauser; 2011.
- 5Bertsekas DP. Nonlinear programming. 3rd ed. Belmont: Athena Scientific; 2016.
- 6Caccioppoli R. Un teorema generale sull’esistenza di elementi uniti in una transformazione funzionale. Rendiconti dell’Academia Nazionale dei Lincei. 1930;11:794–9.
- 7Browder FE. Nonexpansive nonlinear operators in a Banach space. Proc Natl Acad Sci U S A. 1965;54(4):1041–4. doi: 10.1073/pnas.54.4.1041 16578619 PMC 219790 · doi ↗ · pubmed ↗
- 8Göhde D. Zum Prinzip der kontraktiven Abbildung. Mathematische Nachrichten. 1965;30(3–4):251–8. doi: 10.1002/mana.19650300312 · doi ↗
