Trisurf Monte Carlo simulator
Samo Penic
2014-09-11 464443e958f806fdc2b64ae77661a67c199b6204
src/stats.c
@@ -15,20 +15,26 @@
   ts_uint i,j,k;
   ts_double mat[3][3];
   ts_double vec[3];
   for(i = 0; i < 3; i++)
           for(j = 0; j < 3; j++)
                  mat[i][j]=0;
 
   for(k=0;k<vesicle->vlist->n;k++){
      vec[0]=vesicle->vlist->vtx[k]->x;
      vec[1]=vesicle->vlist->vtx[k]->y;
      vec[2]=vesicle->vlist->vtx[k]->z;
      for(i = 0; i < 3; i++)
              for(j = 0; j <= i; j++)
            mat[i][j]+=vec[i]*vec[j];
            mat[i][j]+=vec[i]*vec[j];
   }
// Normalize gyration tensor:
   for(i = 0; i < 3; i++)
           for(j = 0; j <= i; j++)
         mat[i][j]=mat[i][j]/(ts_double)vesicle->vlist->n;
// diagonal elements are copied twice!   
   for(i = 0; i < 3; i++)
           for(j = 0; j <= i; j++){