From dc77e8ef10378188ee0a9657cedd1b7991a0a29f Mon Sep 17 00:00:00 2001 From: Samo Penic <samo@CAE-linux.(none)> Date: Wed, 09 Apr 2014 21:27:38 +0000 Subject: [PATCH] Added spherical harmonics as-is into rbc --- src/sh.c | 25 +++++++++++++++++++++++++ 1 files changed, 25 insertions(+), 0 deletions(-) diff --git a/src/sh.c b/src/sh.c index 0c4b06d..9f28d2f 100644 --- a/src/sh.c +++ b/src/sh.c @@ -379,3 +379,28 @@ sph->N++; return TS_SUCCESS; } + + +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