硬编码双精度数的最低有效字节是一个好的舍入策略吗?

我有一个函数做一些数学计算,然后返回 double . 在windows和android下,由于 std::exp 实现开始不同( Why do I get platform-specific result for std::exp? )。E-17舍入差被传播,最终它不仅仅是我得到的舍入差(结果最终可能会改变2.36到2.47)。当我将结果与一些期望值进行比较时,我希望此函数在所有平台上返回相同的结果。

所以我需要对我的结果进行四舍五入。要做到这一点,最简单的解决办法显然是(就我在网上所能找到的)做到这一点。 std::ceil(d*std::pow<double>(10,precision))/std::pow<double>(10,precision) . 但是,我觉得根据平台的不同,最终还是会有不同的结果(而且,很难决定 precision 应该是)。

我想知道是否硬编码 双重的 可能是一个很好的舍入策略。

这个快速测试似乎表明“是的”:

#include <iostream>
#include <iomanip>

double roundByCast( double d )
{
    double rounded = d;
    unsigned char* temp = (unsigned char*) &rounded;
    // changing least significant byte to be always the same
    temp[0] = 128;
    return rounded;
}

void showRoundInfo( double d, double rounded )
{
    double diff = std::abs(d-rounded);
    std::cout << "cast: " << d << " rounded to " << rounded << " (diff=" << diff << ")" << std::endl;
}

void roundIt( double d )
{
    showRoundInfo( d, roundByCast(d) );
}

int main( int argc, char* argv[] )
{
    roundIt( 7.87234042553191493141184764681 );
    roundIt( 0.000000000000000000000184764681 );
    roundIt( 78723404.2553191493141184764681 );
}

此输出:

cast: 7.87234 rounded to 7.87234 (diff=2.66454e-14)
cast: 1.84765e-22 rounded to 1.84765e-22 (diff=9.87415e-37)
cast: 7.87234e+07 rounded to 7.87234e+07 (diff=4.47035e-07)

我的问题是:

• 是 unsigned char* temp = (unsigned char*) &rounded 这里安全还是有未定义的行为?为什么?
• 如果没有ub(或者如果没有ub有更好的方法可以做到这一点),那么对于所有输入来说,这样一个圆函数是安全和准确的吗?

注意:我知道浮点数不准确。请不要标记为的副本 Is floating point math broken? 或 Why Are Floating Point Numbers Inaccurate? .我理解为什么结果不同,我只是在寻找一种方法,使它们在所有目标平台上都相同。

编辑,我可以重新表述我的问题,因为人们问我为什么我有不同的价值观,为什么我希望他们是一样的。

假设你得到一个 双重的 由于平台特定的实现(如 性病:: )如果你想修复那些不同的 双重的 为了在所有平台上拥有完全相同的内存表示(1),并且您希望尽可能降低最低的精度,那么,固定最低有效字节是一种好方法吗?(因为我觉得四舍五入到任意给定的精度可能会比这个技巧释放更多的信息)。

(1)通过“相同的表示”,我的意思是如果你把它转换成 std::bitset ,您希望看到所有平台的相同位序列。