3

我有以下锈代码：

use std::f64::consts as f64;

fn main() {
    println!("Checking f64 PI...");
    // f64::PI definition: https://github.com/rust-lang/rust/blob/e1fc9ff4a794fb069d670dded1a66f05c86f3555/library/core/src/num/f64.rs#L240
    println!("Definition: pub const PI: f64 = 3.14159265358979323846264338327950288_f64;");
    println!("Print it:                       {:.35}", f64::PI);
    println!("Different after 16 significant digits ----------|                         ");
    println!("##############################################################################");
    println!("Question 1: Why do the digits differ after 16 significant digits when printed?");
    println!("##############################################################################");

    println!("PERFORM ASSERTIONS..."); 
    assert_eq!(f64::PI, 3.14159265358979323846264338327950288_f64); // 36 significant digits definition
    assert_eq!(f64::PI, 3.141592653589793_f64); // 16 significant digits (less then the 36 in definition)
    // compares up to here -------------|
    assert_eq!(f64::PI, 3.14159265358979300000000000000000000_f64); // 36 significant digits (16 used in equality comparison)
    assert_ne!(f64::PI, 3.14159265358979_f64); // 15 significant digits (not equal)

    println!("PERFORM EQUALITY CHECK..."); 
    if 3.14159265358979323846264338327950288_f64 == 3.14159265358979300000000000000000000_f64 {
        println!("BAD: floats considered equal even when they differ past 16 significant digits");
        println!("######################################################################");
        println!("Question 2: Why does equality checking use only 16 significant digits?");
        println!("They are defined using 36 significant digits so why can't we perform");
        println!("an equality check with this accuracy?");
        println!("######################################################################");
    } else {
        println!("GOOD: floats considered different when they differ past 16 significant digits");
        println!("NOTE: This block won't execute :(");
    }
}

我知道浮点运算可能很棘手，但想知道这种棘手是否也会影响 f64 的打印和执行相等性检查。以下是上述代码的输出：

Checking f64 PI...
Definition: pub const PI: f64 = 3.14159265358979323846264338327950288_f64;
Print it:                       3.14159265358979311599796346854418516
Different after 16 significant digits ----------|                         
##############################################################################
Question 1: Why do the digits differ after 16 significant digits when printed?
##############################################################################
PERFORM ASSERTIONS...
PERFORM EQUALITY CHECK...
BAD: floats considered equal even when they differ past 16 significant digits
######################################################################
Question 2: Why does equality checking use only 16 significant digits?
They are defined using 36 significant digits so why can't we perform
an equality check with this accuracy?
######################################################################
4

2 回答 2

5

f64顾名思义，An以 64 位存储。在这个固定数量的存储空间中，我们只能对固定数量的数字进行编码（具体来说，其中 52 个位专用于有效位）。如果您在浮点文字中使用更多数字，则存储在f64变量中的数字将四舍五入为可用位数表示的最接近的数字。因为f64这意味着我们总是可以精确地表示 15 位十进制数字，有时是 16。这解释了为什么有时即使您在源代码中使用不同的浮点文字，数字看起来也相等：这是因为在四舍五入到最接近的可表示数字后，它们是相同的.

打印不同数字的原因是相同的。该数字在存储时四舍五入到最接近的可表示数字，并在打印时再次转换回十进制。额外的数字来源于二进制到十进制的转换，但小数点后 15 或 16 位之后的数字大多没有意义——它们不携带任何关于所表示数字的额外信息。

请注意，这些都不是 Rust 特有的。大多数现代编程语言使用 IEEE 754-1985 标准来表示浮点数，因此它们的行为相同。如果你想要任意精度的算术，你通常需要使用一些库，例如rug crate。

于 2021-01-14T13:03:36.517 回答
2

您从假设传递这些所有这些双文字开始3.14159265358979323846264338327950288_f64，3.141592653589793_f64或者3.14159265358979_f64实际上将这些确切的值分配给您的变量。这个假设是不正确的。

尽管 rust 源代码的作者使用数学常数的前 36 位实际数字来定义f64::PI，但使用 IEEE 754 浮点格式存储的实际 64 位值是不同的。根据在线转换器，最接近的 IEEE 754 64 位浮点值将是0x400921FB54442D18，可以使用转换回十进制3.1415926535897931159979634685时的数字来近似。当您将 IEEE 754 值转换为十进制数时，您会得到相同的值。0x400921FB54442D18

换句话说：

What we wanted to store: 3.14159265358979323846264338327950288 
What is actually stored: 3.14159265358979311599796346854...

也许更简单的可视化方法是想象有一个虚构的数据类型，它可以存储从 0 到 1 的正实数，并在内部使用字符串（字符数组）表示，最大长度为 12 个字符。因此，您采用这种奇异的 96 位类型并创建 5 个变量：

strdouble A = 0.333333;       // internally stored as x = { .raw = "0.333333000" }
strdouble B = 0.333333333;    // internally stored as x = { .raw = "0.333333333" }
strdouble C = 0.333333333333; // internally stored as x = { .raw = "0.333333333" }
strdouble D = 0.333333333444; // internally stored as x = { .raw = "0.333333333" }
strdouble E = 0.333333333555; // internally stored as x = { .raw = "0.333333334" }

您可以看到B,C和D将相等，尽管传递给编译器的文字值完全不同。您还可以看到像 (1/3 + 1/3 + 1/3) 这样的算术将如何返回0.999999999而不是1，因为根本没有办法表示超出最后一个原始数字的任何精度。

于 2021-01-14T13:51:12.487 回答