MatrixMul2

#include <stdio.h>
#include "omp.h"

#define N 10000
#define NUM_THREADS 8

int main(){

    static int a[N][N], b[N][N], c[N][N]={0};
    int i, j, k;

    for(i=0;i<N;i++)
        for(j=0;j<N;j++){
            a[i][j] = i + j;
            b[i][j] = N * N - i - j;
        }


    double start = omp_get_wtime();

    omp_set_num_threads(NUM_THREADS);
    #pragma omp parallel 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];

    double end = omp_get_wtime();

    printf("%f",end-start);
    return 0;
}