#include <omp.h>
#include <stdio.h>
int main()
{
        int id, size;
        double t1,t2,t3;
        int  i, j, k,  N = 4, a[N][ N], b[N][N], c[N][N];
        for (i =0; i < N; i++)
        {
                for (j =0; j < N; j++)
                {
                        a[i][j] = i;
                        b[i][j] = j;

                }

        }
        t1 =  omp_get_wtime();
        #pragma omp parallel private(id, size, i, j)
        {
                id = omp_get_thread_num();
                #pragma omp  for
                for (i =0; i < N;i++)
                {
                        for (j =0; j < N; j++)
                        {
                                for (k =0; k <N ; k++)
                                {
                                        c[i][j]+= a[i][k] * b[k][j];
                                }
                                printf("Parallel c[%i][%i] = %i,  b[%i][%i] = %i,      a[%i][%i] = %i \n", i,j,c[i][j], i, j, b[i][j], i, j,  a[i][j]);
                        }
                }

        }

//      for (i =0; i < N;i++)
  //      {
    //            for (j =0; j < N; j++)                printf("After c[%i][%i]  = %i,      b[%i][%i] = %i,      a[%i][%i] = %i \n", i, j, c[i][j], i, j, b[i][j], i, j,  a[i][j]);
      //  }
        t2 = omp_get_wtime();
        t3 = t2 - t1;
        printf(" Time  t3 = %f \n", t3);
        return 0;
}