关于float奇怪的bug

一直以为浮点数f为NaN的时候，f不会进行+-*/等计算。
偶然间发现一个不知道算不算bug的问题，当f为NaN并且第22位（尾数左一位）为0的时候，对于f的运算会把该位变成1，同时其他位不变。

//f的二进制表示，其中x表示无所谓0或1
//x   1111 1111    0xxxx.....xxxx
//第一部分是符号位，用来表示浮点数正负
//第二部分是阶码，全为1表示浮点数是NaN或者无穷
//第三部分是尾数，当阶码全为1的时候，尾数为0浮点数表示无穷;当尾数不全为0的时候表示NaN
    printf("%x \n", f2u(u2f(0x7f800001) + 2));
    printf("%x \n", f2u(u2f(0x7f900001) - 2));
    printf("%x \n", f2u(u2f(0x7fa00001) × 2));
    printf("%x \n", f2u(u2f(0x7fb00001) / 2));

output：
    0x7fc00001 
    0x7fd00001 
    0x7fe00001 
    0x7ff00001 

奇怪的bug。