projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
ignore valgrind core files
[valse.git]
/
src
/
test
/
test.selectiontotale.c
diff --git
a/src/test/test.selectiontotale.c
b/src/test/test.selectiontotale.c
index
45cb461
..
ebfdf19
100644
(file)
--- a/
src/test/test.selectiontotale.c
+++ b/
src/test/test.selectiontotale.c
@@
-1,72
+1,72
@@
#include "selectiontotale.h"
#include "ioutils.h"
#include "selectiontotale.h"
#include "ioutils.h"
-
mwSize main(mwSize
argc, char** argv)
+
int main(int
argc, char** argv)
{
// read dimensions
{
// read dimensions
- const
I
nt nbDims = 5;
-
I
nt* dimensions = readArray_int("dimensions",&nbDims,1);
-
mwSize
n = dimensions[0];
-
mwSize
p = dimensions[1];
-
mwSize
m = dimensions[2];
-
mwSize
k = dimensions[3];
-
mwSize
L = dimensions[4];
+ const
i
nt nbDims = 5;
+
i
nt* dimensions = readArray_int("dimensions",&nbDims,1);
+
int
n = dimensions[0];
+
int
p = dimensions[1];
+
int
m = dimensions[2];
+
int
k = dimensions[3];
+
int
L = dimensions[4];
free(dimensions);
free(dimensions);
-
mwSize
lengthOne = 1;
+
int
lengthOne = 1;
////////////
// INPUTS //
////////////
// phiInit
////////////
// INPUTS //
////////////
// phiInit
- const
mwSize
dimPhiInit[] = {p, m, k};
-
Real
* phiInit = readArray_real("phiInit",dimPhiInit,3);
+ const
int
dimPhiInit[] = {p, m, k};
+
float
* phiInit = readArray_real("phiInit",dimPhiInit,3);
// rhoInit
// rhoInit
- const
mwSize
dimRhoInit[] = {m, m, k};
-
Real
* rhoInit = readArray_real("rhoInit",dimRhoInit,3);
+ const
int
dimRhoInit[] = {m, m, k};
+
float
* rhoInit = readArray_real("rhoInit",dimRhoInit,3);
// piInit
// piInit
-
Real
* piInit = readArray_real("piInit",&k,1);
+
float
* piInit = readArray_real("piInit",&k,1);
// gamInit
// gamInit
- const
mwSize
dimGamInit[] = {n, k};
-
Real
* gamInit = readArray_real("gamInit",dimGamInit, 2);
+ const
int
dimGamInit[] = {n, k};
+
float
* gamInit = readArray_real("gamInit",dimGamInit, 2);
// min number of iterations
// min number of iterations
-
I
nt* pmini = readArray_int("mini",&lengthOne,1);
-
I
nt mini = *pmini;
+
i
nt* pmini = readArray_int("mini",&lengthOne,1);
+
i
nt mini = *pmini;
free(pmini);
// max number of iterations
free(pmini);
// max number of iterations
-
I
nt* pmaxi = readArray_int("maxi",&lengthOne,1);
-
I
nt maxi = *pmaxi;
+
i
nt* pmaxi = readArray_int("maxi",&lengthOne,1);
+
i
nt maxi = *pmaxi;
free(pmaxi);
// gamma
free(pmaxi);
// gamma
-
Real
* pgamma = readArray_real("gamma",&lengthOne,1);
-
Real
gamma = *pgamma;
+
float
* pgamma = readArray_real("gamma",&lengthOne,1);
+
float
gamma = *pgamma;
free(pgamma);
// lambda
free(pgamma);
// lambda
-
Real
* glambda = readArray_real("glambda",&L,1);
+
float
* glambda = readArray_real("glambda",&L,1);
// X
// X
- const
mwSize
dimX[] = {n, p};
-
Real
* X = readArray_real("X",dimX,2);
+ const
int
dimX[] = {n, p};
+
float
* X = readArray_real("X",dimX,2);
// Y
// Y
- const
mwSize
dimY[] = {n, m};
-
Real
* Y = readArray_real("Y",dimY,2);
+ const
int
dimY[] = {n, m};
+
float
* Y = readArray_real("Y",dimY,2);
// seuil
// seuil
-
Real
* pseuil = readArray_real("seuil",&lengthOne,1);
-
Real
seuil = *pseuil;
+
float
* pseuil = readArray_real("seuil",&lengthOne,1);
+
float
seuil = *pseuil;
free(pseuil);
// tau
free(pseuil);
// tau
-
Real
* ptau = readArray_real("tau",&lengthOne,1);
-
Real
tau = *ptau;
+
float
* ptau = readArray_real("tau",&lengthOne,1);
+
float
tau = *ptau;
free(ptau);
/////////////
free(ptau);
/////////////
@@
-74,19
+74,19
@@
mwSize main(mwSize argc, char** argv)
/////////////
// A1
/////////////
// A1
- const
mwSize
dimA[] = {p, m+1, L};
-
Int* A1 = (Int*)malloc(dimA[0]*dimA[1]*dimA[2]*sizeof(I
nt));
+ const
int
dimA[] = {p, m+1, L};
+
int* A1 = (int*)malloc(dimA[0]*dimA[1]*dimA[2]*sizeof(i
nt));
// A2
// A2
-
Int* A2 = (Int*)malloc(dimA[0]*dimA[1]*dimA[2]*sizeof(I
nt));
+
int* A2 = (int*)malloc(dimA[0]*dimA[1]*dimA[2]*sizeof(i
nt));
// rho
// rho
- const
mwSize
dimRho[] = {m, m, k, L};
-
Real* Rho = (Real*)malloc(dimRho[0]*dimRho[1]*dimRho[2]*dimRho[3]*sizeof(Real
));
+ const
int
dimRho[] = {m, m, k, L};
+
float* Rho = (float*)malloc(dimRho[0]*dimRho[1]*dimRho[2]*dimRho[3]*sizeof(float
));
// pi
// pi
- const
mwSize
dimPi[] = {k, L};
-
Real* Pi = (Real*)malloc(dimPi[0]*dimPi[1]*sizeof(Real
));
+ const
int
dimPi[] = {k, L};
+
float* Pi = (float*)malloc(dimPi[0]*dimPi[1]*sizeof(float
));
//////////////////////////////////////////////
// Main call to constructionModelesLassoMLE //
//////////////////////////////////////////////
// Main call to constructionModelesLassoMLE //
@@
-107,25
+107,25
@@
mwSize main(mwSize argc, char** argv)
free(Y);
// Compare to reference outputs
free(Y);
// Compare to reference outputs
-
I
nt* ref_A1 = readArray_int("A1",dimA, 3);
+
i
nt* ref_A1 = readArray_int("A1",dimA, 3);
compareArray_int("A1", A1, ref_A1, dimA[0]*dimA[1]*dimA[2]);
free(A1);
free(ref_A1);
// A2
compareArray_int("A1", A1, ref_A1, dimA[0]*dimA[1]*dimA[2]);
free(A1);
free(ref_A1);
// A2
-
I
nt* ref_A2 = readArray_int("A2",dimA, 3);
+
i
nt* ref_A2 = readArray_int("A2",dimA, 3);
compareArray_int("A2", A2, ref_A2, dimA[0]*dimA[1]*dimA[2]);
free(A2);
free(ref_A2);
// Rho
compareArray_int("A2", A2, ref_A2, dimA[0]*dimA[1]*dimA[2]);
free(A2);
free(ref_A2);
// Rho
-
Real
* ref_Rho = readArray_real("Rho",dimRho, 4);
+
float
* ref_Rho = readArray_real("Rho",dimRho, 4);
compareArray_real("Rho", Rho, ref_Rho, dimRho[0]*dimRho[1]*dimRho[2]*dimRho[3]);
free(Rho);
free(ref_Rho);
// Pi
compareArray_real("Rho", Rho, ref_Rho, dimRho[0]*dimRho[1]*dimRho[2]*dimRho[3]);
free(Rho);
free(ref_Rho);
// Pi
-
Real
* ref_Pi = readArray_real("Pi",dimPi, 2);
+
float
* ref_Pi = readArray_real("Pi",dimPi, 2);
compareArray_real("Pi", Pi, ref_Pi, dimPi[0]*dimPi[1]);
free(Pi);
free(ref_Pi);
compareArray_real("Pi", Pi, ref_Pi, dimPi[0]*dimPi[1]);
free(Pi);
free(ref_Pi);