reorganize test folder: tests generation with MATLAB almost OK (but suboptimal method)

[valse.git] / src / test / generate_test_data / R_test_helpers / testEMGLLF.R

diff --git a/src/test/generate_test_data/R_test_helpers/testEMGLLF.R b/src/test/generate_test_data/R_test_helpers/testEMGLLF.R
new file mode 100644 (file)
index 0000000..0eb3483
--- /dev/null
+++ b/src/test/generate_test_data/R_test_helpers/testEMGLLF.R
@@ -0,0 +1,46 @@
+testEMGLLF = function(){
+  testFolder = 'data/'
+  
+  #get dimensions
+  dimensions = read.table(strcat(testFolder,'dimensions'), header=FALSE)
+  n = dimensions[1]
+  p = dimensions[2]
+  m = dimensions[3]
+  k = dimensions[4]
+  
+  #get all input arrays
+  phiInit = read.table(strcat(testFolder,'phiInit'), header=FALSE)
+  rhoInit = read.table(strcat(testFolder,'rhoInit'), header=FALSE)
+  piInit = t(read.table(strcat(testFolder,'piInit'), header=FALSE))
+  gamInit = read.table(strcat(testFolder,'gamInit'), header=FALSE)
+  mini = read.table(strcat(testFolder,'mini'), header=FALSE)
+  maxi = read.table(strcat(testFolder,'maxi'), header=FALSE)
+  gamma = read.table(strcat(testFolder,'gamma'), header=FALSE)
+  lambda = read.table(strcat(testFolder,'lambda'), header=FALSE)
+  X = rread.table(strcat(testFolder,'X'), header=FALSE)
+  Y = read.table(strcat(testFolder,'Y'), header=FALSE)
+  tau = read.table(strcat(testFolder,'tau'), header=FALSE)
+  
+  #run EMGLLF.c
+  EMG = .Call("EMGLLF_core",phiInit,rhoInit,piInit1,gamInit,mini,maxi,gamma,lambda,X,Y,tau)
+  phi = EMG$phi
+  rho = EMG$rho
+  pi = EMG$pi
+  LLF = EMG$LLF
+  S = EMG$S
+  
+  #get all stored outputs
+  ref_phi =read.table(strcat(testFolder,'phi'), header=FALSE)
+  ref_rho = read.table(strcat(testFolder,'rho'), header=FALSE)
+  ref_pi = read.table(strcat(testFolder,'pi'), header=FALSE)
+  ref_LLF = read.table(strcat(testFolder,'LLF'), header=FALSE)
+  ref_S = read.table(strcat(testFolder,'S'), header=FALSE)
+  
+  #check that output correspond to stored output
+  tol = 1e-5;
+  checkOutput('phi',phi,ref_phi,tol);
+  checkOutput('rho',rho,ref_rho,tol);
+  checkOutput('pi',pi,ref_pi,tol);
+  checkOutput('LLF',LLF,ref_LLF,tol);
+  checkOutput('S',S,ref_S,tol);
+}