c++ - 加速 Monte Carlo Pi

Question

我编写了一个 C++ 程序，用于通过“将随机点放入四分之一圆中并计算它们等”来计算圆周率。现在我的程序在我看来有点慢，我想过一些改进来加快它的速度(源代码在下面)。
我的第一个想法是使用 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()调用

Answer 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 the time elapsed running a function in C++