Next Article in Journal
On Nonlinear Fractional Difference Equation with Delay and Impulses
Next Article in Special Issue
Direct Comparison between Two Third Convergence Order Schemes for Solving Equations
Previous Article in Journal
Gas of Baby Universes in JT Gravity and Matrix Models
Previous Article in Special Issue
Fast Convergence Methods for Hyperbolic Systems of Balance Laws with Riemann Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems

1
Department of Mathematics, Cameron University, Lawton, OK 73505, USA
2
Department of Theory of Optimal Processes, Ivan Franko National University of Lviv, Universitetska Str. 1, 79000 Lviv, Ukraine
3
Department of Computational Mathematics, Ivan Franko National University of Lviv, Universitetska Str. 1, 79000 Lviv, Ukraine
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(6), 978; https://doi.org/10.3390/sym12060978
Submission received: 25 April 2020 / Revised: 3 June 2020 / Accepted: 3 June 2020 / Published: 8 June 2020
(This article belongs to the Special Issue Iterative Numerical Functional Analysis with Applications)

Abstract

:
A vital role in the dynamics of physical systems is played by symmetries. In fact, these studies require the solution for systems of equations on abstract spaces including on the finite-dimensional Euclidean, Hilbert, or Banach spaces. Methods of iterative nature are commonly used to determinate the solution. In this article, such methods of higher convergence order are studied. In particular, we develop a two-step iterative method to solve large scale systems that does not require finding an inverse operator. Instead of the operator’s inverting, it uses a two-step Schultz approximation. The convergence is investigated using Lipschitz condition on the first-order derivatives. The cubic order of convergence is established and the results of the numerical experiment are given to determine the real benefits of the proposed method.

Graphical Abstract

1. Introduction

