Pasted

// prir lab2.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
#define LOOPS 50
 
double** StworzMacierz(int rows, int columns)
{
     double **macierz;
     macierz = (double**)malloc(rows * sizeof(double*));
     for (int i = 0; i < rows; i++)
         macierz[i] = (double*)malloc(columns * sizeof(double)); 
     return macierz;
}
 
int _tmain(int argc, _TCHAR* argv[])
{
     double **macierz;
     double *v1, *v2;
     long rows = 5000, columns = 5000;
     time_t begin_t, end_t;
     long i, j, k;
 
     v1 = (double*)malloc(columns * sizeof(double));
     v2 = (double*)malloc(rows * sizeof(double));
 
     macierz = StworzMacierz(rows,columns);
 
     for (i = 0; i < columns; i++)
     {
         v1[i] = (double)(i / 1000.0f);
         v2[i] = (double)(i / 1000.0f);
         for ( j = 0; j < rows; j++ )
             macierz[j][i] = (double)(i * j / 1000.0f);
     }
 
     begin_t = time(NULL); 
 
 
     for (k = 0; k < LOOPS; k++)
     {
        #pragma omp parallel for shared(macierz, v1, v2, rows, columns) private(i, j)
         for (i = 0; i < rows; i++)
         {
            // v2[i] = 0.0f;
             for (j = 0; j < columns; j++)
                 v2[i] = macierz[i][j] * v1[j];
         }
     }
 
      end_t = time(NULL);
      printf("Czas obliczen: %f.\n", difftime(end_t, begin_t));
 
      free(v1);
      free(v2);
      for (i = 0; i < rows; i++)
          free(macierz[i]);
      free(macierz);
 
 
    return 0;
}