projects
/
valse.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
| inline |
side by side
(parent:
57f4468
)
fix typo in EMGLLF.c
author
Benjamin Auder
<benjamin.auder@somewhere>
Thu, 23 Mar 2017 10:25:49 +0000
(11:25 +0100)
committer
Benjamin Auder
<benjamin.auder@somewhere>
Thu, 23 Mar 2017 10:25:49 +0000
(11:25 +0100)
pkg/src/sources/EMGLLF.c
patch
|
blob
|
blame
|
history
diff --git
a/pkg/src/sources/EMGLLF.c
b/pkg/src/sources/EMGLLF.c
index
a028919
..
7ce6ded
100644
(file)
--- 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(
n
Y2[,mm,r])
+ //nY2[mm,r] = sum(Y2[,mm,r])
Real sumNy2 = 0.;
for (int u=0; u<n; u++)
- sumNy2 +=
n
Y2[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)]