From 443ba1603235e6cf668d0e24aa794c95941bc597 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Mon, 09 Jul 2012 08:06:53 +0000
Subject: [PATCH] Fix in dipyramid initialization, where we had 0.25*25 as xk
---
src/spherical_trisurf.c | 3 ++-
src/initial_distribution.c | 2 +-
src/sh.c | 2 +-
3 files changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/initial_distribution.c b/src/initial_distribution.c
index 2444e15..382750b 100644
--- a/src/initial_distribution.c
+++ b/src/initial_distribution.c
@@ -20,7 +20,7 @@
ts_vesicle *vesicle=init_vesicle(no_vertices,ncmax1,ncmax2,ncmax3,stepsize);
//TODO: debugging only. Please remove ASAP!
- vesicle->bending_rigidity=25.0*25.0;
+ vesicle->bending_rigidity=25.0;
vesicle->nshell=nshell;
retval = vtx_set_global_values(vesicle);
diff --git a/src/sh.c b/src/sh.c
index 31df74f..5952d39 100644
--- a/src/sh.c
+++ b/src/sh.c
@@ -371,7 +371,7 @@
ts_int i,j;
for(i=0;i<sph->l;i++){
for(j=0;j<2*i+1;j++){
- sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
+ sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
}
}
sph->N++;
diff --git a/src/spherical_trisurf.c b/src/spherical_trisurf.c
index f81e45d..f266c11 100644
--- a/src/spherical_trisurf.c
+++ b/src/spherical_trisurf.c
@@ -43,7 +43,7 @@
-for(i=0;i<10;i++){
+for(i=0;i<10000;i++){
centermass(vesicle);
cell_occupation(vesicle);
for(j=0;j<1000;j++){
@@ -89,6 +89,7 @@
fprintf(fh,"%d,\t%d,\t%e\n", i, j-i, sph->sumUlm2[i][j]/(ts_double)sph->N);
}
+ fprintf(fh,"\n");
}
fclose(fh);
return TS_SUCCESS;
--
Gitblit v1.9.3