vyom_lawda

// CPP Program to multiply two matrix using pthreads
#include <bits/stdc++.h>
using namespace std;

// maximum size of matrix
#define MAX 4

// maximum number of threads
#define MAX_THREAD 4


int MX = 100;

int **matA, **matB, **matC;

int step_i = 0;
int n;

void* multi(void* arg)
{
    int cr = step_i++;

    // Each thread computes 1/4th of matrix multiplication
    for (int i = cr * n / 4; i < (cr + 1) * n / 4; i++)
        for (int j = 0; j < n; j++)
            for (int k = 0; k < n; k++)
                matC[i][j] += matA[i][k] * matB[k][j];
}


// Driver Code
int main()
{
    time_t start, end;
    cin >> n;
    matA = (int**)malloc(n*sizeof(int*));
    matB = (int**)malloc(n*sizeof(int*));
    matC = (int**)malloc(n*sizeof(int*));

    for (int i = 0; i < n; i++) {
        matA[i] = (int*)malloc(n*sizeof(int));
        matB[i] = (int*)malloc(n*sizeof(int));
        matC[i] = (int*)malloc(n*sizeof(int));
        for (int j = 0; j < n; j++) {
            matA[i][j] = rand() % MX;
            matB[i][j] = rand() % MX;
        }
    }
    pthread_t threads[MAX_THREAD];

    time(&start);
    clock_t st, ed;
    st = clock();
    clock_t tic = clock();
    for (int i = 0; i < MAX_THREAD; i++) {
        int* p;
        pthread_create(&threads[i], NULL, multi, (void*)(p));
    }

    for (int i = 0; i < MAX_THREAD; i++)
        pthread_join(threads[i], NULL);
    clock_t toc = clock();
    ed = clock();
    time(&end);
    auto x =double(ed-st)/double(CLOCKS_PER_SEC);
    free(matA);
    free(matB);
    free(matC);
    cout << "n : " << n << " time:" <<(double)(toc - tic) / CLOCKS_PER_SEC << endl;

}