RETRACTED: Computational Analysis of Variational Inequalities Using Mean Extra-Gradient Approach

Abstract

An improved variational inequality strategy for dealing with variational inequality in a Hilbert space is proposed in this article as an alternative; if Hilbert space is used as the domain of interest, the original extra-gradient method is proposed for resolving variational inequality. This improved variational inequality strategy can be used as a substitute for the original extra-gradient method in some situations. Mann’s mean value method, coupled with the widely used sub-gradient extra-gradient strategy, makes it possible to update all of the previous iterations in a single step, thus saving time and effort. All of this is made feasible via the use of Mann’s mean value technique in conjunction with the convex hull of all prior iterations of the algorithm. It is guaranteed that the mean value iteration will result in an acceptable resolution of a variational inequality issue as long as one or more of the criteria for the averaging matrix are fulfilled. Numerous experiments were performed in order to demonstrate the correctness of the theoretical conclusion obtained.

1. Introduction

Suppose that H is the Hilbert space structure with a product existing on the interior of the space structure $\langle .,.\rangle$ in addition to the already defined standard norm $\parallel .\parallel$. To begin, consider C to be a closed convex, non-empty subset of H, and $F:H\to H$ to be a monotone operator with a non-degenerate definition in the space of closed convex subsets of H.

$$\langle \eta -\zeta ,F\left(\eta \right)-F\left(\zeta \right)\rangle \ge 0,$$

The L-Lipschitz operator is formed by combining the $\eta ,\zeta \in H$ and L-Lipschitz operators.

$$\parallel F\left(\eta \right)-F\left(\zeta \right)\parallel \le L\parallel \eta -\zeta \parallel ,$$

for $\eta ,\zeta \in H$. The task was created to help those who wanted to apply the Stampacchia variational inequality to an additive measurement of an array of C objects, such as determining the location of the items, by utilizing the Stampacchia variational inequality [1] as a starting point.

$$\langle F\left({\eta }^{*}\right),z-{\eta }^{*}\rangle \ge 0 \mathrm{for} \mathrm{all} z\in C.$$

(1)

We will represent the solution set of the deliberated variational inequality given above as VIP(F, C). According to the assumption, in the case of the investigated variational inequality, there are unlimited solutions to VIP(F, C). Many iterative techniques have been developed to deal with it, making use of its properties to describe both mathematical and practical issues (see [2] for further discussions). We can use the ${\eta }_1\in H$ equation to derive the answer.

$${\eta }_{k+1}=J_c\left({\eta }_k-\xi F\left({\eta }_k\right)\right), k\in \mathbb{N}.$$

(2)

The metric’s projection onto C is indicated by the letter J_c if the step size is greater than zero. Assume F is η−strongly monotone, L-Lipschitz continuous, and $\tau \in \left(0.2 \eta /L^2\right)$ [3,4] to show that the arrangement produced by (2) meets the unique solution of the problem VIP(F, C).

Korpelevich developed the EM approach [3] in response to the requirement for strong monotonicity, which was needed to aid in the convergence of iterative techniques for F, which was only available in limited quantities at the time of its creation. It is defined as ${\eta }_1 \in H$.

