Bounding the Round-Off Error of the Upwind Scheme for Advection
Louise Ben Salem-Knapp  3, 2, 1@  , Sylvie Boldo  5, 4, 2  , William Weens@
3 : Laboratoire de Recherche en Informatique
CNRS : UMR8623, Université Paris Sud
2 : INRIA
INRIA
1 : Université Paris-Saclay
Université Paris-Saclay, Université Paris Saclay
5 : Laboratoire de Recherche en Informatique  (LRI)  -  Site web
CNRS : UMR8623, Université Paris Sud
LRI - Bâtiments 650-660 Université Paris-Sud 91405 Orsay Cedex -  France
4 : Université Paris-Saclay
Université Paris-Saclay, Université Paris Saclay

Numerical simulations are carefully-written programs, and their correctness is based on mathematical results. Nevertheless, those programs rely on floating-point arithmetic and the corresponding round-off errors are often ignored. This article deals with a specific simple scheme applied to advection, that is a particular equation from hydrodynamics dedicated to the transport of a substance. It shows a tight bound on the round-off error of the 1D and 2D upwind scheme, with an error roughly proportional to the number of steps. The error bounds are generic with respect to the format and exceptional behaviors are taken into account. Some experiments give an insight of the quality of the bounds.



  • Poster
Personnes connectées : 1 Vie privée
Chargement...