An Enhanced Matrix-Free Secant Method via Predictor-Corrector Modified Line Search Strategies for Solving Systems of Nonlinear Equations

Diagonal updating scheme is among the cheapest Newton-like methods for solving system of nonlinear equations. Nevertheless, the method has some shortcomings. In this paper, we proposed an improved matrix-free secant updating scheme via line search strategies, by using the steps of backtracking in the Armijo-type line search as a step length predictor and Wolfe-Like condition as corrector. Our approach aims at improving the overall performance of diagonal secant updating scheme. Under mild assumptions, the global convergence results have been presented. Numerical experiments verify that the proposed approach is very promising.

The famous scheme for finding the solution to (1) is the Newton method.The method is simple to implement, and it produces an iterative sequence {  } from any given initial guess  0 in the neighborhood of  * via  +1 =   − (  (  )) −1  (  ) , where  = 0, 1, 2 . . ..
The attractive features of this method are that it is easy to implement and converges rapidly [1].However, the major difficulty of Newton-type method is matrix storage requirements especially when handling large systems of nonlinear equations [1].To overcome such difficulty, the simple modification on the Newton method is the fixed Newton method.Fixed Newton method for the determination of solution  * is given by  +1 =   − (  ( 0 )) −1  (  ) ,  = 0, 1, 2, . . . .
This method avoids computation and storing the Jacobian in each iterations (except at  = 0).However, it still requires to solve systems of  linear equations and consumes more CPU time as the system's dimension increases [2,3].Quasi-Newton method is another variant of Newtontype methods.It replaces the Jacobian or its inverse with an approximation which can be updated at each iteration and is given as where the matrix   is the approximation of the Jacobian at   .The main idea behind quasi-Newton's method is to International Journal of Mathematics and Mathematical Sciences eliminate the evaluation cost of the Jacobian matrix [4,5].However, the most critical part of such quasi-Newton method is storing the full matrix of the approximate Jacobian, which can be a very expensive task as the dimension of systems increases [6].In contrast, this paper presents an improved variant of quasi-Newton update via the steps of backtracking in the Armijo-type line search and Wolfe-like condition and incorporating restarting strategy whenever the updating matrix is singular or nearly singular.The anticipation has been to improve the overall performance of the diagonal quasi-Newton method.
We organized the rest of this paper as follows.In the next section, we present the proposed method.Section 3 presents convergence results.Numerical experiments are reported in Section 4, and finally, conclusion is given in Section 5.

The Improved Secant Diagonal Updating
It is well known that it is not always feasible to compute the full elements of the Jacobian matrix of the given nonlinear function or it may be very expensive.We often have to approximate the Jacobian matrix by some other approach, and the famous method of doing so is quasi-Newton's method [4].The basic idea underlining this approach has been to reduce the evaluation cost of the Jacobian matrix.This new approach generates a sequence of points {  } via where   is a step length, and   is a diagonal approximation of the inverse Jacobian matrix which can be updated in each iterations.To achieve this, we incorporate some new line search strategies, via the steps of backtracking in the Armijo-type line search as predictor and then improved via Wolfe-like condition as a corrector.The Armijo rule is among the inexact line search methods which guarantees a sufficient degree of accuracy to ensure the algorithm convergence.Nevertheless, the scheme needs high floating points operations and function call.In this paper, we present a simple line search strategy, which is less computational cost, floating points operations and CPU time consumptions respectively compared to classical Armijor line search.The approach is as follows.
In addition, the new strategy is implemented in an iterative way using a fixed initial value of  as follows.
Waziri et al. [7] have set the step length (  = 1, for all ); this approach is mostly used in many Newton-like methods.
Here, we continue in the spirit of diagonal updating, using a new line search strategy to obtain a good step length (  ) in every iteration, anticipating to produce a more accurate approximation of the Jacobian inverse matrix and then employing restating strategy whenever the updating matrix is undefined.To this end,  +1 would be obtained almost similar to the diagonal updating scheme presented in [7] in which   = −    (  ) instead of   = −  (  ).Now, let the deviation between  +1 and   denoted as Φ  =  +1 −   be minimized under some norms; the optimal solution is given as where and tr is the trace operation, respectively.The updated formula for the proposed matrix   is given as [7] To safeguard on the possibilities of generating undefined  +1 , we let  +1 =  0 , whenever this situation happens Now, we can describe the algorithm for our proposed approach as follows.

