From 8431dcc356a53917e94de06b8f08b0a30df406d0 Mon Sep 17 00:00:00 2001 From: sjplimp Date: Sat, 28 Aug 2010 00:01:35 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@4609 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- src/REAX/pair_reax.cpp | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/REAX/pair_reax.cpp b/src/REAX/pair_reax.cpp index 094eee1f97..5191d35f9c 100644 --- a/src/REAX/pair_reax.cpp +++ b/src/REAX/pair_reax.cpp @@ -118,7 +118,8 @@ void PairREAX::compute(int eflag, int vflag) evdwl = ecoul = 0.0; if (eflag || vflag) ev_setup(eflag,vflag); - else evflag = vflag_fdotr = 0; + else evflag = vflag_fdotr = eflag_global = vflag_global = + eflag_atom = vflag_atom = 0; if (vflag_global) FORTRAN(cbkvirial, CBKVIRIAL).Lvirial = 1; else FORTRAN(cbkvirial, CBKVIRIAL).Lvirial = 0; @@ -181,8 +182,7 @@ void PairREAX::compute(int eflag, int vflag) // extract global and per-atom energy from ReaxFF Fortran // compute_charge already contributed to eatom - if (eflag && eflag_global) { - // output_itemized_energy(energy_charge_equilibration); + if (eflag_global) { evdwl += FORTRAN(cbkenergies, CBKENERGIES).eb; evdwl += FORTRAN(cbkenergies, CBKENERGIES).ea; evdwl += FORTRAN(cbkenergies, CBKENERGIES).elp; @@ -203,7 +203,7 @@ void PairREAX::compute(int eflag, int vflag) eng_coul += ecoul; } - if (eflag && eflag_atom) { + if (eflag_atom) { int ntotal = atom->nlocal + atom->nghost; for (i = 0; i < ntotal; i++) eatom[i] += FORTRAN(cbkd,CBKD).estrain[i]; @@ -211,7 +211,7 @@ void PairREAX::compute(int eflag, int vflag) // extract global and per-atom virial from ReaxFF Fortran - if (vflag && vflag_global) { + if (vflag_global) { virial[0] = -FORTRAN(cbkvirial, CBKVIRIAL).virial[0]; virial[1] = -FORTRAN(cbkvirial, CBKVIRIAL).virial[1]; virial[2] = -FORTRAN(cbkvirial, CBKVIRIAL).virial[2]; @@ -220,7 +220,7 @@ void PairREAX::compute(int eflag, int vflag) virial[5] = -FORTRAN(cbkvirial, CBKVIRIAL).virial[5]; } - if (vflag && vflag_atom) { + if (vflag_atom) { int ntotal = atom->nlocal + atom->nghost; j = 0; for (i = 0; i < ntotal; i++) {