projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
R package can now be installed (compilation OK)
[valse.git]
/
src
/
sources
/
constructionModelesLassoMLE.c
diff --git
a/src/sources/constructionModelesLassoMLE.c
b/src/sources/constructionModelesLassoMLE.c
index
272f01b
..
4ab62ad
100644
(file)
--- a/
src/sources/constructionModelesLassoMLE.c
+++ b/
src/sources/constructionModelesLassoMLE.c
@@
-5,7
+5,7
@@
#include <omp.h>
// TODO: comment on constructionModelesLassoMLE purpose
#include <omp.h>
// TODO: comment on constructionModelesLassoMLE purpose
-void constructionModelesLassoMLE(
+void constructionModelesLassoMLE
_core
(
// IN parameters
const double* phiInit, // parametre initial de moyenne renormalisé
const double* rhoInit, // parametre initial de variance renormalisé
// IN parameters
const double* phiInit, // parametre initial de moyenne renormalisé
const double* rhoInit, // parametre initial de variance renormalisé
@@
-83,7
+83,7
@@
void constructionModelesLassoMLE(
double* piLambda = (double*)malloc(k*sizeof(double));
double* LLF = (double*)malloc((maxi+1)*sizeof(double));
double* S = (double*)malloc(lengthA*m*k*sizeof(double));
double* piLambda = (double*)malloc(k*sizeof(double));
double* LLF = (double*)malloc((maxi+1)*sizeof(double));
double* S = (double*)malloc(lengthA*m*k*sizeof(double));
- EMGLLF(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,
+ EMGLLF
_core
(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,
phiLambda,rhoLambda,piLambda,LLF,S,
n,lengthA,m,k);
free(Xa);
phiLambda,rhoLambda,piLambda,LLF,S,
n,lengthA,m,k);
free(Xa);
@@
-138,7
+138,7
@@
void constructionModelesLassoMLE(
for (int mm=0; mm<lengthB; mm++)
{
for (int r=0; r<k; r++)
for (int mm=0; mm<lengthB; mm++)
{
for (int r=0; r<k; r++)
- phi[ai
( A2[ai4(j,0,lambdaIndex,p,m+1,L)]-1, b[mm], r, lambdaIndex, p, m, k, L)] = 0.0
;
+ phi[ai
4( A2[ai(j,0,lambdaIndex,p,m+1,L)]-1, b[mm], r, lambdaIndex, p, m, k, L)] = 0.
;
}
}
}
}