OpenMP

[regina_]

Распараллелить цикл:

For(i=2;i<N;i++) 
 For(j=2;i<N;j++)
   A[i,j] =A[i-2,j] +A[i,j-2];

Не выходит распараллеливание. Результаты неверные.

#include <iostream>
#include <omp.h>
#include <iomanip>
using namespace std;


    
    int main()
    {
    	 int n=1000;
    	 int g = 1;
    	 int **A = new int*[1000];
    	 for (int i = 0; i < n; i++)
        	A[i] = new int[n+1];
    	
    	 for (int i = 0; i < n; i++) {
            for (int j = 0; j < n; j++) {
                A[i][j] = g;
                A[j][i] = g+2;
               // cout << A[i][j] << " ";
                g++;
            }
           // cout << endl;

    }

    cout << endl;
    #pragma omp parallel for num_threads(4)
    	for (int k=0; k<2; ++k) 
    	{
    		int i1 = k/2;
    		int j1 = k%2;
    		for (int i = 2+i1; i < n; ++i)
    		{
    			for (int j = 2+j1; j < n; j+=2)
    			{
    				A[i][j] = A[i - 2][j] + A[i][j - 2];
    			}
    		}
    	}  


  	for (int i = n-2; i < n; i++) {
            for (int j = n-2; j < n; j++) {
                cout << A[i][j] << " ";
            }
            cout << endl;
        }

      for (int i = 0; i < n; i++) {
            for (int j = 0; j < n; j++) {
                A[i][j] = g;
                A[j][i] = g+2;
               // cout << A[i][j] << " ";
                g++;
            }
       }
    for (int i = 2; i<n; i++)
        for (int j = 2; j<n; j++)
            A[i][j] = A[i - 2][j] + A[i][j - 2];  
	for (int i = n-2; i < n; i++) {
            for (int j = n-2; j < n; j++) {
                cout << A[i][j] << " ";
            }
            cout << endl;
        }
    for (int i = 0; i < n; i++)
        delete [] A[i];
    	delete [] A;
}

_____
Код программы нужно выделять (форматировать) тегами [CODE] (читать FAQ)
Модератор