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/spherical_trisurf.c | 24 +----------------------- 1 files changed, 1 insertions(+), 23 deletions(-) diff --git a/src/spherical_trisurf.c b/src/spherical_trisurf.c index 6c5946f..06a84f4 100644 --- a/src/spherical_trisurf.c +++ b/src/spherical_trisurf.c @@ -1,3 +1,4 @@ +/* vim: set ts=4 sts=4 sw=4 noet : */ #include<stdio.h> #include<math.h> #include "general.h" @@ -104,26 +105,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