projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
fix test.constructionModelesLassoMLE; TODO: det(rho[,,r]) ? numerical errors ? ....
[valse.git]
/
src
/
sources
/
EMGrank.c
diff --git
a/src/sources/EMGrank.c
b/src/sources/EMGrank.c
index
2422fc0
..
3a9bf94
100644
(file)
--- a/
src/sources/EMGrank.c
+++ b/
src/sources/EMGrank.c
@@
-132,7
+132,7
@@
void EMGrank_core(
{
Real dotProduct = 0.0;
for (int u=0; u<cardClustR; u++)
{
Real dotProduct = 0.0;
for (int u=0; u<cardClustR; u++)
- dotProduct += Xr[mi(u,j,n,p)] * Yr[mi(u,j,n,m)];
+ dotProduct += Xr[mi(u,j,n,p)] * Yr[mi(u,j
j
,n,m)];
tXrYr[mi(j,jj,p,m)] = dotProduct;
}
}
tXrYr[mi(j,jj,p,m)] = dotProduct;
}
}
@@
-281,7
+281,7
@@
void EMGrank_core(
for (int jj=0; jj<p; jj++)
{
for (int r=0; r<k; r++)
for (int jj=0; jj<p; jj++)
{
for (int r=0; r<k; r++)
- Phi[ai(j
,jj,r,p,m,k)] = phi[ai(j,j
j,r,p,m,k)];
+ Phi[ai(j
j,j,r,p,m,k)] = phi[ai(jj,
j,r,p,m,k)];
}
}
ite++;
}
}
ite++;