Hybrid Monte Carlo is very accurate calculation method of the thermodynamic properties of the substances in the classical molecular mechanics. I would like to announce a significant improvement in the speed of the GPU calculations on the Abalone program.
However, this is the most common Hybrid Monte Carlo algorithm. The only interesting thing is that it uses an VEFRL Omelyan algorithm as MD step. This algorithm is the fourth order, and if you reduce the length of the step only a few times, then you already reach the machine precision.
I do not know for what in the molecular dynamics can be used trajectory with double precision. However, if it is needed, it is easily achieved by Omelyan algorithms.
I assume the relevant reference for the “VEFRL Omelyan algorithm” is this:
[url]http://www.sciencedirect.com/science/article/pii/S0010465502004514[/url]
I.P. Omelyan, I.M. Mryglod, R. Folk
Optimized Forest–Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
Computer Physics Communications
Volume 146, Issue 2, 1 July 2002, Pages 188–202
From the preprint, [url]http://arxiv.org/pdf/cond-mat/0110585.pdf[/url], equation (21) seems to be the relevant computation. If so, it seems very amenable to the application of the FMA (fused multiply-add) operation provided by GPUs. This should help somewhat with maintaining accuracy by reducing round-off errors and giving some protection against subtractive cancellation.