Convergence Analysis
We present the convergence result of EMFM method by proving the existence of the step length   > 0. We will make the following assumptions on nonlinear system .
(ii) There exists  * ∈  such that ( * ) = 0, and   () is continuous for all .(iii) () satisfies Lipschitz condition of order one; that is, there exists a positive constant  such that for all ,  ∈ R  .
(iv) There exist constants To this end, we proceed by given the following result for the step length generated by the proposed strategies.

Numerical Results
In this section, we consider some benchmark problems to illustrate the performance of the method proposed in this paper for solving large-scale systems of nonlinear equations when compared to some Newton-like methods.The computations are performed in MATLAB 7.0 using double precision computer, and the stopping rule used is The identity matrix has been chosen as an initial approximate Jacobian inverse.We further design the codes to terminate whenever one of the following happens: (i) the number of iteration is at least 250, but no point of   that satisfies (34) is obtained; (ii) CPU time in seconds reaches 250; (iii) There is insufficient memory to initiate the run.
The performances of these methods are be compared in terms of number of iterations and CPU time in seconds.In the following, some details on the benchmarks test problems are presented.
Problem 1. System of  nonlinear equations is as follows: Problem 2. Extended Trigonometric function of Spedicator [8] is as follows: ( Problem 3. System of  nonlinear equations is as follows: (A1) Nonadiabatic Stirred Tank Reactors.The first application is the model of [9].The model deals with two continuous nonadiabatic stirred tank reactors.The reactors are in a series, at steady state having a recycle component, and with an exothermic first-order irreversible reaction.By eliminating certain variables, the model results into two nonlinear systems of equations as follows: ) . ( The dimensionless temperatures of the two reactors are represented by  1 and  2 .The parameters , ,  1 ,  2 , and  are given as 1, 1000, 2, 2, and 22, respectively. (A2) Navigation by Range Measurements.Consider two beacons determining position by measuring the distances 1 and 2, where (1, 2) is the position of beacon 1, (1, 2) is the position of beacon 2, and (, ) is an unknown point (see Figure 1).This can be modeled into a two nonlinear systems of equations as follows: The parameters ( 1 ,  1 ), ( 2 ,  2 ),  1 , and  2 are given as (10, 10), (10, −10), 14, and 16, respectively.
Table 1 shows the number of iterations and CPU time for these five methods, respectively.In Table 1, the value "Dim" denotes the dimension of the systems, "NI" denotes number of iterations, while "CPU" is the CPU time in seconds respectively.We analyze the performance of each method via execution time, floating points operations, and storage locations.One can observe that EMFM has the smallest number of iterations compared to the classical diagonal updating (I-VDN) proposed by Waziri et al. [7].This shows that the line search strategies presented in this paper have increased the convergence speed of the classical diagonal updating method.
If we compare the performance of all methods, in terms of CPU time, it is clear that EMFM method consumes less CPU time than the others and still keeping memory requirement and CPU time in seconds to only ().All five methods are able to obtain the solution ( * = −1.1243,1.5001) of A2, but EMFM method consumes less CPU time in second (0.001) compared to the other 4 methods.Moreover, for A1, still proposed method has shown a promising performance with less storage locations (2 locations) whereas NM, FN, and BM, respectively, required 4 locations for each.

Conclusion
From the fact that there is a rapid development in research on solving nonlinear systems, nevertheless, the dimension of the nonlinear system is most of the times so large that it requires a very costly arithmetic operations when using some other Newton-like methods; so we usually employ cheap iterative approach, and the good candidate is classical diagonal updating.Notwithstanding, the classical updating scheme has some lapses.One is that it usually needs high number of iterations.The other is that the method generally has slow convergence behavior.This paper proposes an enhanced diagonal secant updating scheme based on the steps of backtracking in the Armijo-type line search and then improved via Wolfe-like condition.Our approach aims at improving the overall performance of diagonal secant updating scheme.The algorithm only requires to store a row vector while ignoring all the off and low diagonal elements and therefore largely reduces memory locations.In addition, as it uses two line search strategies (predictor and corrector) to obtain a new iterates point, the spectral properties of the diagonal updating scheme is improved, and rapid convergence property is gained.
Computational experiment suggests that it is very vital for diagonal updating scheme to use line search strategy.EMFM method has very good solving speed and the best performance among the Newton-like methods.Finally, it can be concluded that this approach would certainly be quite useful for solving large-scale systems of nonlinear equations.

Theorem 4 .
Assume that  is a strictly convex function.Suppose that the new strategies are employed with   ̸ = 0 and positive   exist, for all sufficiently large .Then, the iterates {  } generated by the line search algorithm have the property that   exist and are positive.We continue to show the convergence of the iterates {  } by recalling that      (  +     )     ≤       (  )     .Let {  } be the sequence generated by (8).