Hi, On 09/09/2015 08:18 AM, Mathieu Malaterre wrote: > 1. removed -ffast-math > 2. add -std=c99 > 3. Change the code for an explicit storage of division (div = dd / dr) > => issue solved ! Not really, no. Al you've done is conceal the bug. I really think you should fix the bug. Otherwise it will keep coming back. -ffloat-store and -fexcess-precision=standard are just workarounds. My problem, of course, is that I don't really know what this code is trying to do, and why it is so sensitive to the last significant bits of floating-point precision. Andrew.