0

我已经编写了一个 C++ 程序，用于通过“将随机点放入四分之一圆并计算它们等”来计算 pi。现在我的程序在我看来有点慢，我已经考虑过一些改进来加快它（源代码如下）。
我的第一个想法是使用 OpenMP 使其成为多线程，即将 (I) 和 (II) 之间的代码分成几个线程，这样我就有了例如近十倍的轮次，而不必等待更长时间（在八核系统上）。
我的另一个想法是使用全局变量和指针，这样我只需要复制指针而不是整数元组。缺点是（idk）？
那么，我还能做些什么来加快程序的速度呢？我主要使用 Windows，但我也可以使用 Unix/Linux。
非常感谢！

代码部分： 

    #include <cstdlib>
    #include <iostream>
    #include <tuple>
    #include <math.h>
    #include <time.h>
    #include <omp.h>
    #include <sys/time.h>

    #define RAND_MAX 32000
    #define LOOPS 1000000

    inline std::tuple<int, int> Throw_points(void)
    {

        int i = 0, j = 0;
        i = rand() % 1000;
        j = rand() % 1000;
        return std::make_tuple(i, j);
    }

    inline bool is_in_circle(std::tuple<int, int> point)
    {
        if ((pow(std::get<0>(point), 2) + pow(std::get<1>(point), 2)) <= pow(1000, 2))
            return true;
        else
            return false;
    }

    inline double pi(void)
    {
        srand(time(NULL));
        long long int in_circle = 0;
        long long int out_circle = 0;
        for (int i = 0; i < LOOPS; i++)
        {
            if (is_in_circle(Throw_points()))
                in_circle++;
            out_circle++;
        }
        return double(in_circle) / double(out_circle) * 4;
    }

呼叫方式pi()

4

3 回答 3

1

我只是玩了一下。实际上，原始帖子（包括我自己的）评论中的所有建议几乎没有任何区别。

但是，摆脱元组

inline void Throw_points(int&i, int&j)
{
    i = rand() % 1000;
    j = rand() % 1000;
}

inline bool is_in_circle(int i, int j)
{
    return i*i + j*j < 1000000;        
}

将程序加快了 5 倍。

顺便说一下，我从这里使用了 boost::progress_timer 解决方案：How to get time elapsed running a function in C++

于 2013-11-11T18:12:27.180 回答
1

观察性能。使用分析工具；这会告诉您代码在哪里花费时间。一般来说，这总是令人惊讶的。

如果你在 gcc 土地使用 gprof

于 2013-11-11T18:23:43.997 回答
0

一些随机观察：

  • 乘法平方可能比调用更快pow。特别是，您不想pow(1000,2)每次都计算常数。
  • intcounters 可能比long long- 您已经将循环数限制为可以由int.
  • 通过引用传递可能会更快。或者它可能会更慢，因为类型很小。或者它可能没有任何区别，因为该函数应该是内联的。
  • if (X) return true; else return false;而不是return X;很奇怪，但可能不会影响性能。
  • rand()对于蒙特卡洛模拟，可能不够随机；它的目的是快速，但质量不高。不幸的是，好的伪随机生成器非常慢。C++11 库有几个选项。

如果您确实将其设为多线程，请确保每个线程都有不同的随机种子；否则，他们只会重复彼此的工作。您将无法使用，rand()因为它不是线程安全的。

于 2013-11-11T18:01:12.087 回答