52   double qtmp,xtmp,ytmp,ztmp,delx,dely,delz,evdwl,ecoul,fpair;  in compute()  local
154         fpair = (forcecoul + factor_lj*forcelj) * r2inv;  in compute()
156         f[i][0] += delx*fpair;  in compute()
157         f[i][1] += dely*fpair;  in compute()
158         f[i][2] += delz*fpair;  in compute()
160           f[j][0] -= delx*fpair;  in compute()
161           f[j][1] -= dely*fpair;  in compute()
162           f[j][2] -= delz*fpair;  in compute()
183                              evdwl,ecoul,fpair,delx,dely,delz);  in compute()
196   double qtmp,xtmp,ytmp,ztmp,delx,dely,delz,fpair;  in compute_inner()  local
260         fpair = (forcecoul + factor_lj*forcelj) * r2inv;  in compute_inner()
263           fpair  *= 1.0 + rsw*rsw*(2.0*rsw-3.0);  in compute_inner()
266         f[i][0] += delx*fpair;  in compute_inner()
267         f[i][1] += dely*fpair;  in compute_inner()
268         f[i][2] += delz*fpair;  in compute_inner()
270           f[j][0] -= delx*fpair;  in compute_inner()
271           f[j][1] -= dely*fpair;  in compute_inner()
272           f[j][2] -= delz*fpair;  in compute_inner()
284   double qtmp,xtmp,ytmp,ztmp,delx,dely,delz,fpair;  in compute_middle()  local
352         fpair = (forcecoul + factor_lj*forcelj) * r2inv;  in compute_middle()
355           fpair *= rsw*rsw*(3.0 - 2.0*rsw);  in compute_middle()
359           fpair *= 1.0 + rsw*rsw*(2.0*rsw - 3.0);  in compute_middle()
362         f[i][0] += delx*fpair;  in compute_middle()
363         f[i][1] += dely*fpair;  in compute_middle()
364         f[i][2] += delz*fpair;  in compute_middle()
366           f[j][0] -= delx*fpair;  in compute_middle()
367           f[j][1] -= dely*fpair;  in compute_middle()
368           f[j][2] -= delz*fpair;  in compute_middle()
380   double qtmp,xtmp,ytmp,ztmp,delx,dely,delz,evdwl,ecoul,fpair;  in compute_outer()  local
490         fpair = (forcecoul + forcelj) * r2inv;  in compute_outer()
492         f[i][0] += delx*fpair;  in compute_outer()
493         f[i][1] += dely*fpair;  in compute_outer()
494         f[i][2] += delz*fpair;  in compute_outer()
496           f[j][0] -= delx*fpair;  in compute_outer()
497           f[j][1] -= dely*fpair;  in compute_outer()
498           f[j][2] -= delz*fpair;  in compute_outer()
554           fpair = (forcecoul + factor_lj*forcelj) * r2inv;  in compute_outer()
558                              evdwl,ecoul,fpair,delx,dely,delz);  in compute_outer()