仅用C语言无符号类型提取48位乘法结果高24位的方法咨询

问题背景 #

x和y是存储IEEE 754单精度浮点数原始位的32位无符号整数，需在C语言中实现浮点数乘法，但仅允许使用unsigned int类型（禁止64位整数）。核心需求是将两个24位尾数（mx和my，已补上隐藏位1）相乘，提取48位乘积的高24位。你尝试通过拆分12位部分模拟乘法，但不确定方法是否正确。

你的实现代码 #

unsigned sx = (x >> 31) & 1;
unsigned sy = (y >> 31) & 1;
unsigned ex = (x >> 23) & 0xFF;
unsigned ey = (y >> 23) & 0xFF;
unsigned mx = (x & 0x7FFFFF) | 0x800000;
unsigned my = (y & 0x7FFFFF) | 0x800000;

unsigned sign = sx ^ sy;
unsigned exp = ex + ey - 127;

unsigned a = mx >> 12;
unsigned b = mx & 0xFFF;
unsigned c = my >> 12;
unsigned d = my & 0xFFF;

unsigned p1 = a * c;
unsigned p2 = a * d;
unsigned p3 = b * c;
unsigned p4 = b * d;

unsigned mid = (p2 + p3) << 12;
unsigned lower = mid + (p4 >> 12);
unsigned upper = (p1 << 8) + (lower >> 4);

代码问题分析 #

拆分24位尾数为12位部分的思路是可行的，但后续计算逻辑错误，无法正确提取48位乘积的高24位：

• 两个24位尾数的乘积公式为：mx*my = (a<<12 + b)*(c<<12 + d) = a*c<<24 + (a*d + b*c)<<12 + b*d
• 我们需要的高24位是(mx*my) >> 24 = a*c + ((a*d + b*c)<<12 + b*d) >> 24
• 你的代码中mid、lower、upper的计算逻辑偏离了这个公式，错误地对中间结果进行移位和截断，导致最终upper并非正确的高24位值。

正确的实现方式 #

通过拆分计算避免64位整数，正确提取高24位的代码如下：

unsigned sx = (x >> 31) & 1;
unsigned sy = (y >> 31) & 1;
unsigned ex = (x >> 23) & 0xFF;
unsigned ey = (y >> 23) & 0xFF;
unsigned mx = (x & 0x7FFFFF) | 0x800000;
unsigned my = (y & 0x7FFFFF) | 0x800000;

unsigned sign = sx ^ sy;
unsigned exp = ex + ey - 127;

// 拆分24位尾数为高12位和低12位
unsigned a = mx >> 12;
unsigned b = mx & 0xFFF;
unsigned c = my >> 12;
unsigned d = my & 0xFFF;

// 计算各部分乘积
unsigned p1 = a * c;          // 12*12=24位，对应a*c<<24的高24位
unsigned p2 = a * d;          // 12*12=24位
unsigned p3 = b * c;          // 12*12=24位
unsigned p4 = b * d;          // 12*12=24位

unsigned temp = p2 + p3;      // (a*d + b*c)，最多22位
unsigned temp_low12 = temp & 0xFFF;
unsigned temp_high = temp >> 12;

// 计算((a*d + b*c)<<12 + p4) >>24的进位值
unsigned carry = temp_high;
// 判断(temp_low12<<12 + p4)是否超过2^24，超过则加1进位
if ((temp_low12 << 12) + p4 >= 0x1000000) {
    carry += 1;
}

// 最终高24位结果
unsigned high24 = p1 + carry;

// 处理尾数溢出（如果high24是25位，说明乘积>=2^24，需要指数+1，尾数右移1位）
if (high24 & 0x1000000) {
    exp += 1;
    high24 >>= 1;
}

// 最终浮点数的位表示（注意：还需处理指数溢出/下溢等边界情况）
unsigned result = (sign << 31) | ((exp & 0xFF) << 23) | (high24 & 0x7FFFFF);

关键逻辑说明 #

1. 进位计算：carry是((a*d + b*c)<<12 + b*d)右移24位的结果，拆分为两部分计算：
   • temp_high是(a*d + b*c)的高12位，直接作为进位的基础部分
   • 检查(temp_low12<<12 + p4)是否达到2^24，若达到则额外加1进位
2. 尾数溢出处理：如果high24最高位为1（即值>=2^24），说明乘积是2.xxxx形式，需将指数加1，尾数右移1位，符合IEEE 754的规格化要求
3. 边界情况：上述代码未处理指数溢出（如exp超过0xFF）、指数下溢（如exp小于0）等极端情况，实际实现时需要补充这些逻辑

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