Next Article in Journal
The g-Good-Neighbor Conditional Diagnosability of Exchanged Crossed Cube under the MM* Model
Next Article in Special Issue
Improvement of Unconstrained Optimization Methods Based on Symmetry Involved in Neutrosophy
Previous Article in Journal
Fabric Surface Defect Detection Using SE-SSDNet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery

by
Aliyu Muhammed Awwal
1,2,
Adamu Ishaku
1,2,
Abubakar Sani Halilu
3,
Predrag S. Stanimirović
4,
Nuttapol Pakkaranang
5,* and
Bancha Panyanak
6,7,*
1
GSU-Mathematics for Innovative Research Group, Gombe State University (GSU), Gombe 760214, Nigeria
2
Department of Mathematics, Faculty of Science, Gombe State University, Gombe 760214, Nigeria
3
Department of Mathematics, School of Chemical Engineering and Physical Sciences, Lovely Professional University, Phagwara 144411, India
4
Faculty of Science and Mathematics, University of Nis, Visegradksa 33, 18000 Nis, Serbia
5
Mathematics and Computing Science Program, Faculty of Science and Technology, Phetchabun Rajabhat University, Phetchabun 67000, Thailand
6
Research Group in Mathematics and Applied Mathematics, Department of Mathematics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
7
Data Science Research Center, Department of Mathematics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
*
Authors to whom correspondence should be addressed.
Symmetry 2022, 14(11), 2375; https://doi.org/10.3390/sym14112375
Submission received: 18 October 2022 / Revised: 5 November 2022 / Accepted: 6 November 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Nonlinear Analysis and Its Applications in Symmetry II)

Abstract

:
Many practical applications in applied sciences such as imaging, signal processing, and motion control can be reformulated into a system of nonlinear equations with or without constraints. In this paper, a new descent projection iterative algorithm for solving a nonlinear system of equations with convex constraints is proposed. The new approach is based on a modified symmetric rank-one updating formula. The search direction of the proposed algorithm mimics the behavior of a spectral conjugate gradient algorithm where the spectral parameter is determined so that the direction is sufficiently descent. Based on the assumption that the underlying function satisfies monotonicity and Lipschitz continuity, the convergence result of the proposed algorithm is discussed. Subsequently, the efficiency of the new method is revealed. As an application, the proposed algorithm is successfully implemented on image deblurring problem.

1. Preliminaries and Motivation

Consider the space R n with a nonempty closed convex subset, A . The symbol · shall mean the Euclidean norm in the entire paper. Recall that a function V is said to be monotone and Lipschitz continuous if it satisfies the following inequalities
0 ( x y ) T ( V ( x ) V ( y ) ) ,
V ( x ) V ( y ) L x y , L > 0 ,
respectively. Let V : R n R n be a monotone function, this paper is concerned with finding a vector x ^ A that solves the following system of monotone nonlinear equations:
V ( x ) = 0 .
Problem (3) is of paramount importance due to its presence in various fields of applied sciences such as physics, engineering, signal processing, robotic motion control, image deblurring problems, etc. (see, [1,2,3,4,5,6,7,8,9]). Among the different classes of iterative algorithms for solving (3) such as Newton-like methods, quasi-Newton methods and trust-region methods (see, [10,11,12,13,14,15]); conjugate gradient (CG) methods (see, [16,17,18,19,20]) stand out. The CG methods are popular, especially for problems with large dimensions, due to the fascinating features like low memory demand, global convergence and simplicity.
Recently, the interesting projection map presented by Solodov and Svaiter [21] has made the CG method for (3) gain more attention. For instance, based on the popular CG_DESCENT method [22], Xiao and Zhu proposed a projection CG method for solving problem (3). The search direction produced by their algorithm is sufficiently descent independent of any line search technique. The sufficient descent property is an important property for any CG method, especially in the proof of its convergence result. However, Liu and Li [17] pointed out the possibility of the CG parameter in [23] being undefined at certain iteration and therefore, proposed another CG projection method. Furthermore, based on the Dai-Yuan CG method [24], Liu and Feng [25] developed an interesting spectral CG method with nice characteristics. Experimental results showed that their method is relatively efficient.
Another set of popular iterative methods for unconstrained optimization problems, min { v ( x ) R : u R n } , and (3) are the quasi-Newton methods. This class of methods compute their directions of search as p k = M k V ( x k ) ,   M k is an approximation of the Hessian inverse of v at x k . The matrix M k is usually computed using different updating formula such as BFGS (Broyden–Fletcher–Goldfarb–Shanno) and SR1 (Symmetric rank-one formula) [26]. Unlike the CG methods, the quasi-Newton method may be unsuitable when the dimension of the problems is large. This inspires Zhang et al. [27] to incorporate the popular Polak–Ribière–Polyak parameter [28,29] into the BFGS updating formula and developed a three-term matrix-free iterative method for minimization problems without constraints. Motivated by the success of the Zhang et al.’s [27] approach, Awwal et al. [30] came up with a spectral CG-like method for (3) where a modified Perry CG parameter [31] is incorporated into a modified BFGS updating formula. The spectral parameter in their method is determined so that the sequence of the search directions is sufficiently descent.
The SR1 quasi-Newton method is a successful iterative method for minimization problems. However, not much has been done to explore its performance on the system of nonlinear Equation (3), perhaps due to the fact that the denominator of the SR1 updating formula could be zero at some point (detail is given in the next section). Motivated by the successes of the methods discussed thus far, in this paper, a new derivative-free spectral projection method for handling problem (3) is proposed with the aid of a modified SR1 updating formula. A few notable contributions of this work are enumerated below:
1.
The search direction of the proposed method is set up based on the modified SR1 updating formula.
2.
The direction of search has some nice features such as boundedness and so on.
3.
The convergence result of the new method is established under some mild conditions.
4.
The numerical performance of the new method is explored by implementing it on a set of test problems.
5.
The proposed method is implemented on image deblurring problem.
The rest part of this paper is arranged as follows. The proposed method and its motivation, as well as its convergence analysis are expressed in Section 2. The efficiency of the new method is shown in Section 3 and the application is illustrated in Section 4. Finally, some remarks are highlighted in Section 5.

2. DFSR1 Method with Its Convergence Results

