6

我尝试搜索有关 long double 的信息，到目前为止，我了解编译器的实现方式不同。

在 Ubuntu (XUbuntu) Linux 12.10 上使用 GCC 时，我得到这个：

double PId = acos(-1);
long double PIl = acos(-1);
std::cout.precision(100);

std::cout << "PId " << sizeof(double) << " : " << PId << std::endl;
std::cout << "PIl " << sizeof(long double)  << " : " << PIl << std::endl;

输出：

PId 8  : 3.141592653589793115997963468544185161590576171875
PIl 16 : 3.141592653589793115997963468544185161590576171875

任何人都明白为什么他们输出（几乎）相同的东西？

4

3 回答 3

8

根据acos 的参考long double，只有将 a 传递给它才会返回long double。您还必须std::acos像狒狒建议的那样使用。这对我有用：

#include <cmath>
#include <iostream>

int main() {

  double PId = acos((double)-1);
  long double PIl = std::acos(-1.0l);
  std::cout.precision(100);

  std::cout << "PId " << sizeof(double) << " :  " << PId << std::endl;
  std::cout << "PIl " << sizeof(long double)  << " : " << PIl << std::endl;
}

输出：

PId 8  : 3.141592653589793115997963468544185161590576171875
PIl 12 : 3.14159265358979323851280895940618620443274267017841339111328125

         3.14159265358979323846264338327950288419716939937510582097494459

最后一行不是输出的一部分，它包含 pi 到此精度的正确数字。

于 2012-12-16T17:44:05.090 回答
5

要获得正确的有效位数，请使用std::numeric_limits. 在 C++11 中，我们有digits10十进制有效数字（与digits给出有效位相反）。

#include <cmath>
#include <iostream>
#include <limits>

int
main()
{
  std::cout.precision(std::numeric_limits<float>::digits10);
  double PIf = acos(-1.0F);
  std::cout << "PIf " << sizeof(float) << " :  " << PIf << std::endl;

  std::cout.precision(std::numeric_limits<double>::digits10);
  double PId = acos(-1.0);
  std::cout << "PId " << sizeof(double) << " :  " << PId << std::endl;

  std::cout.precision(std::numeric_limits<long double>::digits10);
  long double PIl = std::acos(-1.0L);
  std::cout << "PIl " << sizeof(long double)  << " : " << PIl << std::endl;
}

在 x86_64 linux 上，我得到：

PIf 4 :  3.14159
PId 8 :  3.14159265358979
PIl 16 : 3.14159265358979324
于 2012-12-16T18:25:50.210 回答
4

尝试：

long double PIl = std::acos(-1.0L);

这会让你传递一个 long double 而不仅仅是一个被转换的 int。

请注意，无论如何，大多数这些数字都是垃圾。如果您将数字与真实 PI 进行比较，则使用 8 字节双精度数可以获得 15 个精度数

3.1415926535897932384626433

您会看到只有前 15 个数字适合。

如评论中所述，您可能不会获得双倍的精度，因为实现可能只使用 80Bit 表示，然后取决于它为尾数保留多少位。

于 2012-12-16T17:43:31.787 回答