From: Benjamin Auder <benjamin.auder@somewhere>
Date: Sat, 14 Jan 2017 03:50:24 +0000 (+0100)
Subject: revert src/test/utils.c
X-Git-Url: https://git.auder.net/variants/Dynamo/%7B%7B%20asset%28%27mixstore/js/R.css?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 <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+#include <string.h>
+
+// 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<size; i++)
+	{
+		float error = isinteger
+			? fabs(((int*)array)[i] - ((int*)refArray)[i])
+			: fabs(((float*)array)[i] - ((float*)refArray)[i]);
+		if (error >= 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<n; u++)
+	{
+		// read number (as a string)
+		int bufferIndex = 0;
+		while (!feof(file) && curChar!=' ')
+		{
+			curChar = fgetc(file);
+			bufferNum[bufferIndex++] = curChar;
+		}
+		bufferNum[bufferIndex] = 0;
+		// transform string into float, and store it at appropriate location
+		if (isinteger)
+			((int*)array)[u] = atoi(bufferNum);
+		else
+			((float*)array)[u] = atof(bufferNum);
+		// position to next non-separator character
+		curChar = fgetc(file);
+	}
+	fclose(file);
+
+	return array;
+}
+
+int* readArray_int(const char* fileName)
+{
+	return (int*)readArray(fileName, 1);
+}
+
+float* readArray_real(const char* fileName)
+{
+	return (float*)readArray(fileName, 0);
+}
+
+int read_int(const char* fileName)
+{
+	return readArray_int(fileName)[0];
+}
+
+float read_real(const char* fileName)
+{
+	return readArray_real(fileName)[0];
+}
diff --git a/src/test/utils.h b/src/test/utils.h
new file mode 100644
index 0000000..cbffc24
--- /dev/null
+++ b/src/test/utils.h
@@ -0,0 +1,17 @@
+// Check if array == refArray
+void compareArray(const char* ID, const void* array, const void* refArray, int size, int isInteger);
+
+void compareArray_real(const char* ID, const void* array, const void* refArray, int size);
+
+void compareArray_int(const char* ID, const void* array, const void* refArray, int size);
+
+// Read array by columns (as in MATLAB) and return by-rows encoding
+void* readArray(const char* fileName, int isInteger);
+
+int* readArray_int(const char* fileName);
+
+float* readArray_real(const char* fileName);
+
+int read_int(const char* fileName);
+
+float read_real(const char* fileName);