Let us start by stating the following assumption.
Assumption A1.
The point x ^ A satisfies V ( x ^ ) = 0 and the function V satisfies (1) and (2).
We briefly recall the quasi-Newton iterative methods, which produce { x k } using the following recurring formula:
x k + 1 = x k + τ k p k , k = 0 , 1 , 2 , ,
τ k is a step-size obtained via a suitable search line strategy and p k is a search direction that is usually supposed to satisfy
V ( x k ) T p k c V ( x k ) 2 , c > 0 .
One of the popular quasi-Newton matrix updating formula, namely Symmetric rank-one (SR1), is given as follows:
M k = μ k M k 1 + ( s k 1 M k 1 y k 1 ) ( s k 1 M k 1 y k 1 ) T ( s k 1 M k 1 y k 1 ) T y k 1 ,
s k 1 = x k x k 1 ,   y k 1 = V ( x k ) V ( x k 1 ) and the parameter μ k > 0 shall be determined. If μ k = 1 for all k , then the SR1 update (6) reduces to the classical one. By letting M k 1 I (i.e., identity matrix), then (6) becomes
M k = μ k I + ( s k 1 y k 1 ) ( s k 1 y k 1 ) T ( s k 1 y k 1 ) T y k 1 .
The quasi-Newton search direction is defined as follows
p k = M k V ( x k ) , k = 0 , 1 , 2 , .
By substituting (7) into (8), we have
p k = μ k V ( x k ) + β k u k ,
where
β k = u k T V ( x k ) y k 1 T s k 1 y k 1 2 ,
and u k = s k 1 y k 1 . This suggests that (9) imitates the actions of the classical spectral CG method. To ensure that the β k (10) is well-defined, we must ensure that y k 1 2 y k 1 T s k 1 . Now, defining y ¯ k 1 = V ( x k ) V ( x k 1 ) + t s k 1 ,   t > 0 ,   s k 1 = x k x k 1 and using the Lipschitz continuity assumption on V , gives
y ¯ k 1 = V ( x k ) V ( x k 1 ) + t ( x k x k 1 ) V ( x k ) V ( x k 1 ) + t x k x k 1 L x k x k 1 + t x k x k 1 ( L + t ) s k 1 .
Additionally, by the monotonicity of V , we have
y ¯ k 1 2 = ( V ( x k ) V ( x k 1 ) + t ( x k x k 1 ) ) T ( V ( x k ) V ( x k 1 ) + t ( x k x k 1 ) ) = V ( x k ) V ( x k 1 ) 2 + 2 t ( V ( x k ) V ( x k 1 ) ) T ( x k x k 1 ) + t 2 x k x k 1 2 V ( x k ) V ( x k 1 ) 2 + t 2 x k x k 1 2 t 2 s k 1 2 .
This, together with (11), gives
( t + L ) 2 s k 1 2 y ¯ k 1 2 t 2 s k 1 2 > 0 .
However, the monotonicity and Lipschitz continuity of V yield
y ¯ k 1 T s k 1 = ( V ( x k ) V ( x k 1 ) + t ( x k x k 1 ) ) T ( x k x k 1 ) = ( V ( x k ) V ( x k 1 ) ) T ( x k x k 1 ) + t x k x k 1 2 t s k 1 2 > 0 ,
and
y ¯ k 1 T s k 1 = ( V ( x k ) V ( x k 1 ) + t ( x k x k 1 ) ) T ( x k x k 1 ) = ( V ( x k ) V ( x k 1 ) ) T ( x k x k 1 ) + t x k x k 1 2 = V ( x k ) V ( x k 1 ) x k x k 1 + t x k x k 1 2 L x k x k 1 2 + t x k x k 1 2 = ( t + L ) s k 1 2 .
Combining (14) and (15) gives
( t + L ) s k 1 2 y ¯ k 1 T s k 1 t s k 1 2 > 0 .
Therefore, it makes sense to redefine the parameter β k as follows
β ¯ k = u ¯ k T V ( x k ) max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } ,
where u ¯ k = s k 1 y ¯ k 1 .
Next, we determine the parameter μ k so that (9) satisfies (5). Now, taking inner product of (9) and V ( x k ) where β ¯ k is given by (17), we have
V ( x k ) T p k = μ k V ( x k ) 2 + β ¯ k V ( x k ) T u ¯ k = μ k V ( x k ) 2 u ¯ k T V ( x k ) max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) T u ¯ k = μ k + ( u ¯ k T V ( x k ) ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 V ( x k ) 2 .
This means by taking
μ k c ( u ¯ k T V ( x k ) ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 , c > 0 ,
we have
V ( x k ) T p k c V ( x k ) 2 , c > 0 , k = 0 , 1 , 2 , .
Now, without loss of generality, we set the value of μ k as
μ k = c ( u ¯ k T V ( x k ) ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 .
Before we state the steps of the algorithm for the proposed method, the following definition is of great importance.
Definition 1.
Any point x in R n can be projected onto the feasible set A via P A ( x ) = argmin { x y : y A } .
The operator P A ( · ) satisfies the following:
P A ( x ) y x y , y A .
We now state the steps of the proposed Algorithm 1 for solving problem (3).    
Algorithm 1: Derivative-free symmetric rank-one method (DFSR1).
Input:  Given x 0 A , ρ ( 0 , 1 ) , t , σ , κ > 0 ,   0 < < 2 ,   T o l 0 .
Step 1: Determine V ( x k ) . If V ( x k ) T o l , terminate.
Step 2: Calculate p 0 = V ( x 0 ) and
p k = max { μ k , λ k } V ( x k ) + β ¯ k u ¯ k ,
  where u ¯ k = s k 1 y ¯ k 1 ,   β ¯ k and μ k are given by (17) and (21), respectively and
λ k = s k 1 2 y ¯ k 1 T s k 1 .
 Step 3: Set h k = x k + τ k p k ,   τ k = κ ρ i with i being the least non-negative integer satisfying
V ( x k + κ ρ i p k ) T p k σ κ ρ i V ( x k + κ ρ i p k ) 1 / q p k 2 , q 1 .
 Step 4: If V ( h k ) = 0 , terminate the process. If not, compute the new point by
x k + 1 = P A x k V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) .
 Step 5: Increase the counter, i.e., k : = k + 1 and repeat from step 1.
Remark 1.
By the definition of μ k in (21), we cannot guarantee μ k > 0 for all k . Considering the modified Barzilai and Borwein spectral parameter λ k = s k 1 2 y ¯ k 1 s k 1 and using the inequalities (13) and (16), we have
1 t s k 1 2 y ¯ k 1 s k 1 1 t + L > 0 .
Therefore, the use of max { μ ¯ k , λ k } in the definition of the search direction (23) makes sense since max { μ ¯ k , λ k } > 0 , for all k .
Remark 2.
The direction (23) satisfies (20). If max { μ k , λ k } = μ k , then the result follows from (18). Moreover, if max { μ k , λ k } = λ k , then λ k μ k and therefore from (19) we have
λ k μ k c ( u ¯ k T V ( x k ) ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 ,
which means (20) holds.
Remark 3.
The line search (25) is adopted from [32] and its well-definedness has been discussed in [33].
Lemma 1.
Let the sequences { x k } and { h k } be produced by DFSR1. Suppose the function V is Lipschitzian and monotone, then the lim k x k x ^ exists and
p k Γ , for all k 0 , Γ > 0 .
lim k τ k p k = 0 .
Proof. 
Let x ^ A such that V ( x ^ ) = 0 , then V ( x ^ ) T ( x k x ^ ) = 0 . Therefore, it holds that
V ( h k ) T ( x k x ^ ) = V ( h k ) T ( x k h k + h k x ^ ) = V ( h k ) T ( x k h k ) + V ( h k ) T ( h k x ^ ) V ( h k ) T ( x k h k ) + V ( x ^ k ) T ( h k x ^ ) = V ( h k ) T ( x k h k ) ,
Note that the above inequality has been obtained using (1), that is, V ( h k ) T ( h k x ^ ) V ( x ^ ) T ( h k x ^ ) . Next, since 0 < < 2 , then by (22) and (26), we have
x k + 1 x ^ 2 = P A x k V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) x ^ x k V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) x ^ = ( x k x ^ ) V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) 2 x k x ^ 2 2 V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) T ( x k x ^ ) + 2 [ V ( h k ) T ( x k h k ) ] 2 V ( h k ) 2 x k x ^ 2 2 V ( h k ) T ( x k h k ) V ( h k ) 2 V ( h k ) T ( x k h k ) + 2 [ V ( h k ) T ( x k h k ) ] 2 V ( h k ) 2 = x k x ^ 2 ( 2 ) [ V ( h k ) T ( x k h k ) ] 2 V ( h k ) 2
x k x ^ 2 .
Please observe that (32) gives the existence of lim k x k x ^ .
Next, suppose the max { μ k , λ k } = μ k , then from (23), we have
p k = μ k V ( x k ) + β ¯ k u ¯ k | μ k | V ( x k ) + | β ¯ k | u ¯ k = c ( u ¯ k T V ( x k ) ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 V ( x k ) + u ¯ k T V ( x k ) max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } u ¯ k c + u ¯ k 2 V ( x k ) 2 max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } V ( x k ) 2 V ( x k ) + u ¯ k V ( x k ) max { y ¯ k 1 T s k 1 , y ¯ k 1 2 } u ¯ k c + 2 u ¯ k 2 y ¯ k 1 2 V ( x k ) c + 2 ( 1 + t + L ) 2 s k 1 2 t 2 s k 1 2 V ( x k ) c + 2 ( 1 + t + L ) 2 t 2 V ( x k ) .
Since { x k x ^ } is a decreasing sequence, then x k x ^     x 0 x ^ for all k. Setting c 1 = L x 0 x ^ , then by the Lipschitzian continuity of V , we have
V ( x k )   =   V ( x k )     V ( x ^ )     L x k x ^         L x 0 x ^   =   c 1 .
Combining (33) and (34) gives (28) with Γ = c + 2 ( 1 + t + L ) 2 t 2 c 1 . Moreover, if max { μ k , λ k } = λ k , then by (27) the result of (28) equally holds.
Furthermore, since lim k x k x ^ exists, then { x k } is bounded. Therefore by (28) and h k in Step 3 of DFSR1, it means { h k } is bounded. Using similar argument as in (34) gives
V ( h k ) c 2 , c 2 > 0 .
Now, (25) and (31) give
σ 2 τ k 4 V ( h k ) 2 / q p k 4 τ k 2 [ V ( h k ) T p k ] 2 V ( h k ) 2 ( 2 ) x k x ^ 2 x k + 1 x ^ 2 .
(36) together with (35) gives
σ 2 τ k 4 p k 4 V ( h k ) 2 2 / q ( 2 ) x k x ^ 2 x k + 1 x ^ 2 c 2 2 2 / q ( 2 ) x k x ^ 2 x k + 1 x ^ 2 .
Since 0 < < 2 , then taking limits of both sides as k gives
σ 2 lim k τ k 4 p k 4 = 0 .
Hence,
lim k τ k p k = 0 .
Theorem 1.
Let { x k } be a sequence produced by the DFSR1 and x ^ A such that V ( x ^ ) = 0 then { x k } converges to x ^ .
Proof. 
Let the function V be Lipschitz continuous. First of all, we claim that
lim k inf V ( x k ) = 0 .
Assume (39) is false, then
V ( x k ) ω , ω > 0 , holds k 0 .
If τ k κ , by setting h k = x k + τ k p k , then τ k = ρ 1 τ k will not satisfy (25). Therefore,
V ( h k ) T p k + σ τ k V ( h k ) 1 / q p k 2 > 0 .
Since (41) holds, then using (20) yields
c V ( x k ) 2 V ( x k ) T p k < V ( x k ) T p k + V ( h k ) T p k + σ τ k V ( h k ) 1 / q p k 2 = ( V ( h k ) V ( x k ) ) T p k + σ τ k V ( h k ) 1 / q p k 2 L τ k p k 2 + σ τ k V ( h k ) 1 / q p k 2 = ρ 1 τ k ( L + σ V ( h k ) 1 / q ) p k 2 ρ 1 τ k ( L + σ c 2 1 / q ) Γ 2 .
(42) can be re-written as
τ k > ρ c V ( x k ) 2 Γ 2 ( L + σ c 2 1 / q ) ρ c ω 2 Γ 2 ( L + σ c 2 1 / q ) .
Using Cauchy–Schwartz inequality on (20) gives
c V ( x k ) 2 V ( x k ) T p k V ( x k ) p k .
This together with (40) yields
p k c ω .
Combining (43) and (45) gives
τ k p k > ρ c ω 2 Γ 2 ( L + σ c 2 1 / q ) c ω ,
which contradicts (29). Hence, (39) must hold.
Lastly, by the continuity assumption on V, then { x k } possesses some limit point x ^ for which V ( x ^ ) = 0 . However, we have shown { x k x ^ } converges, which means that the conclusion holds. □

