From 3197854e89eb9ffe5ff3137e274d16806c2028a5 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Tue, 27 Jul 2021 12:51:09 +0000
Subject: [PATCH] Found one bug in the code when calculating mprod. Fixed.
---
src/stats.c | 11 +++++++++--
1 files changed, 9 insertions(+), 2 deletions(-)
diff --git a/src/stats.c b/src/stats.c
index 4648d29..b63fecc 100644
--- a/src/stats.c
+++ b/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++){
--
Gitblit v1.9.3