是java。郎。数学。PI等于海湾合作委员会的M_PI?

StackOverflow https://stackoverflow.com/questions/825559

题

我编码的几个基准算法在Java和C/C++。一些算法使用π。我想为两个实现每个算法产生的 相同的 结果，而不舍的不同。要做到这一点的方法之一，工作迄今为止一直是使用定义 pi 恒定是完全相同的两种语言，例如3.14159.然而，它令我愚蠢的定义pi的时候已经有了高精度常数定义在两个Java和海湾合作委员会的图书馆。

我已经花了一些时间编写快速测试程序，在寻找的文件的每一个图书馆和阅读了上浮点类型。但我还没有能够说服自己java。郎。数学。PI(或java。郎。StrictMath.PI)，或者是不是，等于M_PI在数学。h。

海湾合作委员会3.4.4(cygwin)数学。h包括：

#define M_PI            3.14159265358979323846
                                         ^^^^^

但是这个

printf("%.20f", M_PI);

产生

3.14159265358979311600
                 ^^^^^

这表明，最后5个数字不能被信任。

与此同时，Javadocs说，java。郎。数学。PI：

的 double 值，更接近于 任何其他来 pi, 比 圆周的圆其 直径。

和

public static final double PI  3.141592653589793d

其省略了问题的最后五位数字从该定。

System.out.printf("%.20f\n", Math.PI);

产生

3.14159265358979300000
                 ^^^^^

如果你有一些专门知识的浮点的数据类型，你能说服我说这些图书馆常是完全的平等？或者，他们肯定是不平等？

有帮助吗？

解决方案 4

是的，它们是相同的，并且使用它们将确保相同算法的GCC和Java实现处于相同的基础上＆＃8211;至少与使用手工定义的 pi 常量＆＃8224; 一样多。

一个警告，由 S暗示。 Lott ，GCC实现必须在 double 数据类型中保存 M_PI ，而不是 long double ，以确保等效性。 Java和GCC似乎都使用IEEE-754的64位十进制表示来表示各自的 double 数据类型。库值的字节表示（MSB到LSB），表示为 double ，可以按如下方式获得（感谢 JeeBee ）：

pi_bytes.c：

#include <math.h>
#include <stdio.h>
int main()
{
   double pi = M_PI;
   printf("%016llx\n", *((uint64_t*)&pi));
}

pi_bytes.java：

class pi_bytes
{
   public static void main(String[] a)
   {
      System.out.printf("%016x\n", Double.doubleToRawLongBits( Math.PI ) );
   }
}

同时运行：

$ gcc -lm -o pi_bytes pi_bytes.c && ./pi_bytes
400921fb54442d18

$ javac pi_bytes.java && java pi_bytes
400921fb54442d18

M_PI （作为 double ）和 Math.PI 的基础表示是相同的，直到它们的位。

＆＃8224; ＆＃8211;正如 Steve Schnepp 所述，sin，cos，exp等数学函数的输出不能保证是即使这些计算的输入按位相同，也是相同的。

其他提示

请注意以下事项。

两个数字相同，小数点后16位。这几乎是48位是相同的。

在IEEE 64位浮点数中，所有位都不是符号或指数。

#define M_PI 有21位数字;这大约是63位精度，这对于IEEE 80位浮点值很有用。

我认为您看到的是 M_PI 值中的位的普通截断。

什么你想要做的就是打印出来的原位模式的PI的价值观，并对它们进行比较。

在Java使用 http://java.sun.com/j2se/1.5.0/docs/api/java/lang/Double.html#doubleToRawLongBits(double) 法得到长久价值的，你应该打印作为二进制的。

Java5提供了：

  • PI3.141592653589793
  • 原位是4614256656552045848
  • 二是100000000001001001000011111101101010100010001000010110100011000

在C，你可以做的 double pi = M_PI; printf("%lld\n", pi); 获得同样的64位整数：4614256656552045848(谢谢Bruno).

计算相同的值非常困难，即使起始值相同。

