Searching for sinks of Henon map using a multiple-precision GPU arithmetic library

Mioara Joldes (1) , Valentina Popescu (1) , Warwick Tucker (2)

1 LAAS-MAC - Équipe Méthodes et Algorithmes en Commande
2 CAPA - Computer-aided proofs in analysis [Uppsala]

Today, GPUs represent an important hardware development platform for many problems in dynamical systems, where massive parallel computations are needed. Beside that, many numerical studies of chaotic dynamical systems require a computing precision higher than common floating point (FP) formats. One such application is locating invariant sets for chaotic dynamical systems. In particular, we focus on rigorously proving the existence of stable periodic orbits for the Henon map for parameter values close to the classical ones. For that, we present a multiple-precision floating-point arithmetic library in CUDA programming language for the NVIDIA GPU platform. Our library extends the precision using so-called FP expansions, where a number is represented as the unevaluated sum of standard machine precision FP numbers. This format offers the advantage of using directly available and highly optimized hardware FP operations. We generalize algorithms used by multiple-precisions libraries such as Bailey's QD, or the analogue GPU version, GQD.

Mots clés: floating-point arithmetic; multiple precision library; GPGPU computing; error-free transform; Henon map; dynamical systems

Domaines: Arithmétique des ordinateurs; Systèmes dynamiques [math.DS]

Mioara Joldes, Valentina Popescu, Warwick Tucker. Searching for sinks of Henon map using a multiple-precision GPU arithmetic library. 2013, 6p. ⟨hal-00957438⟩