x的平方根等于x

Question

给定一个双精度浮点（非负）数x，其平方的平方根是否总是等于自身？

换句话说，如果执行以下操作，是否会降低精度：

x = <non-negative double>
y = x^2
z = sqrt(y)

这样：

x == z

我对平方变为无穷大或为零的情况不感兴趣，只对适合双精度的数字感兴趣。

Answer 1

对一个数字求平方，产生一个值是原始值位数的两倍。因此，如果 x 太大，那么x^2 中的一些位会丢失，而x 无法从y 中完全恢复[编辑：仍然可以通过适当的舍入从 y 中得到 x]。在 IEEE-754 双精度的情况下，如果 x 的有效数字部分超过 26 位，则 y 的结果将被截断。这是最简单的情况。

如果x 的有效位很少，但指数非常大或非常小，那么x^2 可能对于双精度来说太大，将变成inf 或非正规数，在这种情况下无法恢复x .

如果x 不是太大或太小，那么sqrt(y) 将等于x，因为IEEE-754 标准要求+、-、*、/ 和sqrt 适当的圆角。

Examples:

#include <iostream>
#include <ios>
#include <iomanip>
#include <cmath>

using std::fixed;
using std::hexfloat;
using std::cout;

int main() {
    double x = 1.25e155;
    double y = x*x;
    cout << hexfloat << "x = " << x << ", y = " << y << ", sqrt(y) = " << sqrt(y) << '\n';
    
    x = 1.25e-155;
    y = x*x;
    cout << hexfloat << "x = " << x << ", y = " << y << ", sqrt(y) = " << sqrt(y) << '\n';
    return 0;
}

Answer 2

#include <stdio.h>
#include <math.h>

int main(void) {
  double x = 1.0000000000000001E-160;
  double square = x*x;
  double root = sqrt(square);
  if (root != x) {
    printf("%.20g\n", x);
    printf("%.20g\n", root);
  }
}

输出

1.0000000000000001466e-160
9.9999443357584897793e-161

这里发生的情况是，x 足够大，以至于其平方非零，但又足够小，以至于其平方只能表示为非规范化数字，这会降低可用精度。

我觉得@MarkDickinson 对@LưuVĩnhPhúc 的回答的评论在很大程度上是正确的。如果x 和x*x 都是正规范化数字，那么我无法找到x != sqrt(x*x) 的示例，即使使用快速蛮力（在几个小范围内），尽管这不应被视为证明。