C、C++中多线程积分算pi

最新推荐文章于 2020-10-08 15:14:25 发布

oldbrostable

最新推荐文章于 2020-10-08 15:14:25 发布

阅读量904

点赞数

本文链接：https://blog.csdn.net/oldbrostable/article/details/101160572

版权

$arctanx=\int_{0}^{x} \frac{1}{1+t^{2}}dt$

#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
#include <windows.h>

volatile DWORD dwStart;
volatile int global = 0;

double test2(int num_steps) {
    int i;
    global++;
    double x, pi, sum = 0.0, step;

    step = 1.0 / (double)num_steps;

#pragma omp parallel for reduction(+:sum) private(x)
    for (i = 1; i <= num_steps; i++) {
        x = (i - 0.5) * step;
        sum = sum + 4.0 / (1.0 + x * x);
    }

    pi = step * sum;
    return pi;
}

int main(int argc, char* argv[]) {
    double   d;
    int n = 1000000;

    if (argc > 1)
        n = atoi(argv[1]);

    dwStart = GetTickCount();
    d = test2(n);
    printf_s("For %d steps, pi = %.15f, %d milliseconds\n", n, d, GetTickCount() - dwStart);

    dwStart = GetTickCount();
    d = test2(n);
    printf_s("For %d steps, pi = %.15f, %d milliseconds\n", n, d, GetTickCount() - dwStart);
}