forked from lijiext/lammps
git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@5041 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
parent
6dd3d3e403
commit
93e85b54d6
|
@ -126,12 +126,12 @@ void PairDipoleCut::compute(int eflag, int vflag)
|
||||||
|
|
||||||
// atom can have both a charge and dipole
|
// atom can have both a charge and dipole
|
||||||
// i,j = charge-charge, dipole-dipole, dipole-charge, or charge-dipole
|
// i,j = charge-charge, dipole-dipole, dipole-charge, or charge-dipole
|
||||||
|
|
||||||
|
forcecoulx = forcecouly = forcecoulz = 0.0;
|
||||||
|
tixcoul = tiycoul = tizcoul = 0.0;
|
||||||
|
tjxcoul = tjycoul = tjzcoul = 0.0;
|
||||||
|
|
||||||
if (rsq < cut_coulsq[itype][jtype]) {
|
if (rsq < cut_coulsq[itype][jtype]) {
|
||||||
|
|
||||||
forcecoulx = forcecouly = forcecoulz = 0.0;
|
|
||||||
tixcoul = tiycoul = tizcoul = 0.0;
|
|
||||||
tjxcoul = tjycoul = tjzcoul = 0.0;
|
|
||||||
|
|
||||||
if (qtmp != 0.0 && q[j] != 0.0) {
|
if (qtmp != 0.0 && q[j] != 0.0) {
|
||||||
r3inv = r2inv*rinv;
|
r3inv = r2inv*rinv;
|
||||||
|
|
Loading…
Reference in New Issue