- 当前位置：首页 >域名 >浮点运算潜在的结果不一致问题
浮点运算潜在的结果不一致问题
发布时间：2025-11-04 00:06:45 来源：云智核 作者：域名
昨天阿楠发现了项目中的浮点一个 bug ，是运算因为浮点运算的前后不一致导致的。明明是结果完全相同的 C 代码，参数也严格一致，问题但是浮点计算出了不相同的结果。我对这个现象非常感兴趣，运算仔细研究了一下成因。结果
原始代码比较繁杂。问题在弄清楚原理后，浮点我简化了出问题的运算代码，重现了这个问题：
static void foo(float x) { float xxx = x * 0.01f; printf("%dn",结果 (int)(x * 0.01f)); printf("%dn", (int)xx); } int main() { foo(2000.0f); return 0; }使用 gcc 4.9.2 ，强制使用 x87 浮点运算编译运行，问题你会发现令人诧异的浮点结果。
gcc a.c -mfpmath=387 19 20 前一次的b2b供应网运算输出是 19 ，后一次是结果 20 。这是为什么呢?让我们来看看 gcc 生成的代码，我截取了相关的段落：
flds 16(%rbp) flds .LC0(%rip) fmulp %st, %st(1) fstps -4(%rbp) ; 1. x * 0.01f 结果保存到内存中的 float 变量中 flds 16(%rbp) flds .LC0(%rip) fmulp %st, %st(1) fisttpl -20(%rbp) ; 2. x * 0.01f 结果直接转换为整型 movl -20(%rbp), %eax movl %eax, %edx leaq .LC1(%rip), %rcx call printf flds -4(%rbp) ; 3. 读出 1. 保存的乘法结果 fisttpl -20(%rbp) movl -20(%rbp), %eax movl %eax, %edx leaq .LC1(%rip), %rcx call printf这里我做了三行注释。
首先，0.01 是无法精确表示成 2 进制的，所以 * 0.01 这个操作一定会存在误差。
两次运算都是 x * 0.01f ，虽然按 C 语言的转换规则，表达式中都是 float 时，按 float 精度运算。站群服务器但这里 gcc 生成的代码并没有严格设置 FPU 的精度控制，在注释 2 这个地方，乘法结果是直接从浮点寄存器转换为整数的。而在注释 1 这个地方，把乘法结果通过 fstps 以低精度形式保存到内存，再在注释 3 的地方 flds 读回。
所以在注释 2 和注释 3 的地方，浮点寄存器 st 内的值其实是有差别的，这导致了 fisttpl 转换为整数后结果不同。
原文链接：https://blog.codingnow.com/2017/07/float_inconsistence.html#more
【本文为专栏作者“云风”的原创稿件，转载请通过联系原作者获取授权】
戳这里，看该作者更多好文
企商汇亿华云IT资讯网香港云服务器IT技术网服务器租用源码库源码下载益华科技运维纵横亿华云计算亿华智慧云亿华灵动码上建站益强IT技术网益华IT技术论坛益强数据堂极客编程益强智未来IT资讯网益强科技码力社益强前沿资讯云站无忧益强编程堂益强编程舍思维库云智核益强资讯优选益强智囊团益华科技极客码头创站工坊多维IT资讯益强科技编程之道技术快报全栈开发亿华云亿华科技科技前瞻亿华智造汇智坊益华科技亿华互联
- Copyright © 2025 Powered by 浮点运算潜在的结果不一致问题,云智核 滇ICP备2023006006号-28sitemap