39   double delx1,dely1,delz1,delx2,dely2,delz2,theta,dtheta,dcostheta,tk;  in compute()  local
63     delz1 = x[i1][2] - x[i2][2];  in compute()
65     rsq1 = delx1*delx1 + dely1*dely1 + delz1*delz1;  in compute()
79     c = delx1*delx2 + dely1*dely2 + delz1*delz2;  in compute()
118     f1[2] = (a11*delz1 + a12*delz2)*c0 + (b11*delz1 + b12*delz2)*s0;  in compute()
121     f3[2] = (a22*delz2 + a12*delz1)*c0 + (b22*delz2 + b12*delz1)*s0;  in compute()
144                          delx1,dely1,delz1,delx2,dely2,delz2);  in compute()
156   double delz1 = x[i1][2] - x[i2][2];  in single()  local
157   domain->minimum_image(delx1,dely1,delz1);  in single()
158   double r1 = sqrt(delx1*delx1 + dely1*dely1 + delz1*delz1);  in single()
166   double c = delx1*delx2 + dely1*dely2 + delz1*delz2;  in single()