Readded r=0 check
This commit is contained in:
@ -269,7 +269,7 @@ void ComputeSNAGrid::compute_array()
|
|||||||
int jelem = 0;
|
int jelem = 0;
|
||||||
if (chemflag)
|
if (chemflag)
|
||||||
jelem = map[jtype];
|
jelem = map[jtype];
|
||||||
if (rsq < cutsq[jtype][jtype]) {
|
if (rsq < cutsq[jtype][jtype] && rsq > 1e-20) {
|
||||||
snaptr->rij[ninside][0] = delx;
|
snaptr->rij[ninside][0] = delx;
|
||||||
snaptr->rij[ninside][1] = dely;
|
snaptr->rij[ninside][1] = dely;
|
||||||
snaptr->rij[ninside][2] = delz;
|
snaptr->rij[ninside][2] = delz;
|
||||||
|
|||||||
Reference in New Issue
Block a user