Why doesn't a compiler optimize floating-point *2 into an exponent increment?
For example, 2 * f, might simply increment the exponent of f by 1, saving some cycles.
This simply isn't true.
First you have too many corner cases such as zero, infinity, Nan, and denormals. Then you have the performance issue.
The misunderstanding is that incrementing the exponent is not faster than doing a multiplication.
If you look at the hardware instructions, there is no direct way to increment the exponent. So what you need to do instead is:
- Bitwise convert into integer.
- Increment the exponent.
- Bitwise convert back to floating-point.
There is generally a medium to large latency for moving data between the integer and floating-point execution units. So in the end, this "optimization" becomes much worse than a simple floating-point multiply.
So the reason why the compiler doesn't do this "optimization" is because it isn't any faster.
On modern CPUs, multiplication typically has one-per-cycle throughput and low latency. If the value is already in a floating point register, there's no way you'll beat that by juggling it around to do integer arithmetic on the representation. If it's in memory to begin with, and if you're assuming neither the current value nor the correct result would be zero, denormal, nan, or infinity, then it might be faster to perform something like
addl $0x100000, 4(%eax) # x86 asm example
to multiply by two; the only time I could see this being beneficial is if you're operating on a whole array of floating-point data that's bounded away from zero and infinity, and scaling by a power of two is the only operation you'll be performing (so you don't have any existing reason to be loading the data into floating point registers).
Common floating-point formats, particularly IEEE 754, do not store the exponent as a simple integer, and treating it as an integer will not produce correct results.
In 32-bit float or 64-bit double, the exponent field is 8 or 11 bits, respectively. The exponent codes 1 to 254 (in float) or 1 to 2046 (in double) do act like integers: If you add one to one of these values and the result is one of these values, then the represented value doubles. However, adding one fails in these situations:
- The initial value is 0 or subnormal. In this case, the exponent field starts at zero, and adding one to it adds 2-126 (in float) or 2-1022 (in double) to the number; it does not double the number.
- The initial value exceeds 2127 (in float) or 21023 (in double). In this case, the exponent field starts at 254 or 2046, and adding one to it changes the number to a NaN; it does not double the number.
- The initial value is infinity or a NaN. In this case, the exponent field starts at 255 or 2047, and adding one to it changes it to zero (and is likely to overflow into the sign bit). The result is zero or a subnormal but should be infinity or a NaN, respectively.
(The above is for positive signs. The situation is symmetric with negative signs.)
As others have noted, some processors do not have facilities for manipulating the bits of floating-point values quickly. Even on those that do, the exponent field is not isolated from the other bits, so you typically cannot add one to it without overflowing into the sign bit in the last case above.
Although some applications can tolerate shortcuts such as neglecting subnormals or NaNs or even infinities, it is rare that applications can ignore zero. Since adding one to the exponent fails to handle zero properly, it is not usable.