From fbcbdf4093644e43688d9a44d6553d7153d8a2eb Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Tue, 29 Apr 2014 11:06:35 +0000
Subject: [PATCH] Constant volume in bondflip need lots of debugging
---
src/stats.c | 10 ++++++++--
1 files changed, 8 insertions(+), 2 deletions(-)
diff --git a/src/stats.c b/src/stats.c
index 4648d29..077781e 100644
--- a/src/stats.c
+++ b/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++){
--
Gitblit v1.9.3