code1.c

#include <omp.h>

#define SIZE 1000
double A[SIZE][SIZE];
double B[SIZE][SIZE];
double C[SIZE][SIZE];

int main() {
    // Initialize A and B
    omp_set_num_threads(4);

    // Parallélisation
    #pragma omp parallel for
    for(int i = 0; i < SIZE; i++) {
        for(int j = 0; j < SIZE; j++) {
            C[i][j] = 0.0;
            for(int k = 0; k < SIZE; k++) {
                C[i][j] += A[i][k] * B[k][j];
            }
        }
    }

    return 0;
}