floating point performance, but suffer a huge slowdown when emulating double precision operations. For the
solution of linear systems of equations mixed precision iterative refinement methods are powerful
techniques which allow to concentrate double precision operations in only very few most relevant places
and perform 99% of computations in single precision, while still obtaining the same accuracy as a full
double precision solver. Combining these methods with the fast parallel single precision co-processors
allows to exploit the superior bandwidth and computation power of these devices also in applications with
high accuracy requirements.