2017-04-21 31 views
-1

這是我的代碼,應該用蒙特卡羅方法計算Pi。我們輸入:thread_count - 處理器線程的數量,n - 隨機生成點的數量。 這裏是我的代碼如下蘭特的C++/OpenMP錯誤

using namespace std; 

int main (int argc, char* argv[]) { 
/*sprawdzanie danych: thread_count - liczba wątków, n - liczba punktów*/ 
if (argc != 3) { 
    cout << "Co Ty piszesz!. Ma być: ./pi <thread_count> <n>" << endl; 
    exit(-1); 
} 

/*Sprawiamy by liczby były całkowite*/ 
int thread_count = atoi(argv[1]); 
long n = atoi(argv[2]); 

/*Test wartosci liczb*/ 
if (thread_count <= 0 || n <= 0) { 
    cout << "Co Ty piszesz!. Ma byc większe od 0" << endl; 
    exit(-1); 
} 
unsigned ziarno; 
double x, y; 
long int Ustrzelone = 0; 
double Start = omp_get_wtime(); 
#pragma omp parallel default(none) private(x,y, ziarno) firstprivate(n) reduction(+:Ustrzelone) num_threads(thread_count) 
{ 
    ziarno = 25231 + 16*omp_get_thread_num(); 

    #pragma omp for schedule(dynamic) 
    for(long i = 0; i <= n; i++) { 
     x = (double) rand_r(&ziarno)/RAND_MAX * 2 - 1; 
     y = (double) rand_r(&ziarno)/RAND_MAX * 2 - 1; 

     if ((x*x) + (y*y) <= 1.0) { 
      Ustrzelone += 1; 
     } 
    } 
} 

double Stop = omp_get_wtime(); 

double czas_obliczen = 1000 * (Stop - Start); 

/*Ustawienie ilosci liczb po przecinku i wyswietlenie wyniku*/ 
cout.precision(15); 
//double pi = (double) 4*Ustrzelone/n; 
//cout << "Pi wynosi " << pi << endl; 
cout << czas_obliczen << endl; 

return 0; 

在輸出我得到的錯誤這樣的: [Error] 'rand_r' was not declared in this scope 不會離開這個多少時間所以希望你能帶給我一些結論如何使它更好。

提前致謝。

+1

你缺少一個'#包括' –

+0

都能跟得上,我不會錯過它:'#包括 的#include #include #include ' – panicagaming

+0

YouTube ** are missing it!在你的問題**中提供適當的[mcve] **。 – Zulan

回答

-1

好,因爲你正在使用名稱空間標準; (不推薦),也許嘗試:: rand_r幫助其找到合適的實現......