bgc-c/dev/main.c

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <geometry.h>

#ifdef _WIN64
#include <windows.h>
#else
#include <time.h>
#endif // _WINDOWS_

SPVersor * allocate_versors(const unsigned int amount)
{
    return calloc(amount, sizeof(SPVersor));
}

SPVersor * make_zero_versors(const unsigned int amount)
{
    SPVersor * list = allocate_versors(amount);

    if (list == 0) {
        return 0;
    }

    for (unsigned int i = 0; i < amount; i++) {
        sp_versor_reset(&list[i]);
    }

    return list;
}

SPVersor * make_random_versors(const unsigned int amount)
{
    SPVersor * list = allocate_versors(amount);

    if (list == 0) {
        return 0;
    }

    for (unsigned int i = 0; i < amount; i++) {
        sp_versor_set(
            (2.0f * rand()) / RAND_MAX - 1.0f,
            (2.0f * rand()) / RAND_MAX - 1.0f,
            (2.0f * rand()) / RAND_MAX - 1.0f,
            (2.0f * rand()) / RAND_MAX - 1.0f,
            &list[i]
        );
    }

    return list;
}

void print_versor(const SPVersor* versor)
{
    printf("(%f, %f, %f, %f)\n", versor->_s0, versor->_x1, versor->_x2, versor->_x3);
}

void print_vector(const SPVector3* vector)
{
    printf("(%f, %f, %f) / %f\n", vector->x1,  vector->x2,  vector->x3,  sp_vector3_get_module(vector));
}
/*
int main()
{
    const unsigned int amount = 1000000;

#ifdef _WIN64
    ULONGLONG now;
    now = GetTickCount64();
    srand((unsigned int)(now & 0xfffffff));
#else
    struct timespec now;
    clock_gettime(CLOCK_REALTIME, &now);
    srand((unsigned int)(now.tv_nsec & 0xfffffff));
#endif // _WIN64

    SPVersor * versors = make_random_versors(amount);

    if (versors == 0) {
        printf("Cannot allocate memory for versors");
        return 0;
    }

    SPVector3 initial, result;

    sp_vector3_set_values(1, 2, 3, &initial);
    sp_vector3_copy(&initial, &result);

#ifdef _WIN64
    ULONGLONG start, end;
    start = GetTickCount64();
#else
    struct timespec start, end;
    clock_gettime(CLOCK_REALTIME, &start);
#endif // _WIN64
    for (unsigned int i = 0; i < amount; i++) {
        sp_versor_turn2(&versors[i], &result, &result);
    }

    for (unsigned int i = amount; i > 0; i--) {
        sp_versor_turn_back2(&versors[i - 1], &result, &result);
    }

#ifdef _WIN64
    end = GetTickCount64();

    printf("Time: %lld\n", end - start);
#else
    clock_gettime(CLOCK_REALTIME, &end);

    printf("Time: %lf\n", (end.tv_sec - start.tv_sec) * 1000.0 + (end.tv_nsec - start.tv_nsec) * 0.000001);
#endif // _WIN64


    print_vector(&initial);
    print_vector(&result);

    free(versors);
    return 0;
}
*/


int main()
{
    const unsigned int amount = 1000000;

#ifdef _WIN64
    ULONGLONG now;
    now = GetTickCount64();
    srand((unsigned int)(now & 0xfffffff));
#else
    struct timespec now;
    clock_gettime(CLOCK_REALTIME, &now);
    srand((unsigned int)(now.tv_nsec & 0xfffffff));
#endif // _WIN64

    SPVersor * versors1 = make_random_versors(amount);

    if (versors1 == 0) {
        printf("Cannot allocate memory for versors1");
        return 0;
    }

    SPVersor * versors2 = make_random_versors(amount);

    if (versors2 == 0) {
        printf("Cannot allocate memory for versors2");
        free(versors1);
        return 0;
    }

    SPVersor * results = make_zero_versors(amount);

    if (results == 0) {
        printf("Cannot allocate memory for results");
        free(versors2);
        free(versors1);
        return 0;
    }

#ifdef _WIN64
    ULONGLONG start, end;
    start = GetTickCount64();
#else
    struct timespec start, end;
    clock_gettime(CLOCK_REALTIME, &start);
#endif // _WIN64
    for (int j = 0; j < 1000; j++) {
        for (unsigned int i = 0; i < amount; i++) {
            sp_versor_combine(&versors1[i], &versors2[i], &results[i]);
        }
    }

#ifdef _WIN64
    end = GetTickCount64();

    printf("Time: %lld\n", end - start);
#else
    clock_gettime(CLOCK_REALTIME, &end);

    printf("Time: %lf\n", (end.tv_sec - start.tv_sec) * 1000.0 + (end.tv_nsec - start.tv_nsec) * 0.000001);
#endif // _WIN64

    print_versor(versors1 + 10);
    print_versor(versors2 + 10);
    print_versor(results + 10);

    free(results);
    free(versors2);
    free(versors1);
    return 0;
}