linux多程序程式設計計算圓周率,中值積分定理計算PI值的多線程實作

// Parallel.cpp : 定義控制台應用程式的入口點。

//

#include "stdafx.h"

#include

#include

static long num_steps = 100000;

const int numThreads = 4;

double step, pi;

CRITICAL_SECTION g_cs;

double sum = 0.0;

DWORD WINAPI countFunc(LPVOID pArg) {

double x;

int i;

int temp = *(int *)pArg;

int start = (temp*num_steps) / 4;

int end = start + num_steps / 4;

for (i = start; i < end; i++) {

EnterCriticalSection(&g_cs);

x = (i + 0.5)*step;

sum = sum + 4.0 / (1.0 + x*x);

LeaveCriticalSection(&g_cs);

}

return 0;

}

void main()

{

int i;

HANDLE hThread[numThreads];

step = 1.0 / (double)num_steps;

int tNum[numThreads];

InitializeCriticalSection(&g_cs);

for (int i = 0; i < numThreads; i++) {

tNum[i] = i;

hThread[i] = CreateThread(NULL, 0, countFunc, (LPVOID)&tNum[i], 0, NULL);

}

WaitForMultipleObjects(numThreads, hThread, TRUE, INFINITE);

DeleteCriticalSection(&g_cs);

pi = step * sum;

printf("PI = %12.9f\n", pi);

system("pause");

}