3. Numerical Experiments

We now illustrate the numerical performance of the DFSR1 in comparison with: (i) PDY [25] and (ii) HCGP [34]. In the course of the experiment, we implement the three algorithms (that is, DFSR1, PDY and HCGP) on a collection of test problems (see, Appendix A). We vary the dimension of each test problem as 1000, 5000, 10,000, 50,000, and 100,000 and solve each with the starting points listed in Table 1.
We implement these algorithms in MATLAB R2019b that runs on a PC with intel Core(TM) i5-8250u processor with 4 GB of RAM and CPU 1.60 GHz. We set the termination criterion for the iteration process as V ( x k ) 10 6 . We use the following parameters for Algorithm 1 (DFSR1): ρ = 0.5 ,   c = 0.1 ,   t = 0.01 ,   σ = 0.01 , κ = 1 ,   = 1.99 . The parameters used in the implementation of PDY and HCGP are as given in [25,34].
We use three metrics, namely, ITER (number of iterations), FVAL (number of functions evaluation), and TIME (CPU time), to assess the performance of each algorithm under consideration. We report the ITER, FVAL and TIME recorded by each algorithm in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9. In addition, we report the NORM (final norm of the V ( x k ) as the iteration process stops) to show whether an algorithm obtains a solution. Observing through the NORM in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9, we note that all the algorithms in question satisfied the terminating criteria set. However, the efficiency and performance of each algorithm differ based on the ITER, FVAL and TIME. Perusing through the ITER, FVAL and TIME in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9, we see that DFSR1 recorded smallest ITER, FVAL and TIME in most cases of the entire experiment. This means that DFSR1 was able to meet the terminating criteria faster and with the least ITER and FVAL for most of the problems. We present the summary of these information in Figure 1, Figure 2 and Figure 3 with the aid of the popular Dolan and Moré performance profile [35]. This shows DFSR1 is efficient and outperforms the two existing competitors, that is, PDY [25] and HCGP [34].

4. Application in Image Deblurring

Consider the convex unconstrained minimization problem
min x R n v ( x ) , v ( x ) = 1 2 y Q x 2 2 + η x 1 , η > 0
y R k is an observation and the linear function Q is in R k × n ( k < < n ). The above model (46) has an interesting application in image deblurring.
The 1 regularization model (46) has attracted a great deal of attention and some researchers have developed different iterative algorithms for handling it (see, [36,37,38]). By reformulating (46) into a bound-constrained quadratic program, Figueiredo et al. [39] developed gradient-based projection procedure. The reformulation in [39] is as follows: let the vector x be split into two parts, i.e., x = a b , where a i = ( x i ) + , b i = ( x i ) + for all i = 1 , 2 , , n , and ( . ) + = max { 0 , . } , then (46) is translated into
min w 0 1 2 w T Z w + r T w ,
w = [ a b ] T ,   r = η e 2 n + Q T y Q T y T and Z = Q T Q Q T Q Q T Q Q T Q .
Taking it further, Xiao et al. [40] showed that (47) is equivalent to the following system of nonlinear equation
V ( w ) = min { w , Z w + r } = 0 ,
the “min” means componentwise minimum. Finally, while Pang [41] showed that V satisfies (2), Xiao et al. [40] proved that it also satisfies (1).
One of the immediate advantages associated with the system of monotone nonlinear equations (48) over (46) is that derivative-free algorithms can be employed to efficiently solve it. One of the derivative-free algorithms applied on (48) to recover some disturbed signals as well as restore some blurred images is the SGCS algorithm [40] by Xiao et al. A few years later, based on the popular CG_DESCENT method [42], Xiao and Zhu [23] came up with another derivative-free algorithm (CGD), which performs better than the SGCS algorithm on signal recovery problem. Following this trend, we apply the DFSR1 to recover some blurred color images by solving the model (48). We compare the performance of the DFSR1 with the CGD of Xiao and Zhu [23]. Details of the image deblurring experiment are given as follows.
We coded the two algorithms in MATLAB R2019b and carried out the image deblurring experiment on a PC with an intel Core(TM) i5-8250u processor with 4 GB of RAM and CPU 1.60 GHz. We run each algorithm from the same initial point x 0 = Q T y and terminate it if | v k v k 1 | | v k 1 | < 10 5 , is satisfied, where v k is the merit function evaluation at x k , with v ( x ) = 1 2 y Q x 2 2 + η x 1 . The parameters used for both DFSR1 and CGD in this experiment are taken from [23]. We used four color images for the experiment and the following metrics are used for assessment: (i) ITER, (ii) TIME, (iii) SNR (signal-to-noise-ratio), and (iv) SSIM (structural similarity index measure).
The two algorithms, that is, DFSR1 and CGD, successfully recovered all the four images after being blurred, albeit with different rapidity and level of quality. Original and blurred images, along with images restored by each of the considered algorithms are presented in Figure 4. Furthermore, the data on ITER, TIME, SNR and SSIM are presented in Table 10. In addition, the size of each image used in the experiment are given. The values in Table 10 revealed that DFSR1 recorded less ITER and TIME than CGD in the experiment. As for the quality of the restoration, the SNR and SSIM recorded by DFSR1 reveal that the quality of recoveries by DFSR1 are slightly better than those by the CGD. This shows that the proposed DFSR1 is efficient and applicable.

5. Conclusions

Based on a modified symmetric rank-one matrix updating formula and a quasi-Newton search direction, a new descent direction of searching has been proposed. The proposed algorithm, in this paper, incorporated the popular Solodov and Svaiter projection techniques and a suitable line search used for obtaining the step length. The boundedness of the new direction and the theoretical analysis of the algorithm has been established based on some mild assumptions. The efficiency of the new method as well as its applicability have been illustrated. As future research, the proposed DFSR1 can be modified and implemented on nonlinear least squares problems [43,44].

Author Contributions

Conceptualization, A.M.A. and A.I.; methodology, A.M.A.; software, A.M.A. and N.P.; validation, A.M.A., A.I., A.S.H., P.S.S., N.P. and B.P.; formal analysis, A.M.A., P.S.S. and N.P.; investigation, A.M.A., A.I., A.S.H., P.S.S., N.P. and B.P.; resources, A.M.A., N.P. and B.P.; data curation, A.M.A., A.I., A.S.H. and P.S.S.; writing-original draft preparation, A.M.A. and A.S.H.; writing-review and editing, A.M.A., A.I., A.S.H., P.S.S., N.P. and B.P.; visualization, A.M.A., A.I., A.S.H. and P.S.S.; supervision, A.M.A. and P.S.S.; project administration, A.M.A., N.P. and B.P.; funding acquisition, N.P. and B.P. All authors have read and agreed to the published version of the manuscript.