In this article, F denotes an operator acting on a Banach space X with values inside a Banach space Y and D X is convex. To find the approximate solution p of
F ( x ) = 0
we often utilize the method [1,2,3,4,5]
x k + 1 = x k F ( x k ) 1 F ( x k ) , k 0 ,
due to Newton or its differential or difference modifications. Iterative methods are used since p can be found in closed form only in special cases. That is the first benefit of using iterative methods. However, at each iteration, it is necessary to find one or more inverse operators. This is one setback for using Newton-type methods [1,2,3,4,5]. Since it is not always easy to do, we bypass this obstacle by finding approximations for inversion using only multiplications of linear operators. This is one major benefit. There is the novelty of our paper, since we contribute in this direction.
Methods with approximation of inverse operator are based on different ideas. Some ideas intend to construct approximations to the solution of a nonlinear Equation (1), and others to construct approximations to the inverse operator. There are two approaches to approximate an inverse operator: successive approximation (SA) and parallel (synchronous and asynchronous) approximations. In methods with the successive approximation, the calculations in separate branches are performed alternately. In methods with the parallel approximation requiring the computation of an inverse of a linear operator, the computations in separate branches of the method are performed in parallel. Such methods are effective for numerically solving Equation (1) in the parallel processor system with common memory [6,7,8].
Many authors have investigated methods with SA of inverse operator. For example, the local convergence of modifications of Newton and Steffensen methods are studied in [9]. In [10,11,12], the authors studied the semilocal convergence (SLC) for Ulm method [9]
A k + 1 = A k + ( E A k F ( x k ) ) A k , x k + 1 = x k A k + 1 F ( x k ) , k 0
and its difference analog
A k + 1 = A k + ( E A k F ( x k 1 , x k ) ) A k , x k + 1 = x k A k + 1 F ( x k ) , k 0 .
Here, E is the identity operator in X and F ( u , v ) stands for the first-order divided difference for operator F ( x ) [13]. Moreover, x 0 and A 0 are initial approximations for a solution x * and an inverse operator A * = [ F ( x * ) ] 1 , respectively. Furthermore, the closer x 0 and A 0 are to p, A * the better Equation (4) performs. However, there are no preconditioners. An investigation of the accelerated Newton method and a two-parametric secant-type method with SA of the inverse operator is performed in [14,15] (see also [16,17,18]). It is worth noting that Ulm type methods provide the same order of convergence as Newton’s or Steffensen’s but without using inverses which are very expensive to compute in general. In certain cases, even the ratio of convergence may be smaller (see the numerical section).
A method with synchronous inverse approximation
A k + 1 = A k + ( E A k F ( x k ) ) A k , x k + 1 = x k A k F ( x k ) , k 0
and a method with the asynchronous approximation were considered by A. Rooze [6,7]. Some modifications of a method by Gauss–Newton with SA of the inverse operator to solve a nonlinear least squares problems were proposed by Iakymchuk, Shakhno [8].
In 1983, the authors built third-, fifth-, and sixth-order methods with approximation of inverse operator for solving operator Equations [19]. In particular, the third-order method has the form
y k = x k A k F ( x k ) , x k + 1 = y k A k F ( y k ) , B k = A k ( 2 E F ( x k + 1 , u k + 1 ) A k ) , A k + 1 = B k ( 2 E F ( x k + 1 , u k + 1 ) B k ) , k 0 ,
where u k + 1 = x k + 1 β k + 1 F ( x k + 1 ) , β k + 1 is a real parameter.
To increase the order of convergence, efficiency and applicability of the aforementioned Ulm-type method, by replacing the divided differences by Fréchet derivatives in Equation (6), we develop the method
y k = x k A k F ( x k ) , x k + 1 = y k A k F ( y k ) , B k = A k ( 2 E F ( x k + 1 ) A k ) , A k + 1 = B k ( 2 E F ( x k + 1 ) B k ) , k 0 .
This method is a two-step modification of Newton method [2,20]:
y k = x k F ( x k ) 1 F ( x k ) , x k + 1 = y k F ( x k ) 1 F ( y k ) , k 0 .
It is easy to see that the method in Equation (7) coincides with the third-order iterative process
y k = x k B k F ( x k ) , x k + 1 = y k B k F ( y k ) , B k + 1 = B k + B k ( 2 E F ( x k + 1 ) B k ) ( E F ( x k + 1 ) B k ) , k 0 ,
which was proposed by Esquerro, Hernández [21]. They only investigated the semilocal convergence (SLC) of the method in Equation (9). Obviously, the methods in Equations (7) and (9) can be considered as methods SA for inverse operator.
We examine conditions for local convergence (LC) of the method in Equation (7) including its convergence order. We note that LC results are very important, since they demonstrate the degree of difficulty in selecting starters x 0 that guarantee convergence to p.
Another novelty of the developed method in Equation (7) lies in the fact that the inverse at each step is not calculated compared to other methods; the method is suitable to solve large scale systems; and it has better rate and convergence order than other methods using related information (see also benefits reported in Section 3. Indeed, there are significant developments in the study of Ulm-type methods. Moreover, these ideas can be used to expand the applicability of other methods along the same lines.
Section 2 contains the LC study of method in Equation (7). Large scale systems are solved in Section 3, where X and Y are specialized to be finite-dimensional Euclidean spaces. The study of solving large scale systems is of extreme importance, since most problems from diverse disciplines reduce to solving such systems [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22].

2. LC study for the Method in Equation (7)

The following theorem establishes the conditions under which the iterative process in Equation (7) is convergent.
Theorem 1.
Assume that:
(1) 
Operator A * = [ F ( p ) ] 1 exists;
A * B
and
F ( p ) C ,
where p is the solution of Equation (1).
(2) 
In ball U = { x D : x p R 0 } , the following estimate is satisfied
F ( x ) L ,
where R 0 = max { r 0 , a 1 r 0 2 , a 2 r 0 3 } , r 0 = max { x 0 p , A 0 A * } ,
a 1 = C + 3 2 B L + 3 2 L r 0 ;
a 2 = ( C + 3 2 B L a 1 r 0 + 3 2 L a 1 r 0 2 ) a 1 .
(3) 
Initial approximations x 0 , A 0 are such that
q r 0 < 1 ,
where q = ( max { a 2 , a 3 } ) 1 2 , a 3 = C γ 2 r 0 + L a 2 ( C + γ r 0 2 ) 2 , γ = C + L ( B + r 0 ) 2 a 2 r 0 .
Then, sequences { x k } , { A k } , k 0 , generated by Equation (7), converge to p , A * , respectively, and
r k = max { A k A * , x k p } ( q r 0 ) 3 k 1 r 0 , k 0 .
Proof. 
The proof is performed by mathematical induction.
It follows from
x 0 p r 0 = ( q r 0 ) 3 0 1 r 0 ,
that x 0 U , and the estimate in Equation (14) is true for k = 0 . Suppose that x k U and the estimate in Equation (14) is satisfied for k 0 . It follows that r k r 0 , since q r 0 < 1 is provided by Equation (13). Taking into account Equation (10) and the definition r 0 , we get
A k A * + A k A * B + r 0 .
We obtain from the first equality of Equation (7) and the Taylor’s formula
p y k = p x k + A k ( F ( x k ) F ( p ) ) = p x k + A k F ( x k ) ( x k p ) A k 0 1 F ( x k + t ( p x k ) ) ( p x k ) 2 ( 1 t ) d t = A k 0 1 F ( x k + t ( p x k ) ) ( p x k ) 2 ( 1 t ) d t + ( E A k F ( x k ) ) ( p x k ) .
By the condition in Equation (11) and the estimate in Equation (15), we have since
E A k F ( x k ) = ( A * A k ) F ( p ) + A k ( F ( p ) F ( x k ) ) ,
that
E A k F ( x k ) F ( p ) A * A k + A k F ( p ) F ( x k ) C A * A k + ( B + r 0 ) L p x k r k ( C + L ( B + r 0 ) ) = a r k ,
where a = C + L ( B + r 0 ) .
It follows from Equations (16) and (17) that
p y k A k 0 1 F ( x k + t ( p x k ) ) ( p x k ) 2 ( 1 t ) d t + E A k F ( x k ) p x k ( B + r 0 ) L 2 p x k 2 + ( C + L ( B + r 0 ) ) r k p x k ( B + r 0 ) L 2 r k 2 + ( C + L ( B + r 0 ) ) r k 2 = ( C + 3 2 B L + 3 2 L r 0 ) r k 2 ,
thus
p y k a 1 r k 2 .
In addition, p y k a 1 r 0 2 R 0 , and whence it follows y k U .
We obtain from the second equality of Equation (7) and the Taylor’s formula
x k + 1 p = y k p A k ( F ( y k ) F ( p ) ) = A k 0 1 F ( y k + t ( p y k ) ) ( p y k ) 2 ( 1 t ) d t + ( E A k F ( y k ) ) ( y k p ) .
Hence, we have, given the conditions in Equations (10)–(12) and the estimates in Equations (15), (17), and (18),
x k + 1 p ( B + A * A k ) L 2 p y k 2 + C A * A k + ( B + A * A k ) L p y k p y k C A * A k p y k + 3 2 B L p y k 2 + 3 2 L p y k 2 A * A k p y k ( C r k + 3 2 B L p y k + 3 2 L r k p y k ) ( C + 3 2 B L a 1 r 0 + 3 2 L a 1 r 0 2 ) a 1 r k 3 = a 2 r k 3 .
Thus, x k + 1 U , since x k + 1 p a 2 r 0 3 R 0 .
On the other hand, based on the third formula of Equation (7)
A * B k = ( A * A k ) F ( p ) ( A * A k ) A k ( F ( p ) F ( x k + 1 ) ) A k .
Hence, we have, given Equations (11) and (20),
A * B k F ( p ) A * A k 2 + A k 2 F ( p ) F ( x k + 1 ) C r k 2 + A k 2 L x k + 1 p ( C + L ( B + r 0 ) 2 a 2 r 0 ) r k 2 = γ r k 2 .
In accordance with the fourth equality of Equation (7),
A * A k + 1 = B k ( F ( p ) F ( x k + 1 ) ) B k + ( A * B k ) F ( p ) ( A * B k ) .
From this relationships, based on the conditions in Equations (11) and (12) and estimates in Equations (20) and (21), we get in turn
A * A k + 1 F ( p ) A * B k 2 + B k 2 F ( p ) F ( x k + 1 ) C γ 2 r k 4 + B k 2 L x k + 1 p ( C γ 2 r 0 + L a 2 ( C + γ r 0 2 ) 2 ) r k 3 = a 3 r k 3 .
Given induction assumptions we obtain, based on the estimates in Equations (20) and (22),
r k + 1 = max { x k + 1 p , A k + 1 A * } = max { a 2 , a 3 } r k 3 = q 2 r k 3 = q 2 ( ( q r 0 ) 3 k 1 r 0 ) 3 = ( q r 0 ) 3 k + 1 1 r 0 .
That is, Equation (14) is fulfilled for an iteration k + 1 . The induction is complete.
Moreover, it follows from the estimate in Equation (14) for k the convergence of sequences { x k } and { A k } . ☐
Next, a uniqueness result follows.
Proposition 1.
Suppose: Equation (1) has a solution p D 0 , where D 0 D is convex and
A * ( F ( p ) F ( x ) ) L 0 p x for x D 0 and L 0 > 0 .
Set U 0 = { x D : x p < 1 L 0 } and S = D 0 U 0 .
Then, p uniquely solves (1) on S.
Proof. 
Consider for v S with F ( v ) = 0 , operator T = 0 1 F ( v + a ( p v ) ) d a . Using Equation (23), we obtain
A * ( T F ( p ) ) L 0 0 1 v + a ( p v ) p d a L 0 0 1 ( 1 a ) v p d a L 0 2 v p < 1 ,
thus T is invertible. Hence, p = v follows from
T ( p v ) = F ( p ) F ( v ) = 0 .
 ☐
The method in Equation (7) has a higher convergence rate than the Newton and Steffensen method. Moreover, in contrast, it does not require an inverse operator. Furthermore, the convergence order of method in Equation (7) is larger than for methods with the successive approximation in Equations (3) and (4).

3. Numerical Experiments

We used large scale test problems from Luksan [22] for the numerical study of the methods. Calculations were performed in Octave 5.1.0. Stopping iterative processes occurred under the condition
x k + 1 x k 10 10 .
The initial approximations were calculated by the rules x 0 = w ¯ , y 0 = x 0 + 10 4 . We compared methods by the number of iterations required to obtain an approximate solution. Table 1 shows the results for m = 100 .
Example 1.
Trigonometric-exponential function.
h i = 3 w i 3 + 2 w i + 1 5 + sin w i w i + 1 sin w i + w i + 1 , i = 1 , h i = 3 w i 3 + 2 w i + 1 5 + sin w i w i + 1 sin w i + w i + 1 + + 4 w i w i 1 exp w i 1 w i 3 , 1 < i < m , h i = 4 w i w i 1 exp w i 1 w i 3 , i = m , w ¯ i = 2 s , s I R , i = 1 , m ¯ .
Example 2.
Consider tridiagonal function due to Broyden.
h i = w i ( 0.5 w i 3 ) + 2 w i + 1 1 , i = 1 , h i = w i ( 0.5 w i 3 ) + w i 1 + 2 w i + 1 1 , 1 < i < m , h i = w i ( 0.5 w i 3 ) + w i 1 1 , i = m , w ¯ i = s , s I R , i = 1 , m ¯ .
Example 3.
Counter current reactor problem.
α = 1 / 2 , h i = α ( 1 α ) w i + 2 w i ( 1 + 4 w i + 1 ) , i = 1 , h i = ( 2 α ) w i + 2 w i ( 1 + 4 w i 1 ) , i = 2 , h i = α w i 2 ( 1 α ) w i + 2 w i ( 1 + 4 w i + 1 ) , mod ( i , 2 ) = 0 , 2 < i < m 1 , h i = α w i 2 ( 2 α ) w i + 2 w i ( 1 + 4 w i 1 ) , mod ( i , 2 ) = 1 , 2 < i < m 1 , h i = α w i 2 w i ( 1 + 4 w i + 1 ) , i = m 1 , h i = α w i 2 ( 2 α ) w i ( 1 + 4 w i 1 ) , i = m , w ¯ i = 1000 s , s I R , i = 1 , m ¯ .
The comparison by the number of iterations was performed to confirm the theoretical results about the higher convergence order of the studied method in Equation (7) than for the Ulm method in Equation (3). Obtained results also show that the methods with the approximation of inverse operator are somewhat inferior to the corresponding methods with the calculation of the inverse operator. However, the benefit of these methods is that they can be used when finding the inverse operator is impossible or difficult.
The graphs show the values of error x k p at each iteration (see Figure 1). These are results for s = 5 for Example 1, s = 100 for Example 2, and s = 0.1 for Example 3. We draw graphs using a logarithmic scale for both of the axes.
We see from the results on graphs that the errors of the methods with the approximation of inverse operator in Equations (3) and (7) decrease more slowly than for the basic methods in Equations (2) and (8), respectively. Moreover, obtained numerical results confirm that the methods in Equations (7) and (8) have a higher convergence rate than those in Equations (2) and (3).

4. Conclusions

The LC of method in Equation (7) with approximation of inverse operator for solving large scaled systems of equations is studied. The cubic order and radius of convergence of this method are determined. Numerical results are presented which demonstrate the computational benefits of the method.

Author Contributions

Conceptualization, S.S.; Investigation, I.K.A., S.S., and H.Y.; and Editing, I.K.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Argyros, I.K. Convergence and Applications of Newton-Type Iterations; Springer: New York, NY, USA, 2008; 506p. [Google Scholar]
  2. Amat, S.; Busquier, S.; Plaza, S. Review of some iterative root–finding methods from a dynamical point of view. Ser. A Math. Sci. 2004, 10, 3–35. [Google Scholar]
  3. Amat, S.; Busquier, S.; Plaza, S. Dynamics of King’s and Jarratt iterations. Aequ. Math. 2005, 69, 212–223. [Google Scholar] [CrossRef]
  4. Kantorovich, L.V.; Akilov, G.P. Functional Analysis; Nauka: Moscow, Russia, 1977; 744p. [Google Scholar]
  5. Magrenan, A.A. Different anomalies in a Jarratt family of root-finding methods. Appl. Math. Comput. 2014, 233, 29–38. [Google Scholar]
  6. Rooze, A.F. An iterative method for solving nonlinear equations using parallel inverse operator approximation. Izv. Acad. Nauk Est. SSR Fiz. Mat. 1982, 31, 32–37. (In Russian) [Google Scholar]
  7. Rooze, A.F. An asynchronous iteration method of solving non-linear equations using parallel approximation of an inverse operator. USSR Comput. Math. Math. Phys. 1986, 26, 188–191. (In Russian) [Google Scholar] [CrossRef]
  8. Iakymchuk, R.; Shakhno, S. Methods with successive and parallel approximations of inverse operator for the nonlinear least squares problem. Proc. Appl. Math. Mech. 2015, 15, 569–570. [Google Scholar] [CrossRef] [Green Version]
  9. Ulm, S. On iterative methods with successive approximation of the inverse operator. Izv. Acad. Nauk Est. SSR Fiz. Mat. 1967, 16, 403–411. (In Russian) [Google Scholar]
  10. Argyros, I.K. On Ulm’s method for Fréchet differentiable operators. J. Appl. Math. Comput. 2005, 31, 97–111. [Google Scholar] [CrossRef]
  11. Argyros, I.K. On Ulm’s method using divided differences of order one. Numer. Algorithms 2009, 52, 295–320. [Google Scholar] [CrossRef]
  12. Burmeister, W. Inversionfreie Verfahren zur Lösung nichtlinearer Operatorgleichungen. ZAMM 1972, 52, 101–110. [Google Scholar] [CrossRef]
  13. Ulm, S. On generalized divided differences. Izv. Acad. Nauk. Est. SSR Fiz. Mat. 1967, 16, 13–26. (In Russian) [Google Scholar]
  14. Shakhno, S.M. Convergence of one iterative method with successive approximation of inverse operator. Visnyk Lviv. Univ. Ser. Mech. Math. 1989, 31, 62–66. (In Ukrainian) [Google Scholar]
  15. Shakhno, S.M.; Yarmola, H.P. Two-step secant type method with approximation of the inverse operator. Carpathian Math. Publ. 2013, 5, 150–155. (In Ukrainian) [Google Scholar] [CrossRef] [Green Version]
  16. Vaarmann, O. On some iterative methods with successive approximation of the inverse operator. I. Izv. Acad. Nauk Est. SSR Fiz. Mat. 1968, 17, 379–387. (In Russian) [Google Scholar]
  17. Verzhbitskiy, V.M. Iterative methods with successive approximation of inverse operators. Izv. IMI UdGU 2006, 2, 129–138. (In Russian) [Google Scholar]
  18. Vaarmann, O. On some iterative methods with successive approximation of the inverse operator. II. Izv. Acad. Nauk Est. SSR Fiz. Mat. 1969, 18, 14–21. (In Russian) [Google Scholar]
  19. Vasilyev, F.P.; Khromova, L.M. On high-order methods for solving operator equations. Dokl. Nauk SSSR 1983, 270, 28–31. (In Russian) [Google Scholar]
  20. Esquerro, J.A.; Hernández, M.A.; Salanova, M.A. A discretization scheme for some conservative problems. J. Comp. Appl. Math. 2000, 115, 181–192. [Google Scholar] [CrossRef] [Green Version]
  21. Esquerro, J.A.; Hernández, M.A. An Ulm-type method with R-order of convergence three. Nonlinear Anal. Real World Appl. 2012, 12, 14–26. [Google Scholar] [CrossRef]
  22. Luksan, L. Inexact trust region method for large sparse systems of nonlinear equations. J. Optim. Theory Appl. 1994, 81, 569–590. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The error at each iteration.
Figure 1. The error at each iteration.
Symmetry 12 00978 g001
Table 1. The number of iterations for different values of s.
Table 1. The number of iterations for different values of s.
MethodExample 1Example 2Example 3
0.415151000.010.11
Equation (2)57115711151821
Equation (3)78157915222733
Equation (7)459459131518
Equation (8)458458111316

Share and Cite

MDPI and ACS Style

Argyros, I.K.; Shakhno, S.; Yarmola, H. Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems. Symmetry 2020, 12, 978. https://doi.org/10.3390/sym12060978

AMA Style

Argyros IK, Shakhno S, Yarmola H. Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems. Symmetry. 2020; 12(6):978. https://doi.org/10.3390/sym12060978

Chicago/Turabian Style

Argyros, Ioannis K., Stepan Shakhno, and Halyna Yarmola. 2020. "Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems" Symmetry 12, no. 6: 978. https://doi.org/10.3390/sym12060978

APA Style

Argyros, I. K., Shakhno, S., & Yarmola, H. (2020). Method of Third-Order Convergence with Approximation of Inverse Operator for Large Scale Systems. Symmetry, 12(6), 978. https://doi.org/10.3390/sym12060978

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