From e297565ae34d6cc62f41810a261adbeef576af7b Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Tue, 16 Feb 2016 15:42:49 +0000 Subject: [PATCH] Changed the way the information on the version is inserted into the code --- src/spherical_trisurf_ff.c | 29 +++-------------------------- 1 files changed, 3 insertions(+), 26 deletions(-) diff --git a/src/spherical_trisurf_ff.c b/src/spherical_trisurf_ff.c index 867f37d..157b0e4 100644 --- a/src/spherical_trisurf_ff.c +++ b/src/spherical_trisurf_ff.c @@ -44,10 +44,11 @@ preparationSh(vesicle,r0); calculateYlmi(vesicle); calculateUlm(vesicle); +ts_double vmsr,bfsr; for(i=0;i<500;i++){ cell_occupation(vesicle); for(j=0;j<1000;j++){ - single_timestep(vesicle); + single_timestep(vesicle,&vmsr,&bfsr); } centermass(vesicle); fprintf(stderr, "Preloop %d completed.\n",i+1); @@ -62,7 +63,7 @@ for(i=0;i<10000;i++){ cell_occupation(vesicle); for(j=0;j<1000;j++){ - single_timestep(vesicle); + single_timestep(vesicle,&vmsr,&bfsr); } centermass(vesicle); vesicle_volume(vesicle); @@ -86,27 +87,3 @@ } - -ts_bool saveAvgUlm2(ts_vesicle *vesicle){ - - FILE *fh; - - fh=fopen("sph2out.dat", "w"); - if(fh==NULL){ - err("Cannot open file %s for writing"); - return TS_FAIL; - } - - ts_spharm *sph=vesicle->sphHarmonics; - ts_int i,j; - fprintf(fh,"l,\tm,\tulm^2avg\n"); - for(i=0;i<sph->l;i++){ - for(j=0;j<2*i+1;j++){ - 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