projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
fix typos
[valse.git]
/
src
/
sources
/
EMGLLF.c
diff --git
a/src/sources/EMGLLF.c
b/src/sources/EMGLLF.c
index
6966e9c
..
42419ac
100644
(file)
--- a/
src/sources/EMGLLF.c
+++ b/
src/sources/EMGLLF.c
@@
-3,7
+3,7
@@
#include <gsl/gsl_linalg.h>
// TODO: don't recompute indexes every time......
#include <gsl/gsl_linalg.h>
// TODO: don't recompute indexes every time......
-void EMGLLF(
+void EMGLLF
_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é
@@
-281,7
+281,7
@@
void EMGLLF(
{
YiRhoR[u] = 0.0;
for (int v=0; v<m; v++)
{
YiRhoR[u] = 0.0;
for (int v=0; v<m; v++)
- YiRhoR[u] += Y[
i
mi(i,v,n,m)] * rho[ai(v,u,r,m,m,k)];
+ YiRhoR[u] += Y[mi(i,v,n,m)] * rho[ai(v,u,r,m,m,k)];
}
//compute X(i,:)*phi(:,:,r)
}
//compute X(i,:)*phi(:,:,r)