Fix force calculation time in LAMMPS variant

Signed-off-by: Rafael Ravedutti <rafaelravedutti@gmail.com>
This commit is contained in:
Rafael Ravedutti 2022-03-17 02:53:58 +01:00
parent 887f41871c
commit 5df544637f
2 changed files with 3 additions and 2 deletions

View File

@ -9,3 +9,4 @@ v_out_freq 5
cutforce 0.9 cutforce 0.9
reneigh_every 100 reneigh_every 100
nstat 125000 nstat 125000
half_neigh 0

View File

@ -258,9 +258,9 @@ int main(int argc, char** argv) {
timer[FORCE] += computeForceEam(&eam, &param, &atom, &neighbor, &stats); timer[FORCE] += computeForceEam(&eam, &param, &atom, &neighbor, &stats);
} else { } else {
if(param.half_neigh) { if(param.half_neigh) {
timer[FORCE] = computeForceLJHalfNeigh(&param, &atom, &neighbor, &stats); timer[FORCE] += computeForceLJHalfNeigh(&param, &atom, &neighbor, &stats);
} else { } else {
timer[FORCE] = computeForceLJFullNeigh(&param, &atom, &neighbor, &stats); timer[FORCE] += computeForceLJFullNeigh(&param, &atom, &neighbor, &stats);
} }
} }