Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 283679 | https://doi.org/10.1155/2015/283679

Predrag S. Stanimirović, Gradimir V. Milovanović, Milena J. Petrović, Nataša Z. Kontrec, "A Transformation of Accelerated Double Step Size Method for Unconstrained Optimization", Mathematical Problems in Engineering, vol. 2015, Article ID 283679, 8 pages, 2015. https://doi.org/10.1155/2015/283679

A Transformation of Accelerated Double Step Size Method for Unconstrained Optimization

Academic Editor: Alejandro Ortega-Moñux
Received19 Jan 2015
Accepted24 Feb 2015
Published07 Apr 2015

Abstract

A reduction of the originally double step size iteration into the single step length scheme is derived under the proposed condition that relates two step lengths in the accelerated double step size gradient descent scheme. The proposed transformation is numerically tested. Obtained results confirm the substantial progress in comparison with the single step size accelerated gradient descent method defined in a classical way regarding all analyzed characteristics: number of iterations, CPU time, and number of function evaluations. Linear convergence of derived method has been proved.

1. Introduction and Background

The SM iteration from [1] is defined by the iterative processwhere is a new iterative point, is the previous iterative point, is the gradient vector (search direction), is a step length, and is the acceleration parameter. In [1] it is verified that the accelerated gradient SM iteration (1) outperforms the gradient descent, GD, as well as the Andrei’s accelerated gradient descent AGD method from [2].

Double direction and double step size accelerated methods, denoted by ADD and ADSS methods, respectively, for solving the problems of unconstrained optimization are presented in [3, 4]. These two algorithms can be generally formulated through the next merged expression: where is the previous iterative point and real values and denote two step lengths while the vectors and are two vector directions. The values of step lengths are determined by backtracking line search techniques. The gradient is basically used for defining a search direction, but some new suggestions for deriving a descending vector direction are given in [3, 5]. Taking the substitutionsinto (2) produces the ADD iterative scheme from [3]:where represents the acceleration parameter for the iteration (4). The benefits of the acceleration properties that arise from the usage of the parameter are explained in [3]. The so called nonaccelerated version of ADD method (NADD method shortly) is defined in order to numerically verify the acceleration property of the parameter . Three methods, SM, ADD, and NADD, are numerically compared in [3]. Results show the enormous efficiency of ADD scheme in comparison with its nonaccelerated counterpart NADD. Derivation of the direction vector is explained by the Algorithm 3.2 in [3]. The ADD outperforms its competitive SM method from [1] with respect to the number of iterations.

By replacing the vectors and from (2) by and , respectively, the next iteration is defined as The previous scheme is noted as ADSS model and it is proposed in [4]. In the same paper, a huge improvement in performances of this accelerated gradient descent method when compared to the accelerated gradient descent SM method from [1] is numerically confirmed.

The main contribution of the present paper is a transformation of the double step size iterative scheme (5) for unconstrained optimization into an appropriate accelerated single step size scheme (called TADSS shortly). Convergence properties of the introduced method are investigated. A special contribution is given by the numerical confirmation that the TADSS algorithm developed from the double step size ADSS model (5) is evidently more efficient than the accelerated SM method obtained in a classical way. Surprisingly, numerical experiments show that the TADSS method overcomes the initial ADSS method.

The paper is organized in the following way. The reduction of the double step size ADSS model into the single step size iteration TADSS and the presentation of defined accelerated gradient decent model are given in Section 2. Section 3 contains the convergence analysis of derived algorithm for uniformly convex and strictly convex quadratic functions. The results of numerical experiments as well as their comparative analysis of developed method and its forerunners are illustrated in Section 4.

2. Transformation of ADSS Scheme into a Single Step Size Iteration

Very advanced numerical results obtained in [4] motivated further research on this topic. An idea is to investigate the properties of a single step size method developed as a reduction of the double step size ADSS model. This reduction is defined by an additional assumption which represents a trade-off between two step length parameters and in the ADSS scheme:Taking into account assumption (6) into expression (5), which defines the ADSS iteration, leads to the iterative processThe iteration (7) is noted as transformed ADSS method, or shortly TADSS method. Defined TADSS iteration represents not only a reduction of the double step size ADSS model into the corresponding single step size method, but also a sort of modification of the single step size SM iteration from [1]. This modification can be explained as the substitution of the product , from the SM iteration (1), by the multiplying factor of the gradient from the TADSS iteration (7).

For the sake of simplicity, we use the notation whenever it is possible. The value of the acceleration parameter in th iteration can be derived by using Taylor’s expansion, similarly as described in [1, 3, 4]:The vector in (8) satisfiesFurther, it is reasonable to replace in (8) the Hessian by the diagonal matrix , where is an appropriately chosen real number. This replacement implies The relation (10) allows us to compute the acceleration parameter :

Next, the natural inequality is inevitable. This condition is required in order to fulfill second-order necessary condition and second-order sufficient condition. The choice is reasonable in the case when the inequality appears for some . This choice produces the next iterative point aswhich evidently represents the classical gradient descent step.

We consider now the th iteration, , which is given byExamine the function :defined as the finite part of the Taylor expansion of the functionunder the assumption . This function is convex when , and its derivative is calculated in the following way:

Since the inequality is achieved, the following is valid:Therefore, the function decreases in the case and achieves its own minimum in the case . According to the criteria given by (17), desirable values for are within the interval . Now, (7) is a kind of the gradient descent process in the case . Since , it is easy to verify the following condition for the step length :Since in the case , this fractional number is not appropriate upper bound for in this case. On the other hand, the inequality holds in the case , so that is an appropriate upper bound for in this case.

Further analysis of (10) in the case giveswhich impliesTaking into account , it is not difficult to verify that the criterion (20) restricts desirable values for within the interval

Final conclusion is that the upper bound for is defined in the case , , as the minimum between the upper bounds defined in (18) and (21):

According to the previous discussion, the iterative step is derived by the backtracking line search procedure presented in Algorithm 1.

Algorithm 1 (calculation of the step size by the backtracking line search which starts from the upper bound defined in (22)). Requirement: objective function , the direction of the search at the point , and real numbers and .(1)Set .(2)While take .(3)Return .

Finally, the TADSS algorithm of the defined accelerated gradient descent scheme (7) is presented.

Algorithm 2 (transformed accelerated double step size method (TADSS method)). Requirement: , , , .(1)Set , compute , , and take .(2)If , then go to Step 8; else continue by the next step.(3)Find the step size applying Algorithm 1.(4)Compute using (7).(5)Determine the scalar using (11).(6)If or , then take .(7)Set ; go to Step 2.(8)Return and .

3. Convergence of TADSS Scheme

The content of this section is the convergence analysis of the TADSS method. In the first part of this section a set of uniformly convex functions is considered. The proofs of the following statements can be found in [6, 7] and have been omitted.

Proposition 3 (see [6, 7]). If the function is twice continuously differentiable and uniformly convex on then (1)the function has a lower bound on , where is available;(2)the gradient is Lipschitz continuous in an open convex set which contains ; that is, there exists such that

Lemma 4. Under the assumptions of Proposition 3 there exist real numbers , satisfying such that has a unique minimizer and

The value of decreasing of analyzed function through each iteration is given by the next lemma which is restated and proven in [1]. The same estimation can similarly be found considering iteration (7). Theorem 6 is approved in [1] and confirms a linear convergence of the constructed method.

Lemma 5. For twice continuously differentiable and uniformly convex function on and for the sequence generated by Algorithm (7) the following inequality is valid:where

Theorem 6. If the objective function is twice continuously differentiable as well as uniformly convex on and the sequence is generated by Algorithm 2 thenand the sequence converges to at least linearly.

In the following review the case of strictly convex quadratic functions is analyzed. This set of functions is given byIn the previous expression is a real symmetric positive definite matrix and . It is assumed that the eigenvalues of the matrix are given and lined as . Since the convergence for the most gradient methods is quite difficult to analyze, in many research articles of this profile convergence analysis is reduced on the set of convex quadratics [810]. The convergence of TADSS method is also analyzed under similar presumptions.

Lemma 7. By applying the gradient descent method defined by (7) in which parameters and are given by relation (11) and Algorithm 1 on the strictly convex quadratic function expressed by relation (29) where presents a symmetric positive definite matrix, the next inequalities hold:where and are, respectively, the smallest and the largest eigenvalues of .

Proof. Considering expression (29), the difference between function value in the current and the previous point isApplying expression (7) the following is obtained: Using the facts that the gradient of the function (29) is in conjunction with the equality , one can verify the following: Substituting (33) into (11), the parameter becomesThe last relation confirms that is the Rayleigh quotient of the real symmetric matrix at the vector , so the next inequalities hold:which combined with the fact that prove the right hand side in (30):The estimation proved in [1], is considered in order to prove the left hand side of (30). Using the notation adopted in this paper, expression (37) becomesInequality (38) and the facts that , and lead to the following conclusion: In the last estimation, Lipschitz constant can be replaced by . The conclusion that the eigenvalue of matrix has the property of Lipschitz constant is to be derived from the next analysis. Since matrix is symmetric and the following inequality can be provided:Substituting by , inequality (39) becomes and this proves the left hand side of inequalities (30).

Remark 8. Comparing the estimations resulting from the similarly proposed lemma in [1, 3, 4] with the estimation derived from the previous lemma, considering the TADSS method, it can be concluded that the estimation provided by the TADSS scheme involves only the eigenvalues and and not the parameter from the backtracking procedure.

Theorem 9. Let the additional assumptions for the eigenvalues of matrix be imposed and let be the strictly convex quadratic function given by (29). Assume is the orthonormal eigenvectors of symmetric positive definite matrix and suppose that is the sequence of values constructed by Algorithm 2. The gradients of convex quadratics defined by (29) are and can be expressed as for some real constants and for some integer . Then the application of the gradient descent method (7) on the goal function (29) satisfies the following two statements:

Proof. Taking into account (7) one can verify and by taking (42) we get In order to prove (43), it is enough to show that . So, two cases have to be analyzed. In the first one, it is supposed that . Applying (30) leads to In the other case, it is assumed that . From this condition arrives the following conclusion: Expression (42) impliesThe fact that the parameter , from (43), satisfies confirms expression (44).

4. Numerical Experience

Numerical results provided by applying the implementation of TADSS, ADSS, and SM methods on test functions for unconstrained test problems, proposed in [2, 11], are presented and investigated. We chose most of the functions from the set of test functions presented in [3, 4] and, as proposed in these papers, also investigated the experiments with a large number of variables in each function: 1000, 2000, 3000, 5000, 7000, 8000, 10000, 15000, 20000, and 30000. The stopping criteria are the same as in [1, 3, 4]. Backtracking procedure is developed using the values , of needed parameters. Three main indicators of the efficiency are observed: number of iterations, CPU time, and number of function evaluations. First, we compare the performance of the TADSS scheme with the ADSS method. The reasons for this selection is obvious: the TADSS scheme presents a one-step version of ADSS method. Also, the intention to examine behavior of TADSS and compare it with its forerunner is natural. Obtained numerical values are displayed in Table 1 and refer to the number of iterative steps, the CPU time of executions computed in seconds, and the number of evaluations of the objective function.


Test function Number of iterations CPU time Number of function evaluations
TADSS ADSS TADSS ADSS TADSS ADSS

Extended Penalty 40 50 2 4 1280 1780
Raydan 1 466 34 11 4 8504 4844
Diagonal 1 20 37 0 0 335 1448
Diagonal 3 21 49 0 1 417 1048
Generalized Tridiagonal 1 61 77 0 0 431 719
Extended Tridiagonal 1 60 70 0 0 250 420
Extended Three Expon. Term 40 40 0 0 400 350
Diagonal 4 40 780 0 0 270 2590
Extended Himmelblau 60 70 0 0 300 480
Quadratic QF1 4953 425 15 0 13738 1755
Extended Quad. Penalty QP1 50 60 0 2 571 841
Extended Quad. Penalty QP2 50 60 0 4 569 843
Quadratic QF2 50 60 0 1 583 836
Extended EP1 186 40 1 0 758 487
Extended Tridiagonal 2 638 80 0 0 2052 420
Arwhead 50 64 0 3 601 1082
Engval1 60 70 0 0 290 460
Quartc 10 70 0 0 30 390
Generalized Quartic 60 70 0 0 250 614
Diagonal 7 189 2201 0 33 566 6633
Diagonal 8 160 2213 1 40 586 6709
Diagonal 9 20 43 0 0 352 3646

Obtained numerical results, generally, confirm advantages in favor of TADSS, considering all three tested indicators. More precise, regarding the number of iterative steps TADSS shows better results in 17 out of 22 functions, while ADSS outperforms TADSS in 4 out of 22 experiments and for the extended three exponential terms function both methods require the same number of iterations. Results concerning spanned CPU time confirm that both methods, TADSS and ADSS, are very fast. In 9 out of 22 cases TADSS is faster than ADSS, in 2 out of 22 testings ADSS is faster than TADSS, and even in half of examinations the CPU time of both iterations equals zero. On the issue of the number of evaluations of an objective function, TADSS improves ADSS in 17 out of 22 testings and the opposite case appears in 5 out of 22 cases. Table 2 displays average results of tested values.


