[C++]直接访问float变量内存的类,addEpison和subEpison、almostEqual约等于,及浮点数排序函数。从此忘记0.000001 再也不要使用FLT_EPSILON!

探究浮点数奥秘,这里给出一个直接访问float变量内存的类,二进制兼容float,
并给出addEpison和subEpison两个函数
而这两个函数,是模糊比较所需要的方法的基础。

//负数时随着内存值int变大,浮点值绝对值越大
// 80000000   →→→ bf800000  →  ff7fffff  ff800000  ff800001 ffffffff
//    -0       渐变     -1    渐变   -MAX    -INF       -Nan      -Nan
// 
// 正数时随着内存值int变大,浮点值绝对值越大
// 00000000   →→→ 3f800000  →  ff7fffff  7f7fffff  7f800001 7fffffff
//    +0       渐变     +1    渐变   +MAX    +INF       +Nan      +Nan

//从负到正遍历浮点数的方法:  ff800000 -> 80000000, 00000000 -> 7f7fffff
//                              -INF         -0        +0         +INF

struct Float
{
    union
    {
        struct
        {
            unsigned int Mantissa : 23;
            unsigned int Exponent : 8;
            unsigned int Sign : 1;
        };
        struct 
        {
            unsigned int withoutSign : 31;
            unsigned int : 1;
        };
        int _memoryInt = 0;
        float _float;
    };
    operator float&() { return _float; }
    operator const float&() const { return _float; }
    operator float() const { return _float; }
    Float() = default;
    Float(float f) : _memoryInt((int&)f){}

    bool equalZero() { return withoutSign == 0; }
    bool isNan() { return Exponent == 0xFF && Mantissa != 0; }
    //是正无穷大
    bool isINF_P() { return _memoryInt == 0x7F800000; }
    //是负无穷大
    bool isINF_N() { return _memoryInt == 0xFF800000; }
    //是无穷大
    bool isINF_PN() { return (_memoryInt & 0x7F800000) == 0x7F800000; }
    //非数值
    void makeNan() { _memoryInt = 0x7FFFFFFF; } //正Nan中内存值最大的
    //正无穷大
    void makeINF_P() { _memoryInt = 0x7F800000; }
    //负无穷大
    void makeINF_N() { _memoryInt = 0xFF800000; }

    void makeMax()      { _memoryInt = 0x7F7FFFFF; }
    void makeMax_N()    { _memoryInt = 0xFF7FFFFF; }

    void makeMin()      { _memoryInt = 0x00800000; }    //最小值1.1754943510e-38(0x00800000), 次小值1.175494491e-38(0x00800001)
    void makeMin_N()    { _memoryInt = 0x80800000; }

    void makeTrueMin() { _memoryInt = 0x00000001; }     //最小值1.4012984643e-45F(0x00000001),次小值2.803e-45  2.8025969286496341e-45 (0x00800002)
    void makeTrueMin_N() { _memoryInt = 0x80000001; }


    //增加一个极小值,使得比当前值大;+max,+/-INF,+/-Nan不受影响
    void addEpison()
    {
        if ((unsigned int)_memoryInt >= 0xFF800000u) return; //负无穷大,负Nan
        if (_memoryInt >= 0x7F7FFFFF) return; //正最大值,正无穷大,正Nan
        if (_memoryInt == 0x80000000)   //负0
            _memoryInt = 0x00000001;    //TrueMin
        else
        {
            _memoryInt += 1 - ((_memoryInt<0)<<1);
        }
    }
    //减少一个极小值,是的比当前值小
    void subEpison() 
    {
        if ((unsigned int)_memoryInt >= 0xFF800000u) return; //负最大值,负无穷大,负Nan
        if (_memoryInt >= 0x7F800000) return; //正无穷大,正Nan
        if (_memoryInt == 0x00000000)   //正0
            _memoryInt = 0x80000001;    //TrueNegMin
        else
        {
            _memoryInt -= 1 - ((_memoryInt < 0) << 1);
        }
    }
    bool almostEqual(const Float& other, unsigned int episonCount) const
    {
        if (((unsigned int)_memoryInt >= 0xFF800000u) 
        || (_memoryInt >= 0x7F800000)
        || ((unsigned int)other._memoryInt >= 0xFF800000u) 
        || (other._memoryInt >= 0x7F800000))
            return _float == other._float; //无穷大,Nan

        int absA = _memoryInt < 0 ? 0x80000000 - _memoryInt : _memoryInt;
        int absB = other._memoryInt < 0 ? 0x80000000 - other._memoryInt : other._memoryInt;
        return (episonCount + unsigned int(absA - absB)) <= (episonCount << 1);
    }
};
    Float a = 1.f / 3;
    Float b = 0.5f - 1.f / 6;
    a += a;
    b += b;
    {
        bool test = a == b;
        puts(test ? "true" : "false");  //false
    }
    {
        bool test = a.almostEqual(b, 1); //一个精度差距?
        puts(test ? "true" : "false");  //true
    }
    {
        Float A = a;
        A.addEpison(); //手动增加一个精度差距
        bool test = A.almostEqual(b, 1); //一个精度差距?
        puts(test ? "true" : "false");  //false
    }
    {
        Float A = a;
        A.addEpison(); //手动增加一个精度差距
        bool test = A.almostEqual(b, 2); //二个精度差距?
        puts(test ? "true" : "false");  //true
    }
    {
        Float x1 = a, x2 = a;
        x1.addEpison();
        x2.subEpison();
        bool test2 = a == b || x1 == b || x2 == b;
        puts(test2 ? "true" : "false"); //true
    }

下面给出一个浮点数排序函数,使用他能正确使包含无穷大、Nan的浮点数组成功排序

[](float a, float b) 
    {
        bool isnanA = isnan(a), isnanB = isnan(b);
        if (isnanA || isnanB)
        {
            if (!isnanB)
            {
                return ((int&)a & 0x80000000) != 0;
            }
            if (!isnanA)
            {
                return ((int&)b & 0x80000000) == 0; 
            }
            if (((int&)a < 0) != ((int&)b < 0))
                return (int&)a < (int&)b;
            else //两nan符号相同
                return (((int&)a & 0x7FFFFF) < ((int&)b & 0x7FFFFF)) != ((int&)a < 0);
        }

        if ((int&)a == (int&)b)
            return false;
        else if (a < b)
            return true;
        else if (b < a)
            return false;
        else
            return (int&)a < (int&)b;
    }
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值