Reprint Server

Computation of High-Order Maps to Multiple Machine Precision


The Beam Dynamics simulation package in COSY INFINITY is built upon a differential algebra data type. With it, it is possible to compute transfer maps or arbitrary systems to arbitrary order. However, this data type is limited by the precision of the underlying floating point number model provided by the computer processor.

We will present a method to extend the effective precision of the calculations based purely on standard floating point operations. Those algorithms are then integrated into the differential algebra data type to efficiently extend the available precision, without unnecessarily affecting overall efficiency. To that effect, the precision of each coefficient is adjusted automatically during the calculation.

We will then proceed to show the effectiveness of our implementation by calculating high precision maps of combinations of homogeneous dipole segments, for which the exact results are known, and comparing the high precision coefficients with the results produced by the traditional COSY beam physics package.

A. Wittig, M. Berz, Int. Journal of Modern Physics A 24,5 (2009) 1019-1039


Click on the icon to download the corresponding file.

Download Adobe PDF version (466217 Bytes).

Go Back to the reprint server.
Go Back to the home page.

This page is maintained by Kyoko Makino. Please contact her if there are any problems with it.