Seems that the bugs are fixed. Pegs are self avoidant now

Samo Penic

2016-07-13 82a8abcddd62838bfcd06b6081042caf4d385a88

 src/stats.c

@@ -1,3 +1,4 @@
/* vim: set ts=4 sts=4 sw=4 noet : */
#include <stdio.h>
#include <time.h>
//#include <gsl/gsl_math.h>
@@ -15,20 +16,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++){