$$\left\{\begin{array}{c}{\zeta }_k=J_c\left({\eta }_k-\xi F\left({\eta }_k\right)\right), \\ {\eta }_{k+1}=J_c\left({\eta }_k-\xi F\left({\zeta }_k\right)\right), k\in \mathbb{N}.\end{array}\right.$$

(3)

With the help of EM (3), it is possible to construct a sequence in a finite dimensional space that is controlled by the Lipschitz continuity and monotonicity of F, which can then be used to obtain the VIP(F, C) solution in a finite dimensional space using the EM (3) formula for the VIP(F, C) solution. A number of variants of Korpelevich’s EM have been examined as a consequence of this starting point, e.g., [5,6,7,8,9,10], as well as the sources cited within [5,6,7,8,9,10] and elsewhere. A few of the researchers who have made significant contributions to this area of study include Censor, Gibali, and Reich [11]. Each iteration of EM must be completed in order for the figure to be completed correctly. This is shown by the completion of two metric projections. This means that EM is a suitable method to use if the limited set C is simple enough that a closed-form equation for the metric projection PC onto C exists; otherwise, a hidden minimization sub-issue must be addressed in addition to the main problem, as previously stated. Censor, Gibali, and Reich were the ones who came up with the SEM (sub-gradient extra-gradient technique) to solve this problem, and they were successful in their endeavors. Instead of updating with two metric projections onto C, the SEM only updates with one metric projection onto C when updating the next iteration η_k+1, as opposed to when updating the prior iteration ζ_k. Due to the fact that the SEM was formed during the previous iteration, it includes a half-space containing C that was modified during the previous iteration, which explains why it happened in this case. In order to accomplish the objectives of this research technique, it is necessary to follow the formula below:

$$\left\{\begin{array}{l}{\zeta }_k=J_c\left({\eta }_k-\xi F\left({\eta }_k\right)\right),\\ {\eta }_{k+1}=J_{T_k}\left({\eta }_k-\xi F\left({\zeta }_k\right)\right), k\in \mathbb{N},\end{array}\right.$$

(4)

where

$$T_k=\left\{\omega \in H:\langle \left({\eta }_k-\xi F\left({\zeta }_k\right)\right)-{\zeta }_k,\omega -{\eta }_k\rangle \le 0\right\}.$$

Apart from that, Formula (7) is well documented in the literature and clearly demonstrates the exact formula in an understandable way. The weak convergence outcome is likewise agreed in [9]. For [12,13,14,15,16,17,18,19], other SEM methods, such as electron energy loss spectroscopy, have been explored. When utilizing closed convex simple sets, SEM restricted the performance of the metric projection to a subset of the set’s members, which was not the case when the closed convex simple set was not itself a simple set, as was the case in the absence of such an assumption. A less challenging approach is to predict the intersection of a smaller number of non-empty, convex closed sets first, followed by the intersection of a larger number of such closed sets [20,21,22,23,24,25,26].

Rather than concentrating only on this nonlinear issue, it may be more productive to take a different approach to the problem. The presence of the operator indicates nonlinearity, and vice versa. As a consequence, the equation for the issue is ${\eta }^{*}\in Fix T=\left\{\eta \in H:\eta =T_x\right\}\ne \varnothing$. Starting with the Picard iteration, we find the solution where ${\eta }_{k+1}$ equals the sum of η_k plus an offset, ${\eta }_{k+1}=T{\eta }_k, k\in \mathbb{N}$, and where the solution is defined by the sum of η_k+1 plus an offset and the total of η_k+1 plus an offset, ${\eta }_{k+1}=T{\eta }_k, k\in \mathbb{N}$.

Picard’s iterative technique does not converge, as shown in earlier study studies, indicating that this series of occurrences has no possibility of convergence, as well. T. Mann improved on Picard’s initial method in 1953, but in 2009, he enhanced the process even more by creating an even more complicated iteration. Because of the changes made by T. Mann to this edition, it is frequently referred to as T. Mann’s edition,

$${\eta }_{k+1}=T\overline{{\eta }_k}, k\in \mathbb{N}.$$

In informal conversations about this technique, the phrase “Mann’s mean value iteration” is often used to refer to this method as a whole. It is a widely used technique for resolving optimization difficulties because it helps to avoid numerically unfavorable circumstances such as zigzagging or spiraling behavior in a produced sequence around the solution set, which may occur when using other ways to handle optimization issues [27,28,29,30]. The Mann mean value iteration [24] is useful in a wide range of optimization situations. It is also one of the most extensively studied techniques accessible (see [31] for more information). There has been a great deal of study [32,33,34] that has used the recurrence of Mann’s mean value as a measure of dependability, and it has been shown to be successful. Using the monotone and the Lipschitz continuous operator, as well as concepts from the well-known SEM and Mann’s mean value iteration, this technique is presented here as an iterative approach. It is worth mentioning that some novel schemes given in [35,36] were developed that were used in power control and battery charge planning, resulting in dynamic uncertainties, perturbation of irradiation and temperature, and abrupt faults in output loads.

A new iterative method proposed in the article by using the idea of well-known SEM and Mann’s mean value iteration. A weakly convergent sequence is created at the beginning of the proposed technique, as shown in the illustration. The answer is finally found, and it is both written down in the text and graphically depicted in the picture VIP(F, C). When dealing with a constrained minimization issue, a finite family of non-empty closed convex simple sets is defined as one that is intersected by a constrained set. If certain circumstances are fulfilled, it is conceivable that the new approach will outperform the old one [37].

2. Important Concepts and Preliminaries

References [21,22] may be utilized to acquire more information. The following notations should be considered: When a series is converging, the sign ${\left\{{\eta }_k\right\}}_{k=1}^{\infty }$ indicates whether it is converging strongly or weakly; when a sequence is converging, the symbol ${\left\{{\eta }_k\right\}}_{k=1}^{\infty }$ indicates whether it is converging strongly or weakly; and when a series is converging, the symbol ${\left\{{\eta }_k\right\}}_{k=1}^{\infty }$ indicates whether it is converging strongly or weakly. We represent the strong and weak convergence of the sequence ${\left\{{\eta }_k\right\}}_{k=1}^{\infty }$ to $\eta \in H$ by ${\eta }_k\to \eta$ and ${\eta }_k\to \eta$ correspondingly. As the identifying operator, the letter “I” is utilized to differentiate H from the rest of the alphabet. To answer the question, a closed, convex, and non-empty subset of H and C must be investigated, and this subset must be closed, convex, and non-empty [38,39]. We can obtain an $\eta \in H$ point for any given $\eta \in H$ point in the coordinate system by reversing the direction of the $\eta \in H$ point. $J_c\left(\eta \right)$ is the point in C that is closest to the origin, and it is often referred to as $J_c\left(\eta \right)$.

$$\parallel \eta -J_c\left(\eta \right)\parallel =\underset{\eta \in C}{\inf }\parallel \eta -\zeta \parallel .$$

(5)

As H is projected onto the letter C in this case, it is rendered as J_c(η). It is important to remember that J_c is a non-expansive H to C transformation, and this should be taken into account. $J_c\left(\eta \right):H\to C$. It is tough to understand why this is the case when J_c is non-restrictive and non-expansive H to C mapping.

$$\parallel J_c\left(\eta \right)-J_c\left(\zeta \right)\parallel \le \parallel \eta -\zeta \parallel , \forall \eta ,\zeta \in H.$$

Furthermore, the predictions are based on measurements. J_c fulfils the attribute of variation:

$$\langle \eta -J_c\left(\eta \right),J_c\left(\eta \right)-\zeta \rangle \ge 0, \forall \eta \in H,\zeta \in C.$$

The hyperplane can be defined on the basis of the integer parameters $a\in H\left\{0\right\}$ and $\beta \in ℝ$.

$$H\le \left(\alpha ;\beta \right)=\left\{\eta \in H: a,\eta \le \beta \right\}.$$

The half-space and the hyperplane are both closed and convex sets, and their intersection is likewise a closed set. We can also use the following formula to project the metric onto the half-space $H\le \left(\alpha ;\beta \right)$:

$$J_{H\le \left(\alpha ;\beta \right)}\left(\eta \right)=\left\{\begin{array}{l}\eta -\frac{\langle a,\eta \rangle -\beta }{{\parallel a\parallel }^2}a, \mathrm{if} \langle a,\eta \rangle >\beta , \\ \eta \mathrm{if} \langle a,\eta \rangle \le \beta .\end{array}\right.$$

As illustrated below, we can claim that a point exists. T_η separates C from another point for any non-empty closed convex $C\subset H$, if the point T_η is located on the convex $\eta \ne C$ border (6). An intriguing aspect of the site is that it also provides the following services: When we examine the hyperplane $H\left(\eta -J_c;\langle J_c\left(\eta \right),\eta -J_c\left(\eta \right)\rangle \right)$, we can see that it has two distinct forms, and $\mathcal{H}$ is independent of the value of $\eta .$ It is determined that the first site η is in the first space, and the second site C is in the second space. We know that,

$$C\subset H\le \left(\eta -J_c\left(\eta \right);\langle J_c\left(\eta \right),\eta -J_c\left(\eta \right)\rangle \right).$$

In addition to the hyperplane, $H\le \left(\eta -J_C\left(\eta \right);\langle J_c\left(\eta \right),\eta -J_c\left(\eta \right)\rangle \right)$. If the primary hyperplane fails, another option is to seek the help of a secondary hyperplane to finish the job C at J_c(η).

Let $A:H\to 2^H$. It is capable of executing an operation on a set of values using a set-valued operator, according to the graph.

$$Gr\left(A\right):\left\{\left(\eta ,u\right)\in H\times H:u\in A\eta \right\}.$$

To sum up, all of A’s unmarked papers are marked as A.

$$A^{-1}\left(0\right):\left\{\eta \in H:0\in A\left(\eta \right)\right\}.$$

A monotone operator is defined as follows: Based on the concept of monotonicity, if A is a monotone operator, then B must likewise be a monotone operator.

$$\langle \eta -\zeta ,u-v\rangle \ge 0,$$

$$\forall \left(\eta ,u\right)\left(\zeta ,v\right)\in Gr\left(A\right).$$

Despite the fact that the monotone operator’s graph includes no links to any other monotone operators, it is considered the most monotonous operator [39] that can be found. Furthermore, since A has the greatest degree of monotonicity (even when convex and closed), all of its subsets (including convex and closed) are zeros.

It is conceivable that the set $C\subset H$; furthermore, depending on the circumstances, it can have a concave or convex shape. N_c(η) is the typical daily cone of the same size and form, as seen at $\eta \in C$.

$$N_C\left(\eta \right):\left\{\zeta \in \mathcal{H}: \langle \zeta ,z-\eta \rangle , \forall z\in C\right\}.$$

Allow $F:H\to H$. Assume that $H$ and C are both monotone continuous operators, and that H and C are both sets of the same type. Then, C is a monotone continuous operation that is a closed convex subset of $H$ and is not empty. We can then find out who the operator is. $A:H\to 2^H$ by

$$A\left(\eta \right):\left\{\begin{array}{c}F\left(\eta \right)+N_C\left(\eta \right), \mathrm{for} x\in C,\\ \varnothing , \mathrm{for} x\ne C.\end{array}\right.$$

At that time, A is a maximally monotone operator, and the subsequent significant property is satisfied:

$$\mathrm{VIP}\left(F,C\right)=A^{-1}\left(0\right).$$

3. Methodology of Proposed Scheme

This section is formulated to present an efficient approach, i.e., a mean extra-gradient approach to investigate the solutions of the problems related to the variational inequalities. Before detailing the methodology of the extra-gradient method, we present some preliminaries.

An infinite lower-triangular-row matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ is supposed to be an averaging matrix if the subsequent situations are fulfilled:

 A1.

$a_{l,m}\ge 0, \forall l,m\ge 1;$

 A2.

If $l<m$, then $a_{l,m}=0; \forall l\ge 1;$

 A3.

$a_{l,1}+a_{l,2}+\dots +a_{l,l}=1, \forall l\ge 1;$

 A4.

$\underset{l\to +\infty }{\lim }{a}_{l,m}=0, \forall m\ge 1.$

Considering an averaging matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ and a sequence ${\left\{{\eta }_l\right\}}_{l=1}^{\infty }$ from a real Hilbert space $H$, we represent the mean iterate as;

$${\overline{\eta }}_l=a_{l,1}{\eta }_1+a_{l,2}{\eta }_2+a_{l,3}{\eta }_3+\dots +a_{l,l}{\eta }_l, \forall l\ge 1.$$

The solution procedure of variation inequality by means of Mann’s type mean extra-gradient scheme is given as Algorithm 1.

Algorithm 1: Solution procedure by Mann’s type mean extra-gradient scheme.

1. INITIALIZATION:$\mathrm{Choose} \mathrm{a} \mathrm{point} {\eta }_1 \mathrm{belonging} \mathrm{to} \mathrm{Hilbert} \mathrm{space} H$, a positive 2. $\mathrm{parameter} \xi$$, \mathrm{and} {\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ averaging matrix. 3. STEP 1.$\mathrm{Assumed} \mathrm{a} \mathrm{present} \mathrm{iterate} {\eta }_l\in H$, calculate the mean iterate as; 4. ${\overline{\eta }}_l=a_{l,1}{\eta }_1+a_{l,2}{\eta }_2+a_{l,3}{\eta }_3+\dots +a_{l,l}{\eta }_l,$ 5. also calculate 6. ${\zeta }_l=P_C\left({\overline{\eta }}_l-\xi F\left({\overline{\eta }}_l\right)\right).$ 7. STEP 2.$\mathrm{If} {\zeta }_l={\overline{\eta }}_l,$$\mathrm{then} {\overline{\eta }}_l$$\mathrm{belongs} \mathrm{to} \mathrm{VIP}\left(F,C\right)$ and break the procedure. 8. $\mathrm{Otherwise}, \mathrm{build} \mathrm{half} \mathrm{space} T_l$, which is given by 9. $T_l=\left\{\varsigma \in H:\langle \left({\overline{\eta }}_l-\xi F\left({\overline{\eta }}_l\right)\right)-{\zeta }_l,\varsigma -{\zeta }_l\rangle \le 0\right\},$ 10. and compute the subsequent iterate as; 11. ${\eta }_{l+1}=P_{T_l}\left({\overline{\eta }}_l-\xi F\left({\zeta }_l\right)\right).$ 12. $\mathrm{Update} \mathrm{the} \mathrm{dummy} \mathrm{variable} l$$\mathrm{as} l=l+1$, and perform STEP 1.

Remark 1.

It is important to mention that when ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$, is the identity matrix, and then the above Mann’s type mean extra-gradient scheme becomes the classical sub-gradient extra-gradient scheme given in Ref. [11].

Now, we explain the stopping principles of the proposed scheme in STEP 2.

Proposition 1.

Suppose that the sequences ${\left\{{\overline{\eta }}_l\right\}}_{l=1}^{\infty }$ and ${\left\{{\zeta }_l\right\}}_{l=1}^{\infty }$ are generated by means of the suggested Mann’s type mean extra-gradient scheme. If there exist a constant $l_0\in \mathbb{N}$ so that ${\overline{\eta }}_{l_0}={\zeta }_{l_0}$, then show that ${\overline{\eta }}_{l_0}\in VIP\left(F,C\right).$

 Proof.

Suppose a constant $l_0\in \mathbb{N}$ so that ${\overline{\eta }}_{l_0}={\zeta }_{l_0}$, then by means of the definition ζ_l, we obtain

$${\overline{\eta }}_{l_0}={\zeta }_{l_0}=P_C\left({\overline{\eta }}_{l_0}-\xi F\left({\overline{\eta }}_{l_0}\right)\right),$$

which produces ${\overline{\eta }}_{l_0}\in C.$ For all $z\in C$, we obtain from the following inequality

$$\forall \zeta \in C, \eta \in H;\langle \eta -P_C\left(\eta \right), P_C\left(\eta \right)-\zeta \rangle \ge 0,$$

That

$$\langle z-{\overline{\eta }}_{l_0}, {\overline{\eta }}_{l_0}-\xi F\left({\overline{\eta }}_{l_0}\right)-{\overline{\eta }}_{l_0}\rangle \le 0.$$

This implies that

$$\langle z-{\overline{\eta }}_{l_0}, F\left({\overline{\eta }}_{l_0}\right)\rangle \ge 0,$$

which satisfy that ξ > 0 and this implies that ${\overline{\eta }}_{l_0}\in \mathrm{VIP}\left(F,C\right).$

By the above proposition, for the remaining convergence analysis, we can consider all over this segment that the proposed scheme does not dismiss after some finite number of repetitions; explicitly, we consider that $\forall l\ge 1; {\zeta }_l\ne {\overline{\eta }}_l$. □

Lemma 1.

Suppose the sequence ${\left\{{\overline{\eta }}_l\right\}}_{l=1}^{\infty }$ is obtained by means of Mann’s type mean extra-gradient scheme; then $u\in VIP\left(F,C\right)$ and $\forall l\ge 1$, and the following relation must hold.

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel {\overline{\eta }}_l-{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2, \\ \le {\displaystyle \sum _{m=1}^l}{a}_{l,m}\parallel {\eta }_m-{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2. \end{gathered}$$

 Proof.

Suppose $u\in \mathrm{VIP}\left(F,C\right)$ and l ≥ 1 be fixed. We know that the operator F is monotone, therefore

$$\langle F\left({\zeta }_l\right)-F\left(u\right),{\zeta }_l-u\rangle \ge 0.$$

This implies the following relation

$$0\le \langle F\left(u\right),{\zeta }_l-u\rangle \le \langle F\left({\zeta }_l\right),{\zeta }_l-u\rangle .$$

In the above, the second inequality is true because of $u\in \mathrm{VIP}\left(F,C\right)$ and ${\zeta }_l\in C$. Therefore, we also obtain

$$\langle F\left({\zeta }_l\right),{\eta }_{l+1}-u\rangle \ge \langle F\left({\zeta }_l\right),{\eta }_{l+1}-{\zeta }_l\rangle .$$

By means of the definition of T_l, we obtain the following relation

$$\langle {\eta }_{l+1}-{\zeta }_l,{\overline{\eta }}_l-\xi F\left({\overline{\eta }}_l\right)-{\zeta }_l\rangle \le 0$$

Now, it follows that

$$\begin{gathered} \langle {\eta }_{l+1}-{\zeta }_l,{\overline{\eta }}_l-\xi F\left({\zeta }_l\right)-{\zeta }_l\rangle =\langle {\eta }_{l+1}-{\zeta }_l,{\overline{\eta }}_l-\xi F\left({\overline{\eta }}_l\right)-{\zeta }_l\rangle \\ +\langle {\eta }_{l+1}-{\zeta }_l,{\overline{\eta }}_l-\xi F\left({\zeta }_l\right)+\xi F\left({\overline{\eta }}_l\right)\rangle \le \xi \langle {\eta }_{l+1}-{\zeta }_l,F\left({\overline{\eta }}_l\right)-F\left({\zeta }_l\right)\rangle . \end{gathered}$$

Introducing a parameter z_l as $z_l={\overline{\eta }}_l-\xi F\left({\zeta }_l\right),$ then

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2=\parallel P_{T_l}\left(z_l\right)-{u\parallel }^2=\parallel P_{T_l}\left(z_l\right)-z_l+z_l-{u\parallel }^2 \\ =\parallel P_{T_l}\left(z_l\right)-z_l{\parallel }^2+\parallel z_l-{u\parallel }^2+2\langle P_{T_l}\left(z_l\right)-z_l,z_l-u\rangle . \end{gathered}$$

By means of the property of P_Tl, we have

$$0\ge 2\langle z_l-P_{T_l}\left(z_l\right),u-u-P_{T_l}\left(z_l\right)\rangle =2\parallel z_l-P_{T_l}{\left(z_l\right)\parallel }^2+2\langle P_{T_l}\left(z_l\right)-z_l,z_l-u\rangle ,$$

this implies the following relation

$$\parallel z_l-P_{T_l}{\left(z_l\right)\parallel }^2+2\langle P_{T_l}\left(z_l\right)-z_l,z_l-u\rangle \le -\parallel z_l-P_{T_l}{\left(z_l\right)\parallel }^2.$$

By means of the above relation in $\parallel {\eta }_{l+1}-{u\parallel }^2$ to have

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel z_l-{u\parallel }^2-\parallel z_l-P_{T_l}{\left(z_l\right)\parallel }^2, \\ =\parallel {\overline{\eta }}_l-\xi F\left({\zeta }_l\right)-{u\parallel }^2-\parallel {\overline{\eta }}_l-\xi F\left({\zeta }_l\right)-{\eta }_{l+1}{\parallel }^2, \\ =\parallel {\overline{\eta }}_l-{u\parallel }^2+{\xi }^2\parallel F{\left({\zeta }_l\right)\parallel }^2-2\xi \langle F\left({\zeta }_l\right),{\overline{\eta }}_l-u\rangle -\parallel {\overline{\eta }}_l-{\eta }_{l+1}{\parallel }^2-{\xi }^2\parallel F{\left({\zeta }_l\right)\parallel }^2 \\ +2\xi \langle F\left({\zeta }_l\right),{\overline{\eta }}_l-{\eta }_{l+1}\rangle =\parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\eta }_{l+1}{\parallel }^2+2\xi \langle F\left({\zeta }_l\right),u-{\eta }_{l+1}\rangle . \end{gathered}$$

It can also be rewritten, by means of the above relations, as

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\eta }_{l+1}{\parallel }^2+2\xi \langle F\left({\zeta }_l\right),{\zeta }_l-{\eta }_{l+1}\rangle , \\ \le \parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l+{\zeta }_l-{\eta }_{l+1}{\parallel }^2+2\xi \langle F\left({\zeta }_l\right),{\zeta }_l-{\eta }_{l+1}\rangle , \\ =\parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2-2\langle {\overline{\eta }}_l-{\zeta }_l,{\zeta }_l-{\eta }_{l+1}\rangle +2\xi \langle F\left({\zeta }_l\right),{\zeta }_l-{\eta }_{l+1}\rangle , \\ =\parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2-2\langle {\zeta }_l-{\eta }_{l+1},{\overline{\eta }}_l-{\zeta }_l-\xi F\left({\zeta }_l\right)\rangle , \\ \le \parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2+2\xi \langle {\zeta }_l-{\eta }_{l+1},F\left({\overline{\eta }}_l\right)-F\left({\zeta }_l\right)\rangle , \\ \le \parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2+2\xi \parallel {\zeta }_l-{\eta }_{l+1}\parallel \parallel F\left({\overline{\eta }}_l\right)-F\left({\zeta }_l\right)\parallel . \end{gathered}$$

By means of the L-Lipschitz continuity and using the relation 2xy ≤ x² + y², we obtain the following form

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel {\overline{\eta }}_l-{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2+2\xi L\parallel {\zeta }_l-{\eta }_{l+1}\parallel \parallel {\overline{\eta }}_l-{\zeta }_l\parallel , \\ \le \parallel {\overline{\eta }}_l-{u\parallel }^2-\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2-\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2+{\xi }^2{L}^2\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2+\parallel {\zeta }_l-{\eta }_{l+1}{\parallel }^2, \\ =\parallel {\overline{\eta }}_l-{u\parallel }^2-\left(1-{\tau }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2. \end{gathered}$$

Lastly, by means of the convexity of the norm $\parallel ·{\parallel }^2$ and an averaging matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$, we obtain the following form

$$\begin{gathered} \parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel {\overline{\eta }}_l-{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2, \\ =\parallel {\displaystyle \sum _m^l}{a}_{l,m}{\eta }_m-{\displaystyle \sum _m^l}{a}_{l,m}{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2, \\ =\parallel {\displaystyle \sum _m^l}{a}_{l,m}{\left({\eta }_m-u\right)\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2, \\ \le {\displaystyle \sum _m^l}{a}_{l,m}\parallel {\left({\eta }_m-u\right)\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2. \end{gathered}$$

Now, we discuss a concept which we later use in the convergence analysis of the scheme. □

 Proposition 2

([39]). Consider a real sequence, ${\left\{{\omega }_l\right\}}_{l=1}^{\infty }$, the averaging matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$, and $r\in ℝ.$ If ${\omega }_l\to r$, then ${\overline{\omega }}_l=\sum _{m=1}^l{a}_{l,m}{\omega }_m\to r$.

The averaging matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ is known as M-concentrating, if for all real sequences ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ and ${\left\{{ϵ}_l\right\}}_{l=1}^{\infty }$, so that $\sum _{l=1}^{\infty }{ϵ}_l<+\infty$, and it is satisfied that

$${\omega }_{l+1}\le {\overline{\omega }}_l+{ϵ}_l.$$

In the above, ${\overline{\omega }}_l=\sum _{m=1}^l{a}_{l,m}{\omega }_m, \forall l\ge 1,$ we obtained $\underset{l\to \infty }{\lim }{\omega }_l$.

By means of Lemma 1, if we include an extra previous criterion on ξ, the term on the right-hand side, which is $\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2$, is non-positive. Along with this condition, the ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ averaging matrix is known as M-concentrating.

Theorem 1.

Consider that the matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ is M-concentrating and $\xi \in \left(0,1/L\right)$. Then, ${\left\{{\overline{\eta }}_l\right\}}_{l=1}^{\infty }$ is any sequence produced by means of Mann’s type mean extra-gradient approach and weakly converges to the solution of the problem $VIP\left(F, C\right)$.

 Proof.

Consider an element $u\in \mathrm{VIP}\left(F,C\right)$ and l ≥ 1; then, by means of Lemma 1

$$\parallel {\eta }_{l+1}-{u\parallel }^2\le {\displaystyle \sum _{m=1}^l}{a}_{l,m}\parallel {\eta }_m-{u\parallel }^2-\left(1-{\xi }^2{L}^2\right)\parallel {\overline{\eta }}_l-{\zeta }_l{\parallel }^2.$$

(6)

As we know that $\xi \in \left(0,1/L\right)$, we obtained

$$0<1-{\xi }^2{L}^2<1.$$

The relation ($6$) takes the following form:

$$\parallel {\eta }_{l+1}-{u\parallel }^2\le {\displaystyle \sum _{m=1}^l}{a}_{l,m}\parallel {\eta }_m-{u\parallel }^2.$$

Bearing in mind that $\omega =\parallel {\eta }_l-{u\parallel }^2$ and for all $k\ge 1, {ϵ}_l=0$, and by means of the supposition that the averaging matrix is $M$-concentrating, we determine that the limit $\underset{l\to \infty }{\lim }\parallel {\eta }_l-{u\parallel }^2$ exists and declare $e\left(u\right)\in ℝ$. By means of lemma, we obtain that $\underset{l\to \infty }{\lim }\sum _{m=1}^l{a}_{l,m}\parallel {\eta }_m-{u\parallel }^2$ exists having the limit e(u), and afterwards, it follows from these composed with ($6$) and $0<1-{\xi }^2{L}^2<1$ that

$$\underset{l\to \infty }{\lim }\parallel {\overline{\eta }}_l-{\zeta }_l\parallel =0.$$

(7)

In addition, we observe from Lemma 1

$$\parallel {\eta }_{l+1}-{u\parallel }^2\le \parallel {\overline{\eta }}_l-{u\parallel }^2\le {\displaystyle \sum _{m=1}^l}{a}_{l,m}\parallel {\eta }_m-{u\parallel }^2.$$

We also have the limit $\underset{l\to \infty }{\lim }\parallel {\overline{\eta }}_l-{u\parallel }^2=e\left(u\right).$ As the sequence ${\left\{{\overline{\eta }}_l\right\}}_{l=1}^{\infty }$ is a bounded sequence, there is a weak cluster point, η’, from the Hilbert space H and there is a subset ${\left\{{\overline{\eta }}_{l_i}\right\}}_{i=1}^{\infty }$ so that ${\overline{\eta }}_{l_i}\to {\eta }^{\prime }.$ Therefore, from the relation (7) ${\zeta }_{l_i}\to {\eta }^{\prime }$. We then assume another operator A, which is defined as $A:H\to 2^H$, read as

$$Q\left(\nu \right)=\left\{\begin{array}{c}F\left(\nu \right)+N_C\left(\nu \right), \mathrm{when} \nu \in C,\\ \oslash , \mathrm{otherwise}.\end{array}\right.$$

Now, Q is the operator that is maximally monotone besides $\mathrm{VIP}\left(F,C\right)=Q^{-1}\left(0\right).$ Additionally, as (v, w) belongs to G(Q), this means $w\in Q=F\left(\nu \right)+N_C\left(\nu \right)$, and we obtain $w-F\left(\nu \right)\in N_C\left(\nu \right)$; that is:

$$\langle w-F\left(\nu \right),\nu -\zeta \rangle \ge 0,\hspace{1em}\forall \zeta \in C.$$

(8)

Therefore, by means of the property of ${\zeta }_l$, we obtain

$$\langle {\overline{\eta }}_l-\xi F\left({\overline{\eta }}_l\right)-{\zeta }_l,{\zeta }_l-\nu \rangle \ge 0.$$

This implies that

$$\langle \frac{{\zeta }_l-{\overline{\eta }}_l}{\xi }+F\left({\overline{\eta }}_l\right),\nu -{\zeta }_l\rangle \ge 0,\hspace{1em}\forall l\ge 1.$$

(9)

Hence, by means of the relations ($8$) and ($9$), substituting ζ with ζ_li and ζ_l with ζ_li, respectively, we have

$$\begin{gathered} \langle w,\nu -{\zeta }_{l_i}\rangle \ge \langle F\left(\nu \right),\nu -{\zeta }_{l_i}\rangle \ge \langle F\left(\nu \right),\nu -{\zeta }_{l_i}\rangle -\langle \frac{{\zeta }_{l_i}-{\overline{\eta }}_{l_i}}{\xi }+F\left({\overline{\eta }}_{l_i}\right),\nu -{\zeta }_{l_i}\rangle , \\ =\langle F\left(\nu \right)-F\left({\zeta }_{l_i}\right),\nu -{\zeta }_{l_i}\rangle -\langle F\left({\zeta }_{l_i}\right)+F\left({\overline{\eta }}_{l_i}\right),\nu -{\zeta }_{l_i}\rangle -\langle \frac{{\zeta }_{l_i}-{\overline{\eta }}_{l_i}}{\xi },\nu -{\zeta }_{l_i}\rangle . \end{gathered}$$

Now, taking the limit of the above expression $i\to \infty$, we have

$$w,\nu -{\eta }^{\prime }\ge 0.$$

We know that the operator Q is maximally monotone; we have ${\eta }^{\prime }\in \mathrm{VIP}\left(F,C\right)=Q^{-1}\left(0\right).$ Now, we have to prove that the sequence ${\left\{\overline{\eta }\right\}}_{l=1}^{\infty }$ weakly converges to η’ For this, consider that there is a subsequence ${\left\{{\overline{\eta }}^m\right\}}_{m=1}^{\infty }$ of the sequence ${\left\{\overline{\eta }\right\}}_{l=1}^{\infty }$ so that it converges weakly to ${\zeta }^{\prime }\ne {\eta }^{\prime }$. Considering the above statements, we also have ${\zeta }^{\prime }\in \mathrm{VIP}\left(F,C\right)$ and $\underset{l\to \infty }{\lim }\parallel {\overline{\eta }}_l-{\zeta }^{\prime }\parallel$. Using Opial’s condition, we observe

$$\begin{gathered} \underset{l\to \infty }{\lim }\parallel {\overline{\eta }}_l-{\eta }^{\prime }\parallel =\underset{i\to \infty }{\lim }\inf \parallel {\overline{\eta }}^{l_i}-{\eta }^{\prime }\parallel \le \underset{i\to \infty }{\lim }\inf \parallel {\overline{\eta }}^{l_i}-{\zeta }^{\prime }\parallel =\underset{i\to \infty }{\lim }\parallel {\overline{\eta }}^{l_i}-{\zeta }^{\prime }\parallel =\underset{j\to \infty }{\lim }\parallel {\overline{\eta }}^{l_j}-{\zeta }^{\prime }\parallel , \\ \le \underset{j\to \infty }{\lim }\inf \parallel {\overline{\eta }}^{l_j}-{\eta }^{\prime }\parallel =\underset{l\to \infty }{\lim }\parallel {\overline{\eta }}^l-{\eta }^{\prime }\parallel , \end{gathered}$$

which is a paradox. Thus, η’ = ζ’, and hereafter, we accomplish that ${\left\{\overline{\eta }\right\}}_{l=1}^{\infty }$ converges weakly to η’. □

 Proposition 3

([27]). Suppose the averaging matrix ${\left\{a_{l,m}\right\}}_{l,m=1}^{\alpha }$ fulfills the generalized segmenting condition. Then, the averaging matrix is $M$-concentrating if $lim in{f}_{l\to \infty }{a}_{l,l}>0$.

4. Important Results and Discussion

This section is devoted to the detailed study of the proposed method and its effectiveness by minimizing the distance of assumed point. Suppose that $p,r_i\in {ℝ}^n$ and $s_i\ge 0$ are known data, $\forall i=1,2,3,\dots ,l.$ In this examination, we need to explore the controlled minimization model, which is given as:

$$\min \frac{1}{2}\parallel \eta -{p\parallel }^2, \mathrm{subject} \mathrm{to} \langle r_i,\eta \rangle \le s_i, i=1,2,3,\dots ,l,$$

(10)

It is to be noted that the function $f=0.5\parallel ·-{p\parallel }^2$ is the convex Fréchet differentiable function and Δf is the 1-Lipschitz continuous gradient; besides the constrained set $C_i=\left\{\eta \in {ℝ}^n:\langle r_i,\eta \rangle \le s_i\right\}, i=1\dots m$, is a non-empty set which is closed and convex. Therefore, the considered problem (10) appears as problem (1), with $C={\cap }_{i=1}^m$ and Δf = F. It is noted that the operator F is 1-Lipschitz continuous. In this condition, the attained theoretical solutions satisfy and we can use Mann’s type mean extra-gradient scheme for investigating the problem (10). For simplicity, the classical sub-gradient extra-gradient scheme is denoted as SEM, whereas Mann’s type mean extra-gradient scheme is denoted as Mann-MEM with the general segmenting ${\left\{a_{l,m}\right\}}_{l,m=1}^{\infty }$ defined as:

$$a_{l,m}=\left\{\begin{array}{c}{\left(1-a\right)}^{k-1}, \mathrm{if} j=1 \mathrm{and} k\ge 1,\\ 0, \mathrm{if} j\ge 2 \mathrm{and} k<j,\\ a{\left(1-a\right)}^{k-1}, \mathrm{if} j\ge 2 \mathrm{and} k\ge j.\end{array}\right.$$

(11)

In above, $a\in \left(0, 1\right).$ It is noted that the following set

$$T_k=H_{\le }\left(\left({\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)\right)\right)-{\zeta }_k; \langle {\zeta }_k,{\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)-{\zeta }_k\rangle ,$$

is the Mann-MEM and an auxiliary hyperplane to the constrained set C at the point ζ_k In this condition, J_Tk can be calculated explicitly if the approximation ${\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)-{\zeta }_k\ne 0$. However, if the approximation ${\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)-{\zeta }_k=0$, the half-space T_k becomes the full space H such that the iterate η_k+1 is nothing else but the approximation ${\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)$. In order to investigate the solutions, we first make use the traditional Halpern iteration by accomplishing the inner loop: we choose an arbitrary initial point ${\omega }_1\in {ℝ}^n$ and a sequence ${\left\{{\lambda }_i\right\}}_{i=1}^{\infty }$, and calculate

$${\omega }_{i+1}={\lambda }_i\left({\overline{\eta }}_k-\xi F\left({\overline{\eta }}_k\right)+\left(1-{\lambda }_1\right)\right)P_{C_m}{P}_{C_{m-1}}\dots P_{C_2}{P}_{C_1}{\omega }_i,\hspace{1em}\forall i\ge 1.$$

(12)

We use the following stopping criterion for the inner loop in all the computations to find the numerical value of the point ζ_k

$$\frac{\parallel {\omega }_{i+1}-{\omega }_i\parallel }{\parallel {\omega }_i\parallel +1}\le {10}^{-8}.$$

In the first computation, we deliberated the performance of the proposed scheme in a simple condition. For this, we assume m = 3, n = 2, c = [0.1, 0.1]^T, $a_3={\left[1, -2\right]}^T,a_2={\left[1, -1\right]}^T,a_1={\left[-1.5, 1\right]}^T$, and b₃ = b₂ = b₁ = 0. It can be observed that the sole solution is nothing else than this point [0.1, 0.1]^T. Now, let us begin with the effect of the step size ${\lambda }_k=\lambda /\left(1+k\right)$ for numerous choices of $\lambda \in \left(0, 2\right)$ while applying the suggested Mann-MEM and SEM. We select the initial point ${\eta }_1={\left[-0.2, -0.15\right]}^T,$ step size ξ = 0.5, and α = 0.9. Stopping criteria for both Mann-MEM and SEM are $\parallel {\eta }_k-c\parallel \le {10}^{-5}$ or 100 iterations, whichever comes first.

Table 1. Effects of the step size for numerous parameters λ > 0, with execution sub-gradient extra-gradient and Mann mean extra-gradient methods.

Method	λ	Iterations	Time	Inner. Iter.
SEM	1.3	14	0.1826	47,630
	1.4	15	0.1501	37,476
	1.5	15	0.1177	28,591
	1.6	16	0.0906	23,691
	1.7	>100	>0.2871	>84,555
	1.8	30	0.0899	23,648
	1.9	28	0.0699	17,749
Mann-MEM	1.3	>100	>1.0595	>319,322
	1.4	>100	>0.7589	>224,422
	1.5	18	0.118	33,285
	1.6	18	0.0924	25,846
	1.7	19	0.0906	21,495
	1.8	30	0.0851	23,508
	1.9	23	0.0607	15,925

shows that the significant influence of $\lambda$ belongs to [1.3, 1.9] on the number of iterations, computational time, and total number of inner iterations.

Based on Table 1, both schemes give accurate solutions for enhancing the value of λ. This behavior might be possible because of the larger step size, which is given by the parameter λ, as it can dismiss the inner loop in fewer iterations with the intention of reducing the algorithmic runtime. On the other hand, we can observe that Mann-MEM when λ = 1.3, 1.4 and SEM when λ = 1.7 required >100 iterations to meet the stopping criteria. It is noted that when λ = 1.9, both schemes demonstrate excellent solutions. In addition, when λ = 1.9, the scheme Mann-MEM produced excellent results of algorithm runtime 6.01 × 10⁻² seconds. Figure 1 and Figure 2 are plotted against the step size for the discussed schemes. Taking the same assumption as we considered before and setting the inner-loop step size λ_k = 1.9/(1 + k) for both schemes, we see that for both schemes, the best computation time is attained when ξ = 0.6. In order to learn more about the behavior of convergence analysis of the scheme Mann-MEM, we also assume the effect of a. Figure 3 is plotted against the selection of τ = 0.6 and λ_k = 1.9/(1 + k). It is detected that for the large value of a, we attained the lowest number of iterations and computational time; thus, the superlative algorithm’s performance is attained when a = 0.99.

Table 2. Performance of the SEM and Mann-MEM for various dimensions (n) and number of constraints (m).

m	n	Time		Iteration
		Mann-MEM	SEM	Mann-MEM	SEM
50	500	36.3368	38.7986	51.2	51
100		88.4383	94.3647	51	51
200		239.0405	248.4960	51	50
50	1000	58.6253	61.8089	53	52
100		137.0350	143.5451	53	52
200		344.8198	368.2668	52.7	52
50	2000	118.4089	123.3731	54	53.1
100		245.7529	257.4444	54	53
200		576.3775	604.0555	54	53
50	3000	242.2706	247.8855	55	54
100		440.8821	452.0647	55	54
200		1031.5349	1070.5699	55	54

shows the comparison between SEM and Mann-MEM. It is to be noted that the Mann-MEM is more effective than SEM in that Mann-MEM needs less computational time as compared to SEM. One distinguished performance is that when constraints are quite large, the Mann-MEM needs considerably less computational runtime than average.

5. Conclusions

The aim of this research was to discover a solution to the problem of variational inequality that was controlled by a monotone and Lipschitz continuous operator rather than a single operator. We were able to show that the iteration sequence of Mann’s mean extra-gradient technique was ineffective in addressing the problem at hand. In the case of a specified range of acceptable values, the calculations show that the proposed approach shows better convergence behavior than the traditional sub-gradient extra-gradient method, while the conventional sub-gradient extra-gradient method exhibits poorer convergence behavior. Some conclusions are outlined below.

• In order for the Mann-MEM technique to properly converge, the Lipschitz constant of the operator F must be known. If this knowledge is not accessible, the plan is doomed. The letter F in the formula represents the Lipschitz constant for the element F. Given the difficulties in determining the Lipschitz constant, some may question the validity of the conclusion that Mann-MEM and its convergence properties can be used in real-world situations. However, this is not an unreasonable stance to take. For example, among the many interesting Mann-MEM variants are those that utilize a variable step size rather than a fixed step size ${\left\{{\xi }_k\right\}}_{k=1}^{\infty }$, and those that do not need prior knowledge of the L function, such as the Mann-MEM form that does not require prior knowledge of the L function.
• Another finding that should be noted is that when the average matrix ${\left\{{\alpha }_{l.m}\right\}}_{l,m=1}^{\infty }$ is adjusted to its optimum value, Mann-MEM outperforms SEM, as compared to when it is not altered. Indeed, at this point in time, the search for more examples of average matrices that meet the M-concentration criterion is an interesting alternative to consider.
• It is to be noted that the Mann-MEM is more effective than SEM in that Mann-MEM needs less computation work as compared to SEM. One distinguished performance is that when constraints are quite large, the Mann-MEM needs much less computational runtime than average.