X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=src%2Fsources%2FEMGLLF.c;h=86b606053076242cb6a623861990a4f7517089c4;hp=7e7a3d1d5938d22c607f20e18460a8c67d0b7c02;hb=46a2e676b2d85eef5a1811a6e623b65327fc451d;hpb=ef67d338c7f28ba041abe40ca9a8ab69f8365a90

diff --git a/src/sources/EMGLLF.c b/src/sources/EMGLLF.c
index 7e7a3d1..86b6060 100644
--- a/src/sources/EMGLLF.c
+++ b/src/sources/EMGLLF.c
@@ -92,7 +92,7 @@ void EMGLLF_core(
 				{
 					Real dotProduct = 0.;
 					for (int v=0; v<n; v++)
-						dotProduct += X2[ai(v,u,r,n,m,k)] * Y2[ai(v,mm,r,n,m,k)];
+						dotProduct += X2[ai(v,u,r,n,p,k)] * Y2[ai(v,mm,r,n,m,k)];
 					ps2[ai(u,mm,r,p,m,k)] = dotProduct;
 				}
 			}
@@ -317,7 +317,7 @@ void EMGLLF_core(
 				Real detRhoR = gsl_linalg_LU_det(matrix, signum);
 
 				//FIXME: det(rho[,,r]) too small(?!). See EMGLLF.R
-				Gam[mi(i,r,n,k)] = pi[r] * exp(-0.5*sqNorm2[r] + shift) * detRhoR;
+				Gam[mi(i,r,n,k)] = pi[r] * exp(-0.5*sqNorm2[r] + shift) ; //* detRhoR;
 				sumLLF1 += Gam[mi(i,r,n,k)] / gaussConstM;
 				sumGamI += Gam[mi(i,r,n,k)];
 			}