Funding

The fifth author was partially supported by Phetchabun Rajabhat University. The sixth author was funded by Chiang Mai University and the NSRF via the Program Management Unit for Human Resources & Institutional Development, Research and Innovation (grant number B05F640183).

Acknowledgments

The authors would like to thank the anonymous referees and editor for reading this paper carefully, providing valuable suggestions and comments, which grately improved the final version.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

We use the following monotone nonlinear equation for the second experiments where V ( x ) = ( v 1 ( x ) , v 2 ( x ) , , v n ( x ) ) T , and x = ( x 1 , x 2 , , x n ) T .
Problem A1
([45]).
v 1 ( x ) = e x 1 1 v i ( x ) = e x i + x i 1 1 , i = 1 , 2 , , n 1 ,
where A = R + n .
Problem A2
([46]).
v i ( x ) = 2 x i sin | x i | , i = 1 , 2 , , n ,
where A = R + n .
Problem A3
([47]).
v i ( x ) = e x i 1 , i = 1 , 2 , , n ,
where A = R + n .
Problem A4
([48]).
v i ( x ) = e x i 2 + 3 2 sin ( 2 x i ) 1 , i = 1 , 2 , , n ,
where A = R + n .
Problem A5
([49] ).
v i ( x ) = x i sin ( | x i 1 | ) , i = 1 , 2 , , n 1 ,
where A = { x R n : i = 1 n x i n , x i 1 , i = 1 , 2 , , n } .
Problem A6.
v 1 ( x ) = 2 x 1 x 2 + e x 1 1 , v i ( x ) = x i 1 + 2 x i x i + 1 + e x i 1 , i = 2 , , n 1 , v n ( x ) = x n 1 + 2 x n + e x n 1 ,
where A = R + n .
Problem A7
([25]).
v 1 ( x ) = 5 2 x 1 + x 2 1 , v i ( x ) = x i 1 + 5 2 x i + x i + 1 1 , i = 2 , , n 1 , v n ( x ) = x n 1 + 5 2 x n 1 ,
where A = R + n .
Problem A8.
v 1 ( x ) = x 1 + sin x 1 1 v i ( x ) = x i 1 + 2 x i + sin x i 1 , i = 2 , , n 1 v n ( x ) = x n + sin x n 1 .
where A = R + n .

References

  1. Sulaiman, I.M.; Awwal, A.M.; Malik, M.; Pakkaranang, N.; Panyanak, B. A Derivative-Free MZPRP Projection Method for Convex Constrained Nonlinear Equations and Its Application in Compressive Sensing. Mathematics 2022, 10, 2884. [Google Scholar] [CrossRef]
  2. Al-Sawalha, M.M.; Agarwal, R.P.; Shah, R.; Ababneh, O.Y.; Weera, W. A reliable way to deal with fractional-order equations that describe the unsteady flow of a polytropic gas. Mathematics 2022, 10, 2293. [Google Scholar] [CrossRef]
  3. Mohammad, H.; Awwal, A.M. Globally convergent diagonal Polak–Ribière–Polyak like algorithm for nonlinear equations. Numer. Algorithms 2022, 1–20. [Google Scholar] [CrossRef]
  4. Al-Sawalha, M.M.; Alshehry, A.S.; Nonlaopon, K.; Shah, R.; Ababneh, O.Y. Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm. AIMS Math. 2022, 7, 19739–19757. [Google Scholar] [CrossRef]
  5. Halilu, A.S.; Majumder, A.; Waziri, M.Y.; Ahmed, K.; Awwal, A.M. Motion control of the two joint planar robotic manipulators through accelerated Dai–Liao method for solving system of nonlinear equations. Eng. Comput. 2021, 39, 1802–1840. [Google Scholar] [CrossRef]
  6. Awwal, A.M.; Kumam, P.; Sitthithakerngkiet, K.; Bakoji, A.M.; Halilu, A.S.; Sulaiman, I.M. Derivative-free method based on DFP updating formula for solving convex constrained nonlinear monotone equations and application. AIMS Math. 2021, 6, 8792–8814. [Google Scholar] [CrossRef]
  7. Halilu, A.S.; Majumder, A.; Waziri, M.Y.; Awwal, A.M.; Ahmed, K. On solving double direction methods for convex constrained monotone nonlinear equations with image restoration. Comput. Appl. Math. 2021, 40, 1–27. [Google Scholar] [CrossRef]
  8. Al-Sawalha, M.M.; Khan, A.; Ababneh, O.; Botmart, T. Fractional view analysis of Kersten-Krasil’shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Math. 2022, 7, 18334–18359. [Google Scholar] [CrossRef]
  9. Al-Sawalha, M.M.; Alshehry, A.S.; Nonlaopon, K.; Shah, R.; Ababneh, O.Y. Fractional view analysis of delay differential equations via numerical method. AIMS Math. 2022, 7, 20510–20523. [Google Scholar] [CrossRef]
  10. Wan, Z.; Chen, Y.; Huang, S.; Feng, D. A modified nonmonotone BFGS algorithm for solving smooth nonlinear equations. Optim. Lett. 2014, 8, 1845–1860. [Google Scholar] [CrossRef]
  11. Leong, W.J.; Hassan, M.A.; Waziri, M.Y. A matrix-free quasi-Newton method for solving large-scale nonlinear systems. Comput. Math. Appl. 2011, 62, 2354–2363. [Google Scholar] [CrossRef] [Green Version]
  12. Mohammad, H.; Waziri, M.Y. On Broyden-like update via some quadratures for solving nonlinear systems of equations. Turk. J. Math. 2015, 39, 335–345. [Google Scholar] [CrossRef]
  13. Nocedal, J.; Wright, S. Numerical Optimization; Springer Science & Business Media: New York, NY, USA, 2006. [Google Scholar]
  14. Ahookhosh, M.; Esmaeili, H.; Kimiaei, M. An effective trust-region-based approach for symmetric nonlinear systems. Int. J. Comput. Math. 2013, 90, 671–690. [Google Scholar] [CrossRef]
  15. Yuan, Y.X. Recent advances in trust region algorithms. Math. Program. 2015, 151, 249–281. [Google Scholar] [CrossRef]
  16. Cheng, W.; Xiao, Y.; Hu, Q. A family of derivative-free conjugate gradient methods for large-scale nonlinear systems of equations. J. Comput. Appl. Math. 2009, 224, 11–19. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, J.K.; Li, S.J. A projection method for convex constrained monotone nonlinear equations with applications. Comput. Math. Appl. 2015, 70, 2442–2453. [Google Scholar] [CrossRef]
  18. Waziri, M.Y.; Muhammad, H.U.; Halilu, A.S.; Ahmed, K. Modified matrix-free methods for solving system of nonlinear equations. Optimization 2021, 70, 2321–2340. [Google Scholar] [CrossRef]
  19. Awwal, A.M.; Kumam, P.; Abubakar, A.B. Spectral modified Polak–Ribiére–Polyak projection conjugate gradient method for solving monotone systems of nonlinear equations. Appl. Math. Comput. 2019, 362, 124514. [Google Scholar] [CrossRef]
  20. Abubakar, A.B.; Kumam, P.; Awwal, A.M.; Thounthong, P. A modified self-adaptive conjugate gradient method for solving convex constrained monotone nonlinear equations for signal recovery problems. Mathematics 2019, 7, 693. [Google Scholar] [CrossRef] [Green Version]
  21. Solodov, M.V.; Svaiter, B.F. A globally convergent inexact Newton method for systems of monotone equations. In Reformulation: Nonsmooth, Piecewise Smooth, Semismooth and Smoothing Methods; Springer: Boston, MA, USA, 1998; pp. 355–369. [Google Scholar]
  22. Hager, W.W.; Zhang, H. Algorithm 851: CG_DESCENT, a conjugate gradient method with guaranteed descent. ACM Trans. Math. Softw. (TOMS) 2006, 32, 113–137. [Google Scholar] [CrossRef]
  23. Xiao, Y.; Zhu, H. A conjugate gradient method to solve convex constrained monotone equations with applications in compressive sensing. J. Math. Anal. Appl. 2013, 405, 310–319. [Google Scholar] [CrossRef]
  24. Dai, Y.; Yuan, Y. A nonlinear conjugate gradient method with a strong global convergence property. SIAM J. Optim. 1999, 10, 177–182. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, J.; Feng, Y. A derivative-free iterative method for nonlinear monotone equations with convex constraints. Numer. Algorithms 2019, 82, 245–262. [Google Scholar] [CrossRef]
  26. Sun, W.; Yuan, Y.X. Optimization Theory and Methods: Nonlinear Programming; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006; Volume 1. [Google Scholar]
  27. Zhang, L.; Zhou, W.; Li, D.H. A descent modified Polak–Ribière–Polyak conjugate gradient method and its global convergence. IMA J. Numer. Anal. 2006, 26, 629–640. [Google Scholar] [CrossRef]
  28. Pola, E.; Ribiere, G. Note sur la convergence de methodes de directions conjugées. Rev. Fr. Informat Rech. Oper. Annee 1969, 16, 35–43. [Google Scholar]
  29. Polyak, B.T. The conjugate gradient method in extremal problems. USSR Comput. Math. Math. Phys. 1969, 9, 94–112. [Google Scholar] [CrossRef]
  30. Awwal, A.M.; Kumam, P.; Mohammad, H.; Watthayu, W.; Abubakar, A. A Perry-type derivative-free algorithm for solving nonlinear system of equations and minimizing 1 regularized problem. Optimization 2020, 70, 1231–1259. [Google Scholar] [CrossRef]
  31. Perry, A. A modified conjugate gradient algorithm. Oper. Res. 1978, 26, 1073–1078. [Google Scholar] [CrossRef] [Green Version]
  32. Awwal, A.M.; Wang, L.; Kumam, P.; Mohammad, H.; Watthayu, W. A projection Hestenes–Stiefel method with spectral parameter for nonlinear monotone equations and signal processing. Math. Comput. Appl. 2020, 25, 27. [Google Scholar] [CrossRef]
  33. Awwal, A.M.; Wang, L.; Kumam, P.; Mohammad, H. A two-step spectral gradient projection method for system of nonlinear monotone equations and image deblurring problems. Symmetry 2020, 12, 874. [Google Scholar] [CrossRef]
  34. Koorapetse, M.; Kaelo, P. An efficient hybrid conjugate gradient-based projection method for convex constrained nonlinear monotone equations. J. Interdiscip. Math. 2019, 22, 1031–1050. [Google Scholar] [CrossRef]
  35. Dolan, E.D.; Moré, J.J. Benchmarking optimization software with performance profiles. Math. Program. Ser. 2002, 91, 201–213. [Google Scholar] [CrossRef]
  36. Figueiredo, M.A.; Nowak, R.D. An EM algorithm for wavelet-based image restoration. IEEE Trans. Image Process. 2003, 12, 906–916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. De Mol, C.; Defrise, M. A note on wavelet-based inversion algorithms. Contemp. Math. 2002, 313, 85–96. [Google Scholar]
  38. Yang, J.; Yin, W.; Zhang, Y.; Wang, Y. A fast algorithm for edge-preserving variational multichannel image restoration. SIAM J. Imaging Sci. 2009, 2, 569–592. [Google Scholar] [CrossRef]
  39. Figueiredo, M.; Nowak, R.D.; Wright, S.J. Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems. IEEE J. Sel. Top. Signal Process. 2007, 1, 586–597. [Google Scholar] [CrossRef] [Green Version]
  40. Xiao, Y.; Wang, Q.; Hu, Q. Non-smooth equations based method for 1-norm problems with applications to compressed sensing. Nonlinear Anal. Theory, Methods Appl. 2011, 74, 3570–3577. [Google Scholar] [CrossRef]
  41. Pang, J. Inexact Newton methods for the nonlinear complementarity problem. Math. Program. 1986, 36, 54–71. [Google Scholar] [CrossRef]
  42. Hager, W.W.; Zhang, H. A new conjugate gradient method with guaranteed descent and an efficient line search. SIAM J. Optim. 2005, 16, 170–192. [Google Scholar] [CrossRef] [Green Version]
  43. Yahaya, M.M.; Kumam, P.; Awwal, A.M.; Aji, S. Alternative structured spectral gradient algorithms for solving nonlinear least-squares problems. Heliyon 2021, 7, e07499. [Google Scholar] [CrossRef]
  44. Awwal, A.M.; Kumam, P.; Wang, L.; Yahaya, M.M.; Mohammad, H. On the Barzilai–Borwein gradient methods with structured secant equation for nonlinear least squares problems. Optim. Methods Softw. 2020, 1–20. [Google Scholar] [CrossRef]
  45. La Cruz, W.; Martínez, J.M.; Raydan, M. Spectral residual method without gradient information for solving large-scale nonlinear systems: Theory and experiments. Citeseer 2004. Available online: http://kuainasi.ciens.ucv.ve/mraydan/downloadpapers/TechRep.pdf (accessed on 17 October 2022).
  46. La Cruz, W. A spectral algorithm for large-scale systems of nonlinear monotone equations. Numer. Algorithms 2017, 76, 1109–1130. [Google Scholar] [CrossRef]
  47. Zhou, W.; Shen, D. An Inexact PRP Conjugate Gradient Method for Symmetric Nonlinear Equations. Numer. Funct. Anal. Optim. 2014, 35, 370–388. [Google Scholar] [CrossRef]
  48. Gao, P.; He, C.; Liu, Y. An adaptive family of projection methods for constrained monotone nonlinear equations with applications. Appl. Math. Comput. 2019, 359, 1–16. [Google Scholar] [CrossRef]
  49. Yu, G.; Niu, S.; Ma, J. Multivariate spectral gradient projection method for nonlinear monotone equations with convex constraints. J. Ind. Manag. Optim. 2013, 9, 117–129. [Google Scholar] [CrossRef]
