We propose an analysis of the effects introduced by finite-accuracy and round-off arithmetic on numerical computations of discrete dynamical systems. Our method, which uses the statistical tool of the decay of fidelity, computes the error by directly comparing the numerical orbit with the exact one (or, more precisely, with another numerical orbit computed with a much higher accuracy). Furthermore, as a model of the effects of round-off arithmetic on the map, we also consider a random perturbation of the exact orbit with an additive noise, for which exact results can be obtained for some prototype maps. We investigate the decay laws of fidelity and their relationship with the error probability distribution for regular and chaotic maps, for both additive and numerical noise. In particular, for regular maps we find an exponential decay for additive noise, and a power-law decay for numerical noise. For chaotic maps, numerical noise is equivalent to additive noise, and our method is suitable for identifying a threshold for the reliability of numerical results, i.e., the number of iterations below which global errors can be ignored. This threshold grows linearly with the number of bits used to represent real numbers.

Turchetti G, Vaienti S, Zanlungo F (2010). Asymptotic distribution of global errors in the numerical computations of dynamical systems. PHYSICA. A [10.1016/j.physa.2010.06.060].

Asymptotic distribution of global errors in the numerical computations of dynamical systems

Zanlungo F
2010-01-01

Abstract

We propose an analysis of the effects introduced by finite-accuracy and round-off arithmetic on numerical computations of discrete dynamical systems. Our method, which uses the statistical tool of the decay of fidelity, computes the error by directly comparing the numerical orbit with the exact one (or, more precisely, with another numerical orbit computed with a much higher accuracy). Furthermore, as a model of the effects of round-off arithmetic on the map, we also consider a random perturbation of the exact orbit with an additive noise, for which exact results can be obtained for some prototype maps. We investigate the decay laws of fidelity and their relationship with the error probability distribution for regular and chaotic maps, for both additive and numerical noise. In particular, for regular maps we find an exponential decay for additive noise, and a power-law decay for numerical noise. For chaotic maps, numerical noise is equivalent to additive noise, and our method is suitable for identifying a threshold for the reliability of numerical results, i.e., the number of iterations below which global errors can be ignored. This threshold grows linearly with the number of bits used to represent real numbers.
2010
Settore PHYS-06/A - Fisica per le scienze della vita, l'ambiente e i beni culturali
Turchetti G, Vaienti S, Zanlungo F (2010). Asymptotic distribution of global errors in the numerical computations of dynamical systems. PHYSICA. A [10.1016/j.physa.2010.06.060].
File in questo prodotto:
File Dimensione Formato  
2010_j_physicaa_asymptotic.pdf

accesso aperto

Descrizione: ARTICOLO
Tipologia: Pre-print
Dimensione 1.1 MB
Formato Adobe PDF
1.1 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10447/666284
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact