X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2Fsrc%2Fsources%2FEMGLLF.c;h=7ce6ded7cdcf0bff99786370f9e61dc778a3b608;hp=a028919defff9a15415607c55fb833d01c8531f4;hb=5a190d9d42b0313eace2e85b145b19f74542840f;hpb=57f4468ee9e8395bd78484e73ce467613f5269ef

diff --git a/pkg/src/sources/EMGLLF.c b/pkg/src/sources/EMGLLF.c
index a028919..7ce6ded 100644
--- a/pkg/src/sources/EMGLLF.c
+++ b/pkg/src/sources/EMGLLF.c
@@ -221,10 +221,10 @@ void EMGLLF_core(
 				for (int u=0; u<n; u++)
 					sumPs1 += ps1[ai(u,mm,r,n,m,k)];
 				ps[mi(mm,r,m,k)] = sumPs1;
-				//nY2[mm,r] = sum(nY2[,mm,r])
+				//nY2[mm,r] = sum(Y2[,mm,r])
 				Real sumNy2 = 0.;
 				for (int u=0; u<n; u++)
-					sumNy2 += nY2[ai(u,mm,r,n,m,k)];
+					sumNy2 += Y2[ai(u,mm,r,n,m,k)];
 				nY2[mi(mm,r,m,k)] = sumNy2;
 				//rho[mm,mm,r] = (ps[mm,r]+sqrt(ps[mm,r]^2+4*nY2[mm,r]*(gam2[r]))) / (2*nY2[mm,r])
 				rho[ai(mm,mm,r,m,m,k)] = ( ps[mi(mm,r,m,k)] + sqrt( ps[mi(mm,r,m,k)]*ps[mi(mm,r,m,k)]