Figure 1. Performance profile for ITER.
Figure 1. Performance profile for ITER.
Symmetry 14 02375 g001
Figure 2. Performance profile for FVAL.
Figure 2. Performance profile for FVAL.
Symmetry 14 02375 g002
Figure 3. Performance profile for TIME.
Figure 3. Performance profile for TIME.
Symmetry 14 02375 g003
Figure 4. First row (original images), second row (blurred images), third and last rows (recovered images by methods DFSR1 and CGD, respectively).
Figure 4. First row (original images), second row (blurred images), third and last rows (recovered images by methods DFSR1 and CGD, respectively).
Symmetry 14 02375 g004
Table 1. List of the starting point used for the experiment.
Table 1. List of the starting point used for the experiment.
Starting Points (SP)Values
x 1 ( 1 10 , 1 10 , 1 10 , , 1 10 ) T
x 2 ( 1 2 , 1 2 2 , 1 2 3 , , 1 2 n ) T
x 3 ( 2 , 2 , 2 , , 2 ) T
x 4 ( 1 , 1 2 , 1 3 , , 1 n ) T
x 5 ( 1 1 n , 1 2 n , 1 3 n , , 0 ) T
x 6 rand ( 0 , 1 )
Table 2. Numerical results of DFSR1, PDY and HCGP on Problem 1.
Table 2. Numerical results of DFSR1, PDY and HCGP on Problem 1.
Problem 1 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×16130.051812.2 × 10 9 36740.0316513.49 × 10 7 18380.0401986.82 × 10 7
×2250.0071944.44 × 10 16 42850.0566198.84 × 10 7 20410.0145432.79 × 10 7
×3490.004168046940.0242122.21 × 10 8 20410.0071449.54 × 10 7
×45120.0034284.85 × 10 8 40820.0232133.16 × 10 7 17360.0052646.81 × 10 7
×5130.0039150511040.0310019.89 × 10 7 19390.0088877.54 × 10 7
×6130.001159048980.0313947.21 × 10 7 19390.0089467.44 × 10 7
5000×17150.0158755.38 × 10 8 27560.0553081.3 × 10 8 19390.0429787.61 × 10 7
×2250.0103994.44 × 10 16 42850.114468.84 × 10 7 20410.0186232.79 × 10 7
×3490.010012028580.0528974.05 × 10 7 21430.0408935.3 × 10 7
×45120.0109184.85 × 10 8 36740.0679283.44 × 10 7 17360.020576.82 × 10 7
×5130.0047720511040.140597.98 × 10 7 19400.0307158.34 × 10 7
×6130.0031040501020.112677.91 × 10 7 19400.0341058.28 × 10 7
10,000×16130.0239783.74 × 10 8 29600.108142.14 × 10 8 18380.051385.38 × 10 7
×2250.0080674.44 × 10 16 42850.12698.84 × 10 7 20410.0544152.79 × 10 7
×3490.013192031640.126931.18 × 10 7 20410.0497817.49 × 10 7
×45120.0207754.85 × 10 8 41840.148471.1 × 10 7 17360.0660336.82 × 10 7
×5130.005923047950.195185.62 × 10 7 20410.0540065.89 × 10 7
×6130.0068290541090.16835.66 × 10 7 18370.0507085.87 × 10 7
50,000×1490.078887028580.41788.48 × 10 8 19390.226356.01 × 10 7
×2250.0309884.44 × 10 16 42850.53858.84 × 10 7 20410.208552.79 × 10 7
×3490.088267044900.671252.61 × 10 7 20420.263238.36 × 10 7
×45120.079494.85 × 10 8 45920.669282.84 × 10 7 17360.220426.82 × 10 7
×5130.0229590511040.831046.27 × 10 7 19400.222586.57 × 10 7
×6130.020248046940.722756.27 × 10 7 17360.22416.57 × 10 7
100,000×1490.12576021440.747857.7 × 10 7 18370.385168.5 × 10 7
×2250.057414.44 × 10 16 42851.09368.84 × 10 7 20410.435282.79 × 10 7
×3490.13963037761.15624.43 × 10 8 21430.480125.91 × 10 7
×45120.132494.85 × 10 8 32660.768967.38 × 10 7 17360.347436.82 × 10 7
×5130.0758830511041.44137.08 × 10 7 19400.453439.3 × 10 7
×6130.0367550511041.48017.07 × 10 7 18380.403789.29 × 10 7
Table 3. Numerical results of DFSR1, PDY and HCGP on Problem 2.
Table 3. Numerical results of DFSR1, PDY and HCGP on Problem 2.
Problem 2 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×1250.002032021440.0429257.52 × 10 7 11230.0120527.5 × 10 7
×2130.001357019400.0196415.27 × 10 7 9200.0088229.88 × 10 7
×3130.000939024500.0330428.19 × 10 7 130.0022220
×4130.001431020420.0278115.35 × 10 7 10210.0061968.65 × 10 7
×5130.001253023480.0075259.58 × 10 7 12250.005937.76 × 10 7
×6130.000854023480.0155929.65 × 10 7 12250.0039917.92 × 10 7
5000×1250.00407022460.0315342.34 × 10 22 11240.0164448.38 × 10 7
×2130.003971019400.0324895.27 × 10 7 9200.0160429.88 × 10 7
×3130.003208026540.0469879.45 × 10 7 130.0023450
×4130.007746020420.033225.35 × 10 7 10210.0098368.65 × 10 7
×5130.003147025520.0415255.36 × 10 7 12260.0158028.68 × 10 7
×6130.003627026540.0398635.31 × 10 7 12260.0243568.64 × 10 7
10,000×1250.009787022460.064663.31 × 10 22 12250.0285745.93 × 10 7
×2130.004975019400.107975.27 × 10 7 9200.0229779.88 × 10 7
×3130.00518027560.0789826.68 × 10 7 130.0044610
×4130.0054020420.0665585.35 × 10 7 10210.0215418.65 × 10 7
×5130.005071025520.0741317.58 × 10 7 13270.0413916.14 × 10 7
×6130.009867025520.074157.57 × 10 7 13270.0328966.14 × 10 7
50,000×1250.025063024500.522596.65 × 10 7 12260.191026.63 × 10 7
×2130.04304019400.403025.27 × 10 7 9200.0767669.88 × 10 7
×3130.020766027560.372742.37 × 10 20 130.0131160
×4130.0234020420.205545.35 × 10 7 10210.0909168.65 × 10 7
×5130.016834026540.272158.48 × 10 7 13280.168456.86 × 10 7
×6130.01886026540.559138.47 × 10 7 13280.123616.86 × 10 7
100,000×1250.057742019400.813373.35 × 10 20 12260.213129.37 × 10 7
×2130.038002019400.367345.27 × 10 7 9200.175639.88 × 10 7
×3130.029571034701.04365.62 × 10 7 130.0199060
×4130.033132020420.82765.35 × 10 7 10210.244438.65 × 10 7
×5130.048638027560.709749.24 × 10 7 13280.243689.7 × 10 7
×6130.032773027560.948739.24 × 10 7 13280.311029.72 × 10 7
Table 4. Numerical results of DFSR1, PDY and HCGP on Problem 3.
Table 4. Numerical results of DFSR1, PDY and HCGP on Problem 3.
Problem 3 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×1250.001590130.0238330130.0022520
×2250.003072.22 × 10 16 130.0030509190.0026647.45 × 10 7
×3370.0031730130.0013580250.0024270
×4370.0023151.73 × 10 8 5120.0024987.88 × 10 9 11230.0041187.75 × 10 7
×5490.0038673.2 × 10 12 29590.0101699.22 × 10 7 12250.0035389.66 × 10 7
×6490.0030055.53 × 10 12 14300.0093981.71 × 10 9 12250.0036528.58 × 10 7
5000×1250.0037770130.0021940130.0020150
×2250.0047942.22 × 10 16 130.00209309190.0124837.45 × 10 7
×3370.0057910130.0033940250.0143480
×4370.0073811.85 × 10 8 5120.006866.9 × 10 9 11230.0134587.88 × 10 7
×5490.010287.25 × 10 12 25520.0196556.94 × 10 7 13270.0097095.4 × 10 7
×6490.00954.8 × 10 12 21430.0312174.46 × 10 7 12260.0148549.65 × 10 7
10,000×1250.0077380130.008110130.0031980
×2250.007192.22 × 10 16 130.00888409190.0173697.45 × 10 7
×3370.0107950130.0117140250.0053690
×4370.0095581.87 × 10 8 5120.0266996.79 × 10 9 11230.0224547.9 × 10 7
×5490.0121121.03 × 10 11 27550.0832286.22 × 10 7 13270.0234427.63 × 10 7
×6490.0112631.07 × 10 11 21430.0394028.19 × 10 7 13270.0239767.33 × 10 7
50,000×1250.0219110130.011120130.0095980
×2250.0193172.22 × 10 16 130.00991209190.0565677.45 × 10 7
×3370.073120130.0233740250.0161310
×4370.0339541.88 × 10 8 5120.0494766.69 × 10 9 11230.0816337.92 × 10 7
×5490.0603392.3 × 10 11 130.014165013280.123748.53 × 10 7
×6490.0443743.23 × 10 11 130.028178013280.108538.87 × 10 7
100,000×1250.0440820130.0373950130.0163320
×2250.0499412.22 × 10 16 130.04027409190.113637.45 × 10 7
×3370.0727290130.0836910250.045880
×4370.0710521.88 × 10 8 5120.0923816.68 × 10 9 11230.15117.92 × 10 7
×5490.0796483.25 × 10 11 130.045013014290.247086.03 × 10 7
×6490.146163.22 × 10 11 130.024994014290.183295.74 × 10 7
Table 5. Numerical results of DFSR1, PDY and HCGP on Problem 4.
Table 5. Numerical results of DFSR1, PDY and HCGP on Problem 4.
Problem 4 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×118380.0071328.9 × 10 7 6140.0372132.09 × 10 7 8170.005189.7 × 10 7
×213280.0058948.43 × 10 8 34700.0246626.1 × 10 7 11240.0101932.47 × 10 7
×320420.0077194.95 × 10 7 6140.0063266.55 × 10 7 6130.0029038.28 × 10 9
×410220.0040688.79 × 10 7 33680.0612775.47 × 10 7 10220.0074824.59 × 10 7
×58180.0043276.23 × 10 8 48980.0576618.73 × 10 7 18380.0090346.91 × 10 7
×610220.0053214.02 × 10 7 491000.0787925.53 × 10 7 13270.0065224.87 × 10 9
5000×119400.0422867.94 × 10 7 6140.0278034.68 × 10 7 8180.0249651.8 × 10 7
×215320.0384776.6 × 10 7 35720.137515.59 × 10 7 11240.0231379.68 × 10 8
×321440.051944.42 × 10 7 7160.0383399.35 × 10 8 6130.0133921.85 × 10 8
×412260.0247261.1 × 10 7 36740.132276.96 × 10 7 14300.0340953.04 × 10 7
×58180.0253791.26 × 10 7 34700.15597.37 × 10 7 13270.0178193.58 × 10 7
×68180.017723.49 × 10 8 43880.279227.08 × 10 7 15320.0330456.32 × 10 8
10,000×120420.0960944.48 × 10 7 6140.0603416.62 × 10 7 8180.0316532.55 × 10 7
×215320.0695291 × 10 7 16340.127988.16 × 10 7 11240.0537424.33 × 10 7
×321440.151776.25 × 10 7 7160.0593521.36 × 10 7 6130.0254022.62 × 10 8
×414290.0662883.87 × 10 7 35720.243235.62 × 10 7 15320.0661421.92 × 10 7
×58180.0470451.76 × 10 7 40820.298596.2 × 10 7 13270.059615.05 × 10 7
×69190.0539133.79 × 10 7 48980.397738.52 × 10 7 15310.0994511.68 × 10 7
50,000×121440.366374 × 10 7 7160.183049.46 × 10 8 8180.102055.7 × 10 7
×217350.27328.95 × 10 7 28580.708377.57 × 10 7 19390.238081.91 × 10 8
×322460.414115.57 × 10 7 8180.197572.69 × 10 7 6130.0985795.85 × 10 8
×412260.250778.19 × 10 8 26540.595321.44 × 10 7 15320.268541.64 × 10 7
×58180.138343.89 × 10 7 37760.932977.23 × 10 7 13280.18787.17 × 10 8
×69190.154336.75 × 10 7 48981.34366.36 × 10 7 16340.312173.47 × 10 7
100,000×121440.621115.66 × 10 7 7160.214268.58 × 10 7 8180.196488.06 × 10 7
×217360.604341 × 10 7 35721.03267.61 × 10 7 12260.310845.07 × 10 8
×322460.61027.88 × 10 7 8180.347643.81 × 10 7 6130.12478.28 × 10 8
×414300.429579.3 × 10 8 29600.873755.25 × 10 7 15310.324661.28 × 10 9
×58180.263035.49 × 10 7 35721.20434.43 × 10 7 13280.405341.01 × 10 7
×68180.252331.83 × 10 7 41841.48546.97 × 10 7 13280.412839.35 × 10 8
Table 6. Numerical results of DFSR1, PDY and HCGP on Problem 5.
Table 6. Numerical results of DFSR1, PDY and HCGP on Problem 5.
Problem 5 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×1130.002308016340.0190278.17 × 10 7 6130.0048663.87 × 10 7
×2130.0022962.22 × 10 16 15310.0114577.97 × 10 7 130.0028893.14 × 10 16
×3130.001798018380.0116058.51 × 10 7 130.0015490
×46130.003828017360.0121134.98 × 10 7 9200.0088595.86 × 10 7
×5611230.0757398.25 × 10 7 19400.0327444.77 × 10 7 14300.0148312.7 × 10 7
×6601210.0265419.82 × 10 7 19400.0428734.22 × 10 7 15310.0096474.8 × 10 7
5000×1130.007539017360.128496.85 × 10 7 6130.0230568.66 × 10 7
×2130.0056042.22 × 10 16 15310.0760057.97 × 10 7 130.0047733.14 × 10 16
×3130.005712020420.0584998.59 × 10 7 130.0045120
×46130.010185017360.0484964.99 × 10 7 7150.0154120
×5641290.195848.14 × 10 7 20420.0584354.02 × 10 7 15310.0376371.23 × 10 7
×6641290.150527.9 × 10 7 20420.0517013.98 × 10 7 16330.0403467.64 × 10 7
10,000×1130.020364017360.181649.69 × 10 7 6140.047934.59 × 10 7
×2130.0082062.22 × 10 16 15310.0810927.97 × 10 7 130.0064983.14 × 10 16
×3130.014966021440.289844.56 × 10 7 130.0051540
×46130.033961017360.19584.99 × 10 7 11230.0476134.06 × 10 7
×5651310.406178.75 × 10 7 20420.0978095.69 × 10 7 15310.0685843.8 × 10 7
×6651310.369168.77 × 10 7 20420.095495.5 × 10 7 14290.0586681.06 × 10 7
50,000×1130.030717018380.475868.13 × 10 7 7150.167871.71 × 10 7
×2130.0248482.22 × 10 16 15310.539987.97 × 10 7 130.0176533.14 × 10 16
×3130.039588023480.815559.93 × 10 7 130.0167970
×46130.11358017360.349064.99 × 10 7 11230.17457.52 × 10 7
×5681371.57118.55 × 10 7 20420.468787.76 × 10 7 15310.259827.56 × 10 7
×6681371.44548.59 × 10 7 20420.364117.74 × 10 7 15310.246612.61 × 10 7
100,000×1130.052007019401.2224.31 × 10 7 7150.193572.42 × 10 7
×2130.0466852.22 × 10 16 15310.506667.97 × 10 7 130.0390553.14 × 10 16
×3130.066953027561.78944.72 × 10 7 130.0317850
×46130.29814017360.564454.99 × 10 7 11230.383847.66 × 10 7
×5691392.46129.17 × 10 7 21440.961694.11 × 10 7 15320.625533.75 × 10 7
×6691392.46419.1 × 10 7 21441.26454.12 × 10 7 15310.493388.85 × 10 8
Table 7. Numerical results of DFSR1, PDY and HCGP on Problem 6.
Table 7. Numerical results of DFSR1, PDY and HCGP on Problem 6.
Problem 6 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×144890.0301537.26 × 10 7 1392800.0814339.73 × 10 7 661340.029077.63 × 10 7
×233680.0245837.28 × 10 7 1873760.0872438.84 × 10 7 36740.0155636.67 × 10 7
×347960.0332698.6 × 10 7 2274560.106689.23 × 10 7 551120.0212336.9 × 10 7
×443880.028238.87 × 10 7 2184380.0993539.53 × 10 7 42860.0416438.78 × 10 7
×5521060.0320097.35 × 10 7 2334680.247669.99 × 10 7 43880.0230546.11 × 10 7
×6651320.101198.51 × 10 7 1893800.127499.71 × 10 7 581180.0318027.64 × 10 7
5000×147960.181026.5 × 10 7 1883780.708688.86 × 10 7 811640.172598.38 × 10 7
×233680.116447.28 × 10 7 1873760.665938.84 × 10 7 36740.0787336.67 × 10 7
×347960.173749.37 × 10 7 2775561.06549.68 × 10 7 581180.166636.21 × 10 7
×439800.151589.6 × 10 7 2184380.776479.56 × 10 7 42860.11278.8 × 10 7
×5551120.329667.8 × 10 7 1813640.761639.64 × 10 7 911840.27958.56 × 10 7
×6721460.275989.32 × 10 7 2004020.805189.59 × 10 7 621260.156635.8 × 10 7
10,000×147960.403118.96 × 10 7 1933881.51519.37 × 10 7 821660.441997.28 × 10 7
×233680.253237.28 × 10 7 1873761.31138.84 × 10 7 36740.215376.67 × 10 7
×3541100.403197.71 × 10 7 2204421.64569.39 × 10 7 45920.233669.4 × 10 7
×440820.294238.61 × 10 7 2184381.15529.56 × 10 7 42860.227998.8 × 10 7
×5651320.567449.45 × 10 7 2324662.00159.66 × 10 7 992000.535456.52 × 10 7
×6611240.465677.27 × 10 7 2044101.45899.71 × 10 7 621260.333447.97 × 10 7
50,000×135711.0839.99 × 10 7 1903824.99069.44 × 10 7 941891.81189.9 × 10 7
×233681.00537.28 × 10 7 1873764.7178.84 × 10 7 36740.670276.67 × 10 7
×333680.823187.4 × 10 7 2334685.4539.9 × 10 7 46940.952858.25 × 10 7
×438781.0458.75 × 10 7 2184384.82889.56 × 10 7 42860.781678.8 × 10 7
×5661341.95216.95 × 10 7 2384785.84099.93 × 10 7 1052122.12046.84 × 10 7
×6751522.0318.46 × 10 7 2084185.30679.27 × 10 7 661341.23615.7 × 10 7
100,000×146942.21288.51 × 10 7 1973968.33959.98 × 10 7 921863.34876.48 × 10 7
×233681.63537.28 × 10 7 1873767.42528.84 × 10 7 36741.29836.67 × 10 7
×3531082.24679.03 × 10 7 1843707.59879.73 × 10 7 46941.52856.79 × 10 7
×435721.42498.8 × 10 7 2184389.15829.56 × 10 7 42861.53658.8 × 10 7
×5591202.64047.84 × 10 7 25751611.16289.96 × 10 7 1302625.25846.39 × 10 7
×6841703.62047.37 × 10 7 2064149.65129.17 × 10 7 661342.13158.17 × 10 7
Table 8. Numerical results of DFSR1, PDY and HCGP on Problem 7.
Table 8. Numerical results of DFSR1, PDY and HCGP on Problem 7.
Problem 7 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×1611230.0236739.91 × 10 7 2254520.123589.68 × 10 7 1543100.0637659.18 × 10 7
×2501020.0214147.3 × 10 7 2054120.118539.84 × 10 7 1493000.0409998.68 × 10 7
×3691400.0265747.38 × 10 7 1523060.103569.82 × 10 7 1042100.0284259.8 × 10 7
×4541100.0213719.54 × 10 7 2264540.109599.59 × 10 7 1452920.0404168.83 × 10 7
×545920.0175019.38 × 10 7 2044100.0879279.77 × 10 7 1482980.0420459.49 × 10 7
×6731480.0301568.61 × 10 7 3256520.163619.68 × 10 7 2384780.0807699.94 × 10 7
5000×1571160.174969 × 10 7 2084180.376949.98 × 10 7 1513040.273629 × 10 7
×2651320.142689.6 × 10 7 2084180.448979.66 × 10 7 1493000.223918.85 × 10 7
×3591200.126429.47 × 10 7 1573160.335359.77 × 10 7 1082180.17438.38 × 10 7
×4641300.15618.9 × 10 7 2024060.445259.74 × 10 7 1452920.228779.89 × 10 7
×5661340.149888.6 × 10 7 2054120.370699.48 × 10 7 1452920.399899.42 × 10 7
×6881780.19018.22 × 10 7 3456920.92179.68 × 10 7 2605220.412059.17 × 10 7
10,000×1511030.314615.68 × 10 7 2194401.21469.71 × 10 7 1503020.616389.48 × 10 7
×2551120.291219.62 × 10 7 2004021.15839.98 × 10 7 1452920.636479.38 × 10 7
×3761540.410329.63 × 10 7 1482980.789069.58 × 10 7 1082180.545398.86 × 10 7
×4661340.350869.07 × 10 7 2224461.25429.6 × 10 7 1422860.636589.31 × 10 7
×5561140.309777.82 × 10 7 1973961.13539.86 × 10 7 1452920.644569.89 × 10 7
×6861740.50278.36 × 10 7 3486982.00859.91 × 10 7 2685381.22428.22 × 10 7
50,000×1571161.12889.51 × 10 7 2044103.81619.76 × 10 7 1472962.19169.65 × 10 7
×2551121.08169.99 × 10 7 2014043.62549.87 × 10 7 1452922.12269.5 × 10 7
×3741501.46557.6 × 10 7 1523062.61669.9 × 10 7 1042101.55088.79 × 10 7
×4721461.41349.83 × 10 7 1983982.96679.96 × 10 7 1432882.04979.18 × 10 7
×5681381.35418.71 × 10 7 1963942.94599.62 × 10 7 1422861.82129.97 × 10 7
×6871762.05577.86 × 10 7 3647305.34539.81 × 10 7 2795603.37569.43 × 10 7
100,000×1651322.88449.84 × 10 7 2104226.99919.9 × 10 7 1472963.69069.98 × 10 7
×2651322.67347.82 × 10 7 2154326.91321 × 10 6 1462943.70739.9 × 10 7
×3701422.62978.73 × 10 7 1553125.23519.74 × 10 7 1042102.72559.15 × 10 7
×4611242.21498.42 × 10 7 2114247.07199.6 × 10 7 1422863.58599.94 × 10 7
×5691402.59338 × 10 7 2124266.87229.92 × 10 7 1462943.67648.11 × 10 7
×61022063.95989.32 × 10 7 37475011.95199.55 × 10 7 2885787.24118.77 × 10 7
Table 9. Numerical results of DFSR1, PDY and HCGP on Problem 8.
Table 9. Numerical results of DFSR1, PDY and HCGP on Problem 8.
Problem 8 DFSR1PDYHCGP
DimensionSPITERFVALTIMENORMITERFVALTIMENORMITERFVALTIMENORM
1000×118380.0079589.14 × 10 7 5120.0357156.7 × 10 7 6130.0026377.39 × 10 8
×218370.0074234.41 × 10 7 561140.0226964.91 × 10 7 18380.006692.4 × 10 7
×319400.0120795.7 × 10 7 7160.0061453.73 × 10 7 7150.0025883.49 × 10 8
×416330.0125157.78 × 10 7 531080.0334167.74 × 10 7 18380.0065712.09 × 10 7
×526530.0119848.23 × 10 7 551120.0514539.95 × 10 7 20420.0096473.18 × 10 7
×618370.0083187.92 × 10 7 651320.0388247.53 × 10 7 20420.0073835.83 × 10 7
5000×119400.0262688.15 × 10 7 6140.0170259.57 × 10 8 6130.0164361.65 × 10 7
×218380.0425777.68 × 10 7 44900.104852.79 × 10 7 17360.0391782.35 × 10 7
×320420.0908955.09 × 10 7 7160.0132398.34 × 10 7 7150.0142977.81 × 10 8
×4491000.151769.77 × 10 7 43880.160419.7 × 10 7 20420.0501962.3 × 10 7
×518380.0424965.14 × 10 7 45920.122038.5 × 10 7 23470.0492535.73 × 10 7
×619400.0484089.02 × 10 7 771560.27657.67 × 10 7 21440.0509572.39 × 10 7
10,000×120420.107394.6 × 10 7 6140.0440561.35 × 10 7 6130.0216732.34 × 10 7
×234700.149549.02 × 10 7 511040.202593.85 × 10 7 19400.0774442.08 × 10 7
×320420.102817.2 × 10 7 8180.0363417.53 × 10 8 7150.0228221.1 × 10 7
×4601210.336728.6 × 10 7 37760.291445.71 × 10 7 20420.149966.11 × 10 7
×517360.0856978.71 × 10 7 44900.212189.11 × 10 7 17360.0713953.08 × 10 7
×620420.10035.94 × 10 7 871760.407997.95 × 10 7 22460.120342.88 × 10 7
50,000×121440.528544.11 × 10 7 6140.0959053.03 × 10 7 6130.0815785.22 × 10 7
×232660.758148.53 × 10 7 46941.11089.07 × 10 7 17360.303722.58 × 10 7
×321440.399516.42 × 10 7 10220.216881.43 × 10 7 7150.0930752.47 × 10 7
×424500.454536.22 × 10 7 561141.17699.66 × 10 7 21440.343567.02 × 10 7
×518380.368468.95 × 10 7 46940.843677.11 × 10 7 21440.423532.59 × 10 7
×621440.499059.61 × 10 7 961941.83686.46 × 10 7 21440.381955.7 × 10 7
100,000×121440.814375.81 × 10 7 7160.229413.45 × 10 7 6130.144667.39 × 10 7
×234701.38999.14 × 10 7 531082.19768.8 × 10 7 19400.661132.42 × 10 7
×321440.753369.08 × 10 7 12260.466818.81 × 10 7 7150.208513.49 × 10 7
×418380.689364.98 × 10 7 581182.4075.95 × 10 7 19400.675281.64 × 10 7
×518370.59298.63 × 10 7 45921.64468.52 × 10 7 22460.647098.69 × 10 8
×622460.777847.98 × 10 7 1082183.84863.84 × 10 7 23470.890786.62 × 10 7
Table 10. Test results for DFSR1 and CGD the colored image deblurring experiment.
Table 10. Test results for DFSR1 and CGD the colored image deblurring experiment.
DFSR1CGD
ImageSizeITERTIME(s)SNRSSIMITERTIME(s)SNRSSIM
1st Image256 × 256431.8426.870.91441.9525.550.89
2nd Image256 × 256523.5319.230.81523.5217.940.74
3rd Image256 × 256312.6119.820.86353.0518.450.81
4th Image256 × 256201.7725.420.83232.0124.750.81
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Awwal, A.M.; Ishaku, A.; Halilu, A.S.; Stanimirović, P.S.; Pakkaranang, N.; Panyanak, B. Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery. Symmetry 2022, 14, 2375. https://doi.org/10.3390/sym14112375

AMA Style

Awwal AM, Ishaku A, Halilu AS, Stanimirović PS, Pakkaranang N, Panyanak B. Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery. Symmetry. 2022; 14(11):2375. https://doi.org/10.3390/sym14112375

Chicago/Turabian Style

Awwal, Aliyu Muhammed, Adamu Ishaku, Abubakar Sani Halilu, Predrag S. Stanimirović, Nuttapol Pakkaranang, and Bancha Panyanak. 2022. "Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery" Symmetry 14, no. 11: 2375. https://doi.org/10.3390/sym14112375

APA Style

Awwal, A. M., Ishaku, A., Halilu, A. S., Stanimirović, P. S., Pakkaranang, N., & Panyanak, B. (2022). Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery. Symmetry, 14(11), 2375. https://doi.org/10.3390/sym14112375

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop