Abstract
The expansion of complicated functions of many variables in Taylor polynomials is an important problem for many applications, and in practice can be performed rather conveniently (even to high orders) using polynomial algebras. An important application of these methods is the field of beam physics, where often expansions in about six variables to orders between five and ten are used.
However, often it is necessary to also know bounds for the remainder term of the Taylor formula if the arguments lie within certain intervals. In principle such bounds can be obtained by interval bounding of the (n + 1)-st derivative, which in turn can be obtained with polynomial algebra; but in practice the method is rather inefficient and susceptible to blow-up because of the need of repeated interval evaluations of the derivative. Here we present a new method that allows the computation of sharp remainder intervals in parallel with the accumulation derivatives up to order n.
The method is useful for a variety of numerical problems, including the interval inclusion of very complicated functions prone to blow-up. To this end, the function is represented by a Taylor polynomial with remainder using the above method. Since at least for high orders, the remainder terms have a tendency to be very small, the problem is reduced to an interval evaluation of the Taylor polynomial. The method is used for guaranteed global optimization of blow-up prone functions and compared with some interval-based global optimization schemes.
Similar content being viewed by others
References
Berz, M.: Analyis on a Nonarchimedean Extension of the Real Numbers, Technical Report MSUCL-753, National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI 48824, 1990.
Berz, M.: Arbitrary Order Description of Arbitrary Particle Optical Systems, Nuclear Instruments and Methods A298 (1990), p. 426.
Berz, M.: Automatic Differentiation as Nonarchimedean Analysis, in: Computer Arithmetic and Enclosure Methods, Elsevier Science Publishers, Amsterdam, 1992, p. 439.
Berz, M.: COSY INFINITY Version 7 Reference Manual, Technical Report MSUCL-977, National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI 48824, 1996.
Berz, M.: COSY INFINITY Version 6, in: Berz, M., Martin, S., and Ziegler, K. (eds), Proc. Nonlinear Effects in Accelerators, IOP Publishing, 1992, p. 125.
Berz, M.: Differential Algebraic Formulation of Normal Form Theory, in: Berz, M., Martin, S., and Ziegler, K. (eds), Proc. Nonlinear Effects in Accelerators, IOP Publishing, 1992, p. 77.
Berz, M.: Differential Algebraic Treatment of Beam Dynamics to Very High Orders Including Applications to Spacecharge, AIP Conference Proceedings 177 (1988), p. 275.
Berz, M.: Forward Algorithms for High Orders and Many Variables, Automatic Differentiation of Algorithms: Theory, Implementation and Application, SIAM, 1991.
Berz, M.: High-Order Description of Accelerators Using Differential Algebra and First Applications to the SSC, in: Proceedings, Snowmass Summer Meeting, Snowmass, Co, 1988.
Berz, M.: High-Order Computation and Normal Form Analysis of Repetitive Systems, in: Month, M. (ed.), Physics of Particle Accelerators, volume AIP 249, American Institute of Physics, 1991, p. 456.
Berz, M. and Hoffstätter, G.: Exact Bounds of the Long Term Stability of Weakly Nonlinear Systems Applied to the Design of Large Storage Rings, Interval Computations 2 (1994), pp. 68-89.
Berz, M., Bischof, C., Griewank, A., and Corliss, G. (eds), Computational Differentiation: Techniques, Applications, and Tools, SIAM, Philadelphia, 1996.
Berz, M., Hoffstätter, G., Wan, W., Shamseddine, K., and Makino, K.: COSY INFINITY and Its Applications to Nonlinear Dynamics, in: Computational Differentiation: Techniques, Applications, and Tools, SIAM, 1996.
Csenes, T.: Test Results of Interval Methods for Global Optimization, Computer Arithmetic, Scientific Computation and Mathematical Modelling 12 (1991).
Griewank, A. and Corliss, G. F. (eds), Automatic Differentiation of Algorithms, SIAM, Philadelphia, 1991.
Hansen, E. R. A Generalized Interval Arithmetic, in: Nickel, K. (ed.), Interval Mathematics, 1975.
Hansen, E. R.: Global Optimization Using Interval Analysis—the One-Dimensional Case, J. Optim. Theor. and Appl. 29 (1979), pp. 331-334.
Hansen, E. R.: Global Optimization Using Interval Analysis—the Multidimensional Case, Numerische Mathematik 34 (1980), pp. 247-270.
Jansson, C.: A Global Optimization Method Using Interval Arithmetic, in: Computer Arithmetic and Scientific Computation, Proceedings of the SCAN 91, Amsterdam, North-Holland, Elsevier, 1992.
Makino, K. and Berz, M.: COSY INFINITY Version 7, in: Fourth Computational Accelerator Physics Conference, AIP Conference Proceedings, 1996.
Makino, K. and Berz, M.: Implementation and Applications of Taylor Model Methods.
Makino, K. and Berz, M.: Remainder Differential Algebras and Their Applications, in: Computational Differentiation: Techniques, Applications, and Tools, SIAM, 1996.
Moore, R. E. and Ratschek, H.: Inclusion Functions and Global Optimization II. Mathematical Programming 41 (1988), pp. 341-356.
Nekhoroschev, N. N.: An Exponential Estimate of the Time of Stability of Nearly Integrable Hamiltonian Systems, Uspekhi Mat. Nauk, English translation Russ. Math. Surv. 32(6) (1977).
Rokne, J. and Ratschek, H.: New Computer Methods for Global Optimization, Ellis Horwood Limited, Chichester, England, 1988.
Shamseddine, K. and Berz, M.: Exception Handling in Derivative Computation with Nonarchimedean Calculus, in: Computational Differentiation: Techniques, Applications, and Tools, SIAM, 1996.
Talman, R.: Long Term Prediction and the SSC, in: Nonlinear Problems in Future Accelerators, World Scientific, New York, 1991, pp. 215-231.
Turchetti, G.: Nekhoroshev Stability Estimates for Symplectic Maps and Physical Applications, in: Number Theory and Physics, Springer Proceedings in Physics 47, Springer-Verlag, Berlin, Heidelberg, 1990.
Warnock, R. L.: Close Approximation to Invariant Tori in Nonlinear Mechanics, Physical Review Letters 66(14) (1991), pp. 1803-1806.
Warnock, R. L. and Ruth, R. D.: Long-Term Bounds on Nonlinear Hamiltonian Motion, Physica D 56(14) (1992), pp. 188-215, also SLAC-PUB-5267.
Warnock, R. L. and Ruth, R. D.: Stability of Orbits in Nonlinear Mechanics for Finite but Very Long Times, in: Nonlinear Problems in Future Accelerators, World Scientific, New York, 1991, pp. 67-76, also SLAC-PUB-5304.
Warnock, R. L., Ruth, R. D., Gabella, W., and Ecklund, K.: Methods of Stability Analysis in Nonlinear Mechanics, in: 1987 Accelerator Physics Summer School AIP Conference Proceedings, 1988, also SLAC-PUB-4846, 1989.
Yan, Y.: Applications of Differential Algebra to Single-Particle Dynamics in Storage Rings, Technical Report SSCL-500, Superconducting Super Collider Laboratory, 1991.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Berz, M., Hoffstätter, G. Computation and Application of Taylor Polynomials with Interval Remainder Bounds. Reliable Computing 4, 83–97 (1998). https://doi.org/10.1023/A:1009958918582
Issue Date:
DOI: https://doi.org/10.1023/A:1009958918582