Average performances TADSS ADSS

Number of iterations 331.09 302.86
CPU time (sec) 1.36 4.18
Number of function evaluations 1506.05 1745.23

According to results displayed in Table 2, it can be concluded that although TADSS outperforms ADSS in 17 out of 22 testings with regard to the number of iterations, average results show slight advantage of ADSS on this matter. Considering the average number of evaluations, there is an opposite case in favor of TADSS. Consumed CPU time is averagely three times less in favor to the TADSS comparing to the ADSS. Generally, it can be concluded that the one-step variant of the ADSS method, constructed TADSS scheme, behaves slightly better than the original ADSS iteration, especially when we consider the speed of executions.

Some additional experiments have been carried out in further numerical research. These tests show the comparison between the TADSS and the SM iterations. As mentioned before, both of the schemes, TADSS and SM, are accelerated gradient descent methods with one iterative step size parameter. We choose this additional numerical comparison in order to confirm that the accelerated single step size TADSS algorithm, derived from the accelerated double step size ADSS model, gives better performances with respect to the all three analyzed aspects than the classically defined accelerated single step size SM method. Table 3 with 30 displayed test functions verifies the previous assertion.


Test functions Number of iterations CPU time Number of function evaluations
TADSS SM TADSS SM TADSS SM

Extended Penalty 40 589 2 5 1280 2987
Perturbed Quadratic 4276 111972 13 1868 12017 632724
Raydan 1 466 21125 11 178 8504 116757
Diagonal 1 20 10417 0 116 335 56135
Diagonal 3 21 10574 0 209 417 59425
Generalized Tridiagonal 1 61 278 0 2 431 989
Extended Tridiagonal 1 60 3560 0 35 250 30686
Extended Three Expon. Term 40 164 0 1 400 1224
Diagonal 4 40 80 0 0 270 530
Extended Himmelblau 60 168 0 0 300 566
Quadr. Diag. Perturbed 11542 53133 58 1193 56359 547850
Quadratic QF1 4953 114510 15 2127 13738 649643
Extended Quad. Penalty QP1 50 224 0 12 571 2566
Extended Quad. Penalty QP2 50 162 0 7 569 2057
Quadratic QF2 50 118801 0 2544 583 662486
Extended EP1 186 68 1 1 758 764
Extended Tridiagonal 2 638 584 0 0 2052 2144
Arwhead 50 10 0 0 601 30
Almost Perturbed Quadratic 4202 110121 15 2148 14974 627287
Engval1 60 185 0 7 290 2177
Quartc 10 190 0 0 30 430
Generalized Quartic 60 156 0 0 250 423
Diagonal 7 189 90 0 0 566 220
Diagonal 8 160 96 1 0 586 594
Diagonal 9 20 11235 0 118 352 60566
DIXON3DQ 10 112899 0 1908 30 639957
NONSCOMP 10 10 0 0 30 30
HIMMELH 10 30 0 0 40 80
Power Cute 1455 > 4 > 4567 >
Indef 20 > 0 > 50 >

It can be observed from displayed numerical outcomes that the TADSS method provides better results than the SM method considering the number of iterations in 24 out of 30 testings, while the number of opposite cases is 5 out of 30. For the NONSCOMP test function, both models have the same number of iterations. Concerning the CPU time, both algorithms give the same results for 10 test functions. The TADSS method is faster than SM for 19 test functions, while the SM method is faster than TADSS for one test function only. The greatest progress is obtained with respect to the number of evaluations of the objectives. On this matter, using the TADSS algorithm, better results are obtained in 27 out of 30 test functions, while the opposite case holds for two test functions only. For the NONSCOMP test function both of the compared iterations give the same number of evaluations. From Table 3, we can also notice that for 2 out of 30 test functions testings are lasting more than the time limiter constant defined in [3], while for all 30 test functions when the TADSS algorithm is applied the time of execution is far less than .

The results arranged in Table 4 give even more general view on the benefits provided by applying the TADSS method with regard to the SM method. The average values of 28 test functions, which were possible to test by both methods according to the constant , are presented in the table.


Average performance TADSS SM

Number of iterations 976.21 24336.82
CPU time (in seconds) 4.14 445.68
Number of evaluations 4163.68 146475.96

