Searching for sinks of Henon map using a multiple-precision GPU arithmetic library
LAAS-CNRS, 7 Avenue du Colonel Roche, 31077 Toulouse, France
hal-00957438, (10 March 2014)
@unpublished{joldes:hal-00957438,
hal_id={hal-00957438},
url={http://hal.archives-ouvertes.fr/hal-00957438},
title={Searching for sinks of Henon map using a multiple-precision GPU arithmetic library},
author={Joldes, Mioara and Popescu, Valentina and Tucker, Warwick},
keywords={floating-point arithmetic; multiple precision library; GPGPU computing; error-free transform; Henon map; dynamical systems},
language={Anglais},
affiliation={Laboratoire d’analyse et d’architecture des syst{`e}mes [Toulouse] – LAAS , Department of Mathematics [Uppsala]},
note={7 pages Rapport LAAS ndegre 13492 Rapport LAAS ndegre 13492 },
year={2014},
month={Mar},
pdf={http://hal.archives-ouvertes.fr/hal-00957438/PDF/heart_research_report.pdf}
}
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.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.
March 15, 2014 by hgpu
