ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@9802 f3b2605a-c512-4ea7-a41b-209d697bcdaa

Browse Source
This commit is contained in:
pscrozi
2013-04-29 16:53:14 +00:00
parent bece38546a
commit 2cc636f1aa
2 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/KSPACE/pppm.cpp
View File

@ -1357,7 +1357,7 @@ void PPPM::compute_gf_denom()
gf_b[0] = 4.0 * (gf_b[0]*(l-m)*(l-m-0.5)); gf_b[0] = 4.0 * (gf_b[0]*(l-m)*(l-m-0.5));
} }
int ifact = 1; bigint ifact = 1;
for (k = 1; k < 2*order; k++) ifact *= k; for (k = 1; k < 2*order; k++) ifact *= k;
double gaminv = 1.0/ifact; double gaminv = 1.0/ifact;
for (l = 0; l < order; l++) gf_b[l] *= gaminv; for (l = 0; l < order; l++) gf_b[l] *= gaminv;

2
src/KSPACE/pppm_old.cpp
View File

@ -1194,7 +1194,7 @@ void PPPMOld::compute_gf_denom()
gf_b[0] = 4.0 * (gf_b[0]*(l-m)*(l-m-0.5)); gf_b[0] = 4.0 * (gf_b[0]*(l-m)*(l-m-0.5));
} }
int ifact = 1; bigint ifact = 1;
for (k = 1; k < 2*order; k++) ifact *= k; for (k = 1; k < 2*order; k++) ifact *= k;
double gaminv = 1.0/ifact; double gaminv = 1.0/ifact;
for (l = 0; l < order; l++) gf_b[l] *= gaminv; for (l = 0; l < order; l++) gf_b[l] *= gaminv;