prepare structure for R package

[valse.git] / src / sources / utils / ioutils.h

diff --git a/src/sources/utils/ioutils.h b/src/sources/utils/ioutils.h
new file mode 100644 (file)
index 0000000..31f464a
--- /dev/null
+++ b/src/sources/utils/ioutils.h
@@ -0,0 +1,75 @@
+#ifndef select_ioutils_H
+#define select_ioutils_H
+
+#include <stdlib.h>
+#include <math.h>
+#include <stdint.h>
+#include <uchar.h> //for type wchar16_t
+
+// Include header for mwSize type
+#ifdef Octave
+#include <mex.h>
+#else
+#include <tmwtypes.h>
+#endif
+
+// CHUNK_SIZE = number of lambda values to be treated sequentially by a single core
+#define CHUNK_SIZE 1
+
+// integer type chosen in MATLAB (to be tuned: 32 bits should be enough)
+typedef int64_t Int;
+
+// real number type chosen in MATLAB (default: double)
+typedef double Real;
+
+#ifndef M_PI
+#define M_PI 3.141592653589793 
+#endif
+
+// Fill an array with zeros
+#define zeroArray(array, size)\
+{\
+       for (Int u=0; u<size; u++)\
+               array[u] = 0;\
+}
+
+// Copy an 1D array
+#define copyArray(array, copy, size)\
+{\
+       for (Int u=0; u<size; u++)\
+               copy[u] = array[u];\
+}
+
+// Check if array == refArray
+void compareArray(const char* ID, const void* array, const void* refArray, mwSize size, int isInteger);
+
+#define compareArray_int(ID, array, refArray, size)\
+       compareArray(ID, array, refArray, size, 1)
+#define compareArray_real(ID, array, refArray, size)\
+       compareArray(ID, array, refArray, size, 0)
+
+// Auxiliary to convert from ours ("by-rows") encoding to MATLAB
+void* brToMatlabArray(const void* brArray, const mwSize* dimensions, int nbDims, int isInteger);
+
+#define brToMatlabArray_int(brArray, dimensions, nbDims)\
+       (Int*)brToMatlabArray(brArray, dimensions, nbDims, 1)
+#define brToMatlabArray_real(brArray, dimensions, nbDims)\
+       (Real*)brToMatlabArray(brArray, dimensions, nbDims, 0)
+
+// Auxiliary to convert from MATLAB encoding to ours ("by-rows")
+void* matlabToBrArray(const void* matlabArray, const mwSize* dimensions, int nbDims, int isInteger);
+
+#define matlabToBrArray_int(matlabArray, dimensions, nbDims)\
+       (Int*)matlabToBrArray(matlabArray, dimensions, nbDims, 1)
+#define matlabToBrArray_real(matlabArray, dimensions, nbDims)\
+       (Real*)matlabToBrArray(matlabArray, dimensions, nbDims, 0)
+
+// Read array by columns (as in MATLAB) and return by-rows encoding
+void* readArray(const char* fileName, const mwSize* dimensions, int nbDims, int isInteger);
+
+#define readArray_int(fileName, dimensions, nbDims)\
+       (Int*)readArray(fileName, dimensions, nbDims, 1)
+#define readArray_real(fileName, dimensions, nbDims)\
+       (Real*)readArray(fileName, dimensions, nbDims, 0)
+
+#endif