I am too busy to particitape in this discussion, but you might want to take a look in GMP's mpn/cray directories, in particular in mpn/cray/ieee. It contains some vectorised code for classical Cray machines. Perhaps GPU's can use similar ideas? -- Torbjörn