7 // Check if array == refArray
8 void compareArray(const char* ID
, const void* array
, const void* refArray
, int size
,
11 float EPS
= 1e-5; //precision
12 printf("Checking %s\n",ID
);
14 for (int i
=0; i
<size
; i
++)
16 float error
= isinteger
17 ? fabs(((int*)array
)[i
] - ((int*)refArray
)[i
])
18 : fabs(((float*)array
)[i
] - ((float*)refArray
)[i
]);
19 if (error
>= maxError
)
23 printf(" Inaccuracy: max(abs(error)) = %g >= %g\n",maxError
,EPS
);
28 void compareArray_real(const char* ID
, const void* array
, const void* refArray
, int size
)
30 return compareArray(ID
, array
, refArray
, size
, 0);
33 void compareArray_int(const char* ID
, const void* array
, const void* refArray
, int size
)
35 return compareArray(ID
, array
, refArray
, size
, 1);
38 // Read array by columns (as in MATLAB) and return by-rows encoding
39 void* readArray(const char* fileName
, int isinteger
)
41 // need to prepend 'data/' (not really nice code...)
42 char* fullFileName
= (char*)calloc(5+strlen(fileName
)+1, sizeof(char));
43 strcat(fullFileName
, "data/");
44 strcat(fullFileName
, fileName
);
46 // first pass to know how many elements to allocate
47 char* command
= (char*)calloc(12+strlen(fullFileName
)+8+1, sizeof(char));
48 strcat(command
, "wc -l ");
49 strcat(command
, fullFileName
);
50 FILE *arraySize
= popen(command
, "r");
51 char* bufferNum
= (char*)calloc(64, sizeof(char));
52 fgets(bufferNum
, sizeof(bufferNum
), arraySize
);
53 int n
= atoi(bufferNum
);
56 // open file for reading
57 FILE* arrayFile
= fopen(fullFileName
, "r");
60 // read all values, and convert them to by-rows matrices format
61 size_t elementSize
= isinteger
? sizeof(int) : sizeof(float);
62 void* array
= malloc(n
*elementSize
);
63 for (int i
=0; i
<n
; i
++)
65 fgets(bufferNum
, 64, arrayFile
);
66 // transform buffer content into float or int, and store it at appropriate location
68 ((int*)array
)[i
] = atoi(bufferNum
);
70 ((float*)array
)[i
] = atof(bufferNum
);
78 int* readArray_int(const char* fileName
)
80 return (int*)readArray(fileName
, 1);
83 float* readArray_real(const char* fileName
)
85 return (float*)readArray(fileName
, 0);
88 int read_int(const char* fileName
)
90 int* array
= readArray_int(fileName
);
96 float read_real(const char* fileName
)
98 Real
* array
= readArray_real(fileName
);