66   double delx1,dely1,delz1,delx2,dely2,delz2;  in compute()  local
100     delz2 = x[i3][2] - x[i2][2];  in compute()
102     rsq2 = delx2*delx2 + dely2*dely2 + delz2*delz2;  in compute()
107     c = delx1*delx2 + dely1*dely2 + delz1*delz2;  in compute()
127     f1[2] = a11*delz1 + a12*delz2;  in compute()
130     f3[2] = a22*delz2 + a12*delz1;  in compute()
153                          delx1,dely1,delz1,delx2,dely2,delz2);  in compute()
273   double delz2 = x[i3][2] - x[i2][2];  in single()  local
274   domain->minimum_image(delx2,dely2,delz2);  in single()
275   double r2 = sqrt(delx2*delx2 + dely2*dely2 + delz2*delz2);  in single()
277   double c = delx1*delx2 + dely1*dely2 + delz1*delz2;  in single()