啊,花了一个晚上,大概总结出了一些东西。先说结论吧,这算是 crt 不同而导致的,windows sdk 中的 printf 函数( cout 应该是一致的)调用的是 windows crt 的内容,默认四舍五入; mingw 使用的是自己整的一套 mingw-w64-crt,默认直接截断。
-----
tl;dr:windows crt 使用的是四舍五入,mingw crt 是直接截断。
( v2 排版可能不太好,源代码我有给定位,可以自己开着编辑器去看)
首先是 Windows 部分,进入 Windows SDK ucrt 文件夹(我的版本是 10.0.18362.0 )从 printf() 函数开始追,能够一路追到 convert/_fptostr.cpp 这个 CRT 源码,找到 `__acrt_fp_strflt_to_string` 函数 61 行就能看到四舍五入的策略,就是只要读完 precision 后还有数字,如果这个数字大于等于 5,就往前进一。这个逻辑是完全写进 CRT 的,所以我之前尝试了半天用 `fesetround()` 都没有任何用处(或者说,对 printf 不起作用)。
~~~ cpp
// Do any rounding which may be needed. Note: if digits < 0, we don't do
// any rounding because in this case, the rounding occurs in a digit which
// will not be output because of the precision requested.
if (digits >= 0 && *mantissa_it >= '5')
{
buffer_it--;
while (*buffer_it == '9')
{
*buffer_it-- = '0';
}
*buffer_it += 1;
}
~~~
p.s. 顺带感叹一句,printf 实现原来原来是状态机,DFA 牛逼。而且 windows crt 对 print 这部分似乎比较罗嗦,可能是我见识的太少。打印浮点大致的流程是:设置状态 -> 读取浮点数 -> 根据浮点数转换成高精度小数字符串(这个算法有点看不太懂,我太菜了)-> 根据状态对该字符串进行四舍五入 -> 最后处理该字符串 buffer,输出到 output 设备。不知道为啥要绕这么个弯(可能是算法问题)。
同理,不过 mingw crt 的源码默认没带,要去
https://git.code.sf.net/p/mingw-w64/mingw-w64 克隆一份,我是今天克隆的,不清楚版本(不过 crt 应该不会有太大变化);可以一路追到 mingw-w64/mingw-w64-crt/stdio/mingw_pformat.c,定位到 `__pformat_float_decimal`,处理截断的就是 1796 行的 easy mode 地方。相对 windows crt 的复杂,mingw crt 在这方面倒是比较简单,直接输出 precision 个字符,输出完就行,不做任何其他处理,简而言之就是“截断”。
~~~ cpp
if(decimal_place <= 0){ /* easy mode */
__pformat_putc( '0', stream );
points:
__pformat_emit_radix_point(stream);
for(int32_t written = 0; written < prec; written++){
if(decimal_place < 0){ /* leading 0s */
decimal_place++;
__pformat_putc( '0', stream );
/* significand */
} else if ( sig_written < max_prec ){
__pformat_putc( str_sig[sig_written], stream );
sig_written++;
} else { /* trailing 0s */
__pformat_putc( '0', stream );
}
}
} else // 后面是 hard mode,即小数部分不定长
~~~
p.p.s. mingw 的原理也是状态机,但是相对 windows 简化了许多,没有见到像 windows crt 那样使用的跳转表。大致 printf 流程为:读取 format 根据符号设置状态 -> 根据 IEEE 754 提出指数、小数、符号 -> 直接根据截断进行输出;相比之下没有 windows 那样频繁对 buffer 的操作。
-----
干,搞了一晚上,一开始以为只是单纯的 compiler 的问题,但发现无论怎么调 flag 都没用,才考虑到是 crt 的问题,果然我还是 too young。不过读大厂的代码还是挺舒服的,除了代码定位只能靠 vscode 搜索,其他该有的注释都有,还很详细,不愧是 m$。相比之下 mingw 的代码要逊一点,可能是我不太习惯 c 吧 emm