Nekaj malega. Ne veva več kaj.

mihaf

2014-04-11 611412f383b15904cb0339e77b2f1f6c29453557

 src/stats.c

@@ -15,10 +15,10 @@
   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;
@@ -29,6 +29,12 @@
            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++){