optimizedParallel2.c

#include <stdio.h>
#include <omp.h>
#include <stdlib.h>
#include <time.h>

double **multiplyMatrix(int number, double **matrix1, double **matrix2, double **matrix3);

void fillValues(int number, double **matrix1, double **matrix2);

void runApp(int number);


int main() {

    srand(time(NULL));  // get random number
    for (int i = 200; i <= 2000; i=i+200) {
        runApp(i);
    }

    return 0;
}


void runApp(int number) {
    // Input matrix
    double **matrix1 = (double **) malloc(sizeof(double *) * number);
    double **matrix2 = (double **) malloc(sizeof(double *) * number);
    for (int i = 0; i < number; i++)
        matrix1[i] = (double *) malloc(sizeof(double) * number);
    for (int i = 0; i < number; i++)
        matrix2[i] = (double *) malloc(sizeof(double) * number);

    // Output matrix
    double **matrix3 = (double **) malloc(sizeof(double *) * number);
    for (int i = 0; i < number; i++) {
        matrix3[i] = (double *) malloc(sizeof(double) * number);
    }

    omp_set_num_threads(8);
    fflush(stdout);
    fillValues(number, matrix1, matrix2);


    double sum;
    for (int i=0;i<100;i++) {
        double startTime = omp_get_wtime();
        multiplyMatrix(number, matrix1, matrix2, matrix3);
        double endTime = omp_get_wtime();
        sum = sum + (endTime - startTime);
    }
    printf("Average for n = %d is %f \n",number,sum/100);
}


void fillValues(int number, double **matrix1, double **matrix2) {
    // Fill values
    for (int i = 0; i < number; i++) {
        for (int j = 0; j < number; j++) {
            matrix1[i][j] = (float) rand() / (float) (RAND_MAX / 100);
            matrix2[i][j] = (float) rand() / (float) (RAND_MAX / 100);
        }
    }

}


double **multiplyMatrix(int number, double **matrix1, double **matrix2, double **matrix3) {
    // Multiplication
    double **matrix4 = (double **) malloc(sizeof(double *) * number);
    for (int i = 0; i < number; i++) {
        matrix4[i] = (double *) malloc(sizeof(double) * number);
    }
    #pragma omp parallel for
    for (int i = 0; i < number; i++) {
        #pragma omp parallel for
        for (int j = 0; j < number; j++) {
            matrix4[i][j] = matrix2[j][i];
        }
    }

    #pragma omp parallel for
    for (int i = 0; i < number; i++) {
        for (int j = 0; j < number; j++) {
            double count = 0;

            for (int m = 0; m < number; m++) {
                count += matrix1[i][m] * matrix4[j][m];
            }
            matrix3[i][j] = count;
        }
    }
    return matrix3;
}