Làm thế nào để parallelize một cách chính xác một lồng nhau cho vòng

Question

Tôi đang làm việc với OpenMP parallelize một vô hướng lồng vòng lặp for:

double P[N][N]; 
double x=0.0,y=0.0; 

for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction(x,y); 
     y+=1; 
    } 
    x+=1; 
} 

Trong vòng này, điều quan trọng là ma trận P phải giống nhau ở cả hai vô hướng và phiên bản song song:

Tất cả các thử nghiệm khả năng của tôi đã không thành công ...

Answer 1

vấn đề ở đây là bạn đã thêm lặp-to-lặp phụ thuộc với:

x+=1; 
y+=1; 

Vì vậy, vì mã hiện đang đứng, nó không thể song song. Cố gắng làm như vậy sẽ dẫn đến kết quả không chính xác. (Như bạn có lẽ thấy)

May mắn thay, trong trường hợp của bạn, bạn có thể trực tiếp tính toán chúng mà không cần giới thiệu sự phụ thuộc này:

for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction((double)i, (double)N*i + j); 
    } 
} 

Bây giờ bạn có thể thử ném một OpenMP pragma trên này và xem nếu nó hoạt động:

#pragma omp parallel for 
for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction((double)i, (double)N*i + j); 
    } 
}