The results presented in the previous table confirm that by applying the TADSS method approximately 25 times less iterations and even 35 times less evaluations of the objective function are needed in comparison with the SM method. Finally, when the TADSS is used, the testing is lasting even 107 times shorter than when the SM is applied.

The codes for presented numerical experiments are written in the Visual C++ programming language on a Workstation Intel 2.2 GHz.

5. Conclusion

The accelerated single step size gradient descent algorithm, called TADSS, is defined as a transformation of the accelerated double step size gradient descent model ADSS, proposed in [4]. More precisely, the TADSS scheme is derived from the accelerated double step size gradient descent scheme ADSS by imposing relation (6) between two step parameters in the ADSS iteration. The efficiency of ADSS model regarding all analyzed characteristics in comparison to the accelerated gradient descent single step size SM method has been numerically proved in [4].

The method defined in this way is comparable with its double step size forerunner, ADSS method, as well as with the single step size accelerated gradient descent SM method which is defined in a classical manner.

Results illustrated in Tables 1 and 2 generally indicate that the TADSS method behaves similarly as the ADSS method. From the point of mean values, the ADSS scheme gives slightly better results considering the number of iterations. On the other hand, a certain improvement regarding the number of function evaluations and needed CPU time is obtained by applying the TADSS iterations.

Even greater advantages of derived TADSS method are presented in Tables 3 and 4 where the comparisons between the TADSS and the SM are given. Evidently, the TADSS scheme improves the SM method with respect to all three analyzed characteristics, which was the prime goal in the research presented herein.

Linear convergence of the TADSS method is proved for the uniformly convex and the strictly convex quadratic functions.

Obtained results motivate further investigations of possible accelerated double step size gradient descent models and its transformations into corresponding single step size variants.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors gratefully acknowledge the financial support of the Serbian Ministry of Education, Science and Technological Development.

References

  1. P. S. Stanimirović and M. B. Miladinović, “Accelerated gradient descent methods with line search,” Numerical Algorithms, vol. 54, no. 4, pp. 503–520, 2010. View at: Publisher Site | Google Scholar | MathSciNet
  2. N. Andrei, “An acceleration of gradient descent algorithm with backtracking for unconstrained optimization,” Numerical Algorithms, vol. 42, no. 1, pp. 63–73, 2006. View at: Publisher Site | Google Scholar | MathSciNet
  3. M. J. Petrović and P. S. Stanimirović, “Accelerated double direction method for solving unconstrained optimization problems,” Mathematical Problems in Engineering, vol. 2014, Article ID 965104, 8 pages, 2014. View at: Publisher Site | Google Scholar | MathSciNet
  4. M. J. Petrović, “An accelerated double step size model in unconstrained optimization,” Applied Mathematics and Computation, vol. 250, pp. 309–319, 2015. View at: Publisher Site | Google Scholar | MathSciNet
  5. N. I. Djuranović-Miličić and M. Gardaševic-Filipović, “A multi-step curve search algorithm in nonlinear optimization: nondifferentiable convex case,” Facta Universitatis, Series: Mathematics and Informatics, vol. 25, pp. 11–24, 2010. View at: Google Scholar | MathSciNet
  6. J. M. Ortega and W. C. Rheinboldt, Iterative Solution of Nonlinear Equation in Several Variables, Academic Press, London, UK, 1970. View at: MathSciNet
  7. R. T. Rockafellar, Convex Analysis, Princeton University Press, Princeton, NJ, USA, 1970. View at: MathSciNet
  8. J. Barzilai and J. M. Borwein, “Two-point step size gradient methods,” IMA Journal of Numerical Analysis, vol. 8, no. 1, pp. 141–148, 1988. View at: Publisher Site | Google Scholar | MathSciNet
  9. Y.-H. Dai and L.-Z. Liao, “R-linear convergence of the Barzilai and Borwein gradient method,” IMA Journal of Numerical Analysis, vol. 22, no. 1, pp. 1–10, 2002. View at: Publisher Site | Google Scholar | MathSciNet
  10. B. Molina and M. Raydan, “Preconditioned Barzilai-BORwein method for the numerical solution of partial differential equations,” Numerical Algorithms, vol. 13, no. 1-2, pp. 45–60, 1996. View at: Publisher Site | Google Scholar | MathSciNet
  11. N. Andrei, “An unconstrained optimization test functions collection,” Advanced Modeling and Optimization, vol. 10, pp. 147–161, 2008. View at: Google Scholar

Copyright © 2015 Predrag S. Stanimirović et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

958 Views | 339 Downloads | 4 Citations
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.