Stochastic arithmetic enables one to estimate round-off error propagation using a probabilistic approach. With Stochastic arithmetic, the numerical quality of any simulation program can be controlled. Furthermore by detecting all the instabilities which may occur at run time, a numerical debugging of the user code can be performed. Stochastic arithmetic can be used to dynamically control approximation methods. Such methods provide a result which is affected by a truncation error inherent to the algorithm used and a round-off error due to the finite precision of the computer arithmetic. If the discretization step decreases, the truncation error also decreases, but the round-off error increases. Therefore it can be difficult to control these two errors simultaneously. In order to obtain with an approximation method a result for which the global error (consisting of both the truncation error and the round-off error) is minimal, a strategy, based on a converging sequence computation, has been proposed. Computation is carried out until the difference between two successive iterates has no exact significant digit. Then it is possible to determine which digits of the result obtained are in common with the exact solution. This strategy can apply to the computation of integrals using the trapezoidal rule, Simpson's rule, Romberg's method or the Gauss—Legendre method.
This is a preview of subscription content, log in via an institution.
Buying options
Tax calculation will be finalised at checkout
Purchases are for personal use only
Learn about institutional subscriptionsPreview
Unable to display preview. Download preview PDF.
References
Ahmed, Z. Definitely an integral. American Mathematical Monthly, 109(7):670–671, 2002.
Bailey, D.H. and X.S. Li. A comparison of three high-precision quadrature schemes. In Proceedings of the 5th Real Numbers and Computers Conference, pages 81–95, Lyon, France, September 2003.
Bauer, F.L., H. Rutishauser, and E. Stiefel. New aspects in numerical quadrature. In Proceedings of the Symposia Applied Mathematics, volume XV, pages 199–218, American Mathematical Society, Providence, R.I., 1963.
Burden, R.L. and J.D. Faires. Numerical analysis. Brooks-Cole Publishing, 7th edition, Pacific Grove, CA, USA, 2001.
Chesneaux, J.-M. Study of the computing accuracy by using probabilistic approach. In C. Ullrich, editor, Contribution to Computer Arithmetic and Self-Validating Numerical Methods, pages 19–30, IMACS, New Brunswick, New Jersey, USA, 1990.
Chesneaux, J.-M. The quality relations in scientific computing. Numerical Algorithms, 7:129–143, 1994.
Chesneaux, J.-M. L'arithmétique stochastique et le logiciel CADNA. Habilitation à diriger des recherches, Université Pierre et Marie Curie, Paris, November 1995.
Chesneaux, J.-M. and F. Jézéquel. Dynamical control of computations using the trapezoidal and Simpson's rules. Journal of Universal Computer Science, 4(1):2–10, 1998.
Chesneaux, J.-M. and J. Vignes. Sur la robustesse de la méthode CESTAC. Comptes rendus de l'Académie des sciences, Paris, Série I, Mathématiques, 307:855–860, 1988.
Chesneaux, J.M. and J. Vignes. Les fondements de l'arithmétique stochastique. Comptes Rendus de l'Académie des Sciences de Paris, Série I. Mathématiques, 315:1435–1440, 1992.
Conte, S.D. and C. de Boor. Elementary numerical analysis. McGraw-Hill, International Student edition, 1980.
Engels, H. Numerical quadrature and cubature. Academic Press, London, UK, 1980.
Goldberg, D. What every computer scientist should know about floating-point arithmetic. ACM Computing Surveys, 23(1):5–48, 1991.
Jézéquel, F. Dynamical control of converging sequences computation. Applied Numerical Mathematics, 50(2):147–164, 2004.
Jézéquel, F. Contrôle dynamique de méthodes d'approximation. Habilitation à diriger des recherches, Université Pierre et Marie Curie, Paris, February 2005.
Jézéquel, F. and J.-M. Chesneaux. Computation of an infinite integral using Romberg's method. Numerical Algorithms36(3):265–283, July 2004.
Jézéquel, F., F. Rico, J.-M. Chesneaux, and M. Charikhi. Reliable computation of a multiple integral involved in the neutron star theory. Mathematics and computers in simulation, 71(1):44–61, 2006.
Muller, J.-M. Arithmétique des ordinateurs. Masson, Paris, France, 1989.
Rump, S.M. Reliability in computing. The role of interval methods in scientific computing. Academic Press, London, UK, 1988.
Stoer, J. and R. Bulirsch. Introduction to numerical analysis, volume 12 of Texts in applied mathematics. Springer, 3rd edition, 2002.
Vignes, J. Zéro mathématique et zéro informatique. Comptes rendus de l'Académie des sciences, Paris, Série I, Mathématiques, 303:997–1000, 1986. also: La Vie des Sciences, 4 (1):1–13, 1987.
Vignes, J. Estimation de la précision des résultats de logiciels numériques. La Vie des Sciences, 7(2):93–145, 1990.
Vignes, J. A stochastic arithmetic for reliable scientific computation. Mathematics and Computers in Simulation, 35:233–261, 1993.
Vignes, J. A stochastic approach to the analysis of round-off error propagation. A survey of the CESTAC method. In Proceedings of the 2nd Real Numbers and Computers conference, pages 233–251, Marseille, France, 1996.
Vignes, J. and M. La Porte. Error analysis in computing. In Information Processing 1974, pages 610–614. North-Holland, 1974.
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2009 Springer Science + Business Media B.V
About this paper
Cite this paper
Chesneaux, JM., Jézéquel, F., Lamotte, JL. (2009). Stochastic arithmetic and verification of mathematical models. In: Baveye, P.C., Laba, M., Mysiak, J. (eds) Uncertainties in Environmental Modelling and Consequences for Policy Making. NATO Science for Peace and Security Series C: Environmental Security. Springer, Dordrecht. https://doi.org/10.1007/978-90-481-2636-1_5
Download citation
DOI: https://doi.org/10.1007/978-90-481-2636-1_5
Publisher Name: Springer, Dordrecht
Print ISBN: 978-90-481-2635-4
Online ISBN: 978-90-481-2636-1
eBook Packages: Earth and Environmental ScienceEarth and Environmental Science (R0)