计算Pi不同的方式

时间:2017-05-30 22:12:25

标签: algorithm math trigonometry pi approximation

我知道在代码中计算pi的两种方法。

pi = 4.0 * atan(1.0)

pi = acos(-1.0)

一个与另一个的好处是什么?我知道有些语言有内置的pi表示,但这并不是在询问它们。另外,还有其他常见的计算pi的方法以及它们与其他方法的比较吗?

1 个答案:

答案 0 :(得分:2)

  

一方与另一方的好处是什么?

这些函数的设计不仅仅是为了逼近π的值。在这种情况下,我看不出明显的速度。

事实上,我在我的系统上进行了实验,在那里我看到了两个函数的相同近似值pi,并且或多或少都有相同的速度。

Georgioss-MacBook-Pro:~ gsamaras$ g++ -std=c++0x -Wall -O3 atan.cpp
Georgioss-MacBook-Pro:~ gsamaras$ ./a.out 
It took me on average 1.69e-13 seconds.
3.14159265358979312
Georgioss-MacBook-Pro:~ gsamaras$ g++ -std=c++0x -Wall -O3 acos.cpp
Georgioss-MacBook-Pro:~ gsamaras$ ./a.out 
It took me on average 1.7e-13 seconds.
3.14159265358979312

代码在循环中计算π,并将循环计数器添加到它,确保编译器不会优化它(每次迭代时都是相同的值):

Georgioss-MacBook-Pro:~ gsamaras$ cat acos.cpp
#include <iostream>
#include <ctime>
#include <ratio>
#include <chrono>
#include <cmath>
#include <limits>

#define ITER 1000000

int main ()
{
  using namespace std::chrono;

  high_resolution_clock::time_point t1 = high_resolution_clock::now();

  for(int i = 0; i < ITER; ++i)
  {
    auto pi = acos(-1.0) + i;
    pi += i + i;
  } 
  high_resolution_clock::time_point t2 = high_resolution_clock::now();

  duration<double> time_span = duration_cast<duration<double>>(t2 - t1);

  std::cout << "It took me on average " << time_span.count()/(double)ITER << " seconds.";
  std::cout << std::endl;
  auto pi = acos(-1.0);
  std::cout.precision(std::numeric_limits< double >::max_digits10);
  std::cout << std::fixed << pi << std::endl;

  return 0;
}

基于我的Time measurements (C++)atan()代码是相同的,只是函数更改。

  

此外,还有其他常见的方法来计算pi以及它们与其他人的比较方式吗?

还有许多其他近似π的方法,并且比较它们只是太宽。例如,柏拉图已经接近π,如this