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

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)]