Hello Samuel, Is this some IEEE standard overflow wraparound? For information, Julia (0.6) does not behave exactly in the same way. Both 'print((2^53 + 2^2) + 1 == (2^53 + 2^2))' and 'print((2^53 + 2^1) + 1 == (2^53 + 2^1))' return false. Antoine