On improved three-step schemes with high efficiency index and their dynamics

This paper presents an improvement of the sixth-order method of Chun and Neta as a class of three-step iterations with optimal efficiency index, in the sense of Kung-Traub conjecture. Each member of the presented class reaches the highest possible order using four functional evaluations. Error analysis will be studied and numerical examples are also made to support the theoretical results. We then present results which describe the dynamics of the presented optimal methods for complex polynomials. The basins of attraction of the existing optimal methods and our methods are presented and compared to illustrate their performances.This research was supported by Ministerio de Ciencia y Tecnologia MTM2011-28636-C02-02 and FONDOCYT Republica Dominicana.Babajee, DKR.; Cordero Barbero, A.; Soleymani, F.; Torregrosa Sánchez, JR. (2014). On improved three-step schemes with high efficiency index and their dynamics. Numerical Algorithms. 65(1):153-169. https://doi.org/10.1007/s11075-013-9699-6S153169651Pang, J.S., Chan, D.: Iterative methods for variational and complementary problems. Math. Program. 24(1), 284–313 (1982)Sun, D.: A class of iterative methods for solving nonlinear projection equations. J. Optim. Theory Appl. 91(1), 123–140 (1996)Chun, C., Neta, B.: A new sixth-order scheme for nonlinear equations. Appl. Math. Lett. 25, 185–189 (2012)Kung, H.T., Traub, J.F.: Optimal order of one-point and multipoint iteration. J. ACM 21, 643–651 (1974)Neta, B.: A new family of high-order methods for solving equations. Int. J. Comput. Math. 14, 191–195 (1983)Neta, B.: On Popovski’s method for nonlinear equations. Appl. Math. Comput. 201, 710–715 (2008)Chun, C., Neta, B.: Some modifications of Newton’s method by the method of undeterminate coefficients. Comput. Math. Appl. 56, 2528–2538 (2008)Chun, C., Lee, M.Y., Neta, B., Dzunic, J.: On optimal fourth-order iterative methods free from second derivative and their dynamics. Appl. Math. Comput. 218, 6427–6438 (2012)Cordero, A., Torregrosa, J.R., Vassileva, M.P.: Three-step iterative methods with optimal eighth-order convergence. J. Comput. Appl. Math. 235, 3189–3194 (2011)Cordero, A., Torregrosa, J.R., Vassileva, M.P.: A family of modified Ostrowski’s methods with optimal eighth order of convergence. Appl. Math. Lett. 24, 2082–2086 (2011)Heydari, M., Hosseini, S.M., Loghmani, G.B.: On two new families of iterative methods for solving nonlinear equations with optimal order. Appl. Anal. Dis. Math. 5, 93–109 (2011)Neta, B., Petkovic, M.S.: Construction of optimal order nonlinear solvers using inverse interpolation. Appl. Math. Comput. 217, 2448–2445 (2010)Sharifi, M., Babajee, D.K.R., Soleymani, F.: Finding the solution of nonlinear equations by a class of optimal methods. Comput. Math. Appl. 63, 764–774 (2012)Soleymani, F., Karimi Vanani, S., Khan, M., Sharifi, M.: Some modifications of King’s family with optimal eighth order of convergence. Math. Comput. Model. 55, 1373–1380 (2012)Soleymani, F., Karimi Vanani, S., Jamali Paghaleh, M.: A class of three-step derivative-free root solvers with optimal convergence order. J. Appl. Math. 2012, Article ID 568740, 15 pp. (2012). doi: 10.1155/2012/568740Soleymani, F., Sharifi, M., Mousavi, B.S.: An improvement of Ostrowski’s and King’s techniques with optimal convergence order eight. J. Optim. Theory Appl. 153, 225–236 (2012)Stewart, B.D.: Attractor basins of various root-finding methods. M.S. Thesis, Naval Postgraduate School, Department of Applied Mathematics, Monterey, CA (2001)Amat, S., Busquier, S., Plaza, S.: Review of some iterative root-finding methods from a dynamical point of view. Scientia 10, 3–35 (2004)Amat, S., Busquier, S., Plaza, S.: Dynamics of the King and Jarratt iterations. Aequ. Math. 69, 212–223 (2005)Amat, S., Busquier, S., Plaza, S.: Chaotic dynamics of a third-order Newton type method. J. Math. Anal. Appl. 366, 24–32 (2010)Neta, B., Chun, C., Scott, M.: A note on the modified super-Halley method. Appl. Math. Comput. 218, 9575–9577 (2012)Scott, M., Neta, B., Chun, C.: Basin attractors for various methods. Appl. Math. Comput. 218, 2584–2599 (2011)Ardelean, G.: A comparison between iterative methods by using the basins of attraction. Appl. Math. Comput. 218, 88–95 (2011)Traub, J.F.: Iterative Methods for the Solution of Equations. Prentice Hall, New York (1964)Babajee, D.K.R.: Analysis of higher order variants of Newton’s method and their applications to differential and integral equations and in ocean acidification. Ph.D. Thesis, University of Mauritius (2010

Computing Real Roots of Real Polynomials ... and now For Real!

Very recent work introduces an asymptotically fast subdivision algorithm, denoted ANewDsc, for isolating the real roots of a univariate real polynomial. The method combines Descartes' Rule of Signs to test intervals for the existence of roots, Newton iteration to speed up convergence against clusters of roots, and approximate computation to decrease the required precision. It achieves record bounds on the worst-case complexity for the considered problem, matching the complexity of Pan's method for computing all complex roots and improving upon the complexity of other subdivision methods by several magnitudes. In the article at hand, we report on an implementation of ANewDsc on top of the RS root isolator. RS is a highly efficient realization of the classical Descartes method and currently serves as the default real root solver in Maple. We describe crucial design changes within ANewDsc and RS that led to a high-performance implementation without harming the theoretical complexity of the underlying algorithm. With an excerpt of our extensive collection of benchmarks, available online at http://anewdsc.mpi-inf.mpg.de/, we illustrate that the theoretical gain in performance of ANewDsc over other subdivision methods also transfers into practice. These experiments also show that our new implementation outperforms both RS and mature competitors by magnitudes for notoriously hard instances with clustered roots. For all other instances, we avoid almost any overhead by integrating additional optimizations and heuristics.Comment: Accepted for presentation at the 41st International Symposium on Symbolic and Algebraic Computation (ISSAC), July 19--22, 2016, Waterloo, Ontario, Canad

On the Minimization of Convex Functionals of Probability Distributions Under Band Constraints

The problem of minimizing convex functionals of probability distributions is solved under the assumption that the density of every distribution is bounded from above and below. A system of sufficient and necessary first-order optimality conditions as well as a bound on the optimality gap of feasible candidate solutions are derived. Based on these results, two numerical algorithms are proposed that iteratively solve the system of optimality conditions on a grid of discrete points. Both algorithms use a block coordinate descent strategy and terminate once the optimality gap falls below the desired tolerance. While the first algorithm is conceptually simpler and more efficient, it is not guaranteed to converge for objective functions that are not strictly convex. This shortcoming is overcome in the second algorithm, which uses an additional outer proximal iteration, and, which is proven to converge under mild assumptions. Two examples are given to demonstrate the theoretical usefulness of the optimality conditions as well as the high efficiency and accuracy of the proposed numerical algorithms.Comment: 13 pages, 5 figures, 2 tables, published in the IEEE Transactions on Signal Processing. In previous versions, the example in Section VI.B contained some mistakes and inaccuracies, which have been fixed in this versio