From 4f5ffcfd5ba38b6b7dd6d3b15de8bb8677537b9f Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Thu, 09 Apr 2020 16:15:56 +0000 Subject: [PATCH] Commiting old changes --- src/shcomplex.c | 25 ++++++++++++++++++++++++- 1 files changed, 24 insertions(+), 1 deletions(-) diff --git a/src/shcomplex.c b/src/shcomplex.c index 7d8c9d8..4090346 100644 --- a/src/shcomplex.c +++ b/src/shcomplex.c @@ -1,5 +1,7 @@ +/* vim: set ts=4 sts=4 sw=4 noet : */ #include<math.h> #include<stdlib.h> +#include<string.h> #include<gsl/gsl_complex.h> #include<gsl/gsl_complex_math.h> #include<gsl/gsl_sf_legendre.h> @@ -68,7 +70,7 @@ if(sph->co != NULL) free(sph->co); if(sph->ulm !=NULL) free(sph->ulm); if(sph->ulmComplex !=NULL) free(sph->ulmComplex); - + if(sph->sumUlm2 !=NULL) free(sph->sumUlm2); if(sph->Ylmi!=NULL) { for(i=0;i<sph->l;i++){ if(sph->Ylmi[i]!=NULL){ @@ -116,6 +118,27 @@ return TS_SUCCESS; } +char *Ulm2Complex2String(ts_vesicle *vesicle){ + ts_int i,j; + char *strng=(char *)calloc(5000, sizeof(char)); + char tmpstrng[255]; + for(i=0;i<vesicle->sphHarmonics->l;i++){ + for(j=i;j<2*i+1;j++){ + sprintf(tmpstrng,"%e ", gsl_complex_abs2(vesicle->sphHarmonics->ulmComplex[i][j])); + strcat(strng,tmpstrng); + } + } + //strcat(strng,"\n"); + + return strng; +} + +ts_bool freeUlm2String(char *strng){ + free(strng); + return TS_SUCCESS; +} + + ts_bool storeUlmComplex2(ts_vesicle *vesicle){ ts_spharm *sph=vesicle->sphHarmonics; -- Gitblit v1.9.3