Relaxation to the asymptotic distribution of global errors due to round off

• Published in: Europhysics letters Vol. 89; no. 4; p. 40006
• Main Authors: Turchetti, G, Vaienti, S, Zanlungo, F
• Format: Journal Article
• Language: English
• Published: IOP Publishing 01-02-2010; EDP Sciences
• Subjects: 05.45.-a; 05.45.Pq; Asymptotic properties; Decay; Dynamical systems; Errors; Mathematical models; Orbits; Perturbation methods
• ISSN: 0295-5075; 1286-4854
• DOI: 10.1209/0295-5075/89/40006

We propose an analysis of the effects introduced by finite accuracy and round-off arithmetic on discrete dynamical systems. We investigate, from a statistical viewpoint and using the tool of the decay of fidelity, the error of the numerical orbit with respect to the exact one. As a model we consider a random perturbation of the exact orbit with an additive noise, for which exact results can be obtained for some prototype maps. For regular anysocrounous maps the fidelity has a power law decay, whereas the decay is exponential if a random perturbation is introduced. For chaotic maps the decay is superexponential after an initial plateau and our method is suitable to identify the reliability threshold of numerical results, i.e. a number of iterations below which global errors can be ignored. The same behaviour is observed if a random perturbation is introduced.