v_add.c

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

int main() {
    // Generate input vectors and destination vector
    double *x = gen_array(ARRAY_SIZE);
    double *y = gen_array(ARRAY_SIZE);
    double *z = (double*) malloc(ARRAY_SIZE*sizeof(double));

    // Test framework that sweeps the number of threads and times each run
    double start_time, run_time;
    int num_threads = omp_get_max_threads();

    for (int i = 1; i <= num_threads; i++)
    {
        omp_set_num_threads(i);
        start_time = omp_get_wtime();
        for(int j = 0; j < REPEAT; j++)
        {
            v_add_optimized_adjacent(x,y,z);
        }
        run_time = omp_get_wtime() - start_time;
        if (!verify(x, y, v_add_optimized_adjacent))
        {
            printf("v_add optimized adjacent does not match reference.\n");
            return -1;
        }
        printf("Optimized adjacent: %d thread(s) took %f seconds\n",i,run_time);
    }


    for (int i = 1; i <= num_threads; i++)
    {
        omp_set_num_threads(i);
        start_time = omp_get_wtime();
        for(int j = 0; j < REPEAT; j++)
        {
            v_add_optimized_chunks(x,y,z);
        }
        run_time = omp_get_wtime() - start_time;
        if (!verify(x, y, v_add_optimized_chunks))
        {
            printf("v_add optimized chunks does not match reference.\n");
            return -1;
        }
        printf("Optimized chunks: %d thread(s) took %f seconds\n",i,run_time);
    }

    for (int i = 1; i <= num_threads; i++)
    {
        omp_set_num_threads(i);
        start_time = omp_get_wtime();
        for(int j = 0; j < REPEAT; j++)
        {
            v_add_naive(x,y,z);
        }
        run_time = omp_get_wtime() - start_time;
        printf("Naive: %d thread(s) took %f seconds\n",i,run_time);
    }
    return 0;
}