Search code examples
c#doubleieee-754floating-point-precisionepsilon

Get next smallest Double number


As part of a unit test, I need to test some boundary conditions. One method accepts a System.Double argument.

Is there a way to get the next-smallest double value? (i.e. decrement the mantissa by 1 unit-value)?

I considered using Double.Epsilon but this is unreliable as it's only the smallest delta from zero, and so doesn't work for larger values (i.e. 9999999999 - Double.Epsilon == 9999999999).

So what is the algorithm or code needed such that:

NextSmallest(Double d) < d

...is always true.


Solution

  • If your numbers are finite, you can use a couple of convenient methods in the BitConverter class:

    long bits = BitConverter.DoubleToInt64Bits(value);
    if (value > 0)
        return BitConverter.Int64BitsToDouble(bits - 1);
    else if (value < 0)
        return BitConverter.Int64BitsToDouble(bits + 1);
    else
        return -double.Epsilon;
    

    IEEE-754 formats were designed so that the bits that make up the exponent and mantissa together form an integer that has the same ordering as the floating-point numbers. So, to get the largest smaller number, you can subtract one from this number if the value is positive, and you can add one if the value is negative.

    The key reason why this works is that the leading bit of the mantissa is not stored. If your mantissa is all zeros, then your number is a power of two. If you subtract 1 from the exponent/mantissa combination, you get all ones and you'll have to borrow from the exponent bits. In other words: you have to decrement the exponent, which is exactly what we want.