From: Benjamin Auder <benjamin.auder@somewhere> Date: Thu, 23 Mar 2017 10:25:49 +0000 (+0100) Subject: fix typo in EMGLLF.c X-Git-Url: https://git.auder.net/%7B%7B%20asset%28%27mixstore/images/assets/current/doc/%7B%7B%20pkg.url%20%7D%7D?a=commitdiff_plain;h=5a190d9d42b0313eace2e85b145b19f74542840f;p=valse.git fix typo in EMGLLF.c --- 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)]