1 typedef float rvec[3];
calc_similar_ind(int nind,rvec x[])2 float calc_similar_ind (int nind, rvec x[])
3 {
4   int i, j, d;
5   float m, tm, xd, rd;
6   for(j=0; j<nind; j++) {
7     i = x[j][0];
8     tm += m;
9     for(d=0 ; d<3; d++) {
10       xd = x[i][d] - x[i][d];
11       rd += m * xd;
12     }
13   }
14     return rd/tm;
15 }
16