[gillianseed]

Could you try with -Ofast which enables -ffast-math and post the results?

[nkurz]

I get no significant difference with -Ofast for either icc or gcc. My code is still untested and quite possibly buggy, but I put it up here: http://pastebin.com/qB0DfkXN