From c60a49ec71004bad2af9ca8b3240c035a7ab8d2a Mon Sep 17 00:00:00 2001
From: mihaf <miha.fosnaric@gmail.com>
Date: Tue, 03 Jun 2014 10:04:57 +0000
Subject: [PATCH] Added random seed.
---
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