benchmark_clock_gettime.c

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include "computepi.h"

#define CLOCK_ID CLOCK_MONOTONIC_RAW
#define ONE_SEC 1000000000.0

#define COMPUTE_PI_ERROR(x) (((x > M_PI) ? (x - M_PI) : (M_PI - x)) / M_PI)


int main(int argc, char const *argv[])
{
    struct timespec start = {0, 0};
    struct timespec end = {0, 0};

    if (argc < 2)
        return -1;

    int N = atoi(argv[1]);
    int i, loop = 100;

#if defined(BASELINE)
    // Baseline
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_baseline(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_baseline(N)));

    // OpenMP with 2 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_openmp(N, 2);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_openmp(N, 2)));

    // OpenMP with 4 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_openmp(N, 4);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_openmp(N, 4)));

    // AVX SIMD
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_avx(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_avx(N)));

    // AVX SIMD + Loop unrolling
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_avx_unroll(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf\n", COMPUTE_PI_ERROR(compute_pi_avx_unroll(N)));
#endif

#if defined(LEIBNIZ)
    // Baseline
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_leibniz(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_leibniz(N)));


    // OpenMP with 2 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_leibniz_openmp(N, 2);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_leibniz_openmp(N, 2)));


    // OpenMP with 4 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_leibniz_openmp(N, 4);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_leibniz_openmp(N, 4)));


    // AVX SIMD
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_leibniz_avx(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_leibniz_avx(N)));


    // AVX SIMD + Loop unrolling
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_leibniz_avx_unroll(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf\n", COMPUTE_PI_ERROR(compute_pi_leibniz_avx_unroll(N)));
#endif

#if defined(EULER)
    // Baseline
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_euler(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_euler(N)));


    // OpenMP with 2 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_euler_openmp(N, 2);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_euler_openmp(N, 2)));


    // OpenMP with 4 threads
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_euler_openmp(N, 4);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_euler_openmp(N, 4)));


    // AVX SIMD
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_euler_avx(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_euler_avx(N)));


    // AVX SIMD + Loop unrolling
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_euler_avx_unroll(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf\n", COMPUTE_PI_ERROR(compute_pi_euler_avx_unroll(N)));
#endif

#if defined(CHUD)
    // Baseline
    printf("%lf ", compute_pi_Chud(N));
    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < loop; i++) {
        compute_pi_Chud(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_Chud(N)));

    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < 2; i++) {
        compute_pi_Chud(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_Chud(N)));

    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < 2; i++) {
        compute_pi_Chud(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_Chud(N)));

    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < 2; i++) {
        compute_pi_Chud(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf ", COMPUTE_PI_ERROR(compute_pi_Chud(N)));

    clock_gettime(CLOCK_ID, &start);
    for (i = 0; i < 2; i++) {
        compute_pi_Chud(N);
    }
    clock_gettime(CLOCK_ID, &end);
    printf("%lf ", (double) (end.tv_sec - start.tv_sec) +
                       (end.tv_nsec - start.tv_nsec) / ONE_SEC);
    printf("%lf\n", COMPUTE_PI_ERROR(compute_pi_Chud(N)));
#endif
    return 0;
}