如何实现double到整数尾数与指数转换 计算二进制移位位数

frexp返回的尾数固定是取值落在[0.5, 1)区间的二进制归一化小数，形式等价于二进制的0.b₁b₂b₃…，你需要的整数尾数，本质是把这个二进制小数的所有有效位全部左移到小数点左侧，同时消去末尾无意义的0（满足ASN.1实数编码要求的唯一表示形式，即尾数为奇数），移位过程中同步调整指数值即可。

核心转换关系始终满足：

原始值 = 小数尾数 * 2^原指数 = 整数尾数 * 2^新指数

• 先处理浮点数特殊值：0、无穷、NaN，这三类值ASN.1有单独的编码规则，不需要走尾数指数转换逻辑
• 调用frexp拿到归一化小数尾数和初始指数
• 记录浮点数符号位，将尾数转为正数处理，避免符号干扰移位判断
• 循环将小数尾数乘2（等价于二进制左移1位），每次乘完将指数减1，直到尾数变为纯整数
• 循环将整数尾数右移（除2），每次移位将指数加1，直到尾数为奇数（即二进制最后一位为1，消去所有尾随0）
• 给最终整数尾数补回之前记录的符号位，输出结果即可

#include <math.h>
#include <stdint.h>
#include <stdbool.h>

// 特殊值返回码，对应ASN.1单独编码的场景
#define RET_ZERO    0
#define RET_INF_POS 1
#define RET_INF_NEG 2
#define RET_NAN     3

int double_to_asn1_real(double val, int64_t *out_mantissa, int *out_exp) {
    // 优先处理特殊值
    if (isnan(val)) return RET_NAN;
    if (isinf(val)) return val > 0 ? RET_INF_POS : RET_INF_NEG;
    if (val == 0.0) {
        *out_mantissa = 0;
        *out_exp = 0;
        return RET_ZERO;
    }

    int exp;
    double m = frexp(val, &exp);
    bool neg_flag = false;

    // 分离符号
    if (m < 0) {
        neg_flag = true;
        m = -m;
    }

    // 左移把小数部分全部移到整数位
    while (m != floor(m)) {
        m *= 2;
        exp--;
    }
    int64_t mantissa = (int64_t)m;

    // 右移消去末尾所有冗余的0，保证尾数为奇数，符合ASN.1唯一表示要求
    while (mantissa != 0 && (mantissa & 1) == 0) {
        mantissa >>= 1;
        exp++;
    }

    // 补回符号
    *out_mantissa = neg_flag ? -mantissa : mantissa;
    *out_exp = exp;
    return 0;
}

以测试值0.15625为例，执行流程如下：

1. 调用frexp得到初始小数尾数m=0.625，初始指数exp=-2
2. 第一轮移位：0.625不是整数，乘2得1.25，exp更新为-3
3. 第二轮移位：1.25不是整数，乘2得2.5，exp更新为-4
4. 第三轮移位：2.5不是整数，乘2得5.0，exp更新为-5
5. 此时5是整数，退出左移循环；检查5是奇数，不需要右移消0
6. 最终输出整数尾数5，指数-5，完全匹配预期结果5*2^-5

• 整个移位流程最多执行53次循环（double类型有效二进制位最多为53位，包含隐含的前导1），不会出现死循环
• 实现自动兼容规格化浮点数、非规格化浮点数、负浮点数场景
• 输出结果天然满足ASN.1 BER/DER编码对实数值的格式要求：尾数为整数、无尾随二进制冗余0、表示唯一

内容的提问来源于stack exchange，提问作者FaisalM