matrix, C++ (gcc) - rextester

Language: Layout:

#include <stdio.h>
#include <omp.h>
#include <math.h>
#define N 4096
using namespace std;
int main()
{
   int i, j, k, size;
   double t1, t2, s;
/* Выделение памяти под матрицы */
   double **a=(double**)malloc(sizeof(double*)*N);
   for(i=0; i<N; i++)
      a[i]=(double*)malloc(sizeof(double)*N);
   double **b=(double**)malloc(sizeof(double*)*N);
   for(i=0; i<N; i++)
      b[i]=(double*)malloc(sizeof(double)*N);
   double **c=(double**)malloc(sizeof(double*)*N);
   for(i=0; i<N; i++)
      c[i]=(double*)malloc(sizeof(double)*N);
// Инициализация матриц
   for(i=0; i<N; i++)
      for(j=0; j<N; j++){
         a[i][j]=(i+j)*0.0001e0;
         b[i][j]=(i-j)*0.0001e0;
         c[i][j]=0.0e0;
   }
   t1=omp_get_wtime();
// Основной вычислительный блок
#pragma omp parallel shared(a, b, c) private(i, j, k)
   {
      size=omp_get_num_threads();
#pragma omp for schedule(static)
      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];
   }
   t2=omp_get_wtime();
/* Вычисление контрольной суммы и печать результата */
   for(i=0, s=0.0e0; i<N; i++)
      for(j=0; j<N; j++)
         s+=c[i][j];
   printf("N= %d, Nproc=%d, Sum=%lf, Time=%lf\n",
           N, size, s, t2-t1);
}

Show compiler warnings [ + ] Compiler args [ + ] Show input
Compilation time: 0.33 sec, absolute running time: 1.42 sec, cpu time: 5.38 sec, memory peak: 386 Mb, absolute service time: 1,86 sec