X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=src%2Fsources%2Futils%2Fioutils.h;fp=src%2Fsources%2Futils%2Fioutils.h;h=31f464a75288116b162dd82c83a852cd346c0172;hb=493a35bfea6d1210c94ced8fbfe3e572f0389ea5;hp=0000000000000000000000000000000000000000;hpb=35b42a4bd37b162a3d579693b2b5fa4913a52ed5;p=valse.git

diff --git a/src/sources/utils/ioutils.h b/src/sources/utils/ioutils.h
new file mode 100644
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