浮点计算结果有时会从架构到另一个（例如x86 / PowerPC），从编译器到另一个（想想GCC / MS C ++），甚至使用相同的编译器，但编译选项不同。并非总是，但有时（通常在四舍五入时）。通常只是为了让问题不被注意直到太迟（想想经过多次迭代，以及许多舍入差异）

这使得跨平台多玩家游戏很难同步计算游戏状态的每次迭代（每个节点只接收输入，而不是实际的数据结构）。

因此，即使在相同的语言（C / C ++）中，结果也可能不同，从Java VM到本机主机也可能不同。

更新：

我找不到我读过的来源，但我发现了一篇论文太阳报就此事。

就像你自己回答一样， java.lang.Math.PI和GCC＆＃8217; M_PI可以设法具有相同的值。魔鬼隐藏在这些值的使用中。 IEEE没有指定数学函数的输出（sin，cos，exp，...）。因此，计算的输出不一定相同。

一个double只有52位signficand所以我认为只给你大约15个基数10位数，这可以解释为什么当你要求20位数时你有5个零。

您可以使用BigDecimal获得更高的准确度，如：

private static final BigDecimal PI = new BigDecimal(
"3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679" +
    "8214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196" +
    "4428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273" +
    "7245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094" +
    "3305727036575959195309218611738193261179310511854807446237996274956735188575272489122793818301194912" +
    "9833673362440656643086021394946395224737190702179860943702770539217176293176752384674818467669405132" +
    "0005681271452635608277857713427577896091736371787214684409012249534301465495853710507922796892589235" +
    "4201995611212902196086403441815981362977477130996051870721134999999837297804995105973173281609631859" +
    "5024459455346908302642522308253344685035261931188171010003137838752886587533208381420617177669147303" +
    "5982534904287554687311595628638823537875937519577818577805321712268066130019278766111959092164201989" +
    "3809525720106548586327886593615338182796823030195203530185296899577362259941389124972177528347913151" +
    "5574857242454150695950829533116861727855889075098381754637464939319255060400927701671139009848824012" +
    "8583616035637076601047101819429555961989467678374494482553797747268471040475346462080466842590694912" +
    "9331367702898915210475216205696602405803815019351125338243003558764024749647326391419927260426992279" +
    "6782354781636009341721641219924586315030286182974555706749838505494588586926995690927210797509302955" +
    "3211653449872027559602364806654991198818347977535663698074265425278625518184175746728909777727938000" +
    "8164706001614524919217321721477235014144197356854816136115735255213347574184946843852332390739414333" +
    "4547762416862518983569485562099219222184272550254256887671790494601653466804988627232791786085784383" +
    "8279679766814541009538837863609506800642251252051173929848960841284886269456042419652850222106611863" +
    "0674427862203919494504712371378696095636437191728746776465757396241389086583264599581339047802759009" +
    "9465764078951269468398352595709825822620522489407726719478268482601476990902640136394437455305068203" +
    "4962524517493996514314298091906592509372216964615157098583874105978859597729754989301617539284681382" +
    "6868386894277415599185592524595395943104997252468084598727364469584865383673622262609912460805124388" +
    "4390451244136549762780797715691435997700129616089441694868555848406353422072225828488648158456028506" +
    "0168427394522674676788952521385225499546667278239864565961163548862305774564980355936345681743241125"
);

public static void main(String... args) throws InterruptedException {
    System.out.println("PI to " + PI.scale() + " digits is " + PI);
    System.out.println("PI^2 to " + PI.scale() + " digits is " + 
            PI.multiply(PI).setScale(PI.scale(), BigDecimal.ROUND_HALF_UP));
}

回忆起必须在fortran中获得pi值。

由于没有常量库，我使用过 4 * atan（1。）或acos（-1。）。

不，他们不平等，他们在记忆中有不同的表现。

通常，当您想要比较2个浮点值时，不得使用==（如果是这样，则不能使用终止'等于'进行操作）。你应该使用与epsilon的比较。

double eps = 0.0000001;
if (Math.abs (Java_PI - Another_Pi) <= eps)
  System.out.println ("equals");
许可以下： CC-BY-SA 和 归因
不隶属于 StackOverflow
scroll top