From: Benjamin Auder Date: Sat, 14 Jan 2017 03:50:24 +0000 (+0100) Subject: revert src/test/utils.c X-Git-Url: https://git.auder.net/%7B%7B%20asset%28%27mixstore/images/scripts/current/doc/%7B%7B%20targetUrl%20%7D%7D?a=commitdiff_plain;h=2d22ed81759077a555dd2ceaed41746098dd3c3a;p=valse.git revert src/test/utils.c --- diff --git a/src/test/utils.c b/src/test/utils.c new file mode 100644 index 0000000..a6370b0 --- /dev/null +++ b/src/test/utils.c @@ -0,0 +1,111 @@ +#include +#include +#include +#include + +// Check if array == refArray +void compareArray(const char* ID, const void* array, const void* refArray, int size, + int isinteger) +{ + float EPS = 1e-5; //precision + printf("Checking %s\n",ID); + float maxError = 0.0; + for (int i=0; i= maxError) + maxError = error; + } + if (maxError >= EPS) + printf(" Inaccuracy: max(abs(error)) = %g >= %g\n",maxError,EPS); + else + printf(" OK\n"); +} + +void compareArray_real(const char* ID, const void* array, const void* refArray, int size) +{ + return compareArray(ID, array, refArray, size, 0); +} + +void compareArray_int(const char* ID, const void* array, const void* refArray, int size) +{ + return compareArray(ID, array, refArray, size, 1); +} + +// Read array by columns (as in MATLAB) and return by-rows encoding +void* readArray(const char* fileName, int isinteger) +{ + // need to prepend '../data/' (not really nice code...) + char* fullFileName = (char*)calloc(5+strlen(fileName)+1, sizeof(char)); + strcat(fullFileName, "data/"); + strcat(fullFileName, fileName); + + // first pass to know how many elements to allocate + char* command = (char*)calloc(12+strlen(fullFileName)+8+1, sizeof(char)); + strcat(command, "grep -o ' ' "); + strcat(command, fullFileName); + strcat(command, " | wc -l"); + FILE *countSpaces = popen(command, "r"); + char* buffer = (char*)calloc(32, sizeof(char)); + fgets(buffer, sizeof(buffer), countSpaces); + int n = atoi(buffer) + 1; + free(buffer); + pclose(countSpaces); + + // open file for reading + FILE* file = fopen(fullFileName, "r"); + free(fullFileName); + + int d = 1; + size_t elementSize = isinteger + ? sizeof(int) + : sizeof(float); + + // read all values, and convert them to by-rows matrices format + void* array = malloc(n*elementSize); + char curChar = ' '; + char bufferNum[64]; + for (int u=0; u